[R-br] Modelo de efeito misto com correlacao AR1

[Bernardo Rangel Tura]

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
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: jogov2.csv
Tipo: text/csv
Tamanho: 32795 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110504/8c6a2a29/attachment.csv>