[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)

[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 pda...@gmail.com skrev : On Nov 30, 2010, at 13:58 , Hedberg Peter wrote: aov_data Call: aov(formula = Species1 ~ Site + Obstacle + Treatment + as.factor(Dist_Obstacle) + as.factor

[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 Mainplot

[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 objects. In

[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 the DCA