Diseño 3:Factorial completo completamente al azar(FCCA)

Variable respuesta:Germinación
Factor 1:variedad con niveles:v1,v2,v3/
Factor 2:hormona con 2 niveles: a,p
No bloqueo
Total tratamientos:6 TRATAMIENTOS


 Ecuación del modelo----->Yijk=𝝻+Ti +Bj+(TB)ij+Eijk

 Donde :

y:respuesta en función de número de tratamientos y repeticiones
𝝻:Media global
Ti :Efecto de  los tratamientos/factor 1 --->efecto principal
Bj:Efecto factor 2--------->efecto principal
(TB)ij=Interacción entre factores-------->efecto interacción
Eijk :Error residual

Hipótesis 
Ho=𝝻T1=𝝻T2=𝝻T3=0("Los tratamientos son iguales"por lo tanto,no tienen efecto/efecto nulo sobre la variabilidad datos )
Anova parece ser bueno cuando la variabilidad se debe a los tratamientos
No concluir que no hay diferencias en los tratamientos 

No planteo hipotesis para bloqueo 
mejor tratamiento descriptivo ocn medias aggrgate

Comandos:

germ<-rnorm(120,70,2)
var<-gl(3,40,labels =c("v1","v2","v3"));var
proc<-gl(2,20,120,labels=c("l1","l2"));proc
horm<-gl(2,10,120,labels = c("A","P"));horm
data<-data.frame(germ,var,proc,horm)
attach(data)
View(data)
with(data,aggregate(germ,list(var,horm),mean))
modelo3 <- aov(germ ~var+horm+var:horm, data=data);modelo3
summary(modelo3)

Modelo 2:Fatorial simple en bloques

Factorial simple en bloques

Variable respuesta:Germinación
Factor:variedad con niveles:v1,v2,v3/hormona con 2 niveles a,p
bloqueo:localidades
Total tratamientos:3 TRATAMIENTOS


 Ecuación del modelo----->Yij=𝝻+Ti +Bj+Eij

 Donde :

y:respuesta en función de número de tratamientos y repeticiones
𝝻:Media global
B:Bloqueo
Ti :Efecto de  los tratamientos/factor(i=numero tratamientos)
Ej :Error residual(número de repeticiones)

Hipótesis 
Ho=𝝻T1=𝝻T2=𝝻T3=0("Los tratamientos son iguales"por lo tanto,no tienen efecto/efecto nulo sobre la variabilidad datos )
Anova parece ser bueno cuando la variabilidad se debe a los tratamientos
No concluir que no hay diferencias en los tratamientos 

No planteo hipotesis para bloqueo 
 
Comandos en R
germ<-rnorm(120,70,2)
var<-gl(3,40,labels =c("v1","v2","v3"));var
proc<-gl(2,20,120,labels=c("l1","l2"));proc
horm<-gl(2,10,120,labels = c("A","P"));horm
data<-data.frame(germ,var,proc,horm)
attach(data)
View(data)
modelo2 <- aov(germ ~ horm + proc, data=data);modelo2
summary(modelo2)
modelo2a <- aov(germ ~ var + proc, data=data);modelo2a
summary(modelo2a)

#Profe modelo 2a--->variedad :3,86%<5%----->Rechazo Ho--->Me quedo con Ha en ANOVA=media variedades distinas
#-------->¿Todas las variedades son diferentes o solo una de ellas?---->Utilizo tukey----->Mayor 5% NO RECHAZO HO
#------>Me quedo con Ho en tukey=Todas las variendades son iguales


#COMPARACION DE MEDIAS CON Tukey(Se hace para más de dos tratamientos)

TukeyHSD(modelo2a)
#Explicación porque en ANOVA me dio que las medias son distintas y en tukey que las medias son iguales:No es suficiente utilizar ANOVA,es mejor tukey.
#lsd,duncan

Modelo1:FSCA

MODELO 1:FACTORIAL SIMPLE COMPLETAMENTE AL AZAR


Primer ejemplo 


Variable respuesta:Aceite esencial extraido en limonaria (100g)-->AE
Factor:especie de limonaria con niveles:sp1,sp2,sp3(especies)----->f1
Sin bloqueo
Total tratamientos:3(niveles)


 Ecuación del modelo----->Yij=𝝻+Ti +Ej

 Donde :

y:respuesta en función de número de tratamientos y repeticiones
𝝻:Media global
Ti :Efecto de  los tratamientos
Ej :Error residual

Hipótesis 
Ho=𝝻T1=𝝻T2=𝝻T3=0("Los tratamientos son iguales"por lo tanto,no tienen efecto/efecto nulo sobre la variabilidad datos )
Anova parece ser bueno cuando la variabilidad se debe a los tratamientos
No concluir que no hay diferencias en los tratamientos
Comandos en R


set.seed(1964);options(digits=3)#OJO CON DECIMALES,PORQUE ME DICE LA EXACTITUD DE MEDIDAS
AE=sort (rnorm(120,8.5,0.8),decreasing = T);AE#SORT sirve para ordenar ya sea de menor a mayor o viceversa
f1=gl(4,30,labels = c("sp1","sp2","sp3","sp4"));f1
data=data.frame(AE,f1)
View(data)
MODELO1 <- aov(AE ~ f1, data=data)
summary(MODELO1)
 


 decimales,que equipo use------nuemros cifras apreciaicion de los equipos utilizados


