[R-br] problemas com comparações múltiplas

Maurício Lordêlo mslordelo em gmail.com
Segunda Março 24 21:57:37 BRT 2014


Saudações a todos!
Estou realizando inferência simultânea para um conjunto de dados
correlacionados onde desejo comparar o efeito de quatro tratamentos. Para
isto estou utilizando a função "simplesimint" do pacote "BSagri".
Consigo fazer utilizando Dunnet porém quando mudo para Tukey ocorre um
problema.
Abaixo tenho um CMR onde tenho o vetor de coeficientes estimados (coefs), a
matriz de variância e covariâncias (VarCov), o número de obervações em cada
grupo (n) e a matriz de contrastes (cmat).
O problema principal surge quando  faço utilizando todas as linhas (são
seis) da matriz de contrastes (cmat2). Já para duas, três e até quatro
linhas da matriz de contrastes não há problemas.
Grato caso alguém possa ajudar.
Abraço


library(BSagri)
coefs=c("TratamentoA"=-1.0618,"TratamentoB"=1.1470,"TratamentoC"=0.3068,"TratamentoC"=0.2501)
VarCov=matrix(data = NA, nrow = 4, ncol = 4, byrow = FALSE,dimnames=NULL)
VarCov[1,]=c(0.3133311,-0.1618691,-0.1736806,-0.1712345)
VarCov[2,]=c(-0.1618697,0.4510964,0.1595845,0.1622442)
VarCov[3,]=c(-0.1736809,0.1595843,0.2729564,0.1617069)
VarCov[4,]=c(-0.1712371,0.1622446,0.1617080,0.2757233)
VarCov
n=c("TratamentoA"=40,"TratamentoB"=20,"TratamentoC"=60,"TratamentoD"=60)

cmat1<-contrMat(n=n, type="Dunnett")
cmat1
test1<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat1,
alternative="two.sided")
plotCI(test1)

cmat2<-contrMat(n=n, type="Tukey")
cmat2
test2<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat2,
alternative="two.sided")
##ERRO!!
#Erro em uniroot(pfct, interval = interval) :
#  f() values at end points not of opposite sign



test3<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:2,1:4],
alternative="two.sided")
summary(test3)
plotCI(test3)
test4<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:3,1:4],
alternative="two.sided")
summary(test4)
plotCI(test4)
test5<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:4,1:4],
alternative="two.sided")
summary(test5)
plotCI(test5)
test6<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:5,1:4],
alternative="two.sided")
##ERRO!!
#Erro em uniroot(pfct, interval = interval) :
#  f() values at end points not of opposite sign
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140324/8c9863bb/attachment.html>


Mais detalhes sobre a lista de discussão R-br