En el siguiente experimento se analizaron una serie de datos correspondientes a la germinación de unas semillas plantadas a diferentes niveles de humedad en 8 cajas, las cuales 4 quedaron totalmente cubiertas, mientras que las otras 4 quedaron sin cubrir. Las variables que se analizaron fueron: el porcentaje de germinación, el nivel de humedad y el factor de si se cubrió la caja o no.
require(faraway)
data("seeds")
head(seeds)
Se realizó una exploración univariada de los datos en las que se muestran los siguientes resultados:
require(ggplot2)
seeds$moisture=as.factor(seeds$moisture)
ggplot(seeds, aes(x=moisture,fill=moisture))+geom_bar()+ggtitle("Niveles de humedad")+ylab("Número de Cajas")+xlab("Porcentaje de humedad")

Se logra observar en el gráfico de barras que los diferentes 6 niveles de humedad se realizaron efectivamente en 8 cajas sin ningun tipo de variabilidad.
ggplot(seeds,aes(x=germ))+geom_histogram()+ggtitle("Germinación")+ylab("Número de Cajas")+xlab("Porcentaje de germinación")

Se observa en el histograma una variabilidad muy grande en cuanto al comportamiento de la germinación.
ggplot(seeds,aes(x=covered, fill=covered))+geom_bar()+ggtitle("Cajas cubiertas")+ylab("Cantidad")+xlab("Factor")

Con base a los datos obtenidos se concluye que es un diseño balanceado. Con base a esto, se compáraron las variables del factor de si estaban cubiertas o no y el porcentaje de humedad con la germinación para determinar la eficiencia de cada variable y si determina el porcentaje de germinación.
tapply(seeds$germ, seeds$covered, mean, na.rm=T)
no yes
47.95833 47.65217
ggplot(seeds, aes(y=germ, x=covered, fill=covered))+geom_boxplot()

Como se logra observar en los datos y en el diagrama de cajas, es una muy baja variabilidad en las cajas que estan cubiertas y las que no lo estan, por lo que no influye en el porcentaje de germinación y se recomienda no cubrirlas para evitar el gasto de plástico innecesario.
tapply(seeds$germ, seeds$moisture, mean, na.rm=T)
1 3 5 7 9 11
33.50000 60.62500 71.37500 65.00000 57.57143 0.00000
ggplot(seeds, aes(y=germ, x=moisture, fill=moisture))+geom_boxplot()+ggtitle("Niveles de humedad")+ylab("Germinación")+xlab("Porcentaje de humedad")

