Comandos Proyecto

Comandos para la elaboración de las gráficas 1 y 2

read.table("clipboard",h=T)->datos
datos
head(datos)
attach(datos)
aov(Conteo~Trata)->anova1
summary(anova1)
plot(anova1)
library(Rcmdr)
par(mfrow=c(2,2))
plot(anova1,xlab="Numero Celulas")
lm(Conteo~Trata)->lm1
shapiro.test(residuals(lm1))
read.table("clipboard",h=T)->a
a}

Comandos para la elaboración de los boxplots (Gráficas 3 y 4)

library(ggplot2)
p<-ggplot(a,aes(Tratamiento,Conteo))+
  geom_boxplot()
p
read.table("clipboard",h=T)->b
b
k<-ggplot(b, aes(Tratamiento,Conteo))+
  geom_boxplot()
k}

Calculo de las medias

ct<-c(3175000,4100000,4125000,2500000,2650000,2875000,1700000,900000,1450000,1200000,1175000,2850000)
ct
mean(ct)
c<-c(1900000, 2075000, 2025000, 2225000, 2150000, 2050000, 700000, 1850000, 1150000, 1050000, 925000, 925000)
c
mean(c)
st<-c(1575000, 1925000, 1925000, 1725000, 1675000, 1700000, 1525000, 1050000, 1400000, 2550000, 925000)
st
mean(st)
cn<-c(3175000, 2125000, 1725000, 2100000, 1875000, 2900000, 900000, 825000, 1450000, 1150000, 1000000, 325000)
cn
mean(cn)
vn<-c(1600000, 2875000, 2650000, 3225000, 1675000, 1700000, 2000000, 1650000, 1275000, 925000, 1175000, 1175000, 1400000)
vn
mean(vn)
dz<-c(3375000, 3250000, 2525000, 2500000, 2975000, 3925000, 1125000, 1550000, 800000, 850000, 525000, 2400000)
dz
mean(dz)}