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

Question

我在Stata中有一个glm-formula，我需要将其“翻译”为R.我对它的外观有一个粗略的想法，但我不确定，因为我从未使用过Stata。

这是一个带有不合理术语的glm，其中交互项ridit_survey是两个术语ridit*survey的乘积（即在Stata中：ridit_survey=ridit*survey。

这是Stata中的模型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))

第二个模型：

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

我会翻译成

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

这是否正确（如果没有，它会如何正确）？