[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