. * praccum:      version 1.6.2 1/5/01
. * verified:     2/23/01
. which praccum
c:\stata6\ado\stbplus\p\praccum.ado
*! version 1.6.3 3/17/01 backwards compatibility with v6 bug

. 
. * Note: this is to check if praccum matches results of prvalue.
. * It does not check that prvalue is correct.
. 
. * cloglog
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. quietly cloglog lfp k5 k618 age age2 wc hc lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=inLF|x):       0.7014   95% ci: (0.3840,0.9510)
  Pr(y=NotInLF|x):    0.2986   95% ci: (0.0490,0.6160)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=inLF|x):       0.7123   95% ci: (0.5126,0.8846)
  Pr(y=NotInLF|x):    0.2877   95% ci: (0.1154,0.4874)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=inLF|x):       0.7027   95% ci: (0.5975,0.8015)
  Pr(y=NotInLF|x):    0.2973   95% ci: (0.1985,0.4025)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=inLF|x):       0.6725   95% ci: (0.6165,0.7274)
  Pr(y=NotInLF|x):    0.3275   95% ci: (0.2726,0.3835)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=inLF|x):       0.6210   95% ci: (0.5672,0.6751)
  Pr(y=NotInLF|x):    0.3790   95% ci: (0.3249,0.4328)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=inLF|x):       0.5494   95% ci: (0.4941,0.6064)
  Pr(y=NotInLF|x):    0.4506   95% ci: (0.3936,0.5059)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=inLF|x):       0.4613   95% ci: (0.4073,0.5190)
  Pr(y=NotInLF|x):    0.5387   95% ci: (0.4810,0.5927)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=inLF|x):       0.3647   95% ci: (0.2904,0.4513)
  Pr(y=NotInLF|x):    0.6353   95% ci: (0.5487,0.7096)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=inLF|x):       0.2698   95% ci: (0.1706,0.4105)
  Pr(y=NotInLF|x):    0.7302   95% ci: (0.5895,0.8294)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqp0 |       9    .4383192   .1628331   .2877253    .730227  
  agsqp1 |       9    .5616808   .1628331    .269773   .7122747  

. list agsq* in 1/9

         agsqx     agsqp0     agsqp1 
  1.        20   .2985594   .7014406  
  2.        25   .2877253   .7122747  
  3.        30   .2972831   .7027169  
  4.        35    .327544    .672456  
  5.        40   .3789573   .6210427  
  6.        45   .4506453   .5493547  
  7.        50   .5386655   .4613345  
  8.        55   .6352658   .3647342  
  9.        60    .730227    .269773  

. 
. * cnreg
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. gen lwgcens = lwg

. replace lwgcens = 0 if lwg < 0
(18 real changes made)

. replace lwgcens = 2 if lwg > 2
(45 real changes made)

. gen censor = 0

. replace censor = -1 if lwgcens == 0
(22 real changes made)

. replace censor = 1 if lwgcens == 2
(45 real changes made)

. quietly cnreg lwgcens lfp k5 k618 age age2 wc hc inc, censor(censor)

. prvalue, x(age 20 age2 400) rest(mean) brief

Predicted value of y*:  .7988954   95% ci: ( .5062962, 1.091494)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

Predicted value of y*:  .9361635   95% ci: ( .7580404, 1.114287)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

Predicted value of y*:  1.039235   95% ci: ( .9439359, 1.134533)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

Predicted value of y*:  1.108108   95% ci: ( 1.056405, 1.159812)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

Predicted value of y*:  1.142785   95% ci: ( 1.094899, 1.190671)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

Predicted value of y*:  1.143265   95% ci: ( 1.093008, 1.193521)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

Predicted value of y*:  1.109547   95% ci: ( 1.057112, 1.161982)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

Predicted value of y*:  1.041633   95% ci: ( .9571576, 1.126107)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

Predicted value of y*:  .9395207   95% ci: ( .7825143, 1.096527)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqxb |       9    1.028795   .1165047   .7988954   1.143265  

. list agsq* in 1/9

         agsqx     agsqxb 
  1.        20   .7988954  
  2.        25   .9361635  
  3.        30   1.039235  
  4.        35   1.108108  
  5.        40   1.142785  
  6.        45   1.143265  
  7.        50   1.109547  
  8.        55   1.041633  
  9.        60   .9395207  

. 
. /*
> * gologit
> use binlfp2, clear
> gen age2 = age*age
> gen colboth = ((wc==1)*2) + (hc==1)
> recode colboth 0=-6 1=3 2=15 3=59
> quietly gologit colboth lfp k5 k618 age age2 lwg inc
> prvalue, x(age 20 age2 400) rest(mean) brief
> praccum, saving(mage) xis(20)
> prvalue, x(age 25 age2 625) rest(mean) brief
> praccum, using(mage) xis(25)
> prvalue, x(age 30 age2 900) rest(mean) brief
> praccum, using(mage) xis(30)
> prvalue, x(age 35 age2 1225) rest(mean) brief
> praccum, using(mage) xis(35)
> prvalue, x(age 40 age2 1600) rest(mean) brief
> praccum, using(mage) xis(40)
> prvalue, x(age 45 age2 2025) rest(mean) brief
> praccum, using(mage) xis(45)
> prvalue, x(age 50 age2 2500) rest(mean) brief
> praccum, using(mage) xis(50)
> prvalue, x(age 55 age2 3025) rest(mean) brief
> praccum, using(mage) xis(55)
> prvalue, x(age 60 age2 3600) rest(mean) brief
> praccum, using(mage) xis(60) gen(agsq)
> list agsq* in 1/9
> */
. 
. * logit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. quietly logit lfp k5 k618 age age2 wc hc lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=inLF|x):       0.7324   95% ci: (0.3972,0.9191)
  Pr(y=NotInLF|x):    0.2676   95% ci: (0.0809,0.6028)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=inLF|x):       0.7318   95% ci: (0.5320,0.8675)
  Pr(y=NotInLF|x):    0.2682   95% ci: (0.1325,0.4680)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=inLF|x):       0.7164   95% ci: (0.6111,0.8024)
  Pr(y=NotInLF|x):    0.2836   95% ci: (0.1976,0.3889)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=inLF|x):       0.6847   95% ci: (0.6276,0.7368)
  Pr(y=NotInLF|x):    0.3153   95% ci: (0.2632,0.3724)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=inLF|x):       0.6343   95% ci: (0.5810,0.6846)
  Pr(y=NotInLF|x):    0.3657   95% ci: (0.3154,0.4190)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=inLF|x):       0.5627   95% ci: (0.5041,0.6196)
  Pr(y=NotInLF|x):    0.4373   95% ci: (0.3804,0.4959)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=inLF|x):       0.4699   95% ci: (0.4098,0.5308)
  Pr(y=NotInLF|x):    0.5301   95% ci: (0.4692,0.5902)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=inLF|x):       0.3619   95% ci: (0.2787,0.4543)
  Pr(y=NotInLF|x):    0.6381   95% ci: (0.5457,0.7213)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=inLF|x):       0.2520   95% ci: (0.1416,0.4078)
  Pr(y=NotInLF|x):    0.7480   95% ci: (0.5922,0.8584)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqp0 |       9    .4282142   .1752595   .2676314   .7479599  
  agsqp1 |       9    .5717858   .1752595   .2520402   .7323686  

. list agsq* in 1/9

         agsqx     agsqp0     agsqp1 
  1.        20   .2676314   .7323686  
  2.        25   .2682353   .7317647  
  3.        30   .2836163   .7163837  
  4.        35   .3152536   .6847464  
  5.        40   .3656723   .6343277  
  6.        45   .4373158   .5626842  
  7.        50   .5301194   .4698806  
  8.        55   .6381241   .3618759  
  9.        60   .7479599   .2520402  

. 
. * logit with a long name
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. quietly logit lfp k5 k618 age age2 wc hc lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=inLF|x):       0.7324   95% ci: (0.3972,0.9191)
  Pr(y=NotInLF|x):    0.2676   95% ci: (0.0809,0.6028)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=inLF|x):       0.7318   95% ci: (0.5320,0.8675)
  Pr(y=NotInLF|x):    0.2682   95% ci: (0.1325,0.4680)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=inLF|x):       0.7164   95% ci: (0.6111,0.8024)
  Pr(y=NotInLF|x):    0.2836   95% ci: (0.1976,0.3889)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=inLF|x):       0.6847   95% ci: (0.6276,0.7368)
  Pr(y=NotInLF|x):    0.3153   95% ci: (0.2632,0.3724)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=inLF|x):       0.6343   95% ci: (0.5810,0.6846)
  Pr(y=NotInLF|x):    0.3657   95% ci: (0.3154,0.4190)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=inLF|x):       0.5627   95% ci: (0.5041,0.6196)
  Pr(y=NotInLF|x):    0.4373   95% ci: (0.3804,0.4959)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=inLF|x):       0.4699   95% ci: (0.4098,0.5308)
  Pr(y=NotInLF|x):    0.5301   95% ci: (0.4692,0.5902)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=inLF|x):       0.3619   95% ci: (0.2787,0.4543)
  Pr(y=NotInLF|x):    0.6381   95% ci: (0.5457,0.7213)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=inLF|x):       0.2520   95% ci: (0.1416,0.4078)
  Pr(y=NotInLF|x):    0.7480   95% ci: (0.5922,0.8584)

