我正在尝试使用数据集在R中进行比例测试,该数据集报告全球队列以及来自东西方人口的亚群体。在示例代码中,6/15 pts是EAST,5/15是WEST,4/15是GLOBAL,因为它是未知的,但15/15应报告为全局,因为它们是整体的子集。
set.seed(1)
dataset<-data.frame("sex"=sample(c("M","F"),15,replace=TRUE),"cohort"=sample(c("GLOBAL","EAST","WEST"),15,replace=TRUE),"grade"=sample(c(1,2,3),15,replace=TRUE))
table(dataset$sex,dataset$cohort)
EAST GLOBAL WEST
F 4 2 2
M 2 2 3
prop.table(table(dataset$sex,dataset$cohort))
EAST GLOBAL WEST
F 0.2666667 0.1333333 0.1333333
M 0.1333333 0.1333333 0.2000000
如何强制它反映:
EAST GLOBAL WEST
F 4 8 2
M 2 7 3
然后对这些“3”人群进行比例测试?对于数据集$ grade,需要完成相同的操作,该等级有3个类别。