set.seed(7234)
rendimiento = c(rnorm(250, 3, 0.5), runif(150, 2.5, 2.8))
genotipo = gl(4, 100, 400, c('Criolla', 'Pastusa', 'Tuquerreña', 'Suprema'))
fertilizacion = gl(2, 50, 400, c('Fertilizacion', 'Control'))
df=data.frame(rendimiento, genotipo, fertilizacion)

Librerias usadas

library(collapsibleTree)
library(ggplot2)

Histogrma

hist(rendimiento)

Tratamientos Factor 1: Cualitativo; 4 niveles Factor 2: Cualitativo; 2 niveles Tratamientos: 8

collapsibleTree(df, hierarchy = c('fertilizacion','genotipo'))

Revision de datos:

summary(df)
##   rendimiento           genotipo         fertilizacion
##  Min.   :0.7697   Criolla   :100   Fertilizacion:200  
##  1st Qu.:2.5738   Pastusa   :100   Control      :200  
##  Median :2.7171   Tuquerreña:100                      
##  Mean   :2.8618   Suprema   :100                      
##  3rd Qu.:3.1187                                       
##  Max.   :4.9150
med_genotipo= tapply(df$rendimiento,df$genotipo,mean)
boxplot(df$rendimiento~df$genotipo)
points(c(1,2,3,4), med_genotipo, col = 'red', pch = '16')

med_fertilizacion= tapply(df$rendimiento,df$fertilizacion,mean)
boxplot(df$rendimiento~df$fertilizacion)
points(c(1,2), med_fertilizacion, col = 'blue', pch = 16)

media1= tapply(df$rendimiento, list(df$genotipo,df$fertilizacion),mean)
media1
##            Fertilizacion  Control
## Criolla         3.057334 2.958199
## Pastusa         2.990311 3.010931
## Tuquerreña      2.957725 2.636749
## Suprema         2.640784 2.642420

los prometiods mas altos de rendimientos los tiene el gonotipo criolla-fertilizada y el genotipo pastusa-control

df$trt= interaction(df$genotipo,df$fertilizacion)
ggplot(df, aes(x= trt, y = rendimiento)) + geom_bar(stat = 'identity', position = 'dodge')

tapply(df$rendimiento, list(df$genotipo,df$fertilizacion), median)
##            Fertilizacion  Control
## Criolla         2.945848 3.021637
## Pastusa         2.948578 3.073258
## Tuquerreña      3.001324 2.634660
## Suprema         2.649275 2.632671
dv= tapply(df$rendimiento, list(df$genotipo,df$fertilizacion), sd)
dv
##            Fertilizacion    Control
## Criolla       0.61006359 0.60201138
## Pastusa       0.47351646 0.56638278
## Tuquerreña    0.56236323 0.08503767
## Suprema       0.08273656 0.08907774
cv= 100*(dv/media1);cv
##            Fertilizacion   Control
## Criolla         19.95410 20.350605
## Pastusa         15.83502 18.810885
## Tuquerreña      19.01337  3.225095
## Suprema          3.13303  3.371067

El mejor tratamiento es el pastusa-control, ya que tiene el segundo mejor pormedio, la mejor mediana, y su coeficiente de variacion es inferior al 20%.