<div dir="ltr"><div><br></div>Saudações a todos!<div>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".</div>
<div>Consigo fazer utilizando Dunnet porém quando mudo para Tukey ocorre um problema.</div><div>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).</div>
<div>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.</div><div>Grato caso alguém possa ajudar.</div>
<div>Abraço</div><div><br></div><div><br><div><div>library(BSagri)</div><div>coefs=c("TratamentoA"=-1.0618,"TratamentoB"=1.1470,"TratamentoC"=0.3068,"TratamentoC"=0.2501)</div><div>
VarCov=matrix(data = NA, nrow = 4, ncol = 4, byrow = FALSE,dimnames=NULL)</div><div>VarCov[1,]=c(0.3133311,-0.1618691,-0.1736806,-0.1712345)</div><div>VarCov[2,]=c(-0.1618697,0.4510964,0.1595845,0.1622442)</div><div>VarCov[3,]=c(-0.1736809,0.1595843,0.2729564,0.1617069)</div>
<div>VarCov[4,]=c(-0.1712371,0.1622446,0.1617080,0.2757233)</div><div>VarCov</div><div>n=c("TratamentoA"=40,"TratamentoB"=20,"TratamentoC"=60,"TratamentoD"=60)</div><div><br></div><div>
cmat1<-contrMat(n=n, type="Dunnett")</div><div>cmat1</div><div>test1<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat1, alternative="two.sided")</div><div>plotCI(test1)</div><div><br></div>
<div>cmat2<-contrMat(n=n, type="Tukey")</div><div>cmat2</div><div>test2<-simplesimint(coef=coefs, vcov=VarCov, df=NULL, cmat=cmat2, alternative="two.sided")</div><div>##ERRO!!</div><div>#Erro em uniroot(pfct, interval = interval) : </div>
<div>#  f() values at end points not of opposite sign</div><div><br></div><div><br></div></div></div><div><br></div><div><div>test3<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:2,1:4], alternative="two.sided")</div>
<div>summary(test3)</div><div>plotCI(test3)</div><div>test4<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:3,1:4], alternative="two.sided")</div><div>summary(test4)</div><div>plotCI(test4)</div><div>
test5<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:4,1:4], alternative="two.sided")</div><div>summary(test5)</div><div>plotCI(test5)</div><div>test6<-simplesimint(coef=coefs, vcov=VarCov, df=114, cmat=cmat2[1:5,1:4], alternative="two.sided")</div>
<div>##ERRO!!</div><div>#Erro em uniroot(pfct, interval = interval) : </div><div>#  f() values at end points not of opposite sign</div></div><div><br></div></div>