[R-br] função glht{multcomp} comparaçãos função lme
Fernando Souza
nandodesouza em gmail.com
Segunda Janeiro 30 14:05:58 BRST 2017
Caros colegas,
Estou realizando um exercício de parcela subdividida utilizando a
função lme{nlme}. Acontece que quando rodo a função summary.glht as
comparações pairwise são feitas e eu recebo as seguintes mensagens
warnings:
Warning messages:
1: In RET$pfunction("adjusted", ...) : Completion with error > abseps
2: In RET$pfunction("adjusted", ...) : Completion with error > abseps
3: In RET$pfunction("adjusted", ...) : Completion with error > abseps
4: In RET$pfunction("adjusted", ...) : Completion with error > abseps
Essas mensagens no entanto não aparecem quando defino outro teste na
função summary.glht
summary.glht(compara, test=adjusted("fdr"))
Porque isso ocorre?
Segue o CMR.
att
#install.packages(c("nlme","multcomp"))
library(nlme)
library(multcomp)
PARCSUB<-read.csv("https://www.dropbox.com/s/y8qmhwau0i8euva/ExemploSub.csv?raw=1",head=TRUE)
PARCSUB<-transform(PARCSUB,PARCELA=factor(PARCELA),BLOCO=factor(BLOCO),PASTO=factor(PASTO),MINERAL=factor(MINERAL))
modeloPARSUBlme<-lme(PRODLEIT~BLOCO+PASTO*MINERAL,random=~1|PARCELA,data=PARCSUB)
compara<- glht(modeloPARSUBlme, linfct=mcp(PASTO="Tukey"))
summary.glht(compara)
==================================
Fernando Souza
Zootecnista, DSc. Produção e Alimentação Animal
celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim)
e-mail:nandodesouza em gmail.com
Lattes: http://lattes.cnpq.br/6519538815038307
blog: https://producaoanimalcomr.wordpress.com/
=================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170130/b7a01d21/attachment.html>
Mais detalhes sobre a lista de discussão R-br