我收到了一个包含多个虚拟变量和其他变量的数据集。基本上我想做的是用rms的summary.formula创建汇总表。但是,我不知道如何从多个虚拟变量中创建单个变量,并且它们不是互斥的。这是可能吗。当然我可以创建一个表格等,但是我不能使用 summary.formula 并且我希望 summary.formula 输出只包含虚拟变量的各个级别。
编辑:澄清:a & b 需要总结,但它们并不相互排斥。由于每一行都记录了年龄,因此我需要将 a 和 b 汇总为一个变量,以便在 summary.formula 中使用。我已经编辑了下面的代码,以便将 0 和 1 分别更改为 NA 或 a、b。
我希望 summary.formula 输出是这样的:
h<-data.frame(a=sample(c("A",NA),100,replace=T),b=sample(c("B",NA),100,replace=T),age=rnorm(100,50,25),epo=sample(c("Y","N"),100,T))
library(rms)
summary.formula(epo~age####+summary variable of a & b######,method="reverse",data=h)
#-----------------
Descriptive Statistics by epo
+---------+--------------------------+--------------------------+
| |N |Y |
| |(N=56) |(N=44) |
+---------+--------------------------+--------------------------+
|age |31.53434/48.90788/67.69096|28.63689/43.93502/57.81834|
+---------+--------------------------+--------------------------+
|sab : A | 25% (14) | 16% ( 7) |
+---------+--------------------------+--------------------------+
| B | 27% (15) | 32% (14) |
+---------+--------------------------+--------------------------+