Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Forecasting
- /
- Inverse mills ratio or heckman question

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 03-30-2015 09:23 AM
(4733 views)

I am trying to determine if there is possibility unmeasured bias in my model. I have a continuous dep variable and binary ind variable which is TREATED (0/1). I want to determine if I have unmeasured bias. I don't have proc qlim so I need to create an inverse Mills ratio and run it through a GLM.

I first create a PROBIT model and output the estimated probabilities (prob) of being treated.

Next, I calculate the Inverse Mills Ratio:

IMR = pdf('NORMAL', prob ) / cdf('NORMAL', prob ); /*inverse mills ratio*/

Then run my GLM:

proc glm data = weighted_PS;

class RHS;

model LHS = RHS IMR/ ss3 solution;

weight weights;

run;

Is this correct?

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

You are correct that PROC QLIM is the most efficient, accurate and simplest way to accomplish Heckman's 2-step estimator.

You will need to correct your standard errors in your second stage. See the QLIM documentation. SAS/ETS(R) 13.2 User's Guide

You might just want to use PROC REG with some HCCME= options to correct your standard errors. SAS/STAT(R) 9.2 User's Guide, Second Edition

But you are close as is.

1 REPLY 1

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

You are correct that PROC QLIM is the most efficient, accurate and simplest way to accomplish Heckman's 2-step estimator.

You will need to correct your standard errors in your second stage. See the QLIM documentation. SAS/ETS(R) 13.2 User's Guide

You might just want to use PROC REG with some HCCME= options to correct your standard errors. SAS/STAT(R) 9.2 User's Guide, Second Edition

But you are close as is.

**Don't miss out on SAS Innovate - Register now for the FREE Livestream!**

Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.

Multiple Linear Regression in SAS

Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.

Find more tutorials on the SAS Users YouTube channel.