. praccum, using(mage) xis(60) gen(agsq_long_name_here)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
  agsq_x |       9          40   13.69306         20         60  
 agsq_p0 |       9    .4282142   .1752595   .2676314   .7479599  
 agsq_p1 |       9    .5717858   .1752595   .2520402   .7323686  

. list agsq* in 1/9

        agsq_x    agsq_p0    agsq_p1 
  1.        20   .2676314   .7323686  
  2.        25   .2682353   .7317647  
  3.        30   .2836163   .7163837  
  4.        35   .3152536   .6847464  
  5.        40   .3656723   .6343277  
  6.        45   .4373158   .5626842  
  7.        50   .5301194   .4698806  
  8.        55   .6381241   .3618759  
  9.        60   .7479599   .2520402  

. 
. * mlogit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. gen colboth = ((wc==1)*2) + (hc==1)

. recode colboth 0=-6 1=3 2=15 3=59
(753 changes made)

. quietly mlogit colboth lfp k5 k618 age age2 lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.4871
  Pr(y=15|x):         0.0145
  Pr(y=59|x):         0.4077
  Pr(y=-6|x):         0.0906

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.4179
  Pr(y=15|x):         0.0273
  Pr(y=59|x):         0.3493
  Pr(y=-6|x):         0.2055

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.3241
  Pr(y=15|x):         0.0404
  Pr(y=59|x):         0.2787
  Pr(y=-6|x):         0.3568

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.2379
  Pr(y=15|x):         0.0492
  Pr(y=59|x):         0.2167
  Pr(y=-6|x):         0.4962

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.1775
  Pr(y=15|x):         0.0528
  Pr(y=59|x):         0.1763
  Pr(y=-6|x):         0.5935

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.1423
  Pr(y=15|x):         0.0529
  Pr(y=59|x):         0.1588
  Pr(y=-6|x):         0.6461

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.1264
  Pr(y=15|x):         0.0510
  Pr(y=59|x):         0.1631
  Pr(y=-6|x):         0.6595

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.1251
  Pr(y=15|x):         0.0476
  Pr(y=59|x):         0.1923
  Pr(y=-6|x):         0.6351

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=3|x):          0.1355
  Pr(y=15|x):         0.0422
  Pr(y=59|x):         0.2556
  Pr(y=-6|x):         0.5667

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqp3 |       9    .2415206   .1370607   .1250517   .4871354  
 agsqp15 |       9    .0419887   .0130994   .0145019   .0529001  
 agsqp59 |       9    .2442806   .0873539   .1587618   .4077388  
 agsqp_6 |       9      .47221   .2079981   .0906239   .6595201  

. list agsq* in 1/9

         agsqx     agsqp3    agsqp15    agsqp59    agsqp_6 
  1.        20   .4871354   .0145019   .4077388   .0906239  
  2.        25   .4178609   .0273289   .3493013   .2055089  
  3.        30   .3241276   .0404455   .2786702   .3567568  
  4.        35   .2379427   .0491967    .216682   .4961786  
  5.        40   .1774519    .052796   .1762696   .5934825  
  6.        45   .1422765   .0529001   .1587618   .6460616  
  7.        50   .1263636   .0509911   .1631252   .6595201  
  8.        55   .1250517   .0475618    .192329   .6350575  
  9.        60   .1354754   .0421767   .2556478   .5667002  

. 
. * mlogit with base specified
. drop agsq*

. quietly mlogit colboth lfp k5 k618 age age2 lwg inc, b(3)

. prvalue, x(age 20 age2 400) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.0906
  Pr(y=15|x):         0.0145
  Pr(y=59|x):         0.4077
  Pr(y=3|x):          0.4871

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.2055
  Pr(y=15|x):         0.0273
  Pr(y=59|x):         0.3493
  Pr(y=3|x):          0.4179

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.3568
  Pr(y=15|x):         0.0404
  Pr(y=59|x):         0.2787
  Pr(y=3|x):          0.3241

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.4962
  Pr(y=15|x):         0.0492
  Pr(y=59|x):         0.2167
  Pr(y=3|x):          0.2379

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.5935
  Pr(y=15|x):         0.0528
  Pr(y=59|x):         0.1763
  Pr(y=3|x):          0.1775

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.6461
  Pr(y=15|x):         0.0529
  Pr(y=59|x):         0.1588
  Pr(y=3|x):          0.1423

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.6595
  Pr(y=15|x):         0.0510
  Pr(y=59|x):         0.1631
  Pr(y=3|x):          0.1264

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.6351
  Pr(y=15|x):         0.0476
  Pr(y=59|x):         0.1923
  Pr(y=3|x):          0.1251

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

Predicted probabilities for each category:
  Pr(y=-6|x):         0.5667
  Pr(y=15|x):         0.0422
  Pr(y=59|x):         0.2556
  Pr(y=3|x):          0.1355

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
 agsqp_6 |       9      .47221   .2079981   .0906239   .6595201  
 agsqp15 |       9    .0419887   .0130994   .0145019   .0529001  
 agsqp59 |       9    .2442806   .0873539   .1587618   .4077388  
  agsqp3 |       9    .2415206   .1370607   .1250517   .4871354  

. list agsq* in 1/9

         agsqx    agsqp_6    agsqp15    agsqp59     agsqp3 
  1.        20   .0906239   .0145019   .4077388   .4871354  
  2.        25   .2055089   .0273289   .3493013   .4178609  
  3.        30   .3567568   .0404455   .2786702   .3241276  
  4.        35   .4961786   .0491967    .216682   .2379427  
  5.        40   .5934825    .052796   .1762696   .1774519  
  6.        45   .6460616   .0529001   .1587618   .1422765  
  7.        50   .6595201   .0509911   .1631252   .1263636  
  8.        55   .6350575   .0475618    .192329   .1250517  
  9.        60   .5667002   .0421767   .2556478   .1354754  

. 
. * nbreg
. use couart2, clear
(Academic Biochemists / S Long)

. gen phd2 = phd*phd

. quietly nbreg art fem mar kid5 phd phd2 ment

. prvalue, x(phd 1 phd2 1) rest(mean) brief

Predicted rate: 1.3502 
Predicted probabilities:

  Pr(y=0|x):   0.3465  Pr(y=1|x):   0.2935
  Pr(y=2|x):   0.1790  Pr(y=3|x):   0.0950
  Pr(y=4|x):   0.0467  Pr(y=5|x):   0.0218
  Pr(y=6|x):   0.0099  Pr(y=7|x):   0.0043
  Pr(y=8|x):   0.0019  Pr(y=9|x):   0.0008

. praccum, saving(mphd) xis(1)

. prvalue, x(phd 1.5 phd2 2.25) rest(mean) brief

Predicted rate: 1.4681 
Predicted probabilities:

  Pr(y=0|x):   0.3222  Pr(y=1|x):   0.2874
  Pr(y=2|x):   0.1846  Pr(y=3|x):   0.1032
  Pr(y=4|x):   0.0534  Pr(y=5|x):   0.0263
  Pr(y=6|x):   0.0125  Pr(y=7|x):   0.0058
  Pr(y=8|x):   0.0026  Pr(y=9|x):   0.0012

. praccum, using(mphd) xis(1.5)

. prvalue, x(phd 2 phd2 4) rest(mean) brief

Predicted rate: 1.5632 
Predicted probabilities:

  Pr(y=0|x):   0.3043  Pr(y=1|x):   0.2819
  Pr(y=2|x):   0.1880  Pr(y=3|x):   0.1091
  Pr(y=4|x):   0.0586  Pr(y=5|x):   0.0300
  Pr(y=6|x):   0.0148  Pr(y=7|x):   0.0071
  Pr(y=8|x):   0.0034  Pr(y=9|x):   0.0016

. praccum, using(mphd) xis(2)

. prvalue, x(phd 2.5 phd2 6.25) rest(mean) brief

Predicted rate: 1.6299 
Predicted probabilities:

  Pr(y=0|x):   0.2926  Pr(y=1|x):   0.2778
  Pr(y=2|x):   0.1899  Pr(y=3|x):   0.1129
  Pr(y=4|x):   0.0622  Pr(y=5|x):   0.0326
  Pr(y=6|x):   0.0165  Pr(y=7|x):   0.0081
  Pr(y=8|x):   0.0039  Pr(y=9|x):   0.0019

. praccum, using(mphd) xis(2.5)

. prvalue, x(phd 3 phd2 9) rest(mean) brief

Predicted rate: 1.6642 
Predicted probabilities:

  Pr(y=0|x):   0.2868  Pr(y=1|x):   0.2756
  Pr(y=2|x):   0.1907  Pr(y=3|x):   0.1148
  Pr(y=4|x):   0.0640  Pr(y=5|x):   0.0339
  Pr(y=6|x):   0.0174  Pr(y=7|x):   0.0087
  Pr(y=8|x):   0.0043  Pr(y=9|x):   0.0021

