DCA

La primera imagen es la representación gráfica generada por el software estadístico R, de comparación de medias, empleando la prueba de Tukey. La segunda imagen, muestra los numeros en tabla, generados por SAS, con medias agrupadas.

El Diseño Completamente Aleatorio (DCA) es el primer cuadrúpedo de batalla de los cursos de Diseños Experimentales. Los datos de estos resultados los tomamos de la Tabla 6.1, página 87, de Experimentación en Agricultura. Recomendamos lean el Capítulo 6. Agregamos códigos a los análisis que aparecen en esa sección, para dos programas más: R y SAS.

Explicación en página 91. 

Códigos de R

#***********************************************************
# EXPERIMENTACIÓN EN AGRICULTURA                              
# Fernández Escobar, R.; Trapero, A.; Domínguez, J. 2010       
# CAPÍTULO 6 - DISEÑO COMPLETAMENTE ALEATORIO                  
#************************************************************

# Importando datos
trigo <- read.csv("C:\\Users\\Administrator\\Desktop\\tabla 6.1.csv", header=TRUE)
attach(trigo)
trigo

# Estadísticas descriptivas (instalar "psych" y "doBy")
library(psych)

describeBy(trigo$rendimiento, variedad)

library(doBy)

summaryBy(rendimiento ~ variedad, data = trigo, 
          FUN = function(x){c(sum = sum(x), m = mean(x), min = min(x), max = max(x),
                                var = var(x), s = sd(x), SEM = sd(x)/sqrt(length(x)))})

# Prueba de Bartlett
bartlett.test(rendimiento ~ variedad, data = trigo)

# Análisis de varianza
anova <- aov(rendimiento ~ variedad, data=trigo)
summary(anova)

# Comparación de medias utilizando LSD y Tukey
library(agricolae)

LSD <- LSD.test(anova, "variedad", group = T)
LSD

TukeyHSD(anova)

plot(TukeyHSD(anova))

Códigos de SAS

/***********************************************************

 EXPERIMENTACIÓN EN AGRICULTURA                             

 Fernández Escobar, R.; Trapero, A.; Domínguez, J. 2010      

 CAPÍTULO 6 - DISEÑO COMPLETAMENTE ALEATORIO                 

************************************************************/

* Lectura de datos

  Tabla 6.1 Rendimiento (t/ha) de cuatro variedades de  trigo

            en cada una de las cuatro parcelas elementales;

DATA trigo;

  INPUT variedad $ rendimiento;

  LABEL variedad = 'Variedades de trigo'

        rendimiento = 'Rendimiento (t/ha)';

CARDS;

VAR1   4.58

VAR1   4.33

VAR1   4.6

VAR1   3.26

VAR2   3.57

VAR2   3.21

VAR2   3.59

VAR2   2.35

VAR3   5.68

VAR3   5.43

VAR3   5.1

VAR3   4.57

VAR4   4.68

VAR4   4.21

VAR4   4.36

VAR4   3.79

;

ODS HTML;

PROC PRINT;

PROC MEANS SUM MEAN MIN MAX VAR STD STDERR;

CLASS variedad;

RUN;

* ANOVA, separación de medias y test de igual de varianzas;

PROC ANOVA;

  CLASS variedad;

  MODEL rendimiento = variedad;

  MEANS variedad / HOVTEST = BARTLETT;

  MEANS variedad / LSD TUKEY;

  QUIT;

RUN;

ODS HTML CLOSE;

-------------------------------------------------------

Tabla 6.1, capítulo 6.R, capítulo 6.sas son los archivos que necesitan para correr los análisis en sus programas.