Hello, I am wondering how to test for the relative effects of different species in an abundance matrix that I am analyzing with manyglm in the mvabund package.
e.g., in the package example,
data(spider)
$abund
Alopacce Alopcune Alopfabr Arctlute Arctperi Auloalbi Pardlugu Pardmont
Pardnigr Pardpull Trocterr Zoraspin
1 25 10 0 0 0 4 0 60
12 45 57 4
2 0 2 0 0 0 30 1 1
15 37 65 9
3 15 20 2 2 0 9 1 29
18 45 66 1
4 2 6 0 1 0 24 1 7
29 94 86 25
5 1 20 0 2 0 9 1 2
135 76 91 17
6 0 6 0 6 0 6 0 11
27 24 63 34
7 2 7 0 12 0 16 1 30
89 105 118 16
8 0 11 0 0 0 7 55 2
2 1 30 3
9 1 1 0 0 0 0 0 26
1 1 2 0
10 3 0 1 0 0 0 0 22
0 0 1 0
11 15 1 2 0 0 1 0 95
0 1 4 0
12 16 13 0 0 0 0 0 96
1 8 13 0
13 3 43 1 2 0 18 1 24
53 72 97 22
14 0 2 0 1 0 4 3 14
15 72 94 32
15 0 0 0 0 0 0 6 0
0 0 25 3
16 0 3 0 0 0 0 6 0
2 0 28 4
17 0 0 0 0 0 0 2 0
0 0 23 2
18 0 1 0 0 0 0 5 0
0 0 25 0
19 0 1 0 0 0 0 12 0
1 0 22 3
20 0 2 0 0 0 0 13 0
0 0 22 2
21 0 1 0 0 0 0 16 1
0 1 18 2
22 7 0 16 0 4 0 0 2
0 0 1 0
23 17 0 15 0 7 0 2 6
0 0 1 0
24 11 0 20 0 5 0 0 3
0 0 0 0
25 9 1 9 0 0 2 1 11
6 0 16 6
26 3 0 6 0 18 0 0 0
0 0 1 0
27 29 0 11 0 4 0 0 1
0 0 0 0
28 15 0 14 0 1 0 0 6
0 0 2 0
is there a way analogous to the simper method for distance matrices to test
whether Alopacce is important in determining a significant effect of a
treatment factor
in a model that looked at manyglm(abund~treatment)?
Many thanks,
Beth
--
****************************************************
Elizabeth G. Pringle
Michigan Society of Fellows
University of Michigan
****************************************************
[[alternative HTML version deleted]]
_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
