Random ou Repeated

[Mag35]

Bonjour,

Je souhaiterai savoir quand les deux options RANDOM et REPEATED sont équivalentes dans un modèle mixte. Je travaille sur des données avec le temps en mesures répétées. Mon modèle est le suivant :

Code :

1
2
3
4
5
6
7
proc mixed DATA=donnees  METHOD=REML;
WHERE time ^=0;
class dog treat TIME;
model ECG=treat | TIME / ddfm=KR;
repeated TIME / type=CS subject=dog;
parms /nobound ;
run;

En remplaçant l'option REPEATED par les options suivantes :

Code :

random intercept / subject=dog;

Code :

random dog;

Je retrouve le même résultat mais je n'arrive pas à l'expliquer. Y a t-il d'autres équivalences entre RANDOM et REPEATED?

De plus l'option "parms /nobound ;" modifie la contrainte de positivité des composantes de la variance. Mais faut-il la spécifier? N'y a-t-il pas de risque d'avoir un modèle non ajusté?

En vous remerciant d'avance,

Magali

[Manoutz]

Je souhaiterai savoir quand les deux options RANDOM et REPEATED sont équivalentes dans un modèle mixte.

Jamais, normalement. je n'ai pas sufisamment d'informations sur ton modele, mais si tu as des bases en statistiques tu dois savoir que ces deux statement ne modelisent pas du tout la meme chose. Chaque statement correspond a la specification de 'un modele particulier.

Jouer avec l'option parms/nobound implique l'utilisation de modeles bien specifiques. Dans quel contexte s'inscrit ta question?

[Mag35]

Merci pour la rapidité de la réponse!

Un modèle à effets sujet aléatoire et un modèle ayant un structure de type CS pour la matrice des résidus donnent les mêmes résultats. Il n'est donc pas possible de parler d'équivalence?

En ce qui concerne l'option params je voulais savoir si cette option était à spécifier dans le cas d'un problème de converge du modèle ou s'il y a d'autres raisons pour l'utiliser? je travaille sur un modèle simple parametre = traitement + temps + traitement * temps avec le temps en mesures répétées.

[Manoutz]

J'ai l'impression que tu te compliques un peu la tache:

je travaille sur un modèle simple parametre = traitement + temps + traitement * temps avec le temps en mesures répétées.

Alors pourquoi inclues tu un statement random?

En ce qui concerne l'option params je voulais savoir si cette option était à spécifier dans le cas d'un problème de converge du modèle ou s'il y a d'autres raisons pour l'utiliser?

cette option permet de relaxer des contraintes, donc c'est UNE alternative en cas de non convergence. Il y a peut etre d'autres raisons, a voir suivant le contexte et la problematique. Mais je dirais que la convergence doit regrouper la majorite des situations de son utilisation.

[Mag35]

Merci encore pour la rapidité.

Pour l'option RANDOM c'est une question que l'on se pose pour savoir pourquoi nous obtenons les mêmes résultats qu'avec REPEATED avec CS comme structure mais comme il s'agit de mesures répétées cette dernière option est la plus adéquate.

En ce qui concerne parms je vais continuer à regarder de plus près.

Merci encore