是否有R等效的Stata包装Grinter?我想用R中的逻辑回归方程来描绘交互效应,我似乎无法使我的任何图形起作用。如果R中有一个命令可以完成Grinter在Stata中的操作,那将非常有帮助。
这是我的等式:
svyglm(onset.new~Gini+preci.ma.30.CRU+lag.polity+
pop.growth+lpop+lag.lGDP.pc+Oil+ethfrac+lmtnest+Gini*preci.ma.30.CRU,
design=NEWD,na.action=na.omit,family=quasibinomial)
我感兴趣的是绘制交互项Gini * preci.ma.30.CRU(不等式乘以降水量变化)以查看这两个变量组合如何影响冲突起始(onset.new)。所有其他变量都只是控件。
谢谢!
答案 0 :(得分:0)
您可以使用回归方程来计算预测的发病率。新的基尼系数和降水量(新的基尼系数和低/高降水量)。然后,您可以创建具有这些预测值的数据集,并将其绘制为条形图或折线图。数据集如下所示:
predict.onset.new基尼降水量 .34低低 .23高低 .81低高 .28高高