[R] Extracting part of anova.cca()
Nevil Amos
nevil.amos at gmail.com
Mon Jan 23 12:31:06 CET 2012
I have an anova.cca() output.
If it is possible I would like to extract the table it contains as a data frame
ie
> > myAnova
> Permutation test for rda under reduced model
> Terms added sequentially (first to last)
>
> Model: rda(formula = mygenind at tab ~ mds3dCS_NULL + mds3dTRE_25_2_CS25 + mydata$LAT.x + mydata$LONG.x + mydata$Hab_Config + mds3dTRE_25_100_CS25 + mydata$TreeCov + mydata$Site_No + mds3dSFW_EO_100_CS25 + mds3dSFW_EO_5000_CS25 + mydata$Landscape + mds3dSFW_TH_10_CS25, scale = T, na.action = "na.omit")
> Df Var F N.Perm Pr(>F)
> mds3dCS_NULL 3 6.527 2.6380 99 0.01 **
> mds3dTRE_25_2_CS25 3 5.596 2.2616 99 0.01 **
> mydata$LAT.x 1 1.892 2.2939 99 0.01 **
> mydata$LONG.x 1 1.891 2.2929 99 0.01 **
> mydata$Hab_Config 2 2.772 1.6807 99 0.01 **
> mds3dTRE_25_100_CS25 3 3.952 1.5970 99 0.01 **
> mydata$TreeCov 1 1.272 1.5418 99 0.01 **
> mydata$Site_No 1 1.345 1.6306 99 0.01 **
> mds3dSFW_EO_100_CS25 3 3.741 1.5120 99 0.01 **
> mds3dSFW_EO_5000_CS25 3 3.410 1.3783 99 0.01 **
> mydata$Landscape 8 9.300 1.4094 99 0.01 **
> mds3dSFW_TH_10_CS25 3 3.658 1.4785 99 0.01 **
> Residual 139 114.644
> ---
> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> >
I would like to extract as a table:
> Df Var F N.Perm Pr(>F)
> mds3dCS_NULL 3 6.527 2.6380 99 0.01 **
> mds3dTRE_25_2_CS25 3 5.596 2.2616 99 0.01 **
> mydata$LAT.x 1 1.892 2.2939 99 0.01 **
> mydata$LONG.x 1 1.891 2.2929 99 0.01 **
> mydata$Hab_Config 2 2.772 1.6807 99 0.01 **
> mds3dTRE_25_100_CS25 3 3.952 1.5970 99 0.01 **
> mydata$TreeCov 1 1.272 1.5418 99 0.01 **
> mydata$Site_No 1 1.345 1.6306 99 0.01 **
> mds3dSFW_EO_100_CS25 3 3.741 1.5120 99 0.01 **
> mds3dSFW_EO_5000_CS25 3 3.410 1.3783 99 0.01 **
> mydata$Landscape 8 9.300 1.4094 99 0.01 **
> mds3dSFW_TH_10_CS25 3 3.658 1.4785 99 0.01 **
> Residual 139 114.644
Is this possible? what is the syntax?
I have looked at str(myAnova) but cannot see how to get the table with variable names included
Many thanks
Nevil Amos
More information about the R-help
mailing list