Interpretación 



Segundo ejemplo


Variable respuesta:Germinación
f1=Hormona/variedad
 

Más conceptos(pendiente)

Conceptos

Datos crudos:Datos que no han sido transformados
Datos no crudos:Datos transformados,por lo general provienen de sumatorias o promedios.Como por ejemplo la precipitación,la radiación.

Las medidas más importantes en agronomía son el promedio(media,medida de tendencia central) y la varianza(medida de dispersión),normalmente,se le calculan a datos crudos.El promedio suele ser muy sensible a datos atipícos y por lo tanto muchas veces su valor no nos representa adecuadamente nuestros datos;por otra parte la varianza es importante ya que mis datos deben tener la menor variabilidad posible dentro de un mismo conjunto datos, ya que estos provienen de un mismo lugar y por lo tanto tienen que ser muy similares,allí,es donde cobra la importancia de la repetición de toma de datos(la variabilidad de datos solo puede ser posible en diferentes tratamientos).Cuando hay valores que no son similares estos reciben el nombre de datos atipícos,los cuales pueden ser producidos por errores intrumentales o humano.

Si hay muchos datos atipícos  se pueden eliminar o puede repetirse la toma de datos.

En R la varianza se calcula con 

 for

 como se que r me alculo varianza muestral o poblacional

sd---desviacion estandar

var---varianza

en R los anetriores son meustrales

varianza=sum cuadrados/grados libertad

Imágenes tomadas de :https://es.slideshare.net/helena0912/stats-muestreo

http://www.monografias.com/trabajos88/dispersion-relativa/image001.png 

ANOVA----->Análisis de varianza

El ANOVA me permite saber si mi varianza se debe a los tratamientos(ideal) o a las repeticiones.

Dicho ANOVA es diferente para cada diseño.

El ANOVA realiza una razón entre VarY/VarX,donde Y son los tratamientos y X el error.Dicha razón al graficarse se comporta como una curva de f .En el ANOVA se encuentra el f value y me dice cuanto son los tratamientos más variables que el error,donde lo que se desea es que dicho número sea grande .Sin embargo nuestro pr (f) es el valor con el que tomamos la decisión de la variabilidad de los datos .

 

Crear una tabla en R

¿Cómo crear una tabla en R?

A continuación se presentan los comandos para crear una tabla de R,en los numerales se explica cada línea del programa.

¿Cómo se ve en R?

Se debe tener en cuenta que en la tablas tanto horizontales como verticales no se ven la totalidad de datos.


¿Cuáles son los comandos?


#SE TIENEN 3 VARIEDADES,SE PRETENDE GENERAR UNA TABLA QUE UNIFIQUE LOS DATOS DE CADA VARIEDAD
#Y POSTERIORMENTE GENERAR EN UNA MISMA PANTALLA LAS GRAFICAS PARA CADA VARIEDAD

criolla=c(5.5,5.8,6,6.4,6.2)#Crear vector con función
criolla=rnorm(50,6,1) #RNORM(CANTIDAD DATOS,MEDIA,VARIANZA),esta función nos sirve para normalizar nuestros datos del vector inicial
pastusa=c(4.3,4.4,4.5,4.1,4.9)
pastusa=rnorm(50,5.5,0.8 )
R12=c(6.4,6.7,6.9,6.8,7)
R12=rnorm(50,6.2,0.4)
TABLAVERTICAL=cbind(criolla,pastusa,R12)#organizar datos en columnas
View(TABLAVERTICAL)
TABLAHORIZONTAL=rbind(criolla,pastusa,R12)#organizar datos en filas
View(TABLAHORIZONTAL)
boxplot(TABLAVERTICAL)# EN UN MISMA PANTALLA PUEDO OBTENER LAS GRÀFICAS DE CADA VARIEDAD
boxplot(TABLAHORIZONTAL)

¿Cómo generar niveles de un factor en R y como cálcular el coeficiente de variación?


Se quiere crear la siguiente tabla:

Donde :

y:Rendimiento (aceite)
factor 1:especies con 2 niveles
factor 2:dosis con 4 niveles
bloqueo:localidad

¿Cómo se ve en R?

 CORREGIR

Tipos de diseños experimentales

Darghan,Universidad Nacional de Colombia

Factorial Simple(FS):De un solo factor

Factorial Completo(FC):Más de un factor,los niveles son iguales o se prueban por igual.No mirar solamente que los números coincidan,sino,si estan anidados o no.

Factorial Incompleto(FI):Más de un factor,los niveles no se prueban por igual.

Aparte de mirar si los niveles se prueban por igual o no ,también,debemos mirar si estan anidados los factores y una manera fácil de entenderlo es con el siguiente cuadro sinóptico:

Completamente al azar(CA):No hay bloqueo,por lo general se da en ambientes controlados.

Bloques al azar(BA):Bloqueo debido a una nueva variable.A su vez los bloques pueden ser completos(En ambos bloques se prueban los mismos factores)/incompletos(lo contrario).

Parcelas Divididas(PD):El bloqueo se da entre los factores

Ejemplos de diseños experimentales

1.Factorial simple en arreglo completamente al azar(Se caracteriza por tener un solo factor)(FSCA ó CA)

y:Rendimiento(materia seca)

factor:Riego

Niveles del factor:riego,goteo,gravedad.

 2.Factorial completo completamente al azar

3 clase de diseño-Más conceptos