Caros colegas, Estou precisando de um ajuda para um modelo de efeito misto com correlação AR1 A auto correlação se justifica pois os pacientes devem apreender e obter melhores resultado de um bloco para outro então usei os seguinte comando ------------------------------------------------------- base <- read.csv('jogov2.csv') require(nlme) modelo <- lme(valor~Gr+Idd+Escol,random=~1|Bloco,correlation=corAR1(, form=~1|Bloco) ,data=base) summary(modelo) ----------------------------------------------------- Que resultou em ----------------------------------------------------- Linear mixed-effects model fit by REML Data: base AIC BIC logLik 3928.017 3958.571 -1957.009 Random effects: Formula: ~1 | Bloco (Intercept) Residual StdDev: 2.480127 6.846336 Correlation Structure: AR(1) Formula: ~1 | Bloco Parameter estimate(s): Phi 0.1407996 Fixed effects: valor ~ Gr + Idd + Escol Value Std.Error DF t-value p-value (Intercept) -1.7057456 1.8249072 577 -0.934703 0.3503 GrB 1.7242603 0.6606029 577 2.610131 0.0093 Idd -0.0479666 0.0178547 577 -2.686502 0.0074 Escol 0.3441736 0.0926290 577 3.715615 0.0002 Correlation: (Intr) GrB Idd GrB -0.068 Idd -0.415 -0.101 Escol -0.602 -0.227 0.001 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -3.16689925 -0.56838036 -0.09641692 0.55523377 3.66357685 Number of Observations: 585 Number of Groups: 5 ------------------------------------------------------- Agora vem o problema eu preciso saber se o efeito aleatório tem um média diferente diferente entre o grupo B e o grupo A e além disso se o valor de Phi é diferente para os dois grupos Alguém pode me ajudar ? Ps. Base em anexo -- []s Tura