[R-br] Criando um fatorial completo com dados desbalanceados

[ASANTOS]

Prezados Listeiros,

     Gostaria de criar  um fatorial completo entre 4 tratamentos com 
dados desbalanceados usando a função gen.factorial() com o pacote 
AlgDesign. Alguém poderia me ajudar, segue CRM abaixo:

#Dados artificiais
tratamentos_b <- gl(4, 15, labels = paste("t", 1:3, sep=""))## Dados 
balanceados
tratamentos_nb <- gl(1, 30, labels = paste("t", 4, sep="")) ### Dados 
desbalanceados
tratamentos<-as.character(tratamentos_b,tratamentos_nb)

sp  <- cbind(c(rnorm(10,  5, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 
2.5, 0.25),
              c(rnorm(10, 12, 0.25), rnorm(50, 2.5, 0.25)), rnorm(60, 
2.5, 0.25))
colnames(sp) <- c("sp1", "sp2", "sp3", "sp4")

dados<-data.frame(sp,tratamentos)

## Criando o fatorial completo
require(AlgDesign)
contr2df2 <-AlgDesign::gen.factorial(dados$tratamentos, TRUE, 
varNames=c("T1_2", "T1_3", "T1_4", "T2_3", "T2_4", "T3_4"))
contr2df2

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
OrcID: orcid.org/0000-0001-8232-6722
Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635
======================================================================