671409_ja-JP

liliconcom

JSL フィット モデルでの相互作用効果の使用

コミュニティさん、こんにちは。

主効果と相互作用を持つモデルを動的に作成するスクリプトを作成しようとしています。

 

たとえば、次のように定義された効果列のサブセットがあるとします。

 

 

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

 

次に、effect_cols を削減する一連のモデルを構築します。

続いて相互作用項を考えます。 ただし、相互作用用語の列名または列参照のリストを作成するのに問題があります。

私は次のようなことを試しました：

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

そしてこれをモデルに構築します:

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

ただし、これでは一貫してエラーが発生し、parse() と eval() を使用する他の試みは機能しません。

 

ご協力をよろしくお願いいたします。

 

ありがとう

高度な統計モデリング自動化とスクリプト作成マックOSRe: JSL フィット モデルでの相互作用効果の使用

こんにちは@mkiesz 、

ステップワイズ プラットフォームのスクリプトを作成しようとしましたか?モデル選択基準の一部が組み込まれている場合があります。

 

Re: JSL フィット モデルでの相互作用効果の使用

私の記憶が正しければ、このようなものを試しました

 

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