Question

我有类似的数据集:

val<-c("Y","N")
test<-data.frame(age=rnorm(n=100,mean=50,sd=10),var1=sample(val,100,T),var2=sample(val,100,T),var3=sample(val,100,T),sex=sample(c("F","M"),100,T))

I´d like to create a summary reporting the mean age for each category using Hmisc.

library(Hmisc)
summary.formula(age~sex+var1+var2+var3,data=test)

然而,瓦尔1-瓦尔3实际上属于相同的分类变量,其等级为瓦尔1、瓦尔和瓦尔3,而不是Y/N。此外,这些并非相互排斥。因此,有可能产生一种变数的变量,即不同的等级,不相互排斥和类型

summary.formula(age~sex+var4,data=test)

and have an output like:

+-------+-+---+----+
|       | |N  |age |
+-------+-+---+----+
|sex    |F| 44|48.0|
|       |M| 56|50.8|
+-------+-+---+----+
|var4   |var1| xx|yy|
|       |var2| xx|yy|
        |var3| xx|yy|
+-------+-+---+----+
|Overall| |100|49.6|
+-------+-+---+----+

任何帮助都将受到高度赞赏。

Answer 1

如何填写<代码>paste0? (或过去(......,ep=......)

> test$var4 <- paste0(test$var1, test$var2, test$var3)
> summary.formula(age~sex+var4, data=test)
age    N=100

+-------+---+---+--------+
|       |   |  N|     age|
+-------+---+---+--------+
|    sex|  F| 50|50.25440|
|       |  M| 50|51.32134|
+-------+---+---+--------+
|   var4|NNN| 13|46.64417|
|       |NNY| 17|51.34456|
|       |NYN| 15|52.92185|
|       |NYY| 17|47.35685|
|       |YNN|  9|50.91647|
|       |YNY|  7|48.04489|
|       |YYN| 10|53.23713|
|       |YYY| 12|56.14394|
+-------+---+---+--------+
|Overall|   |100|50.78787|
+-------+---+---+--------+
>

Answer 2

I think the problem lies in that you are trying to combine statistics for two different data sets:

按个人分列的数据:

summary.formula(age~sex, test)

# age    N=100
# 
# +-------+-+---+--------+
# |       | |N  |age     |
# +-------+-+---+--------+
# |sex    |F| 35|49.99930|
# |       |M| 65|48.96266|
# +-------+-+---+--------+
# |Overall| |100|49.32548|
# +-------+-+---+--------+

按汽车分列的数据

在这里,你需要每车一行;这里是创造数据的一种方式,但我确信,必须采取更恶毒的方式:

    var1 <- subset(test, var1 == "Y", c("age", "sex"))
    var2 <- subset(test, var2 == "Y", c("age", "sex"))
    var3 <- subset(test, var3 == "Y", c("age", "sex"))
    var1$var <- "var1"
    var2$var <- "var2"
    var3$var <- "var3"
    vars <- rbind(var1, var2, var3)

然后是汇总统计数据:

    summary.formula(age~var, data=vars)
    # age    N=147
    # 
    # +-------+----+---+--------+
    # |       |    |N  |age     |
    # +-------+----+---+--------+
    # |var    |var1| 47|48.91983|
    # |       |var2| 43|46.31811|
    # |       |var3| 57|49.35292|
    # +-------+----+---+--------+
    # |Overall|    |147|48.32672|
    # +-------+----+---+--------+

如你所知,两份摘要的<代码>Overall各节并不匹配,因为它们来自两个不同的数据集。 (不可能把你要求的方式结合起来。)

友情链接