假设我有一个m x n数据矩阵M,其中第一列是我的二元响应变量,其余列是预测变量。
我使用
创建了一个多元逻辑回归模型my.model = glm(M[,1] ~ M[,2] + M[,3] + ... M[,n], family=binomial("logit"))
我的问题是如何在不必写出M [,2] + ... + M [,n]项的情况下如何做到这一点,例如:有些事情只是:
my.model = glm(M[,1] ~M[2:n], family = binomial("logit"))
(我知道后者不会工作,但我正在寻找一些可以使用矩阵的列表/块来做类似的事情。)
答案 0 :(得分:0)
你为什么不这样做:
glm(M[,1] ~., family = binomial("logit))
这将调用所有列,或者如果您需要特定列,请事先从数据框中过滤掉。