昨天找到了一个例子可以解决这个问题
https://www.rdocumentation.org/packages/sampleSelection/versions/1.2-0/topics/invMillsRatio

library("sampleSelection")
data(Mroz87)
myProbit <- glm( lfp ~ nwifeinc + educ + exper + I( exper^2 ) + age +
   kids5 + kids618, family = binomial( link = "probit" ), data=Mroz87 )
Mroz87$IMR <- invMillsRatio( myProbit )$IMR1
myHeckit <- lm( log( wage ) ~ educ + exper + I( exper^2 ) + IMR,
   data = Mroz87[ Mroz87$lfp == 1, ] )

谢谢楼主的分享，受用了