Comment effectuer un test du carré en T de Hotelling sur 2 échantillons dans JMP ?J'essaie de comparer plusieurs spécifications ou variables de produits différentes (k = 8) entre 2 groupes différents. Je ne trouve pas actuellement de moyen de comparer facilement ces multiples facteurs entre les 2 groupes dans JMP en utilisant le test du carré en T de Hotelling à 2 échantillons. J'aimerais également savoir si JMP peut ensuite créer les CI si nous rejetons la statistique du test T2 de Hotelling et examinons également les tracés de profil ?
Analyse et modélisation des données de baseRe : Comment effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?J'ai trouvé sur Internet du code SAS qui effectuait le calcul de 2 échantillons Hotelling T2.J'ai converti le code en JSL.Je l'ai testé avec quelques exemples de données et tout s'est bien passé.C'est un peu chiant mais ça marche.Jetez un œil à l'annotation dans le JSL pour voir ce qui doit être configuré pour fonctionner avec vos données.
Names Default To Here( 1 );
// Open the data table
dt = open("$SAMPLE_DATA/Prostate Cancer.jmp");
groupVar = "Status"; // Enter the grouping column name
groupList = {"CCD", "Nor"}; // Set the values for the 2 groups
startColumn = 660; // Set the starting column number for the data to be included
// Create the group 1 and group 2 matrices for analysis
gp1Rows = dt << get rows where( as column( groupVar) == groupList[1] );
group1 = dt[gp1Rows, startColumn :: N Cols( dt )];
gp1Rows = dt << get rows where( as column( groupVar ) == groupList[2] );
group2 = dt[gp1Rows, startColumn :: N Cols( dt )];
// Get the dimensions of the matricies
nRows1 = N Row( group1 );
nRows2 = N Row( group2 );
numCols = N Col( group1 );
// Create a filled vector for the analysis
j1 = J( nRows1, 1, 1 );
j2 = J( nRows2, 1, 1 );
// Create an identity matrix for each group
identity1 = Identity( nRows1 );
identity2 = Identity( nRows2 );
// Calculate the means for each of the columns by each group
MeanY1 = group1` * j1 / nRows1;
MeanY2 = group2` * j2 / nRows2;
// Create an output window to display the results in
nw = New Window("Hotelling T2", outlinebox("Hotelling 2 Sample T2 Analysis",spacerbox(size(1,15)),lub = vlistbox()));
// Add the means table to the output
obMean = Outline Box("Group Means",
mTB = Table Box(
string Col box("Groups", groupList),
number col box("# Rows", matrix(nRows1)||matrix(nRows2))
)
);
for(i=startColumn, i<=NCols(dt),i++,
mTB << append(number col box(column(dt,i)<<GET name=""></GET>
Re : Comment effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?Merci pour votre réponse. Je n'ai pas pensé à essayer de trouver une solution en utilisant SAS.Maintenant que j'y pense, j'ai peut-être du code SAS pour faire une comparaison multivariée entre 2 groupes qui peuvent également fonctionner. Cela dit, j'espérais que quelqu'un savait si : analyser > ajuster le menu du modèle > Manova , pouvait effectuer une comparaison multivariée de seulement 2 groupes au lieu de plusieurs ?
Re : Comment puis-je effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?Depuis l'aide JMP :
Re : Comment puis-je effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?Merci Marc. J'ai vu cette option de test T2 à 1 échantillon dans le menu multivarié, mais je n'ai rien trouvé pour un test à 2 échantillons - c'est ce que je recherche pour comparer les spécifications de produits entre 2 groupes. Je pense que je peux jouer avec l'option MANOVA, mais utiliser 2 groupes au lieu des 3+ groupes qu'elle appelle habituellement et voir si je peux obtenir les sorties correctes.
Re : Comment effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?Pour ceux qui sont intéressés - JE PENSE (quelqu'un devra peut-être vérifier) que vous pouvez utiliser l'option MANOVA pour effectuer une analyse T2 de Hotelling à 2 échantillons comme suit :
- Analyser > Ajuster le modèle > sélectionnez toutes vos variables de réponse que vous essayez de comparer entre vos 2 groupes et déplacez-les vers "Y", sélectionnez la colonne contenant vos étiquettes pour vos 2 groupes, sélectionnez "Ajouter" sous "Construire les effets du modèle" > sélectionnez " Manova" dans le menu déroulant "Personnalité" > Exécuter > sous le menu " Choisir la réponse ", sélectionnez " Identité " > Exécuter.
- Cela devrait afficher votre valeur F et votre probabilité pour votre test.
- Je ne sais pas comment faire les intervalles simultanés ou ajustés par Bonferroni à partir d'ici.
Re : Comment effectuer un test du carré en T de Hotelling à 2 échantillons dans JMP ?cher Jim
Puis-je vous persuader de publier également le code SAS ? Ce serait très utile.
Merci
Ré: Comment effectuer un test T-carré de Hotelling à 2 échantillons dans JMP ?Voici un lien vers un exemple SAS. Je crois que c'est l'exemple que j'ai utilisé, mais je n'en suis pas sûr.
7.1.15 - Statistique du test en T de Hotelling à deux échantillons | STAT 505 (psu.edu)