. praccum, using(mphd) xis(3)

. prvalue, x(phd 3.5 phd2 12.25) rest(mean) brief

Predicted rate: 1.6639 
Predicted probabilities:

  Pr(y=0|x):   0.2869  Pr(y=1|x):   0.2757
  Pr(y=2|x):   0.1907  Pr(y=3|x):   0.1148
  Pr(y=4|x):   0.0639  Pr(y=5|x):   0.0339
  Pr(y=6|x):   0.0174  Pr(y=7|x):   0.0087
  Pr(y=8|x):   0.0042  Pr(y=9|x):   0.0020

. praccum, using(mphd) xis(3.5)

. prvalue, x(phd 4 phd2 16) rest(mean) brief

Predicted rate: 1.6291 
Predicted probabilities:

  Pr(y=0|x):   0.2927  Pr(y=1|x):   0.2778
  Pr(y=2|x):   0.1898  Pr(y=3|x):   0.1129
  Pr(y=4|x):   0.0621  Pr(y=5|x):   0.0325
  Pr(y=6|x):   0.0165  Pr(y=7|x):   0.0081
  Pr(y=8|x):   0.0039  Pr(y=9|x):   0.0019

. praccum, using(mphd) xis(4) gen(phsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   phsqx |       7         2.5   1.080123          1          4  
  phsqmu |       7    1.566925   .1179342   1.350161   1.664171  
  phsqp0 |       7    .3045746   .0222781   .2868401   .3464993  
  phsqp1 |       7    .2813983   .0067789   .2756362    .293544  
  phsqp2 |       7    .1875197   .0043167   .1790195   .1906729  
  phsqp3 |       7     .108953   .0073824   .0950148   .1147904  
  phsqp4 |       7    .0586899   .0064806   .0466711    .063957  
  phsqp5 |       7    .0301376   .0045512   .0218172   .0339129  
  phsqp6 |       7    .0149758   .0028394   .0098537   .0173736  
  phsqp7 |       7    .0072664   .0016439    .004339   .0086778  
  phsqp8 |       7    .0034628   .0009037   .0018739    .004251  
  phsqp9 |       7    .0016272   .0004783   .0007969   .0020506  
  phsqs0 |       7    .3045746   .0222781   .2868401   .3464993  
  phsqs1 |       7    .5859729   .0290436   .5624763   .6400433  
  phsqs2 |       7    .7734926   .0247501   .7531492   .8190628  
  phsqs3 |       7    .8824456   .0173757   .8679396   .9140776  
  phsqs4 |       7    .9411355   .0108991   .9318966   .9607487  
  phsqs5 |       7    .9712731   .0063499   .9658095   .9825659  
  phsqs6 |       7    .9862489   .0035116   .9831831   .9924196  
  phsqs7 |       7    .9935154   .0018683   .9918609   .9967586  
  phsqs8 |       7    .9969782   .0009648   .9961118   .9986326  
  phsqs9 |       7    .9986054   .0004866   .9981624   .9994295  

. list phsq* in 1/7

Observation 1

       phsqx            1      phsqmu     1.350161      phsqp0     .3464993
      phsqp1      .293544      phsqp2     .1790195      phsqp3     .0950148
      phsqp4     .0466711      phsqp5     .0218172      phsqp6     .0098537
      phsqp7      .004339      phsqp8     .0018739      phsqp9     .0007969
      phsqs0     .3464993      phsqs1     .6400433      phsqs2     .8190628
      phsqs3     .9140776      phsqs4     .9607487      phsqs5     .9825659
      phsqs6     .9924196      phsqs7     .9967586      phsqs8     .9986326
      phsqs9     .9994295


Observation 2

       phsqx          1.5      phsqmu     1.468101      phsqp0     .3221629
      phsqp1     .2874145      phsqp2     .1845858      phsqp3     .1031695
      phsqp4     .0533668      phsqp5     .0262715      phsqp6     .0124952
      phsqp7     .0057943      phsqp8     .0026352      phsqp9     .0011802
      phsqs0     .3221629      phsqs1     .6095774      phsqs2     .7941632
      phsqs3     .8973327      phsqs4     .9506995      phsqs5      .976971
      phsqs6     .9894662      phsqs7     .9952605      phsqs8     .9978957
      phsqs9     .9990759


Observation 3

       phsqx            2      phsqmu       1.5632      phsqp0     .3042924
      phsqp1     .2818929      phsqp2     .1879895      phsqp3     .1091055
      phsqp4     .0586038      phsqp5     .0299571      phsqp6     .0147952
      phsqp7     .0071242      phsqp8     .0033644      phsqp9     .0015646
      phsqs0     .3042924      phsqs1     .5861852      phsqs2     .7741747
      phsqs3     .8832802      phsqs4      .941884      phsqs5      .971841
      phsqs6     .9866362      phsqs7     .9937604      phsqs8     .9971249
      phsqs9     .9986895


Observation 4

       phsqx          2.5      phsqmu     1.629904      phsqp0     .2925977
      phsqp1     .2777964      phsqp2     .1898624      phsqp3     .1129314
      phsqp4     .0621666      phsqp5     .0325681      phsqp6     .0164845
      phsqp7     .0081349      phsqp8     .0039372      phsqp9     .0018765
      phsqs0     .2925977      phsqs1     .5703942      phsqs2     .7602565
      phsqs3     .8731879      phsqs4     .9353545      phsqs5     .9679226
      phsqs6     .9844071      phsqs7     .9925421      phsqs8     .9964793
      phsqs9     .9983559


Observation 5

       phsqx            3      phsqmu     1.664171      phsqp0     .2868401
      phsqp1     .2756362      phsqp2     .1906729      phsqp3     .1147904
      phsqp4      .063957      phsqp5     .0339129      phsqp6     .0173736
      phsqp7     .0086778      phsqp8      .004251      phsqp9     .0020506
      phsqs0     .2868401      phsqs1     .5624763      phsqs2     .7531492
      phsqs3     .8679396      phsqs4     .9318966      phsqs5     .9658095
      phsqs6     .9831831      phsqs7     .9918609      phsqs8     .9961118
      phsqs9     .9981624


Observation 6

       phsqx          3.5      phsqmu     1.663883      phsqp0     .2868879
      phsqp1     .2756545      phsqp2     .1906665      phsqp3      .114775
      phsqp4     .0639421      phsqp5     .0339016      phsqp6      .017366
      phsqp7     .0086731      phsqp8     .0042483      phsqp9     .0020491
      phsqs0     .2868879      phsqs1     .5625424      phsqs2      .753209
      phsqs3      .867984      phsqs4     .9319261      phsqs5     .9658276
      phsqs6     .9831936      phsqs7     .9918668      phsqs8      .996115
      phsqs9     .9981641


Observation 7

       phsqx            4      phsqmu     1.629056      phsqp0     .2927422
      phsqp1     .2778494      phsqp2      .189841      phsqp3     .1128845
      phsqp4     .0621219      phsqp5     .0325349      phsqp6     .0164627
      phsqp7     .0081217      phsqp8     .0039297      phsqp9     .0018723
      phsqs0     .2927422      phsqs1     .5705917      phsqs2     .7604327
      phsqs3     .8733172      phsqs4     .9354392      phsqs5     .9679741
      phsqs6     .9844368      phsqs7     .9925585      phsqs8     .9964882
      phsqs9     .9983605


. 
. * ologit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. gen colboth = ((wc==1)*2) + (hc==1)

. recode colboth 0=-6 1=3 2=15 3=59
(753 changes made)

. quietly ologit colboth lfp k5 k618 age age2 lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=-6|x):         0.1355
  Pr(y=3|x):          0.1444
  Pr(y=15|x):         0.0782
  Pr(y=59|x):         0.6418

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=-6|x):         0.2508
  Pr(y=3|x):          0.2028
  Pr(y=15|x):         0.0901
  Pr(y=59|x):         0.4564

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=-6|x):         0.3815
  Pr(y=3|x):          0.2232
  Pr(y=15|x):         0.0824
  Pr(y=59|x):         0.3130

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=-6|x):         0.4952
  Pr(y=3|x):          0.2135
  Pr(y=15|x):         0.0687
  Pr(y=59|x):         0.2226

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=-6|x):         0.5740
  Pr(y=3|x):          0.1957
  Pr(y=15|x):         0.0578
  Pr(y=59|x):         0.1726

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=-6|x):         0.6150
  Pr(y=3|x):          0.1834
  Pr(y=15|x):         0.0520
  Pr(y=59|x):         0.1496

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=-6|x):         0.6204
  Pr(y=3|x):          0.1817
  Pr(y=15|x):         0.0512
  Pr(y=59|x):         0.1467

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=-6|x):         0.5908
  Pr(y=3|x):          0.1909
  Pr(y=15|x):         0.0554
  Pr(y=59|x):         0.1629

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=-6|x):         0.5241
  Pr(y=3|x):          0.2079
  Pr(y=15|x):         0.0648
  Pr(y=59|x):         0.2033

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
 agsqp_6 |       9    .4652543   .1732685   .1355445   .6204205  
  agsqp3 |       9     .193701   .0230012   .1444141   .2231658  
 agsqp15 |       9    .0667188   .0141198   .0512065   .0900796  
 agsqp59 |       9    .2743259   .1703102   .1466944   .6418417  
 agsqs_6 |       9    .4652543   .1732685   .1355445   .6204205  
  agsqs3 |       9    .6589552   .1816045   .2799585   .8020991  
 agsqs15 |       9    .7256741   .1703102   .3581583   .8533056  
 agsqs59 |       9           1          0          1          1  

