Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Ich versuche, mehrere verschiedene Produktspezifikationen oder Variablen zu vergleichen (k = 8). zwischen 2 verschiedenen Gruppen. Ich kann derzeit in JMP keine Möglichkeit finden, diese mehreren Faktoren zwischen den beiden Gruppen einfach zu vergleichen, indem ich den Hotelling-T-Quadrat-Test für zwei Stichproben verwende. Ich würde auch gerne wissen, ob JMP dann die CIs erstellen kann, wenn wir die T2-Teststatistik von Hotelling ablehnen und auch die Profildiagramme überprüfen?
Grundlegende Datenanalyse und ModellierungBetreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Ich habe im Internet einen SAS-Code gefunden, der die 2 Sample Hotelling T2-Berechnung durchgeführt hat.Ich habe den Code in JSL konvertiert.Ich habe es anhand einiger Beispieldaten getestet und es war in Ordnung.Es klappert zwar etwas, aber es funktioniert.Sehen Sie sich die Anmerkung in der JSL an, um zu sehen, was eingestellt werden muss, damit es mit Ihren Daten funktioniert.
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>
Betreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Danke für die Antwort. Ich habe nicht daran gedacht, mit SAS eine Lösung zu finden.Wenn ich jetzt darüber nachdenke, habe ich möglicherweise SAS-Code für einen multivariaten Vergleich zwischen zwei Gruppen, der möglicherweise auch funktioniert. Trotzdem hatte ich gehofft, dass jemand wüsste, ob mit „analysieren“ > „Modell anpassen“ > „Manova“ ein multivariater Vergleich von nur zwei Gruppen anstelle mehrerer Gruppen möglich ist?
Betreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Aus der JMP-Hilfe:
Betreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Danke Mark. Ich habe die Option für den T2-Test mit einer Stichprobe im multivariaten Menü gesehen, konnte aber nichts für einen Test mit zwei Stichproben finden – was ich suche, um Produktspezifikationen zwischen zwei Gruppen zu vergleichen. Ich glaube, ich probiere vielleicht mit der MANOVA-Option herum, verwende aber 2 Gruppen statt der 3+ Gruppen, die normalerweise erforderlich sind, und schaue, ob ich die richtigen Ausgaben bekomme.
Betreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Für diejenigen, die Interesse haben – ich denke (jemand muss dies möglicherweise überprüfen) – Sie können die MANOVA-Option verwenden, um eine Hotelling-T2-Analyse für zwei Stichproben wie folgt durchzuführen:
- Analysieren > Modell anpassen > Wählen Sie alle Ihre Antwortvariablen aus, die Sie zwischen Ihren beiden Gruppen vergleichen möchten, und verschieben Sie sie nach „Y“, wählen Sie die Spalte mit Ihren Beschriftungen für Ihre beiden Gruppen aus, wählen Sie „Hinzufügen“ unter „Modelleffekte konstruieren“ > wählen Sie „ Manova“ aus dem Dropdown-Menü „Persönlichkeit“ > Ausführen > im Menü „Antwort wählen“ wählen Sie „Identität“ > Ausführen.
- Dies sollte Ihren F-Wert und Ihre Wahrscheinlichkeit für Ihren Test ausgeben.
- Ich bin mir nicht sicher, wie ich die simultanen oder Bonferroni-angepassten Intervalle von hier aus durchführen soll.
Betreff: Wie führe ich in JMP einen Hotelling-T-Quadrat-Test mit zwei Stichproben durch?Lieber Jim
Könnte ich Sie überreden, auch den SAS-Code zu veröffentlichen? Das wäre am hilfreichsten.
Danke
Betreff: Wie führe ich in JMP einen Hotelling T-Quadrat-Test mit zwei Stichproben durch?Hier ist ein Link zu einem SAS-Beispiel. Ich glaube, das ist das Beispiel, das ich verwendet habe, aber ich bin nicht sicher.
7.1.15 – Die Hotelling-T-Quadrat-Teststatistik für zwei Stichproben | STAT 505 (psu.edu)
