<![CDATA[
<p>사회과학 논문에 Probit으로 추정을 하고 이것을 바탕으로 PSM(성향점수매칭)을 한다음 ols를 하려고 합니다.</p><p> </p><p>일반적인 서양매칭방법 수순인듯 합니다.</p><p> </p><p>물론 앞에 로짓분석을 할수도 있다고 합니다만...</p><p> </p><p>아무튼 먼저 순서형 probit분석을 R로 할 수 있어야 되는데요. 대체 어디서 이런걸 공부할 수 있을까요?ㅠㅠ</p><p> </p><p>공부할 수 있는 교재나, 인터넷 자료 구하기가 너무 어렵네요.</p><p> </p><p>R3.2.4최신 버전을 사용하는데 인터넷 영문자료는 옛날게 많아서 그대로 되질 않네요ㅠㅠ</p><p> </p><p>예컨데 인터넷에는 다음과 같이 설명되어있는데요.(</p><div class="f kv _SWb" style="font: small/16px arial, sans-serif; height: 17px; text-align: left; color: rgb(128, 128, 128); text-transform: none; text-indent: 0px; letter-spacing: normal; word-spacing: 0px; display: block; white-space: nowrap; widows: 1; font-size-adjust: none; font-stretch: normal; background-color: rgb(255, 255, 255); -webkit-text-stroke-width: 0px;"><cite class="_Rm" style="color: rgb(0, 102, 33); font-size: 14px; font-style: normal;"><a href="http://www.upa.pdx.edu/IOA/newsom/da2/ho_ordinal%20examples.pdf">www.upa.pdx.edu/IOA/newsom/da2/ho_<b>ordinal</b>%20examples.pdf</a></cite><div class="action-menu ab_ctl" style="margin: 0px 3px; vertical-align: middle; display: inline; position: relative; -webkit-user-select: none;"><a class="_Fmb ab_button" id="am-b18" role="button" aria-expanded="false" aria-haspopup="true" style="margin: 0px; padding: 0px; border-radius: 0px; border: 0px currentColor; transition:none; width: 13px; height: 12px; text-align: center; color: rgb(128, 128, 128); line-height: 27px; font-family: arial, sans-serif; font-size: 11px; font-weight: bold; text-decoration: none; display: inline-block; cursor: default; min-width: 0px; box-shadow: none; background-image: none; background-color: white; -webkit-user-select: none;" aria-label="Result details" href="https://www.google.co.kr/#" data-ved="0ahUKEwifj6eG__vLAhUG6KYKHSZ2Ay04ChDsHQhMMAg" jsaction="m.tdd;keydown:m.hbke;keypress:m.mskpe"><span class="mn-dwn-arw" style="border-width: 5px 4px 0px !important; border-style: solid; border-color: rgb(0, 102, 33) transparent; left: 0px; top: 7.5px; width: 0px; height: 0px; margin-top: -4px; margin-left: 3px; position: absolute;"></span></a></div></div><p> </p><p>Probit Model in R
(Same precautions regarding missing data and nonnumeric variables apply to p</p><p>robit models).
</p><p>> library(MASS)
</p><p>> model <-polr(cutmeal ~ mosmed + depress1 + educat +
marital,data=mydata,contrasts=NULL,method=c("probit"))</p><p>
> summary(model,digits = 3)
</p><p> </p><p>Re-fitting to get Hessian
</p><p> </p><p>Call:
polr(formula = cutmeal ~ mosmed + depress1 + educat + marital,
data = mydata, contrasts = NULL, method = c("probit"))
</p><p> </p><p>Coefficients:
Value Std. Error t value
mosmed 0.00667 0.0136 0.488
depress1 0.12133 0.0230 5.281
educat -0.01374 0.0649 -0.212
marital 0.26103 0.1359 1.921
</p><p> </p><p>Intercepts:
Value Std. Error t value
0|1 1.152 0.198 5.829
1|2 1.587 0.203 7.833</p><p> </p><p>
Residual Deviance: 717.29
AIC: 729.29 </p><p>---------------------------------------------------------------------------------------------------------------</p><p>위에거 안되구요. 교재 189page에 있는 프로빗 분석을 하려고 하면 reponse의 범위가 0~1사이여야 하기 때문에 안된다고 합니다.</p><p> </p><p>제가 하려는건 다항 순서형 프로빗이기 때문에 당연히 response가 더 많습니다.</p>
<!-- -->
]]>
카페 게시글
묻고 답하기
성향점수매칭 전 순서형 프로빗 분석법 문의
마가렛
추천 0
조회 445
16.04.07 16:16
댓글 2
북마크 번역하기 공유하기 기능 더보기
다음검색
첫댓글 공부는... 어렵습니다. logistic/probit regression 자체가 GLMs(generalized linear models)라고 어려운 이론인데, 이건 종속변수가 ordinal로 그걸 확장한 개념입니다. Agresti가 쓴 Analysis of Ordinal Categorical Data가 있는데 어려울겁니다.
polr()을 사용하면 됩니다. method="probit"도 맞는데요. 10여년 전에 polr()을 사용했었는데 그때 잘됐었는데 지금 안될 것 같지는 않은데요.