我想从书中复制一些计算(logit回归)。这本书给出了列联表和结果。
表格如下:
.
example <- matrix(c(21,22,6,51), nrow = 2, byrow = TRUE)
#Labels:
rownames(example) <- c("Present","Absent")
colnames(example) <- c(">= 55", "<55")
它给了我这个:
>= 55 <55
Present 21 22
Absent 6 51
但是要使用glm()函数,数据必须采用以下方式:
(两个列,一个带有“年龄”,另一个带有“当前”,填充为0/1)
age <- c(rep(c(0),27), rep(c(1),73))
present <- c(rep(c(0),21), rep(c(1),6), rep(c(0),22), rep(c(1),51))
data <- data.frame(present, age)
> data
present age
1 0 0
2 0 0
3 0 0
. . .
. . .
. . .
100 1 1
是否有一种简单的方法可以从表/矩阵中获取此结构?
reshape2::melt(example)
这会给你,
Var1 Var2 value
1 Present >= 55 21
2 Absent >= 55 6
3 Present <55 22
4 Absent <55 51
您可以轻松使用 glm
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句