Displaying 6 results from an estimated 6 matches for "anovamodel".
2012 Feb 06
1
multiple comparisons in nested design
...ution of a example:
#> Mesures <- as.data.frame(matrix(rnorm(45*1, mean=10, sd=1), ncol=1))
#> colnames(Mesures) <- "VR"
#> Mesures$trat <- as.factor(rep(1:5,rep(9,5)))
#> Mesures$patient <- as.factor(rep(1:15,rep(3,15)))
The Analysis of Variance could be:
#> AnovaModel.1 <-aov(VR ~ trat + trat/patient, data=Mesures)
#> summary(AnovaModel.1)
Df Sum Sq Mean Sq F value Pr(>F)
trat 4 1.580 0.3950 0.751 0.565
trat:patient 10 4.811 0.4811 0.915 0.533
Residuals 30 15.778 0.5259
But this methods causes pseudoreplication bec...
2012 Feb 04
5
Comparaciones múltiples en ANOVA anidadp
...> Medidas <- as.data.frame(matrix(rnorm(45*1, mean=10, sd=1), ncol=1))
> colnames(Medidas) <- "VR"
> Medidas$Tratamiento <- as.factor(rep(1:5,rep(9,5)))
> Medidas$Individuo <- as.factor(rep(1:15,rep(3,15)))
El análisis de la varianza sugerido podría ser:
> AnovaModel.1 <-aov(VR ~ Tratamiento + Tratamiento/Individuo,
data=Medidas)
> summary(AnovaModel.1)
Df Sum Sq Mean Sq F value Pr(>F)
Tratamiento 4 1.15 0.2871 0.25 0.907
Tratamiento:Individuo 10 12.39 1.2395 1.08 0.407
Residuals 30 34.44 1...
2012 Apr 20
4
Problem with Tukey test
I'm new using R im trying to do a tukey test, but when i see the results the
p value results in NA im guessing its because i have missing values but im
not sure how to fix it
AnovaModel.2 <- aov(area ~ trat, data=apilados)
> summary(AnovaModel.2)
Df Sum Sq Mean Sq F value Pr(>F)
trat 4 11847 2961.76 9.9905 1.500e-06 ***
Residuals 76 22531 296.46
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*'...
2005 Jan 26
1
Specification of factorial random-effects model
...77.02,81.62,
59.32,62.29,60.7,55.42,55.5,51.17,50.54,53.54,51.85,
64.5,63.6,65.19,55.07,50.26,53.73,54.57,47.8,48.8,91.56,
94.49,92.17,82.14,83.16,81.31,83.58,78.63,77.08,60.53,
60.79,58.57,51.28,52.9,51.54,49.15,48.97,51.61,59.44,
60.07,60.07,51.94,52.2,50.2,49.45,50.75,49.56)
anovamodel <- aov(y ~ c + Error(a*b))
summary(anovamodel)
lmemodel <- lme(y ~ c, random = ~ 1|a*b)
anova(lmemodel)
N.B. You have to load the package 'nlme' before lme() will run at all.
When it does run, it gives the error
Error in getGroups.data.frame(dataMix, groups) :
Inv...
2009 Oct 19
1
access elements of summary.aov?
Dear all
How do I access individual elements of a "summary.aov" object?
> data(iris)
> AnovaModel.1 <- aov(Sepal.Length ~ Species, data=iris)
> tmp <- summary(AnovaModel.1)
> tmp
Df Sum Sq Mean Sq F value Pr(>F)
Species 2 63.2 31.6 119 <2e-16 ***
Residuals 147 39.0 0.3
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0...
2010 Aug 31
1
anova and lm results differ
...4 24 6.03 0.4987 0.7367
Residuals 3205 38754 12.09
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
2. (I run this from Rcmdr, this is what my students would normally do)
>AnovaModel.2 <- (lm(eduyrs ~ cntry*edf, data=ESS1))
> Anova(AnovaModel.2)
Anova Table (Type II tests)
Response: eduyrs
Sum Sq Df F value Pr(>F)
cntry 38 1 3.1158 0.07763 .
edf 11010 4 227.6296 < 2e-16 ***
cntry:edf 24 4...