671409_fr-FR

liliconcom

Utilisation des effets d'interaction dans le modèle d'ajustement JSL

Salut la communauté,

J'essaie de créer un script qui crée dynamiquement un modèle avec les principaux effets et interactions.

 

Par exemple, je pourrais avoir un sous-ensemble de colonnes d'effets définies comme suit :

 

 

all_cols = Current Data Table () << Get Column Names(String);
effect_cols = all_cols[3::8];

 

Et puis je construis une série de modèles réduisant les effect_cols.

Par la suite, je considère les termes d'interaction. Cependant, j'ai du mal à créer une liste de noms de colonnes ou de références de colonnes pour les termes d'interaction.

J'ai essayé quelque chose comme ceci:

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

et le construire dans un modèle avec ceci :

Fit Model(
  Y( :"Result"n ),
  Effects(
   concat(effect_cols,interaction_effects),
  Personality( "Standard Least Squares" ),
  Emphasis( "Effect Screening" ),
  Run( ),
)

Cependant, cela donne systématiquement une erreur et mes autres tentatives d'utilisation de parse() et eval() ne fonctionnent pas.

 

Toute aide est très appréciée.

 

Merci

Modélisation statistique avancéeAutomatisation et scriptsmacOSRe : Utilisation des effets d'interaction dans le modèle d'ajustement JSL

Salut @mkiesz ,

avez-vous essayé de scripter la plateforme pas à pas ? certains de vos critères de sélection de modèle peuvent être intégrés.

 

Re : Utilisation des effets d'interaction dans le modèle d'ajustement JSL

Si je me souviens bien, j'ai essayé avec quelque chose comme ça

 

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