2018-10-24
閱讀量:
1611
怎么計算inverse mills ratio(基于R)
R語言做probit回歸之后,如何計算出inverse mills ratio?
在經(jīng)歷多方幫助(此處真的心累)算是找到了一個解決的辦法:
library("sampleSelection")
data(Mroz87)
myProbit <- glm( lfp ~ nwifeinc + educ + exper + I( exper^2 ) + age +
? ?kids5 + kids618, family = binomial( link = "probit" ), data=Mroz87 )
Mroz87IMR<-invMillsRatio(myProbit)IMR<?invMillsRatio(myProbit)IMR1
myHeckit <- lm( log( wage ) ~ educ + exper + I( exper^2 ) + IMR,
? ?data = Mroz87[ Mroz87$lfp == 1, ] )
寫貼不易,你確定不點個贊嗎??






評論(0)


暫無數(shù)據(jù)
推薦帖子
0條評論
0條評論
0條評論
0條評論