glmfitのインジケータ関数やハンドルを作成するには？

Question

ハンドルまたはインジケータ機能の作成に関する質問があります。私は、4つの説明変数と1つの列（最後の1列、5）の1と2の観測値がグループ1に属しているか2にグループ2に属しているかを示すXを含むX行列を持っています。グループ1に属する観測値に対して1つ、グループ2に属する観測値に対して1つです。したがって、glmfitが特定のグループの観測値のみを計算するように、ある種の指示関数が必要です。誰かが私にこれをどうすればできるのか教えてくれますか？私は次のglmfitを使用します：

[B1, dev, stats1] = glmfit(X(:,1:4), Y, 'binomial', 'link', 'logit'); 

Answer 1

次は、この作業を行いますか？上記で

indicator = X(:, 5) 
[B1, dev, stats1] = glmfit(X(indicator==1,1:4), Y, 'binomial', 'link', 'logit'); 
[B2, dev, stats2] = glmfit(X(indicator==2,1:4), Y, 'binomial', 'link', 'logit'); 

、X(indicator==1, 1:4)は論理的インデックスと呼ばれるもの採用し、hereも参照。 indicatorのエントリが1である行のみからなるXの部分行列を提供します。