6772_fr-FR

liliconcom

Régression logistique - situation compliquée<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />

Salut ,

Avant, je décris mon problème, je dois mentionner que je travaille avec JMP7.

J'ai un tableau avec 12 colonnes, les 6 premières colonnes sont des paramètres descriptifs pour chacune des 6 prochaines colonnes. Les 6 premières colonnes sont toutes des caractères et les 6 dernières colonnes sont numériques, binomiales.

Je dois exécuter une régression logistique pour chacune des 6 dernières colonnes et utiliser les 6 premières colonnes comme prédicteurs pour extraire la valeur "Prob> ChiSq" du test de "Pearson".

Maintenant, j'ai créé le script suivant pour l'une des colonnes "Y, Réponse":

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);

Quand je le lance, cela fonctionne parfaitement. Maintenant, si j'ajoute la colonne suivante au script en tant que:

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);

Le script cesse de fonctionner. J'ai vérifié la structure arborescente et les deux NumberColBox font référence à des valeurs correctes, mais aucun journal n'est créé et le script cesse de fonctionner.

Maintenant, mes questions:

1) Comment puis-je résoudre mon problème lorsque plusieurs valeurs de Y sont statiques?

2) Est-il possible de définir la valeur dynamique Y à travers une boucle? Cela devient compliqué pour moi d'avoir deux boucles.

J'apprécie pour toutes les suggestions.