将Stata代码转换为R代码：glm

我在Stata中有一个glm公式，我需要将其“翻译”为R。我大致知道它是什么样子，但我不确定，因为我从未使用过Stata

这是一个具有惯性项的glm，其中交互项

ridit_survey

是两个项

ridit*survey

的产物（即，在Stata:

ridit_survey=ridit*survey

）

这是斯塔塔的模型1

xi: glm smoking ridit survey ridit_survey age if men==1, fam(bin) link(log) nolog eform

我会把它翻译成R，像这样：

fit <- glm(xi ~ smoking + age + ridit * survey, family = binomial("log"),
           data = subset(dataframe, subset = men == 1))

我会翻译成

fit <- glm(xi ~ smoking + age + ridit * survey, family = binomial("identity"),
           data = subset(dataframe, subset = men == 1))

---

fit我手头没有STATA，因此很遗憾，我无法安装和检查。您的STATA中的
xi:
[NB拼写]当您通过附加变量指定交互时，调用是多余的。R人们可以注意到那里没有什么可翻译的。
nolog
只是调整您得到的输出量。
eform
指定指数系数。帮助在Ok公开，所以结果是什么？谢谢，我想是这样的回答我的问题！如果你愿意，你可以将你的评论转换成答案，这样我就可以接受。绝对：
xi:
是一个命令前缀，而不是变量名。我不能确定使用什么R代码，因为我根本不是一个经常使用的用户。我手头没有STATA，所以很遗憾，我无法安装和检查。
xi:
你的状态[注意拼写]当您通过附加变量指定交互时，调用是多余的。R人们可以注意到那里没有什么可翻译的。
nolog
只是调整您得到的输出量。
eform
指定指数系数。帮助在Ok公开，所以结果是什么？谢谢，我想是这样的回答我的问题！如果你愿意，你可以将你的评论转换为答案，这样我就可以接受。绝对：
xi:
是命令前缀而不是变量名。我不能确定使用什么R代码，因为我根本不是一个经常使用的用户。
fit <- glm(xi ~ smoking + age + ridit * survey, family = binomial("identity"),
           data = subset(dataframe, subset = men == 1))