Arbol diseño factorial incompleto en bloques al azar

Javier Stiven Bustos Sepulveda

rto=rnorm(18,3,0.3)
lote=gl(2,9,18, c("L1","L2"))
Genotipo=gl(3,3,9, c("G1","G2","G3"))
caudal1=gl(3,1,3, c("C0","C10","C20"))
caudal2=gl(3,1,3, c("C0","C5","C15"))
q=c(caudal1,caudal2,caudal1)

dt=data.frame(lote,Genotipo,Caudal=q,Gravedad_espc=rto)

library(collapsibleTree)
collapsibleTree(dt,hierarchy = c("lote","Genotipo","Caudal","Gravedad_espc"))