OpenBugs/R范畴变量“是”，“否”作为1,0

Question

我有一个简单的后勤回归

logit(p[i]) <- b0 + b1*x[i]

但是圆周率的数据是“是”，“不是”。如何将(是/否)编码到(1,0)？

像这样的smt

logit(encode(p[i])) <- b0 + b1*x[i]

Answer 1

试试旧的时尚方式！我假设p和x是数据列，我将它们合并在一个名为df的数据框架中。

library(tidyverse)
library(stats) # for glm
df <- as.dataframe(cbind(p,x))
df %>% mutate(p=ifelse(p=="yes",1,0)
logit.model <- glm(p~x, data=df, family=binomial(link="logit"))
summary(logit.model)