Se a única variável manipulada ("independente") foram os tratamentos, e a resposta é unicamente uma única variável também, matematicamente você não acaba com a mesma análise, apenas por caminhos diferentes?
On Wed, Nov 4, 2020 at 12:28 PM Fernando Souza por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Boa tarde pessoal, > > Estou analisando um experimento em Switchback com vacas leiteiras e estou > com duvidas sobre a montagem do modelo. Gostaria da ajuda de vocês. > O experimento comparou dois tratamentos A e B. Utilizou 30 vacas, porém as > vacas não foram avaliadas simultaneamente. 10 vacas entraram no momento 1, > 10 entraram no momento 2 e 10 entraram no momento 3. As vacas dentro de > cada momento foram sorteadas para uma das duas sequencias possíveis A-B-A > ou B-A-B, sendo a sequencia aplicadas em 3 períodos diferentes de igual > duração. > Entendo que o desenho experimental desse trabalho é semelhante ao descrito > nesse trabalho: > https://www.sciencedirect.com/science/article/pii/S0022030287802734 > <https://link.getmailspring.com/link/6aac746e-76fd-40c3-9bb7-28232aac3...@getmailspring.com/0?redirect=https%3A%2F%2Fwww.sciencedirect.com%2Fscience%2Farticle%2Fpii%2FS0022030287802734&recipient=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy> > sendo o block os momentos (1,2 e 3) > > Porém nesse trabalho os autores utilizam o PROC GLM dos SAS. Eu no entanto > desejo analisá-lo através da abordagem dos modelos mistos, com a função > lme{nlme} > Fiz o seguinte modelo e gostaria de uma opinião para saber se o modelo > está correto? Se não como posso melhorá-lo > > dados<-structure(list(Bloco = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > 3L, 3L, 3L), .Label = c("1", "2", "3"), class = "factor"), Ordem = > structure(c(1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = > "factor"), > Periodo = c(1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, > 1L, 2L, 3L), Tratamento = structure(c(1L, 2L, 1L, 1L, 2L, > 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, > 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, > 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, > 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, > 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, > 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("A", > "B"), class = "factor"), Sequencia = structure(c(1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L), .Label = c("1", > "2"), class = "factor"), Vacas = structure(c(1L, 1L, 1L, > 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, > 7L, 7L, 7L, 8L, 8L, 8L, 9L, 9L, 9L, 10L, 10L, 10L, 11L, 11L, > 11L, 12L, 12L, 12L, 13L, 13L, 13L, 14L, 14L, 14L, 15L, 15L, > 15L, 16L, 16L, 16L, 17L, 17L, 17L, 18L, 18L, 18L, 19L, 19L, > 19L, 20L, 20L, 20L, 21L, 21L, 21L, 22L, 22L, 22L, 23L, 23L, > 23L, 24L, 24L, 24L, 25L, 25L, 25L, 26L, 26L, 26L, 27L, 27L, > 27L, 28L, 28L, 28L, 29L, 29L, 29L, 30L, 30L, 30L), .Label = c("1", > "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", > "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", > "23", "24", "25", "26", "27", "28", "29", "30"), class = "factor"), > Leite = c(18.28, 15.73, 14.96, 16.08, 15.64, 15.41, 16.21, > 15.52, 14.93, 17.29, 14.52, 12.98, 13.28, 12.57, 12.31, 21.82, > 19.28, 18.75, 16.69, 15.15, 14, 17.11, 15.51, 14.73, 15.38, > 14.77, 13.37, 15.21, 14.46, 13.5, 12.67, 11.91, 11.55, 15.78, > 14.68, 14.61, 13.29, 12.12, 11.77, 17.32, 14.81, 14.66, 16.13, > 14.81, 14.29, 14.94, 13.45, 13.09, 15.27, 12.36, 11.29, 14.02, > 13, 12.06, 15.25, 14.9, 14.59, 12.44, 11.8, 11.56, 12.56, > 12.22, 7.68, 15.15, 14.56, 10.08, 14.26, 13.34, 10.21, 13.56, > 12.85, 8.9, 14.09, 14.63, 11.52, 17.7, 17.24, 14.89, 16.69, > 16.28, 12.25, 12.66, 12.14, 9, 18.74, 18.23, 14.98, 17.62, > 17.4, 15.67), Per = structure(c(1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, > 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L), .Label = c("1", "2", "3" > ), class = "factor")), class = "data.frame", row.names = c(NA, > -90L)) > > #install.packages("nlme") > library(nlme) > options(contrasts=c("contr.sum","contr.poly")) > controle<-lmeControl(maxIter = 200, msMaxIter = 200,opt = 'optim') > mlme <- lme(Leite ~Tratamento,random=~1|Bloco/Vacas/Per,control=controle, > data=dados) > anova(mlme,type='marginal') > > _______________________________________________ > 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. >
_______________________________________________ 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.