Question

我有关于创建句柄或指标功能的问题。我有一个X矩阵，包含4个解释变量和一个列（最后一个，第5列）的一个和两个，表示观察是否属于组1乘1或组2乘2.我想执行2 glmfit。一个用于属于组1的观察，一个用于属于组2的观察。因此，我需要某种指示函数，以便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。它提供X的子矩阵，仅包含indicator条目为1的行。

如何为glmfit制作指示功能或处理？

1 个答案: