[R-sig-eco] Understanding the output from lme anova

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Mon Dec 6 15:57:27 CET 2010


Dear Peter,

A * B is just a shorthand notation for A + B + A:B The model will
auto-expand those term to the ":" notation.

Best regards,

Thierry

PS Creating the factors in your dataset and the use of spaces will
increase the readability of your code!

------------------------------------------------------------------------
----
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek
team Biometrie & Kwaliteitszorg
Gaverstraat 4
9500 Geraardsbergen
Belgium

Research Institute for Nature and Forest
team Biometrics & Quality Assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium

tel. + 32 54/436 185
Thierry.Onkelinx at inbo.be
www.inbo.be

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to
say what the experiment died of.
~ Sir Ronald Aylmer Fisher

The plural of anecdote is not data.
~ Roger Brinner

The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of
data.
~ John Tukey
  

> -----Oorspronkelijk bericht-----
> Van: r-sig-ecology-bounces at r-project.org 
> [mailto:r-sig-ecology-bounces at r-project.org] Namens Hedberg Peter
> Verzonden: maandag 6 december 2010 15:46
> Aan: r-sig-ecology at r-project.org
> Onderwerp: [R-sig-eco] Understanding the output from lme anova
> 
> Dear all.
> 
> I have a question regarding how to interpret the output from 
> my lme.anova.
> 
> My model looks like this
> 
> lme_Cyperaceae<-lme(Cyperaceae~Site+Ditch_block+Tree_cut+as.fa
> ctor(Dist_start)+((as.factor(Dist_start))*Tree_cut)+((as.facto
> r(Dist_ditch))/Ditch_block),random=~1|Mainplot/Transect/Obsplot)
> 
> My question concerns that I use the divission sign to 
> indicate nesting in ((as.factor(Dist_ditch))/Ditch_block), 
> and I use * to indicate two-way interaction in 
> ((as.factor(Dist_start))*Tree_cut), but in the output both 
> terms are presented with a colon (:) sign.
> 
> Tree_cut:as.factor(Dist_start)        3     6  0.259899  0.8520
> Ditch_block:as.factor(Dist_ditch)     4    56  2.695738  0.0399
> 
> Is there something wrong in how I wrote the model, or is the 
> colon sign  used for all interactions weather nested or the * type.
> 
> Would be greatly thankful for any help in understanding this.
> 
> The output mentioned above is from an analysis of lme.anova 
> on the result of the lme model.
> 
> Best regards,
> 
> Petter Hedberg
> 
> 
> 
> ___________________________________
> NOCC, http://nocc.sourceforge.net
> 
> 
> 
> 
> --
> This email was Anti Virus checked by Astaro Security Gateway. 
> http://www.astaro.com
> 
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
> 



More information about the R-sig-ecology mailing list