¿Cómo realizo una prueba T-cuadrado de Hotelling de 2 muestras en JMP?Estoy intentando comparar varias especificaciones o variables de productos diferentes (k = 8) entre 2 grupos diferentes. No puedo encontrar una manera de comparar fácilmente esos múltiples factores entre los 2 grupos en JMP en este momento usando la prueba T cuadrada de Hotelling de 2 muestras. También me gustaría saber si JMP puede crear los CI si rechazamos la estadística de prueba T2 de Hotelling y también revisamos los gráficos del perfil.
Análisis y modelado de datos básicosRe: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?Encontré un código SAS en Internet que realizó el cálculo de 2 muestras de Hotelling T2.Convertí el código a JSL.Lo probé con algunos datos de muestra y salió bien.Es un poco ruidoso pero funciona.Eche un vistazo a la anotación en JSL para ver qué se debe configurar para que funcione con sus datos.
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: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?Gracias por la respuesta. No pensé en intentar encontrar una solución usando SAS.Ahora que lo pienso, es posible que tenga un código SAS para realizar una comparación multivariante entre 2 grupos que también puede funcionar. Dicho esto, esperaba que alguien supiera si: analizar > ajustar menú de modelo > Manova, ¿podría funcionar una comparación multivariada de solo 2 grupos en lugar de varios?
Re: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?De la ayuda de JMP:
Re: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?Gracias Marcos. Vi la opción de prueba T2 de 1 muestra en el menú multivariado, pero no pude encontrar nada para una prueba de 2 muestras, que es lo que estoy buscando para comparar las especificaciones de productos entre 2 grupos. Creo que puedo jugar con la opción MANOVA, pero uso 2 grupos en lugar de los 3+ grupos que normalmente requiere y veré si puedo obtener los resultados correctos.
Re: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?Para aquellos que estén interesados, CREO (alguien puede necesitar verificar) que pueden usar la opción MANOVA para realizar un análisis T2 de Hotelling de 2 muestras de la siguiente manera:
- Analizar > Ajustar modelo > seleccione todas las variables de respuesta que está intentando comparar entre sus 2 grupos y muévalas a "Y", seleccione la columna que contiene las etiquetas para sus 2 grupos, seleccione "Agregar" en 'Construir efectos de modelo' > seleccione " Manova" en el menú desplegable 'Personalidad' > Ejecutar > en el menú 'Elegir respuesta', seleccione "Identidad" > Ejecutar.
- Esto debería generar su valor F y probabilidad para su prueba.
- No estoy seguro de cómo hacer los intervalos simultáneos o ajustados por Bonferroni desde aquí.
Re: ¿Cómo realizo una prueba T cuadrada de Hotelling de 2 muestras en JMP?querido Jim
¿Puedo convencerle de que publique también el código SAS? Eso sería de gran ayuda.
Gracias
Re: ¿Cómo realizo una prueba T-cuadrado de Hotelling de 2 muestras en JMP?A continuación se incluye un enlace a un ejemplo de SAS. Creo que este es el ejemplo que utilicé, pero no estoy seguro.
7.1.15 - Estadística de la prueba T-cuadrada de Hotelling para dos muestras | STAT 505 (psu.edu)
