2014-10-27
## Loading required package: splines
## Loading required package: RcmdrMisc
## Loading required package: car
## Loading required package: sandwich
## La interfaz R-Commander sólo funciona en sesiones interactivas
> melon2 <- read.table("/home/arn/Tesis/Analisis 2do muestreo/muestreo2",
+ header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE)
> library(agricolae, pos=16)
> summary(aov(area ~ concentracion*tipo, data=melon2))
Df Sum Sq Mean Sq F value Pr(>F)
concentracion 2 59061516 29530758 3.206 0.06441 .
tipo 1 113317517 113317517 12.303 0.00251 **
concentracion:tipo 2 67839247 33919623 3.683 0.04563 *
Residuals 18 165785460 9210303
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> modelo<-(aov(area ~ concentracion*tipo, data=melon2))
> HSD.test(modelo, "tipo", console=TRUE)
Study: modelo ~ "tipo"
HSD Test for area
Mean Square Error: 9210303
tipo, means
area std r Min Max
Foliar 12369.336 3268.095 12 8890.90 19116.59
Suelo 8023.504 3990.914 12 3384.62 14992.43
alpha: 0.05 ; Df Error: 18
Critical Value of Studentized Range: 2.971152
Honestly Significant Difference: 2602.983
Means with the same letter are not significantly different.
Groups, Treatments and means
a Foliar 12370
b Suelo 8024
> HSD.test(modelo, "concentracion", console=TRUE)
Study: modelo ~ "concentracion"
HSD Test for area
Mean Square Error: 9210303
concentracion, means
area std r Min Max
1x10-3molar 10180.86 4494.303 8 4907.17 19116.59
2x10-3molar 8282.96 4847.487 8 3384.62 17576.58
Agua 12125.44 2422.042 8 8890.90 14992.43
alpha: 0.05 ; Df Error: 18
Critical Value of Studentized Range: 3.609304
Honestly Significant Difference: 3872.714
Means with the same letter are not significantly different.
Groups, Treatments and means
a Agua 12130
a 1x10-3molar 10180
a 2x10-3molar 8283
> m<-HSD.test(modelo, "concentracion",console=TRUE)
Study: modelo ~ "concentracion"
HSD Test for area
Mean Square Error: 9210303
concentracion, means
area std r Min Max
1x10-3molar 10180.86 4494.303 8 4907.17 19116.59
2x10-3molar 8282.96 4847.487 8 3384.62 17576.58
Agua 12125.44 2422.042 8 8890.90 14992.43
alpha: 0.05 ; Df Error: 18
Critical Value of Studentized Range: 3.609304
Honestly Significant Difference: 3872.714
Means with the same letter are not significantly different.
Groups, Treatments and means
a Agua 12130
a 1x10-3molar 10180
a 2x10-3molar 8283
> m
$statistics
Mean CV MSerror HSD
10196.42 29.76386 9210303 3872.714
$parameters
Df ntr StudentizedRange
18 3 3.609304
$means
area std r Min Max
1x10-3molar 10180.86 4494.303 8 4907.17 19116.59
2x10-3molar 8282.96 4847.487 8 3384.62 17576.58
Agua 12125.44 2422.042 8 8890.90 14992.43
$comparison
NULL
$groups
trt means M
1 Agua 12125.44 a
2 1x10-3molar 10180.86 a
3 2x10-3molar 8282.96 a
> W<-HSD.test(modelo, "tipo", console=TRUE)
Study: modelo ~ "tipo"
HSD Test for area
Mean Square Error: 9210303
tipo, means
area std r Min Max
Foliar 12369.336 3268.095 12 8890.90 19116.59
Suelo 8023.504 3990.914 12 3384.62 14992.43
alpha: 0.05 ; Df Error: 18
Critical Value of Studentized Range: 2.971152
Honestly Significant Difference: 2602.983
Means with the same letter are not significantly different.
Groups, Treatments and means
a Foliar 12370
b Suelo 8024
> W
$statistics
Mean CV MSerror HSD
10196.42 29.76386 9210303 2602.983
$parameters
Df ntr StudentizedRange
18 2 2.971152
$means
area std r Min Max
Foliar 12369.336 3268.095 12 8890.90 19116.59
Suelo 8023.504 3990.914 12 3384.62 14992.43
$comparison
NULL
$groups
trt means M
1 Foliar 12369.336 a
2 Suelo 8023.504 b
> with(melon2, plotMeans(area, concentracion, tipo, error.bars="se",
+ xlab="Concentración de yoduro de potasio", ylab="Area (cm2)", main=""))
