Teste glht:multcomp com lista de contrastes
Estou tentando usar a função glht com uma lista de contrastes. Tem solução? require(multcomp) x=structure(list(t = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L ), .Label = c("t1", "t2", "t3", "t4"), class = "factor"), r = c(5L, 6L, 23L, 29L, 12L, 16L, 45L, 49L)), .Names = c("t", "r"), class = "data.frame", row.names = c(NA, -8L)) m=lm(r~t, data=x) teste3=glht (m, linfct = mcp(t = "Tukey" )) mc=teste3$linfct y=row.names(mc) dd=list(y[c(1,2)],y[3]) c=lapply(dd, function(ddi) { mc[ddi,] }) # ?????lapply(c, function(ci) {summary(glht(m, linfct =ci))})
Você tem que garantir que tá passando matrizes para linfct=, a mensagem de erro é um elemento da lista era vetor e não matriz. Isso pode resolver. lapply(c, class) c2 <- lapply(c, matrix, ncol=length(coef(m))) lapply(c2, class) lapply(c2, function(ci) { summary(glht(m, linfct =ci)) }) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Emmanuel Arnhold -
Walmes Zeviani