Criando um fatorial completo com dados desbalanceados
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@yahoo.com.br alexandre.santos@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 ======================================================================
Não sei se entendi bem, mas para gerar fatores com níveis com número diferente de ocorrência você pode começar usando rep(), se o número de repetições é fixo/conhecido, ou sample() se o número de repetições é aleatório, enfim. da <- data.frame(F1 = rep(LETTERS[1:4], c(4, 6, 10, 20)), F2 = sample(1:3, size = 40, replace = TRUE)) str(da) xtabs(~F1 + F2, data = da) À disposição. Walmes.
Muito obrigado Walmes!! Era exatamente a dica que eu queria, da <- data.frame(F1 = rep(LETTERS[1:4], c(4, 6, 10, 20)), F2 = sample(1:3, size = 40, replace = TRUE)) str(da) length(da[,1]) xtabs(~F1 + F2, data = da) -- ====================================================================== 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@yahoo.com.br alexandre.santos@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 ====================================================================== Em 29/03/2016 08:57, Walmes Zeviani escreveu:
Não sei se entendi bem, mas para gerar fatores com níveis com número diferente de ocorrência você pode começar usando rep(), se o número de repetições é fixo/conhecido, ou sample() se o número de repetições é aleatório, enfim.
da <- data.frame(F1 = rep(LETTERS[1:4], c(4, 6, 10, 20)), F2 = sample(1:3, size = 40, replace = TRUE)) str(da)
xtabs(~F1 + F2, data = da)
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
participantes (2)
-
ASANTOS -
Walmes Zeviani