data=read.csv('F:/SRC/GPSM.csv',stringsAsFactors = F)
data=data[,-1]
data$treat=as.factor(data$group)
mnps.AOD <- mnps(treat ~ Gender + Age + Location + stage + size+Pre_chemo+resection,
data = data,
estimand = "ATE",
verbose = FALSE,
stop.method = c("es.mean", "ks.mean"),
n.trees = 3000)
library(survey)
data$w <- get.weights(mnps.AOD, stop.method = "es.mean")
design.mnps <- svydesign(ids=~1, weights=~w, data=data)
st=strata(data, stratanames=c('treat'), size=c(300,300,300), method=c("srswor","srswr","poisson",
"systematic")[4], pik=1/data$w,description=T)
sda=data[st$ID_unit,]
write.csv(sda,'f:/SRC/SDA.CSV')
NULL
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- Given a binary treestruct TreeLinkNode {TreeLinkNode *lef...
- [NSNullnull]-- nil ---:(null)NULL--(null) Nil--(null)
- 默认情况下，总是有三个文件处于打开状态，标准输入(键盘输入)、标准输出（输出到屏幕）、标准错误（也是输出到屏幕），...
