[R] Understanding the output from lme anova

2010-12-06 Thread Hedberg Peter
I have a question regarding understanding the output from a nested linnear mixed effect model. My model looks like the following lme(Poaceae~Site+Ditch_block+Tree_cut+as.factor(Dist_start)+(as.factor(Dist_start)*Tree_cut)+((as.factor(Dist_ditch))/Ditch_block),random=~1|Mainplot/Transect/Obsplot) I

[R] Understanding the output from lme.anova

2010-12-06 Thread Hedberg Peter
I have a question regarding understanding the output from a nested linnear mixed effect model. My model looks like the following lme(Poaceae~Site+Ditch_block+Tree_cut+as.factor(Dist_start)+(as.factor(Dist_start)*Tree_cut)+((as.factor(Dist_ditch))/Ditch_block),random=~1|Mainplot/Transect/Obsplot)

Re: [R] Why do I get (not estimable) in this aov data

2010-12-02 Thread Hedberg Peter
variables for Transect is better. Thank you for the help. Best regards, Petter Peter Dalgaard skrev : > On Nov 30, 2010, at 13:58 , Hedberg Peter wrote: > > >> aov_data > > Call: > > aov(formula = Species1 ~ Site + Obstacle + Treatment + > > as.factor(D

[R] Why do I get (not estimable) in this aov data

2010-11-30 Thread Hedberg Peter
Dear all. I am trying to analyze my data with an anova, using aov, and the result looks like this > my_data<-read.csv(file.choose(),header=T,dec=".",sep=";") > attach(my_data) > names(my_data) [1] "Site" "Obstacle" "Treatment" "Dist_Obstacle" "Dist_start" "Transect"

[R] How do I get "total inertia"+"R-squared" in R?

2010-11-17 Thread Hedberg Peter
Dear people at the R-help-list. I have run into a problem in the package "vegan", while using the function "vegan", and would greatly appreciate any advice how to solve the problem. I have put in my data R, using the "decorana" function in the package "vegan" and received nice coordinates for th

[R] Total inertia in package Vegan?

2010-11-17 Thread Hedberg Peter
Dear all at the R-project help list. I have run into a problem when it comes to getting values for "total inertia" and "R-squared" for my DCA using "decorana" in the package "vegan". I have tried the "goodness" function, but the reply indicates that it does not work with "decorana" class object