R语言做倾向性评分的代码谁有？求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人最近在研究倾向性评分PSM，打算用R语言完成，想知道具体代码怎么写！求大神指点！

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：R语言 倾向性 PSM

看看郭申阳著的《倾向值分析：统计方法与应用》一书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)

好的，太感谢您了，我会根据你上面的代码自己去研究一下。谢啦1

pkl 发表于 2016-7-21 17:32
看看郭申阳著的《倾向值分析：统计方法与应用》一书P114-115中有用R的optmatch函数进行全匹配分析的例子。
...

好的，太感谢您了，我会根据你上面的代码自己去研究