671409_es-ES

liliconcom

Uso de efectos de interacción en el modelo de ajuste JSL

Hola comunidad,

Intento crear un script que cree dinámicamente un modelo con efectos e interacciones principales.

 

Por ejemplo, podría tener un subconjunto de columnas de efectos definido por lo siguiente:

 

 

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

 

Y luego construyo una serie de modelos reduciendo los effect_cols.

Posteriormente, considero términos de interacción. Sin embargo, tengo problemas para crear una lista de nombres de columnas o referencias de columnas para términos de interacción.

He probado algo como esto:

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

y construyéndolo en un modelo con esto:

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

Sin embargo, esto constantemente da un error y mis otros intentos de usar parse() y eval() no funcionan.

 

Cualquier ayuda es muy apreciada.

 

Gracias

Modelado estadístico avanzadoAutomatización y secuencias de comandosMac OSRe: Uso de efectos de interacción en el modelo de ajuste JSL

Hola @mkiesz ,

¿Intentaste escribir un script para la plataforma paso a paso? es posible que tenga incorporados algunos de sus criterios de selección de modelo.

 

Re: Uso de efectos de interacción en el modelo de ajuste JSL

Si no recuerdo mal he probado con algo como esto

 

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