Modérateur : Groupe des modérateurs
Code : Tout sélectionner
> gls(Y~t, data=df, na.action = (na.omit), method = "REML", correlation=corCAR1(form=~as.factor(td)|id))
Code : Tout sélectionner
Generalized least squares fit by REML
Model: Y ~ t
Data: df
Log-restricted-likelihood: -139
Coefficients:
(Intercept) t
4.509 0.134
Correlation Structure: Continuous AR(1)
Formula: ~as.factor(td) | id
Parameter estimate(s):
Phi
0.725
Degrees of freedom: 950 total; 948 residual
Residual standard error: 0.378
Florian Heinze a écrit :CA Y EST !!! En remplaçant corrAR1 par corrCAR1, on obteint les mêmes résultats que SAS
Retourner vers « Questions en cours »
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité