全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5138 3
2016-07-21
本人最近在研究倾向性评分PSM,打算用R语言完成,想知道具体代码怎么写!求大神指点!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-7-21 17:32:54
看看郭申阳著的《倾向值分析:统计方法与应用》一书P114-115中有用R的optmatch函数进行全匹配分析的例子。
#optmatch using glm to create propensity scores and then optmatch
set.seed(10)
library(foreign)
cds <- read.dta("chpt5_2.dta")
attach(cds)
#logistic regression
lcds <- glm(kuse ~ pcg_adc + age97 + mratio96 + pcged97 + black, family=binomial, data=cds)
summary(lcds)

library(optmatch)
#create propensity scores based on the logistic regression
pdist <- mdist(lcds)
#run full match
fm <- fullmatch(pdist)
(fm.d <- matched.distances(fm, pdist, pres=TRUE))
unlist(fm.d, max)
mean(unlist(fm.d))
sum(unlist(fm.d))
stratumStructure(fm)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-7-22 17:43:09
好的,太感谢您了,我会根据你上面的代码自己去研究一下。谢啦1
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-7-24 17:43:35
pkl 发表于 2016-7-21 17:32
看看郭申阳著的《倾向值分析:统计方法与应用》一书P114-115中有用R的optmatch函数进行全匹配分析的例子。
...
好的,太感谢您了,我会根据你上面的代码自己去研究
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群