Por último se concluye que a diferentes niveles de humedad hay mayor variación en el porcentaje de germinación y el nivel cuyo porcentaje de germinación es mayor y menos variable, por ende, más eficiente es al 5% de humedad y el menos eficiente es al 11% de humedad.
LS0tDQp0aXRsZTogIkluZm9ybWUgZGUgZXhwZXJpbWVudG8gZGUgc2VtaWxsYXMiDQphdXRob3I6ICJKZWFuIEthcmxvIERlbGdhZG8iDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KRW4gZWwgc2lndWllbnRlIGV4cGVyaW1lbnRvIHNlIGFuYWxpemFyb24gdW5hIHNlcmllIGRlIGRhdG9zIGNvcnJlc3BvbmRpZW50ZXMgYSBsYSBnZXJtaW5hY2nDs24gZGUgdW5hcyBzZW1pbGxhcyBwbGFudGFkYXMgYSBkaWZlcmVudGVzIG5pdmVsZXMgZGUgaHVtZWRhZCBlbiA4IGNhamFzLCBsYXMgY3VhbGVzIDQgcXVlZGFyb24gdG90YWxtZW50ZSBjdWJpZXJ0YXMsIG1pZW50cmFzIHF1ZSBsYXMgb3RyYXMgNCBxdWVkYXJvbiBzaW4gY3VicmlyLg0KTGFzIHZhcmlhYmxlcyBxdWUgc2UgYW5hbGl6YXJvbiBmdWVyb246IGVsIHBvcmNlbnRhamUgZGUgZ2VybWluYWNpw7NuLCBlbCBuaXZlbCBkZSBodW1lZGFkIHkgZWwgZmFjdG9yIGRlIHNpIHNlIGN1YnJpw7MgbGEgY2FqYSBvIG5vLg0KDQoNCmBgYHtyIENhcmdhciBsb3MgZGF0b3MgY29ycmVzcG9uZGllbnRlcyBhbCBleHBlcmltZW50bywgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0NCnJlcXVpcmUoZmFyYXdheSkNCmRhdGEoInNlZWRzIikNCmhlYWQoc2VlZHMpDQpgYGANClNlIHJlYWxpesOzIHVuYSBleHBsb3JhY2nDs24gdW5pdmFyaWFkYSBkZSBsb3MgZGF0b3MgZW4gbGFzIHF1ZSBzZSBtdWVzdHJhbiBsb3Mgc2lndWllbnRlcyByZXN1bHRhZG9zOg0KYGBge3IgRXhwbG9yYWNpw7NuIFVuaXZhcmlhZG8gZGUgbG9zIGRhdG9zLCBmaWcuaGVpZ2h0PTMsIGZpZy53aWR0aD0zLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQ0KcmVxdWlyZShnZ3Bsb3QyKQ0Kc2VlZHMkbW9pc3R1cmU9YXMuZmFjdG9yKHNlZWRzJG1vaXN0dXJlKQ0KZ2dwbG90KHNlZWRzLCBhZXMoeD1tb2lzdHVyZSxmaWxsPW1vaXN0dXJlKSkrZ2VvbV9iYXIoKStnZ3RpdGxlKCJOaXZlbGVzIGRlIGh1bWVkYWQiKSt5bGFiKCJOw7ptZXJvIGRlIENhamFzIikreGxhYigiUG9yY2VudGFqZSBkZSBodW1lZGFkIikNCmBgYA0KU2UgbG9ncmEgb2JzZXJ2YXIgZW4gZWwgZ3LDoWZpY28gZGUgYmFycmFzIHF1ZSBsb3MgZGlmZXJlbnRlcyA2IG5pdmVsZXMgZGUgaHVtZWRhZCBzZSByZWFsaXphcm9uIGVmZWN0aXZhbWVudGUgZW4gOCBjYWphcyBzaW4gbmluZ3VuIHRpcG8gZGUgdmFyaWFiaWxpZGFkLg0KYGBge3IgZmlnLmhlaWdodD0zLCBmaWcud2lkdGg9M30NCmdncGxvdChzZWVkcyxhZXMoeD1nZXJtKSkrZ2VvbV9oaXN0b2dyYW0oKStnZ3RpdGxlKCJHZXJtaW5hY2nDs24iKSt5bGFiKCJOw7ptZXJvIGRlIENhamFzIikreGxhYigiUG9yY2VudGFqZSBkZSBnZXJtaW5hY2nDs24iKQ0KYGBgDQpTZSBvYnNlcnZhIGVuIGVsIGhpc3RvZ3JhbWEgdW5hIHZhcmlhYmlsaWRhZCBtdXkgZ3JhbmRlIGVuIGN1YW50byBhbCBjb21wb3J0YW1pZW50byBkZSBsYSBnZXJtaW5hY2nDs24uDQoNCmBgYHtyIGZpZy5oZWlnaHQ9MywgZmlnLndpZHRoPTN9DQpnZ3Bsb3Qoc2VlZHMsYWVzKHg9Y292ZXJlZCwgZmlsbD1jb3ZlcmVkKSkrZ2VvbV9iYXIoKStnZ3RpdGxlKCJDYWphcyBjdWJpZXJ0YXMiKSt5bGFiKCJDYW50aWRhZCIpK3hsYWIoIkZhY3RvciIpDQpgYGANCkNvbiBiYXNlIGEgbG9zIGRhdG9zIG9idGVuaWRvcyBzZSBjb25jbHV5ZSBxdWUgZXMgdW4gZGlzZcOxbyBiYWxhbmNlYWRvLiBDb24gYmFzZSBhIGVzdG8sIHNlIGNvbXDDoXJhcm9uIGxhcyB2YXJpYWJsZXMgZGVsIGZhY3RvciBkZSBzaSBlc3RhYmFuIGN1YmllcnRhcyBvIG5vIHkgZWwgcG9yY2VudGFqZSBkZSBodW1lZGFkIGNvbiBsYSBnZXJtaW5hY2nDs24gcGFyYSBkZXRlcm1pbmFyIGxhIGVmaWNpZW5jaWEgZGUgY2FkYSB2YXJpYWJsZSB5IHNpIGRldGVybWluYSBlbCBwb3JjZW50YWplIGRlIGdlcm1pbmFjacOzbi4NCg0KYGBge3IgZmlnLmhlaWdodD00LCBmaWcud2lkdGg9NH0NCnRhcHBseShzZWVkcyRnZXJtLCBzZWVkcyRjb3ZlcmVkLCBtZWFuLCBuYS5ybT1UKQ0KZ2dwbG90KHNlZWRzLCBhZXMoeT1nZXJtLCB4PWNvdmVyZWQsIGZpbGw9Y292ZXJlZCkpK2dlb21fYm94cGxvdCgpDQoNCmBgYA0KQ29tbyBzZSBsb2dyYSBvYnNlcnZhciBlbiBsb3MgZGF0b3MgeSBlbiBlbCBkaWFncmFtYSBkZSBjYWphcywgZXMgdW5hIG11eSBiYWphIHZhcmlhYmlsaWRhZCBlbiBsYXMgY2FqYXMgcXVlIGVzdGFuIGN1YmllcnRhcyB5IGxhcyBxdWUgbm8gbG8gZXN0YW4sIHBvciBsbyBxdWUgbm8gaW5mbHV5ZSBlbiBlbCBwb3JjZW50YWplIGRlIGdlcm1pbmFjacOzbiB5IHNlIHJlY29taWVuZGEgbm8gY3VicmlybGFzIHBhcmEgZXZpdGFyIGVsIGdhc3RvIGRlIHBsw6FzdGljbyBpbm5lY2VzYXJpby4NCg0KYGBge3J9DQp0YXBwbHkoc2VlZHMkZ2VybSwgc2VlZHMkbW9pc3R1cmUsIG1lYW4sIG5hLnJtPVQpDQpnZ3Bsb3Qoc2VlZHMsIGFlcyh5PWdlcm0sIHg9bW9pc3R1cmUsIGZpbGw9bW9pc3R1cmUpKStnZW9tX2JveHBsb3QoKStnZ3RpdGxlKCJOaXZlbGVzIGRlIGh1bWVkYWQiKSt5bGFiKCJHZXJtaW5hY2nDs24iKSt4bGFiKCJQb3JjZW50YWplIGRlIGh1bWVkYWQiKQ0KYGBgDQpQb3Igw7psdGltbyBzZSBjb25jbHV5ZSBxdWUgYSBkaWZlcmVudGVzIG5pdmVsZXMgZGUgaHVtZWRhZCBoYXkgbWF5b3IgdmFyaWFjacOzbiBlbiBlbCBwb3JjZW50YWplIGRlIGdlcm1pbmFjacOzbiB5IGVsIG5pdmVsIGN1eW8gcG9yY2VudGFqZSBkZSBnZXJtaW5hY2nDs24gZXMgbWF5b3IgeSBtZW5vcyB2YXJpYWJsZSwgcG9yIGVuZGUsIG3DoXMgZWZpY2llbnRlIGVzIGFsIDUlIGRlIGh1bWVkYWQgeSBlbCBtZW5vcyBlZmljaWVudGUgZXMgYWwgMTElIGRlIGh1bWVkYWQuDQo=