Ciao Comunità,
Sto provando a creare uno script che crei dinamicamente un modello con effetti e interazioni principali.
Ad esempio, potrei avere un sottoinsieme di colonne di effetti definite da quanto segue:
all_cols = Current Data Table () << Get Column Names(String);
effect_cols = all_cols[3::8];
E poi creo una serie di modelli riducendo l'effetto_cols.
Successivamente, considererò i termini di interazione. Tuttavia, ho difficoltà a creare un elenco di nomi di colonna o riferimenti di colonna per i termini di interazione.
Ho provato qualcosa del genere:
if(length(effect_cols) > 1,
interaction_effects = {};
for(i = 1, i < length(effect_cols), i++,
for(j=i+1,j<= length(effect_cols), j ++,
inter_tmp = effect_cols[i] || "*" || effect_cols[j];
interaction_effects = insert(interaction_effects,inter_tmp);
)
)
);e costruendolo in un modello con questo:
Fit Model(
Y( :"Result"n ),
Effects(
concat(effect_cols,interaction_effects),
Personality( "Standard Least Squares" ),
Emphasis( "Effect Screening" ),
Run( ),
)
Tuttavia, questo genera costantemente un errore e gli altri miei tentativi di utilizzare parse() ed eval() non funzionano.
Ogni aiuto è molto apprezzato.
Grazie
CIAO @mkiesz ,
hai provato a creare script per la piattaforma graduale? potrebbe avere alcuni dei criteri di selezione del modello integrati.
Se ricordo bene ho provato con qualcosa del genere
Names Default To Here(1);
dt = Open("$SAMPLE_DATA/Tiretread.jmp");
effect_cols = {"SILICA", "HARDNESS", "SILANE", "SULFUR"};
terms = {};
For Each({first_term, idx}, effect_cols[1::N Items(effect_cols) - 1],
For Each({second_term}, effect_cols[idx+1::N Items(effect_cols)],
cur_interaction = Substitute(Expr(_firstterm_*_secondterm_),
Expr(_firstterm_), NameExpr(AsColumn(first_term)),
Expr(_secondterm_), NameExpr(AsColumn(second_term))
);
Insert Into(terms, Name Expr(cur_interaction));
);
);
fitm = dt << Fit Model(
Y(:ABRASION),
Effects(Eval(terms)),
Personality("Standard Least Squares"),
Emphasis("Effect Screening")
);