. list agsq* in 1/9

Observation 1

       agsqx           20     agsqp_6     .1355445      agsqp3     .1444141
     agsqp15     .0781998     agsqp59     .6418417     agsqs_6     .1355445
      agsqs3     .2799585     agsqs15     .3581583     agsqs59            1


Observation 2

       agsqx           25     agsqp_6     .2507621      agsqp3     .2027676
     agsqp15     .0900796     agsqp59     .4563907     agsqs_6     .2507621
      agsqs3     .4535297     agsqs15     .5436093     agsqs59            1


Observation 3

       agsqx           30     agsqp_6     .3814667      agsqp3     .2231658
     agsqp15     .0823621     agsqp59     .3130054     agsqs_6     .3814667
      agsqs3     .6046325     agsqs15     .6869946     agsqs59            1


Observation 4

       agsqx           35     agsqp_6     .4952184      agsqp3     .2134665
     agsqp15      .068668     agsqp59     .2226471     agsqs_6     .4952184
      agsqs3     .7086849     agsqs15     .7773529     agsqs59            1


Observation 5

       agsqx           40     agsqp_6     .5739658      agsqp3      .195657
     agsqp15     .0578017     agsqp59     .1725755     agsqs_6     .5739658
      agsqs3     .7696228     agsqs15     .8274245     agsqs59            1


Observation 6

       agsqx           45     agsqp_6     .6149585      agsqp3     .1834435
     agsqp15     .0519847     agsqp59     .1496134     agsqs_6     .6149585
      agsqs3      .798402     agsqs15     .8503867     agsqs59            1


Observation 7

       agsqx           50     agsqp_6     .6204205      agsqp3     .1816786
     agsqp15     .0512065     agsqp59     .1466944     agsqs_6     .6204205
      agsqs3     .8020991     agsqs15     .8533056     agsqs59            1


Observation 8

       agsqx           55     agsqp_6     .5908401      agsqp3      .190855
     agsqp15     .0554139     agsqp59     .1628911     agsqs_6     .5908401
      agsqs3     .7816951     agsqs15      .837109     agsqs59            1


Observation 9

       agsqx           60     agsqp_6     .5241117      agsqp3     .2078609
     agsqp15     .0647532     agsqp59     .2032741     agsqs_6     .5241117
      agsqs3     .7319727     agsqs15     .7967259     agsqs59            1


. 
. * oprobit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. gen colboth = ((wc==1)*2) + (hc==1)

. recode colboth 0=-6 1=3 2=15 3=59
(753 changes made)

. quietly oprobit colboth lfp k5 k618 age age2 lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=-6|x):         0.1408
  Pr(y=3|x):          0.1532
  Pr(y=15|x):         0.0751
  Pr(y=59|x):         0.6309

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=-6|x):         0.2647
  Pr(y=3|x):          0.1979
  Pr(y=15|x):         0.0826
  Pr(y=59|x):         0.4547

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=-6|x):         0.3943
  Pr(y=3|x):          0.2109
  Pr(y=15|x):         0.0771
  Pr(y=59|x):         0.3176

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=-6|x):         0.5022
  Pr(y=3|x):          0.2034
  Pr(y=15|x):         0.0672
  Pr(y=59|x):         0.2272

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=-6|x):         0.5761
  Pr(y=3|x):          0.1903
  Pr(y=15|x):         0.0586
  Pr(y=59|x):         0.1751

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=-6|x):         0.6144
  Pr(y=3|x):          0.1812
  Pr(y=15|x):         0.0537
  Pr(y=59|x):         0.1507

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=-6|x):         0.6189
  Pr(y=3|x):          0.1800
  Pr(y=15|x):         0.0531
  Pr(y=59|x):         0.1480

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=-6|x):         0.5898
  Pr(y=3|x):          0.1872
  Pr(y=15|x):         0.0569
  Pr(y=59|x):         0.1661

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=-6|x):         0.5256
  Pr(y=3|x):          0.1999
  Pr(y=15|x):         0.0646
  Pr(y=59|x):         0.2099

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
 agsqp_6 |       9    .4696523   .1691886   .1407753   .6188918  
  agsqp3 |       9    .1893253   .0170091   .1532007   .2109297  
 agsqp15 |       9    .0654346   .0108294   .0531302   .0826132  
 agsqp59 |       9    .2755878   .1661103   .1479953     .63094  
 agsqs_6 |       9    .4696523   .1691886   .1407753   .6188918  
  agsqs3 |       9    .6589776   .1749774    .293976   .7988745  
 agsqs15 |       9    .7244122   .1661103     .36906   .8520047  
 agsqs59 |       9           1          0          1          1  

. list agsq* in 1/9

Observation 1

       agsqx           20     agsqp_6     .1407753      agsqp3     .1532007
     agsqp15      .075084     agsqp59       .63094     agsqs_6     .1407753
      agsqs3      .293976     agsqs15       .36906     agsqs59            1


Observation 2

       agsqx           25     agsqp_6      .264742      agsqp3     .1979088
     agsqp15     .0826132     agsqp59     .4547359     agsqs_6      .264742
      agsqs3     .4626509     agsqs15     .5452641     agsqs59            1


Observation 3

       agsqx           30     agsqp_6     .3943483      agsqp3     .2109297
     agsqp15     .0771485     agsqp59     .3175735     agsqs_6     .3943483
      agsqs3      .605278     agsqs15     .6824265     agsqs59            1


Observation 4

       agsqx           35     agsqp_6     .5022054      agsqp3     .2033907
     agsqp15     .0671819     agsqp59     .2272221     agsqs_6     .5022054
      agsqs3      .705596     agsqs15     .7727779     agsqs59            1


Observation 5

       agsqx           40     agsqp_6      .576051      agsqp3     .1902862
     agsqp15     .0585881     agsqp59     .1750747     agsqs_6      .576051
      agsqs3     .7663372     agsqs15     .8249253     agsqs59            1


Observation 6

       agsqx           45     agsqp_6     .6144038      agsqp3     .1811521
     agsqp15     .0537156     agsqp59     .1507284     agsqs_6     .6144038
      agsqs3     .7955559     agsqs15     .8492715     agsqs59            1


Observation 7

       agsqx           50     agsqp_6     .6188918      agsqp3     .1799827
     agsqp15     .0531302     agsqp59     .1479953     agsqs_6     .6188918
      agsqs3     .7988745     agsqs15     .8520047     agsqs59            1


Observation 8

       agsqx           55     agsqp_6      .589815      agsqp3     .1871864
     agsqp15     .0568675     agsqp59      .166131     agsqs_6      .589815
      agsqs3     .7770014     agsqs15      .833869     agsqs59            1


Observation 9

       agsqx           60     agsqp_6     .5256376      agsqp3     .1998908
     agsqp15     .0645827     agsqp59     .2098889     agsqs_6     .5256376
      agsqs3     .7255284     agsqs15     .7901111     agsqs59            1


. 
. * poisson
. use couart2, clear
(Academic Biochemists / S Long)

. gen phd2 = phd*phd

. quietly poisson art fem mar kid5 phd phd2 ment

. prvalue, x(phd 1 phd2 1) rest(mean) brief

Predicted rate: 1.3507 
Predicted probabilities:

  Pr(y=0|x):   0.2590  Pr(y=1|x):   0.3499
  Pr(y=2|x):   0.2363  Pr(y=3|x):   0.1064
  Pr(y=4|x):   0.0359  Pr(y=5|x):   0.0097
  Pr(y=6|x):   0.0022  Pr(y=7|x):   0.0004
  Pr(y=8|x):   0.0001  Pr(y=9|x):   0.0000

. praccum, saving(mphd) xis(1)

. prvalue, x(phd 1.5 phd2 2.25) rest(mean) brief

Predicted rate: 1.4733 
Predicted probabilities:

  Pr(y=0|x):   0.2292  Pr(y=1|x):   0.3376
  Pr(y=2|x):   0.2487  Pr(y=3|x):   0.1221
  Pr(y=4|x):   0.0450  Pr(y=5|x):   0.0133
  Pr(y=6|x):   0.0033  Pr(y=7|x):   0.0007
  Pr(y=8|x):   0.0001  Pr(y=9|x):   0.0000

. praccum, using(mphd) xis(1.5)

. prvalue, x(phd 2 phd2 4) rest(mean) brief

Predicted rate: 1.5721 
Predicted probabilities:

  Pr(y=0|x):   0.2076  Pr(y=1|x):   0.3264
  Pr(y=2|x):   0.2565  Pr(y=3|x):   0.1344
  Pr(y=4|x):   0.0528  Pr(y=5|x):   0.0166
  Pr(y=6|x):   0.0044  Pr(y=7|x):   0.0010
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

