[R-br] Ajuda para realizar post-hoc em Anova de Medidas repetidas.
Augusto Ribas
ribas.aca em gmail.com
Sexta Novembro 25 15:57:48 BRST 2011
Omg. Vida difícil.
Não consigo instalar o pacote 'contrast"
Como foi comentado em outro e-mail, ele necessita do pacote Design, que foi
descontinuado.
Ai eu fui la e baixei a ultima versão no link do R que o Benilton comentou
no outro e-mail
Ai quando vou instalar a partir do zip da o seguinte erro:
##########################################################################
> utils:::menuInstallLocal()
Erro em read.dcf(file.path(pkgname, "DESCRIPTION"), c("Package", "Type")) :
não é possível abrir a conexão
Além disso: Mensagens de aviso perdidas:
1: In unzip(zipname, exdir = dest) :
erro 1 na extração a partir de arquivo zip
2: In read.dcf(file.path(pkgname, "DESCRIPTION"), c("Package", "Type")) :
não foi possível abrir o arquivo comprimido
'Design_2.3-0.tar.gz/DESCRIPTION', motivo provável 'No such file or
directory'
#######################################################################################################
Mas esse arquivo DESCRIPTION ta la dentro do arquivo compactado do pacote,
eu fui la olhar, não sei o que acontece.
Dei um sys.info() pra ver se ajudo.
###
> Sys.info()
sysname release
"Windows" "7"
version nodename
"build 7601, Service Pack 1" "AUGUSTO_NOTE"
machine login
"x86" "Augusto Notebook"
user effective_user
"Augusto Notebook" "Augusto Notebook"
###
To vendo que ta chegando a hora de mudar pro linux...
Em 25 de novembro de 2011 10:56, Walmes Zeviani
<walmeszeviani em gmail.com>escreveu:
> Augusto,
>
> Você deve crias os contrastes de Tukey. Não há função para isso. Você deve
> criar exatamente da maneira como você quer comparar, ou seja, pelas suas
> explicações, comparar níveis de B (pairwise) fixando níveis de A.
> Combinando as funções que o R possui (combn, apply, contrast, multcomp) a
> tarefa fica simples.
>
> da <- expand.grid(A=gl(3,1), B=gl(4,4))
> da$y <- rnorm(nrow(da))
>
> require(nlme)
>
> m0 <- gls(y~A*B, da)
> anova(m0)
>
> # comparar níveis (pairwise) de B fixando nível de A (contrastes de Tukey)
> choose(nlevels(da$B), 2) # número de comparações
> cpr <- combn(nlevels(da$B), 2); cpr # as comparações
>
> Blev <- levels(da$B)
> Alev <- levels(da$A)
> cpr.names <- paste("B", cpr[1,], "-", "B", cpr[2,], sep="")
>
> require(contrast)
>
> cpr.list <- lapply(cpr.list, t)
> cpr.list # matriz de contrastes de Tukey entre B fixando A
>
> c0 <- apply(cpr, 2,
> function(i){
> c1 <- contrast(m0, list(A=Alev[1], B=Blev[i[1]]),
> list(A=Alev[1], B=Blev[i[2]]))
> c1$X
> })
> colnames(c0) <- paste("A", Alev[1], "/", cpr.names, sep="")
> c0
>
> require(multcomp)
> glht(m0, linfct=t(c0)) # repetir para os demais níveis de A
>
> À 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 em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
--
Grato
Augusto C. A. Ribas
Site Pessoal: http://augustoribas.heliohost.org
Lattes: http://lattes.cnpq.br/7355685961127056
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111125/93c01603/attachment-0001.html>
Mais detalhes sobre a lista de discussão R-br