Extracting inverse mills ratios from Proc QULIM

[markutrgv]

I am trying to download inverse mills ratios from PROC QULIM, but can't seem to get the output statement right.  Below is the code I have written in order to simply look at the various new data that QULIM is supposed to generate.  Anyone know what I'm doing wrong?

proc qlim;
model lf = fceoedu;
endogenous lf ~ discrete;
output out=entrfndr2;
data entrfndr3; set work.entrfndr2;
proc print;
var ipofirms mills;

[StatDave]

Use the MARGINAL option in the OUTPUT statement. See the example in this note.