. praccum, using(mphd) xis(2)

. prvalue, x(phd 2.5 phd2 6.25) rest(mean) brief

Predicted rate: 1.641  
Predicted probabilities:

  Pr(y=0|x):   0.1938  Pr(y=1|x):   0.3180
  Pr(y=2|x):   0.2609  Pr(y=3|x):   0.1427
  Pr(y=4|x):   0.0585  Pr(y=5|x):   0.0192
  Pr(y=6|x):   0.0053  Pr(y=7|x):   0.0012
  Pr(y=8|x):   0.0003  Pr(y=9|x):   0.0000

. praccum, using(mphd) xis(2.5)

. prvalue, x(phd 3 phd2 9) rest(mean) brief

Predicted rate: 1.6757 
Predicted probabilities:

  Pr(y=0|x):   0.1872  Pr(y=1|x):   0.3137
  Pr(y=2|x):   0.2628  Pr(y=3|x):   0.1468
  Pr(y=4|x):   0.0615  Pr(y=5|x):   0.0206
  Pr(y=6|x):   0.0058  Pr(y=7|x):   0.0014
  Pr(y=8|x):   0.0003  Pr(y=9|x):   0.0001

. praccum, using(mphd) xis(3)

. prvalue, x(phd 3.5 phd2 12.25) rest(mean) brief

Predicted rate: 1.6739 
Predicted probabilities:

  Pr(y=0|x):   0.1875  Pr(y=1|x):   0.3139
  Pr(y=2|x):   0.2627  Pr(y=3|x):   0.1466
  Pr(y=4|x):   0.0613  Pr(y=5|x):   0.0205
  Pr(y=6|x):   0.0057  Pr(y=7|x):   0.0014
  Pr(y=8|x):   0.0003  Pr(y=9|x):   0.0001

. praccum, using(mphd) xis(3.5)

. prvalue, x(phd 4 phd2 16) rest(mean) brief

Predicted rate: 1.6358 
Predicted probabilities:

  Pr(y=0|x):   0.1948  Pr(y=1|x):   0.3186
  Pr(y=2|x):   0.2606  Pr(y=3|x):   0.1421
  Pr(y=4|x):   0.0581  Pr(y=5|x):   0.0190
  Pr(y=6|x):   0.0052  Pr(y=7|x):   0.0012
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

