63200_fr-FR

liliconcom

Option Enregistrer les colonnes avec un modèle mixte

Bonjour tous le monde,

 

Lorsque nous exécutons un modèle mixte à l'aide de la plateforme Fit Model, nous pouvons enregistrer les valeurs résiduelles et prédites, mais nous avons également la possibilité d'enregistrer les valeurs résiduelles conditionnelles et prédites conditionnelles (qui prennent en compte les effets fixes et aléatoires). Cependant, nous n'avons pas les mêmes options pour les résidus étudiés et les colonnes StdErr Pred Formula.

Quel est le meilleur moyen d’obtenir ou de recalculer la "version conditionnelle" de ces deux colonnes?

 

Notez que je ne travaille pas avec JMP Pro.

 

Merci d'avance pour votre aide!

 

 

Option de sauvegarde des colonnes avec modèle mixte

Oui c'est vrai je vais faire ça!

Option de sauvegarde des colonnes avec modèle mixte

Vous pouvez suggérer cette fonctionnalité dans la liste de souhaits de la communauté JMP.

Option de sauvegarde des colonnes avec modèle mixte

Bonjour GM,

Oui, vous avez raison, nous pouvons extraire les erreurs types d’effets aléatoires de ce type, mais ce que je recherche, c’est

les erreurs types des résidus et des valeurs prédites. Je ne pense pas que ce soit la même information ici.

Option de sauvegarde des colonnes avec modèle mixte

Merci Cameron pour cette réponse détaillée, je vais la tester!

Mais en effet, ce n’est pas un moyen facile d’obtenir ce que nous voulons. Je me demande pourquoi il n'est pas directement disponible dans le menu Enregistrer les colonnes.

Option de sauvegarde des colonnes avec modèle mixte

Si vous estimez les effets aléatoires à l'aide de la plateforme Fit Model avec Personality = Standard Least Square et Method = REML, vous pouvez extraire les erreurs standard des effets aléatoires en créant un tableau de données à partir du tableau de prédiction d'effets aléatoires. Les erreurs standard sont incluses.   La personnalité de modèle mixte ne génère pas les erreurs standard avec les coefficients Aléatoire du rapport.

 

Option de sauvegarde des colonnes avec modèle mixte

Tu as raison.  Duh.  Je n'ai pas accès à JMP Pro, je ne peux donc pas le vérifier, mais je ne pense pas que vous pourrez obtenir tous les composants dont vous auriez besoin pour calculer les résidus conditionnels normalisés directement à partir du rapport.  La colonne PredSE donne la matrice X comme second argument de Vec Quadratic (), et le rapport vous donne le MSE.  Il faudrait construire G et Z. 

 

Vous pouvez obtenir Z de quelque chose comme:

Design Nom( :RandomEffectCol << Get Values, {list of levels})

 Si vous avez plusieurs effets aléatoires, vous pouvez simplement les concaténer.

Pour obtenir G, cela dépend de la structure de corrélation choisie.  Si vous avez une structure de corrélation simple avec tous les 0 en diagonale, vous pouvez construire G en multipliant votre vecteur de composantes de variance avec une matrice d'identité.

Une fois que vous avez X, G, Z et le MSE, vous pouvez alors calculer la variance de la matrice d'erreur de la manière suivante:

n = N Row(dt);

V = Z*G*Z` + MSE*Identity(n);
K = Identity(n) - Z*G*Z`*Inv(V);
Q = X*Inv(X`*Inv(V)*X)*X`;

Var_e = K*(V-Q)*K`;

A partir de là, prenez le carré des éléments diagonaux de Var_e pour obtenir les erreurs standard conditionnelles de vos résidus conditionnels.  Je ne pense pas que JMP ait une fonction pour retourner uniquement les éléments diagonaux de la matrice, mais en voici une simple qui fera l'affaire:

getDiag = function({X},
 n = ncol(X);
 diags = {};
 for(i=1, i<=N Col(X),i++,
  Insert Into(diags,X[i,i])
 )
);

C'est beaucoup d'effort! Désolé, j'aurais aimé connaître un moyen plus simple.

Option de sauvegarde des colonnes avec modèle mixte

Merci pour votre réponse Cameron!

 

Cependant, si vous avez sauvegardé le Std Err de résidus, je pense que cela est lié à ceux qui n'incluent pas d'effet aléatoire et non aux résidus conditionnels. Il ne semble pas possible de sauvegarder le Conditional Std Err des résidus. Je pense que nous pouvons refaire le calcul basé sur les valeurs de Hats mais le problème est le même: elles sont probablement liées aux "résidus résiduels".

Option de sauvegarde des colonnes avec modèle mixte

Vous pouvez étudier les résidus conditionnels en les divisant par leurs erreurs-types.  Enregistrez les deux colonnes dans la table, puis calculez les valeurs studentisées avec une formule de colonne.  Je ne sais pas tout de suite comment vous aider avec la formule Std Error Pred.