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