<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Hola ,
Antes, describo mi problema, debo mencionar que trabajo con JMP7.
Tengo una tabla con 12 columnas, las primeras 6 columnas son parámetros descriptivos para cada una de las 6 columnas siguientes. Las primeras 6 columnas son todas de caracteres y las últimas 6 columnas son todas numéricas, binomiales.
Necesito ejecutar una regresión logística para cada una de las últimas 6 columnas y usar las primeras 6 columnas como predictores para extraer el valor "Prob> ChiSq" de la prueba de "Pearson".
Ahora, hice el siguiente script para una de las columnas "Y, Respuesta":
dt = Current Data Table();
column_names = dt << Get Column Names(Character, String);
results = New Table("Results", Invisible,
New Column("Parameter", character),
New Column("GEN", numeric),
Add Rows(NItems(column_names))
);
for(i=1,i<=NItems(column_names),i++,
colName = column_names;
cmd = Expr(
ow = dt << Contingency(
Y( :GEN),
X( :cName ),
Contingency Table( 0 ),
Mosaic Plot( 0 )
);
Wait(0.1);
owr = ow<
pValueBox= owr[NumberColBox(6)];
pValue = pValueBox[2];
Column(results,1) = colName;
Column(results,2) = pValue;
ow << Close Window;
);
Substitute Into(cmd,Expr(cName),colName);
Eval(cmd);
);
results << journal;
Close(results);
Cuando lo ejecuto, funciona perfectamente. Ahora, si agrego la siguiente columna al script como:
dt = Current Data Table();
column_names = dt << Get Column Names(Character, String);
results = New Table("Results", Invisible,
New Column("Parameter", character),
New Column("GEN", numeric),
New Column("Ne", numeric),
Add Rows(NItems(column_names))
);
for(i=1,i<=NItems(column_names),i++,
colName = column_names;
cmd = Expr(
ow = dt << Contingency(
Y( :GEN, :Ne),
X( :cName ),
Contingency Table( 0 ),
Mosaic Plot( 0 )
);
Wait(0.1);
owr = ow<
pValueBox= owr[NumberColBox(6)];
pValue = pValueBox[2];
pValueBox1= owr[NumberColBox(12)];
pValue1 = pValueBox1[2];
Column(results,1) = colName;
Column(results,2) = pValue;
Column(results,3) = pValue1;
ow << Close Window;
);
Substitute Into(cmd,Expr(cName),colName);
Eval(cmd);
);
results << journal;
Close(results);
El guión deja de funcionar. Revisé la estructura de árbol y NumberColBoxes se refiere a los valores correctos, pero no se crea ningún diario y el script deja de funcionar.
Ahora, mis preguntas:
1) ¿Cómo puedo resolver mi situación cuando tengo varios valores estáticos de Y?
2) ¿Es posible definir el valor dinámico de Y a través de un bucle? Se me complica tener dos bucles.
Agradezco cualquier sugerencia.