. praccum, using(mphd) xis(4) gen(phsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   phsqx |       7         2.5   1.080123          1          4  
  phsqmu |       7    1.574642   .1215293   1.350741   1.675674  
  phsqp0 |       7    .2084441   .0267157   .1871819   .2590483  
  phsqp1 |       7    .3254445    .013676   .3136559   .3499071  
  phsqp2 |       7    .2555176   .0097962   .2363169   .2627926  
  phsqp3 |       7    .1344546   .0150809   .1064009   .1467849  
  phsqp4 |       7    .0533221    .009605     .03593   .0614909  
  phsqp5 |       7    .0169927     .00412   .0097064   .0206078  
  phsqp6 |       7     .004531    .001357   .0021851   .0057553  
  phsqp7 |       7    .0010394   .0003655   .0004217   .0013777  
  phsqp8 |       7    .0002093   .0000836   .0000712   .0002886  
  phsqp9 |       7    .0000376   .0000166   .0000107   .0000537  
  phsqs0 |       7    .2084441   .0267157   .1871819   .2590483  
  phsqs1 |       7    .5338887   .0403593   .5008378   .6089554  
  phsqs2 |       7    .7894063   .0306144   .7636304   .8452723  
  phsqs3 |       7    .9238608   .0155462   .9104154   .9516732  
  phsqs4 |       7     .977183   .0059449   .9719063   .9876032  
  phsqs5 |       7    .9941757   .0018259   .9925141   .9973097  
  phsqs6 |       7    .9987067   .0004691   .9982694   .9994949  
  phsqs7 |       7     .999746   .0001037   .9996471   .9999165  
  phsqs8 |       7    .9999553   .0000201   .9999357   .9999877  
  phsqs9 |       7    .9999929   3.46e-06   .9999894   .9999983  

. list phsq* in 1/7

Observation 1

       phsqx            1      phsqmu     1.350741      phsqp0     .2590483
      phsqp1     .3499071      phsqp2     .2363169      phsqp3     .1064009
      phsqp4       .03593      phsqp5     .0097064      phsqp6     .0021851
      phsqp7     .0004217      phsqp8     .0000712      phsqp9     .0000107
      phsqs0     .2590483      phsqs1     .6089554      phsqs2     .8452723
      phsqs3     .9516732      phsqs4     .9876032      phsqs5     .9973097
      phsqs6     .9994949      phsqs7     .9999165      phsqs8     .9999877
      phsqs9     .9999983


Observation 2

       phsqx          1.5      phsqmu     1.473301      phsqp0     .2291678
      phsqp1     .3376331      phsqp2     .2487176      phsqp3     .1221453
      phsqp4     .0449892      phsqp5     .0132565      phsqp6     .0032551
      phsqp7     .0006851      phsqp8     .0001262      phsqp9     .0000207
      phsqs0     .2291678      phsqs1     .5668008      phsqs2     .8155184
      phsqs3     .9376637      phsqs4     .9826528      phsqs5     .9959094
      phsqs6     .9991645      phsqs7     .9998496      phsqs8     .9999758
      phsqs9     .9999965


Observation 3

       phsqx            2      phsqmu     1.572054      phsqp0     .2076182
      phsqp1     .3263871      phsqp2     .2565491      phsqp3     .1344364
      phsqp4     .0528353      phsqp5      .016612      phsqp6     .0043525
      phsqp7     .0009775      phsqp8     .0001921      phsqp9     .0000336
      phsqs0     .2076182      phsqs1     .5340053      phsqs2     .7905545
      phsqs3      .924991      phsqs4     .9778263      phsqs5     .9944383
      phsqs6     .9987908      phsqs7     .9997683      phsqs8     .9999604
      phsqs9     .9999939


Observation 4

       phsqx          2.5      phsqmu     1.640969      phsqp0     .1937922
      phsqp1      .318007      phsqp2     .2609197      phsqp3     .1427204
      phsqp4     .0585499      phsqp5     .0192157      phsqp6     .0052554
      phsqp7      .001232      phsqp8     .0002527      phsqp9     .0000461
      phsqs0     .1937922      phsqs1     .5117992      phsqs2     .7727189
      phsqs3     .9154393      phsqs4     .9739892      phsqs5      .993205
      phsqs6     .9984604      phsqs7     .9996923      phsqs8      .999945
      phsqs9     .9999911


Observation 5

       phsqx            3      phsqmu     1.675674      phsqp0     .1871819
      phsqp1     .3136559      phsqp2     .2627926      phsqp3     .1467849
      phsqp4     .0614909      phsqp5     .0206078      phsqp6     .0057553
      phsqp7     .0013777      phsqp8     .0002886      phsqp9     .0000537
      phsqs0     .1871819      phsqs1     .5008378      phsqs2     .7636304
      phsqs3     .9104154      phsqs4     .9719063      phsqs5     .9925141
      phsqs6     .9982694      phsqs7     .9996471      phsqs8     .9999357
      phsqs9     .9999894


Observation 6

       phsqx          3.5      phsqmu     1.673923      phsqp0     .1875099
      phsqp1     .3138773      phsqp2     .2627033      phsqp3     .1465817
      phsqp4     .0613416      phsqp5     .0205362      phsqp6     .0057294
      phsqp7     .0013701      phsqp8     .0002867      phsqp9     .0000533
      phsqs0     .1875099      phsqs1     .5013872      phsqs2     .7640905
      phsqs3     .9106722      phsqs4     .9720139      phsqs5     .9925501
      phsqs6     .9982795      phsqs7     .9996496      phsqs8     .9999363
      phsqs9     .9999896


Observation 7

       phsqx            4      phsqmu      1.63583      phsqp0     .1947906
      phsqp1     .3186443      phsqp2      .260624      phsqp3     .1421122
      phsqp4     .0581178      phsqp5     .0190142      phsqp6      .005184
      phsqp7     .0012114      phsqp8     .0002477      phsqp9      .000045
      phsqs0     .1947906      phsqs1     .5134349      phsqs2     .7740589
      phsqs3     .9161711      phsqs4      .974289      phsqs5     .9933032
      phsqs6     .9984872      phsqs7     .9996986      phsqs8     .9999464
      phsqs9     .9999914


. 
. * probit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. quietly probit lfp k5 k618 age age2 wc hc lwg inc

. prvalue, x(age 20 age2 400) rest(mean) brief

  Pr(y=inLF|x):       0.7198   95% ci: (0.3937,0.9242)
  Pr(y=NotInLF|x):    0.2802   95% ci: (0.0758,0.6063)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

  Pr(y=inLF|x):       0.7242   95% ci: (0.5294,0.8681)
  Pr(y=NotInLF|x):    0.2758   95% ci: (0.1319,0.4706)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

  Pr(y=inLF|x):       0.7124   95% ci: (0.6099,0.8000)
  Pr(y=NotInLF|x):    0.2876   95% ci: (0.2000,0.3901)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

  Pr(y=inLF|x):       0.6833   95% ci: (0.6278,0.7349)
  Pr(y=NotInLF|x):    0.3167   95% ci: (0.2651,0.3722)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

  Pr(y=inLF|x):       0.6349   95% ci: (0.5829,0.6846)
  Pr(y=NotInLF|x):    0.3651   95% ci: (0.3154,0.4171)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

  Pr(y=inLF|x):       0.5654   95% ci: (0.5088,0.6206)
  Pr(y=NotInLF|x):    0.4346   95% ci: (0.3794,0.4912)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

  Pr(y=inLF|x):       0.4744   95% ci: (0.4164,0.5330)
  Pr(y=NotInLF|x):    0.5256   95% ci: (0.4670,0.5836)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

  Pr(y=inLF|x):       0.3665   95% ci: (0.2828,0.4570)
  Pr(y=NotInLF|x):    0.6335   95% ci: (0.5430,0.7172)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

  Pr(y=inLF|x):       0.2525   95% ci: (0.1352,0.4086)
  Pr(y=NotInLF|x):    0.7475   95% ci: (0.5914,0.8648)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqp0 |       9    .4296299   .1713803   .2757552   .7475491  
  agsqp1 |       9    .5703701   .1713803   .2524509   .7242448  

. list agsq* in 1/9

         agsqx     agsqp0     agsqp1 
  1.        20   .2802111   .7197889  
  2.        25   .2757552   .7242448  
  3.        30   .2876104   .7123896  
  4.        35   .3167318   .6832682  
  5.        40   .3650725   .6349275  
  6.        45   .4346277   .5653723  
  7.        50   .5255651   .4744348  
  8.        55   .6335461   .3664539  
  9.        60   .7475491   .2524509  

. 
. * regress
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. quietly regress lwg lfp k5 k618 age age2 wc hc inc

. prvalue, x(age 20 age2 400) rest(mean) brief

Predicted value of y:  .8231874   95% ci: ( .4817171, 1.164658)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

Predicted value of y:  .9444012   95% ci: ( .7366136, 1.152189)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

Predicted value of y:  1.036023   95% ci: ( .9249351,  1.14711)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

Predicted value of y:  1.098052   95% ci: ( 1.037786, 1.158319)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

Predicted value of y:   1.13049   95% ci: ( 1.074569,  1.18641)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

Predicted value of y:  1.133335   95% ci: ( 1.074658, 1.192012)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

Predicted value of y:  1.106588   95% ci: ( 1.045447, 1.167729)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

Predicted value of y:  1.050249   95% ci: ( .9516996, 1.148798)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

Predicted value of y:  .9643174   95% ci: ( .7810022, 1.147632)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqxb |       9    1.031849   .1037433   .8231874   1.133335  

. list agsq* in 1/9

         agsqx     agsqxb 
  1.        20   .8231874  
  2.        25   .9444012  
  3.        30   1.036023  
  4.        35   1.098052  
  5.        40    1.13049  
  6.        45   1.133335  
  7.        50   1.106588  
  8.        55   1.050249  
  9.        60   .9643174  

. 
. * tobit
. use binlfp2, clear
(PSID 1976 / T Mroz)

. gen age2 = age*age

. gen lwgcens = lwg

. replace lwgcens = 0 if lwg < 0
(18 real changes made)

. replace lwgcens = 2 if lwg > 2
(45 real changes made)

. gen censor = 0

. replace censor = -1 if lwgcens == 0
(22 real changes made)

. replace censor = 1 if lwgcens == 2
(45 real changes made)

. quietly tobit lwgcens lfp k5 k618 age age2 wc hc inc, ll(0) ul(2)

. prvalue, x(age 20 age2 400) rest(mean) brief

Predicted value of y*:  .7988954   95% ci: ( .5062962, 1.091494)

. praccum, saving(mage) xis(20)

. prvalue, x(age 25 age2 625) rest(mean) brief

Predicted value of y*:  .9361635   95% ci: ( .7580404, 1.114287)

. praccum, using(mage) xis(25)

. prvalue, x(age 30 age2 900) rest(mean) brief

Predicted value of y*:  1.039235   95% ci: ( .9439359, 1.134533)

. praccum, using(mage) xis(30)

. prvalue, x(age 35 age2 1225) rest(mean) brief

Predicted value of y*:  1.108108   95% ci: ( 1.056405, 1.159812)

. praccum, using(mage) xis(35)

. prvalue, x(age 40 age2 1600) rest(mean) brief

Predicted value of y*:  1.142785   95% ci: ( 1.094899, 1.190671)

. praccum, using(mage) xis(40)

. prvalue, x(age 45 age2 2025) rest(mean) brief

Predicted value of y*:  1.143265   95% ci: ( 1.093008, 1.193521)

. praccum, using(mage) xis(45)

. prvalue, x(age 50 age2 2500) rest(mean) brief

Predicted value of y*:  1.109547   95% ci: ( 1.057112, 1.161982)

. praccum, using(mage) xis(50)

. prvalue, x(age 55 age2 3025) rest(mean) brief

Predicted value of y*:  1.041633   95% ci: ( .9571576, 1.126107)

. praccum, using(mage) xis(55)

. prvalue, x(age 60 age2 3600) rest(mean) brief

Predicted value of y*:  .9395207   95% ci: ( .7825143, 1.096527)

. praccum, using(mage) xis(60) gen(agsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   agsqx |       9          40   13.69306         20         60  
  agsqxb |       9    1.028795   .1165047   .7988954   1.143265  

. list agsq* in 1/9

         agsqx     agsqxb 
  1.        20   .7988954  
  2.        25   .9361635  
  3.        30   1.039235  
  4.        35   1.108108  
  5.        40   1.142785  
  6.        45   1.143265  
  7.        50   1.109547  
  8.        55   1.041633  
  9.        60   .9395207  

. 
. * zip
. use couart2, clear
(Academic Biochemists / S Long)

. gen phd2 = phd*phd

. quietly zip art fem mar kid5 phd phd2 ment, inf(fem phd phd2)

. prvalue, x(phd 1 phd2 1) rest(mean) brief

Predicted rate: 1.375  
Predicted probabilities:

  Pr(y=0|x,z): 0.3679  Pr(y=1|x):   0.2941
  Pr(y=2|x):   0.2022  Pr(y=3|x):   0.0927
  Pr(y=4|x):   0.0319  Pr(y=5|x):   0.0088
  Pr(y=6|x):   0.0020  Pr(y=7|x):   0.0004
  Pr(y=8|x):   0.0001  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1540

. praccum, saving(mphd) xis(1)

. prvalue, x(phd 1.5 phd2 2.25) rest(mean) brief

Predicted rate: 1.4719 
Predicted probabilities:

  Pr(y=0|x,z): 0.3697  Pr(y=1|x):   0.2763
  Pr(y=2|x):   0.2034  Pr(y=3|x):   0.0998
  Pr(y=4|x):   0.0367  Pr(y=5|x):   0.0108
  Pr(y=6|x):   0.0027  Pr(y=7|x):   0.0006
  Pr(y=8|x):   0.0001  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1820

. praccum, using(mphd) xis(1.5)

. prvalue, x(phd 2 phd2 4) rest(mean) brief

Predicted rate: 1.5499 
Predicted probabilities:

  Pr(y=0|x,z): 0.3682  Pr(y=1|x):   0.2639
  Pr(y=2|x):   0.2045  Pr(y=3|x):   0.1057
  Pr(y=4|x):   0.0409  Pr(y=5|x):   0.0127
  Pr(y=6|x):   0.0033  Pr(y=7|x):   0.0007
  Pr(y=8|x):   0.0001  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1979

. praccum, using(mphd) xis(2)

. prvalue, x(phd 2.5 phd2 6.25) rest(mean) brief

Predicted rate: 1.6125 
Predicted probabilities:

  Pr(y=0|x,z): 0.3588  Pr(y=1|x):   0.2575
  Pr(y=2|x):   0.2076  Pr(y=3|x):   0.1116
  Pr(y=4|x):   0.0450  Pr(y=5|x):   0.0145
  Pr(y=6|x):   0.0039  Pr(y=7|x):   0.0009
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1991

. praccum, using(mphd) xis(2.5)

. prvalue, x(phd 3 phd2 9) rest(mean) brief

Predicted rate: 1.6581 
Predicted probabilities:

  Pr(y=0|x,z): 0.3405  Pr(y=1|x):   0.2574
  Pr(y=2|x):   0.2134  Pr(y=3|x):   0.1179
  Pr(y=4|x):   0.0489  Pr(y=5|x):   0.0162
  Pr(y=6|x):   0.0045  Pr(y=7|x):   0.0011
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1852

. praccum, using(mphd) xis(3)

. prvalue, x(phd 3.5 phd2 12.25) rest(mean) brief

Predicted rate: 1.6786 
Predicted probabilities:

  Pr(y=0|x,z): 0.3159  Pr(y=1|x):   0.2635
  Pr(y=2|x):   0.2212  Pr(y=3|x):   0.1237
  Pr(y=4|x):   0.0519  Pr(y=5|x):   0.0174
  Pr(y=6|x):   0.0049  Pr(y=7|x):   0.0012
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1589

. praccum, using(mphd) xis(3.5)

. prvalue, x(phd 4 phd2 16) rest(mean) brief

Predicted rate: 1.6617 
Predicted probabilities:

  Pr(y=0|x,z): 0.2909  Pr(y=1|x):   0.2761
  Pr(y=2|x):   0.2294  Pr(y=3|x):   0.1271
  Pr(y=4|x):   0.0528  Pr(y=5|x):   0.0175
  Pr(y=6|x):   0.0049  Pr(y=7|x):   0.0012
  Pr(y=8|x):   0.0002  Pr(y=9|x):   0.0000

Pr(Always0|z): 0.1247

. praccum, using(mphd) xis(4) gen(phsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   phsqx |       7         2.5   1.080123          1          4  
  phsqmu |       7    1.572518   .1138613   1.374963   1.678554  
  phsqp0 |       7    .3445365   .0307107   .2908717   .3696901  
  phsqp1 |       7    .2698241   .0132464   .2573594   .2941026  
  phsqp2 |       7    .2116526   .0102729     .20219   .2293765  
  phsqp3 |       7    .1112016   .0126411   .0926679   .1270538  
  phsqp4 |       7     .044012   .0079031   .0318537   .0527822  
  phsqp5 |       7    .0139929   .0033871   .0087595    .017542  
  phsqp6 |       7    .0037215    .001119   .0020073    .004877  
  phsqp7 |       7    .0008514   .0003026   .0003943   .0011695  
  phsqp8 |       7     .000171   .0000695   .0000678   .0002454  
  phsqp9 |       7    .0000306   .0000139   .0000104   .0000458  
 phsqinf |       7    .1716897    .027096   .1247377   .1990634  
  phsqs0 |       7    .3445365   .0307107   .2908717   .3696901  
  phsqs1 |       7    .6143606   .0349703   .5669414   .6620474  
  phsqs2 |       7    .8260132   .0254049   .7963178   .8642374  
  phsqs3 |       7    .9372148   .0127955   .9233716   .9569054  
  phsqs4 |       7    .9812268   .0048947   .9761539   .9887591  
  phsqs5 |       7    .9952197   .0015079   .9936534   .9975187  
  phsqs6 |       7    .9989412    .000389   .9985304    .999526  
  phsqs7 |       7    .9997926   .0000864   .9996998   .9999203  
  phsqs8 |       7    .9999636   .0000169   .9999452   .9999881  
  phsqs9 |       7    .9999942   2.96e-06    .999991   .9999985  

. list phsq* in 1/7

Observation 1

       phsqx            1      phsqmu     1.374963      phsqp0     .3679448
      phsqp1     .2941026      phsqp2       .20219      phsqp3     .0926679
      phsqp4     .0318537      phsqp5     .0087595      phsqp6     .0020073
      phsqp7     .0003943      phsqp8     .0000678      phsqp9     .0000104
     phsqinf     .1540462      phsqs0     .3679448      phsqs1     .6620474
      phsqs2     .8642374      phsqs3     .9569054      phsqs4     .9887591
      phsqs5     .9975187      phsqs6      .999526      phsqs7     .9999203
      phsqs8     .9999881      phsqs9     .9999985


Observation 2

       phsqx          1.5      phsqmu     1.471877      phsqp0     .3696901
      phsqp1      .276326      phsqp2     .2033589      phsqp3     .0997731
      phsqp4     .0367134      phsqp5     .0108075      phsqp6     .0026512
      phsqp7     .0005575      phsqp8     .0001026      phsqp9     .0000168
     phsqinf      .181953      phsqs0     .3696901      phsqs1     .6460161
      phsqs2     .8493751      phsqs3     .9491482      phsqs4     .9858616
      phsqs5     .9966691      phsqs6     .9993203      phsqs7     .9998778
      phsqs8     .9999804      phsqs9     .9999971


Observation 3

       phsqx            2      phsqmu     1.549926      phsqp0     .3681547
      phsqp1     .2638861      phsqp2      .204502      phsqp3     .1056543
      phsqp4     .0409391      phsqp5     .0126905      phsqp6     .0032782
      phsqp7     .0007259      phsqp8     .0001406      phsqp9     .0000242
     phsqinf     .1978975      phsqs0     .3681547      phsqs1     .6320407
      phsqs2     .8365427      phsqs3      .942197      phsqs4     .9831361
      phsqs5     .9958267      phsqs6     .9991049      phsqs7     .9998308
      phsqs8     .9999714      phsqs9     .9999956


Observation 4

       phsqx          2.5      phsqmu     1.612524      phsqp0     .3587571
      phsqp1       .25751      phsqp2     .2076205      phsqp3     .1115977
      phsqp4     .0449885      phsqp5      .014509      phsqp6     .0038994
      phsqp7     .0008983      phsqp8     .0001811      phsqp9     .0000324
     phsqinf     .1990634      phsqs0     .3587571      phsqs1     .6162671
      phsqs2     .8238876      phsqs3     .9354853      phsqs4     .9804738
      phsqs5     .9949828      phsqs6     .9988821      phsqs7     .9997804
      phsqs8     .9999614      phsqs9     .9999939


Observation 5

       phsqx            3      phsqmu     1.658055      phsqp0     .3404659
      phsqp1     .2573594      phsqp2     .2133581      phsqp3     .1179198
      phsqp4     .0488794      phsqp5     .0162089      phsqp6     .0044792
      phsqp7      .001061      phsqp8     .0002199      phsqp9     .0000405
     phsqinf     .1852483      phsqs0     .3404659      phsqs1     .5978253
      phsqs2     .8111833      phsqs3     .9291031      phsqs4     .9779825
      phsqs5     .9941915      phsqs6     .9986707      phsqs7     .9997317
      phsqs8     .9999515      phsqs9     .9999921


Observation 6

       phsqx          3.5      phsqmu     1.678554      phsqp0     .3158713
      phsqp1     .2635151      phsqp2     .2211622      phsqp3     .1237442
      phsqp4     .0519278      phsqp5     .0174327      phsqp6      .004877
      phsqp7     .0011695      phsqp8     .0002454      phsqp9     .0000458
     phsqinf      .158882      phsqs0     .3158713      phsqs1     .5793865
      phsqs2     .8005487      phsqs3     .9242929      phsqs4     .9762207
      phsqs5     .9936534      phsqs6     .9985304      phsqs7     .9996998
      phsqs8     .9999452      phsqs9      .999991


Observation 7

       phsqx            4      phsqmu     1.661729      phsqp0     .2908717
      phsqp1     .2760697      phsqp2     .2293765      phsqp3     .1270538
      phsqp4     .0527822      phsqp5      .017542      phsqp6     .0048583
      phsqp7     .0011533      phsqp8     .0002396      phsqp9     .0000442
     phsqinf     .1247377      phsqs0     .2908717      phsqs1     .5669414
      phsqs2     .7963178      phsqs3     .9233716      phsqs4     .9761539
      phsqs5     .9936958      phsqs6     .9985541      phsqs7     .9997074
      phsqs8      .999947      phsqs9     .9999912


. 
. * zinb
. use couart2, clear
(Academic Biochemists / S Long)

. gen phd2 = phd*phd

. quietly zinb art fem mar kid5 phd phd2 ment, inf(fem mar)

. prvalue, x(phd 1 phd2 1) rest(mean) brief

Predicted rate: 1.3502 
Predicted probabilities:

  Pr(y=0|x,z): 0.3465  Pr(y=1|x):   0.2935
  Pr(y=2|x):   0.1790  Pr(y=3|x):   0.0950
  Pr(y=4|x):   0.0467  Pr(y=5|x):   0.0218
  Pr(y=6|x):   0.0099  Pr(y=7|x):   0.0043
  Pr(y=8|x):   0.0019  Pr(y=9|x):   0.0008

Pr(Always0|z): 0.0000

. praccum, saving(mphd) xis(1)

. prvalue, x(phd 1.5 phd2 2.25) rest(mean) brief

Predicted rate: 1.4681 
Predicted probabilities:

  Pr(y=0|x,z): 0.3222  Pr(y=1|x):   0.2874
  Pr(y=2|x):   0.1846  Pr(y=3|x):   0.1032
  Pr(y=4|x):   0.0534  Pr(y=5|x):   0.0263
  Pr(y=6|x):   0.0125  Pr(y=7|x):   0.0058
  Pr(y=8|x):   0.0026  Pr(y=9|x):   0.0012

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(1.5)

. prvalue, x(phd 2 phd2 4) rest(mean) brief

Predicted rate: 1.5632 
Predicted probabilities:

  Pr(y=0|x,z): 0.3043  Pr(y=1|x):   0.2819
  Pr(y=2|x):   0.1880  Pr(y=3|x):   0.1091
  Pr(y=4|x):   0.0586  Pr(y=5|x):   0.0300
  Pr(y=6|x):   0.0148  Pr(y=7|x):   0.0071
  Pr(y=8|x):   0.0034  Pr(y=9|x):   0.0016

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(2)

. prvalue, x(phd 2.5 phd2 6.25) rest(mean) brief

Predicted rate: 1.6299 
Predicted probabilities:

  Pr(y=0|x,z): 0.2926  Pr(y=1|x):   0.2778
  Pr(y=2|x):   0.1899  Pr(y=3|x):   0.1129
  Pr(y=4|x):   0.0622  Pr(y=5|x):   0.0326
  Pr(y=6|x):   0.0165  Pr(y=7|x):   0.0081
  Pr(y=8|x):   0.0039  Pr(y=9|x):   0.0019

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(2.5)

. prvalue, x(phd 3 phd2 9) rest(mean) brief

Predicted rate: 1.6641 
Predicted probabilities:

  Pr(y=0|x,z): 0.2869  Pr(y=1|x):   0.2756
  Pr(y=2|x):   0.1907  Pr(y=3|x):   0.1148
  Pr(y=4|x):   0.0640  Pr(y=5|x):   0.0339
  Pr(y=6|x):   0.0174  Pr(y=7|x):   0.0087
  Pr(y=8|x):   0.0043  Pr(y=9|x):   0.0021

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(3)

. prvalue, x(phd 3.5 phd2 12.25) rest(mean) brief

Predicted rate: 1.6639 
Predicted probabilities:

  Pr(y=0|x,z): 0.2869  Pr(y=1|x):   0.2756
  Pr(y=2|x):   0.1907  Pr(y=3|x):   0.1148
  Pr(y=4|x):   0.0639  Pr(y=5|x):   0.0339
  Pr(y=6|x):   0.0174  Pr(y=7|x):   0.0087
  Pr(y=8|x):   0.0042  Pr(y=9|x):   0.0020

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(3.5)

. prvalue, x(phd 4 phd2 16) rest(mean) brief

Predicted rate: 1.629  
Predicted probabilities:

  Pr(y=0|x,z): 0.2928  Pr(y=1|x):   0.2778
  Pr(y=2|x):   0.1898  Pr(y=3|x):   0.1129
  Pr(y=4|x):   0.0621  Pr(y=5|x):   0.0325
  Pr(y=6|x):   0.0165  Pr(y=7|x):   0.0081
  Pr(y=8|x):   0.0039  Pr(y=9|x):   0.0019

Pr(Always0|z): 0.0000

. praccum, using(mphd) xis(4) gen(phsq)

New variables created by praccum:

Variable |     Obs        Mean   Std. Dev.       Min        Max
---------+-----------------------------------------------------
   phsqx |       7         2.5   1.080123          1          4  
  phsqmu |       7     1.56693   .1178974   1.350232   1.664144  
  phsqp0 |       7    .3045891   .0222701    .286861   .3464985  
  phsqp1 |       7    .2813902   .0067762   .2756306   .2935313  
  phsqp2 |       7    .1875111   .0043145   .1790153   .1906628  
  phsqp3 |       7    .1089495   .0073791   .0950175   .1147842  
  phsqp4 |       7    .0586899    .006478   .0466758   .0639548  
  phsqp5 |       7    .0301389   .0045497   .0218213   .0339129  
  phsqp6 |       7    .0149773   .0028386   .0098564   .0173743  
  phsqp7 |       7    .0072676   .0016435   .0043407   .0086786  
  phsqp8 |       7    .0034636   .0009036   .0018748   .0042516  
  phsqp9 |       7    .0016277   .0004783   .0007974   .0020511  
 phsqinf |       7    4.47e-10          0   4.47e-10   4.47e-10  
  phsqs0 |       7    .3045891   .0222701    .286861   .3464985  
  phsqs1 |       7    .5859793   .0290329   .5624915   .6400297  
  phsqs2 |       7    .7734903   .0247415   .7531543    .819045  
  phsqs3 |       7    .8824398   .0173704   .8679386   .9140625  
  phsqs4 |       7    .9411297   .0108964   .9318934   .9607383  
  phsqs5 |       7    .9712686   .0063487   .9658063   .9825596  
  phsqs6 |       7    .9862459   .0035112   .9831806    .992416  
  phsqs7 |       7    .9935135   .0018682   .9918591   .9967567  
  phsqs8 |       7    .9969771   .0009648   .9961107   .9986315  
  phsqs9 |       7    .9986047   .0004866   .9981618   .9994289  

. list phsq* in 1/7

Observation 1

       phsqx            1      phsqmu     1.350232      phsqp0     .3464985
      phsqp1     .2935313      phsqp2     .1790153      phsqp3     .0950175
      phsqp4     .0466758      phsqp5     .0218213      phsqp6     .0098564
      phsqp7     .0043407      phsqp8     .0018748      phsqp9     .0007974
     phsqinf     4.47e-10      phsqs0     .3464985      phsqs1     .6400297
      phsqs2      .819045      phsqs3     .9140625      phsqs4     .9607383
      phsqs5     .9825596      phsqs6      .992416      phsqs7     .9967567
      phsqs8     .9986315      phsqs9     .9994289


Observation 2

       phsqx          1.5      phsqmu     1.468137      phsqp0     .3221712
      phsqp1     .2874038      phsqp2     .1845785      phsqp3     .1031686
      phsqp4     .0533689      phsqp5     .0262742      phsqp6     .0124975
      phsqp7     .0057958      phsqp8     .0026361      phsqp9     .0011807
     phsqinf     4.47e-10      phsqs0     .3221712      phsqs1      .609575
      phsqs2     .7941535      phsqs3      .897322      phsqs4     .9506909
      phsqs5     .9769651      phsqs6     .9894626      phsqs7     .9952584
      phsqs8     .9978945      phsqs9     .9990752


Observation 3

       phsqx            2      phsqmu     1.563207      phsqp0     .3043071
      phsqp1     .2818844      phsqp2     .1879805      phsqp3     .1091018
      phsqp4     .0586039      phsqp5     .0299585      phsqp6     .0147967
      phsqp7     .0071254      phsqp8     .0033652      phsqp9     .0015651
     phsqinf     4.47e-10      phsqs0     .3043071      phsqs1     .5861915
      phsqs2      .774172      phsqs3     .8832738      phsqs4     .9418777
      phsqs5     .9718363      phsqs6      .986633      phsqs7     .9937584
      phsqs8     .9971237      phsqs9     .9986888


Observation 4

       phsqx          2.5      phsqmu     1.629888      phsqp0     .2926165
      phsqp1     .2777898      phsqp2     .1898526      phsqp3     .1129261
      phsqp4     .0621651      phsqp5     .0325686      phsqp6     .0164856
      phsqp7     .0081359      phsqp8      .003938      phsqp9      .001877
     phsqinf     4.47e-10      phsqs0     .2926165      phsqs1     .5704063
      phsqs2     .7602589      phsqs3      .873185      phsqs4     .9353501
      phsqs5     .9679188      phsqs6     .9844043      phsqs7     .9925402
      phsqs8     .9964782      phsqs9     .9983552


Observation 5

       phsqx            3      phsqmu     1.664144      phsqp0      .286861
      phsqp1     .2756306      phsqp2     .1906628      phsqp3     .1147842
      phsqp4     .0639548      phsqp5     .0339129      phsqp6     .0173743
      phsqp7     .0086786      phsqp8     .0042516      phsqp9     .0020511
     phsqinf     4.47e-10      phsqs0      .286861      phsqs1     .5624915
      phsqs2     .7531543      phsqs3     .8679386      phsqs4     .9318934
      phsqs5     .9658063      phsqs6     .9831806      phsqs7     .9918591
      phsqs8     .9961107      phsqs9     .9981618


Observation 6

       phsqx          3.5      phsqmu     1.663856      phsqp0     .2869087
      phsqp1     .2756488      phsqp2     .1906564      phsqp3     .1147689
      phsqp4     .0639399      phsqp5     .0339016      phsqp6     .0173668
      phsqp7     .0086739      phsqp8     .0042489      phsqp9     .0020495
     phsqinf     4.47e-10      phsqs0     .2869087      phsqs1     .5625575
      phsqs2     .7532139      phsqs3     .8679828      phsqs4     .9319227
      phsqs5     .9658243      phsqs6     .9831911      phsqs7      .991865
      phsqs8      .996114      phsqs9     .9981635


Observation 7

       phsqx            4      phsqmu     1.629042      phsqp0     .2927608
      phsqp1     .2778427      phsqp2     .1898313      phsqp3     .1128793
      phsqp4     .0621206      phsqp5     .0325355      phsqp6     .0164638
      phsqp7     .0081227      phsqp8     .0039304      phsqp9     .0018728
     phsqinf     4.47e-10      phsqs0     .2927608      phsqs1     .5706035
      phsqs2     .7604347      phsqs3      .873314      phsqs4     .9354346
      phsqs5     .9679701      phsqs6     .9844338      phsqs7     .9925565
      phsqs8     .9964869      phsqs9     .9983597


. 
. log close