Code SAS dans JMP ?J'aimerais savoir s'il existe un moyen d'exécuter un code SAS (ci-dessous) dans JMP. Ou comment je peux émuler le code SAS ci-dessous dans JMP.
proc mixed
Class year week treatment rep;
model stdmisummax=year|week|treatment/ddfm=kr residual;
random rep(year) treatment*rep(year);
repeated week/type=ar(1) subject=treatment*rep(year);
run;
Mes facteurs sont : Traitement, semaine (mesure répétée), année, rep (répétition).
Merci d'avance.
PS: j'ai essayé d'exécuter quelque chose de similaire dans JMP en utilisant le complément de mesures répétées, mais cela ne fonctionnera pas, car je ne peux pas croiser trois facteurs pour créer l'identifiant du sujet.
Je pense que ma version JMP n'est pas pro ou c'est une sorte de version pro pour l'université.
Je n'ai pas de logiciel SAS.
Re : Code SAS dans JMP ?Oui, l'effet aléatoire, Année(Sujet), ou quel que soit le premier, est la variation estimée entre les sujets .
Re : Code SAS dans JMP ?@jiancao Merci pour ton explication.
Chaque fois que j'inclus plus d'un effet aléatoire, j'obtiens "0" pour Estimation, Erreur Std (voir capture d'écran) pour l'un des effets aléatoires. Pourquoi serait-ce? Est-ce à cause de la nidification ?
Merci pour votre aide,
Isadora
Re : Code SAS dans JMP ?Vous pouvez ajouter plusieurs effets aléatoires dans l'analyse des mesures répétées. Cependant, vous recevrez une note disant " La colonne sujet répété a été imbriquée dans tous les effets aléatoires .” Voir la capture d'écran ci-dessous. Je pense que c'est pourquoi Adam a suggéré de laisser tomber Rép(Année) parce qu'il nichera dans Traitement.
Avez-vous essayé d'inclure deux effets aléatoires ? Qu'as-tu vu?
Re : Code SAS dans JMP ?Salut @AdamMorris et @DonMcCormack
J'aimerais savoir si vous pouviez m'aider à comprendre un peu mieux pourquoi vous n'avez utilisé qu'un seul terme comme effet aléatoire. Est-ce parce que je suis autorisé à utiliser un seul terme aléatoire sous Random Effects lorsque j'utilise AR1 ?
Merci,
Isadora
Re : Code SAS dans JMP ?@DonMcCormack C'est très utile !
Merci de l'avoir partagé avec moi et de m'avoir donné votre avis sur la meilleure solution. J'apprécie vraiment ça.
Isadora
Re : Code SAS dans JMP ?Ma solution est ci-dessous:
Fit Model(
Y( :Response ),
Effects(
:Treatment,
:Week,
:Treatment * :Week,
:Year,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Rep[:Treatment, :Year] ),
NoBounds( 1 ),
Personality( "Mixed Model" ),
Subject( :Treatment, :Year, :Rep ),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Run( Variogram( AR1( 1 ) ) )
);
Il n'y a que deux différences. Pour l'effet aléatoire supplémentaire, j'ai imbriqué Rep dans Année. J'ai utilisé Traitement, Année et Rep pour définir implicitement Sujet. Cependant, je trouve l'approche d'Adam consistant à définir explicitement le sujet plus intuitivement attrayante.
Re : Code SAS dans JMP ?@AdamMorris , vous avez été si utile!
Merci beaucoup pour le partage de vos connaissances et toutes ces bonnes informations ! J'apprécie cela.
Re : Code SAS dans JMP ?@DonMcCormack Merci beaucoup pour votre réponse.
Accepteriez-vous également de partager votre solution ? J'aimerais les comparer, car je souhaite en savoir plus sur le codage dans JMP et d'autres approches pour mes données.
Merci!
Isadora
Re : Code SAS dans JMP ?Bonjour Isadora,
Pour partager avec d'autres lecteurs de ce fil, la solution proposée (du support technique de JMP) est de créer une nouvelle colonne Sujet qui est la concaténation des 3 variables identifiant les sujets (Traitement, Année et Rep).
Il est également nécessaire de créer une copie continue de la colonne Semaine à utiliser comme variable Répétée.
JMP Pro ne peut pas correspondre exactement au même modèle que celui spécifié par le code SAS (pour les modèles AR(1), JMP imbrique les effets aléatoires dans la variable sujet).Cependant, JMP peut s'adapter au modèle en omettant le terme Rep(Year).Le script ci-dessous montre comment on peut configurer cela.
Merci,
Adam
New Column("Subject", Character, Formula( :Treatment || ", " || Char( :Year ) || ", " || Char( :Rep ) ));
New Column("Week2", numeric, continuous, formula(:Week));
Fit Model(
Y( :response ),
Effects(
:Treatment,
:Week,
:Year,
:Treatment * :Week,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Treatment*:Year*:Rep ),
Subject(:Subject),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Personality( "Mixed Model" )
);Re : Code SAS dans JMP ?Isadora,
J'ai parlé avec Adam Morris de votre question. Indépendamment, j'ai obtenu les mêmes résultats que lui, même si j'ai configuré le modèle d'une manière légèrement différente. J'aimais un peu mieux sa configuration, alors j'utiliserais celle-là. S'il vous plaît laissez-moi savoir si vous avez des questions supplémentaires.
Enfiler
Re : Code SAS dans JMP ?
Salut @DonMcCormack , j'ai joint une capture d'écran montrant à quoi ressemble le tableau de données pour chaque sujet.
J'ai 4 répétitions, 6 traitements et 2 ans. Des échantillons ont été prélevés à 3 jours différents après l'application des traitements. J'appelle ces jours d'échantillonnage du numéro d'échantillonnage dans la table de données et ils sont appelés "semaine" dans le script SAS.
J'ai JMPPro 13.
Merci d'avance,
Isadora
Re : Code SAS dans JMP ?Il serait utile d'en savoir un peu plus sur la structure du tableau de données. À quoi cela ressemblerait-il pour un sujet ? Les répétitions se produisent-elles dans la semaine ? Vous aurez besoin de JMP Pro pour utiliser la structure de covariance AR(1).
Re : Code SAS dans JMP ?Je ne suis pas un expert des modèles d'effets mixtes en général ni de cette plateforme JMP Pro en particulier. Peut-être que quelqu'un qui est un tel expert se joindra à la discussion. Si ce n'est pas le cas, vous pouvez l'envoyer au support technique de JMP pour plus de suggestions.
Re : Code SAS dans JMP ?@Mark_Bailey Merci pour le partage de cette note. Je l'ai déjà lu.
J'ai la plate-forme mixte dans ma version JMP. Le problème est que je ne peux pas stipuler correctement l'identifiant du sujet dans cette plateforme, car mon sujet est une interaction de traitement*rep(année). Cette plate-forme ne me permet pas de croiser et d'imbriquer pour ajouter le sujet "traitement * représentant (année)" dans la zone ID du sujet. Il vous permet d'ajouter un seul facteur.
Merci
Re : Code SAS dans JMP ?Désolé pour la confusion. Non, il n'est pas possible d'exécuter des programmes SAS dans JMP.
JMP est capable d'effectuer des analyses simples de mesures répétées. Regarde ça note .Les situations plus complexes (c'est-à-dire les structures de covariance non simples) nécessitent la plate-forme Modèle mixte dans JMP Pro.
Re : Code SAS dans JMP ?@Mark_Bailey Merci pour votre réponse. Le problème est que je ne souhaite pas exécuter les analyses dans SAS. Je veux savoir s'il est possible d'exécuter ce code dans JMP d'une manière ou d'une autre ou comment je peux émuler ce code dans un script JMP ou dans Fit Model.
Merci.
Re : Code SAS dans JMP ?Si cela s'applique à votre situation, sachez que JMP est un client complet de SAS. Ainsi, tout comme SAS s'exécute en tant que client, JMP le peut également. Il nécessite une connexion SAS Server. Ce serveur peut fonctionner sur le même ordinateur, comme c'est le cas de mon ordinateur portable professionnel, ou il peut être distant via une connexion réseau. Avez-vous accès à un serveur SAS distant (entreprise, université) ? Si tel est le cas, collez votre code SAS dans un éditeur de script JMP. Faites un clic droit dans l'éditeur et sélectionnez Soumettre SAS.
Sélectionnez Fichier > Préférences et configurez d'abord la connexion au serveur SAS.
Re : Code SAS dans JMP ?@uday_guntupalli Merci pour votre réponse. Oui, je le sais, mais je ne veux pas entrer dans SAS. J'ai reçu ce code d'un collègue, car je n'ai pas pu ou pas pu savoir comment exécuter les mêmes analyses avec des facteurs stipulés de cette façon dans JMP.
Merci 
Re : Code SAS dans JMP ?@bordini ,
Je sais que vous dites que vous n'avez pas SAS, mais saviez-vous que vous pouviez télécharger gratuitement l'édition SAS University ? Si vous avez vraiment besoin de tester le code SAS, vous pouvez télécharger l'édition universitaire SAS et l'essayer.
https://www.sas.com/en_us/software/university-edition/download-software.html
Je ne sais pas comment effectuer l'étape équivalente dans JMP, donc je surveillerai également d'autres réponses
