我想要一个平滑的代码,以便我可以将带有预测变量列表的函数作为矩阵的一部分来调用,然后再对预测变量使用逻辑回归。当我想在我的逻辑回归中使用一个选择的主题时,我的问题就出现了,这是我需要为一个估计量做的事情,因为逻辑回归中系数的名称变成了矩阵名称+变量名称。因此,当我想基于逻辑回归的系数获得预测值时,会出现错误。因此,我想找到一个解决方案,使系数的标签不包括矩阵的名称。当我不减少数据集时,代码将起作用,因为它将认为两个数据集具有相同数量的观测值。
我希望有人可以帮助我。
true.conf <- cbind(var1=simdata$var1,var2=simdata$var2,simdata$var3=var3)
true.conf = true.conf[simdata$treatment==0,]
glm0.out=glm(Y~true.conf,family=binomial,data=simdata)
t.0=predict(model=glm0.out,newdata=simdata,type="response")
这是现在我输入glm0.out时的样子:
"Coefficients:
(Intercept) true.confvar1 true.confvar2 true.confvar3"
这是我希望输入glm0.out时的样子:
"Coefficients:
(Intercept) var1 var2 var3"
当然,如果有一种方法可以使用预测而不改变键入时系数的名称,那对我也很好。