[R-sig-eco] rarefied data output
Jari Oksanen
jari.oksanen at oulu.fi
Wed Mar 28 17:36:48 CEST 2012
On 28 Mar 2012, at 12:59, cristabel.duran wrote:
> Dear list,
>
> I have a species matrix with species abundances (spp) and in other
> matrix with environmental information (env) with several factors. The
> factor I am interested in (mhab5) has five levels.
> The spp matrix looks like that:
>
> Aba_sp Alc_gla Ann Ann_1uc Apo Asp_dec Asp_sp Bur Bur_1 Bur_1uc Bur_2
> P1 0 0 0 0 0 0 0 0 0 0 0
> P2 0 0 0 0 0 0 0 0 0 0 0
> P3 0 0 0 0 0 0 0 1 1 0 0
> P4 0 0 0 0 0 0 0 0 0 0 0
> P5 0 0 0 0 0 0 0 1 0 0 2
> P6 0 0 0 0 0 0 0 0 0 0 0
> P7 0 0 0 0 1 0 0 0 0 0 0
>
>
> Beside of testing for differences in species richness between the levels
> of the factor mhab5, I want to use the rarefied data for a model. In
> other words I need the rarefy richness for P1, P2, P3, and so on.....
>
> I used the function accumcomp (BiodiversityR) which internally uses the
> function rarefy (vegan)
> acumhab5 <- accumcomp(x=spp, y=env, factor="mhab5",
> method="rarefaction", conditioned=TRUE, legend=T)
>
> I am able to extract the rarefied species richness with
> richnessmhab5 <- t(acumhab5[ , , 2])
>
> I get this output:
>
> mhab5
> obs Border Branches Crown Trunk Undist.
> 1 15.66701 11.88617 8.582444 19.37374 16.85696
> 2 23.31066 18.17601 13.536749 27.53000 23.01568
> 3 28.83948 22.61587 17.677169 33.35105 26.53217
> 4 32.98145 26.06969 20.730801 37.93189 29.04925
> 5 36.39266 28.90257 23.480157 41.76465 30.97731
> 6 39.39713 31.30142 25.594925 45.01781 32.62629
> 7 41.92287 33.47224 27.560948 47.98189 34.00000
>
>
> With this output I am not able to know which results correspond to which
> plot (see in spp matrix: P1, P2, P3, etc)
> any clue how can I get the information in the order I need?
Cristabel,
I am not quite sure what you try to get. However, you used function that gives you the richness of *accumulated* plots (=several plots combined together). There is no way of knowing the rarefied richness of single plots if you ask for accumulated results. Now I really don't know what you want to get, but here two alternatives:
- if you want to know the rarefied richness of single plots, use rarefy().
- if you want to know which plots belong to which mhab5 data, look at your data.
Cheers, Jari Oksanen
More information about the R-sig-ecology
mailing list