I have the following code to generate predicted probabilities of child death when education=edu1 at different origin of births (SADC, ROA, and SA). However, I want to set the continuous regressors to their mean, some categorical regressors to 1, and some categorical regressors to 0. Can you help me code these "adjustments"?   What I have:   proc qlim data=sa.Ipumsi_sa_reg;     model CHDIED=SADC ROA AGE EDU2 EDU3 EDU4 EDU5 EDU6 MARRIED PERSONS EMPLOY OWNER URBAN_I URBAN_F IMOVER RECENT / DISCRETE(d=probit);     weight WTPER;     where OSEAS=0 and 20 le AGE le 35 and EDU1=1;     output out=outpp_edu1 proball; run; quit; proc means data=outpp_edu1; var Prob2_CHDIED; where sadc=1; run; proc means data=outpp_edu1; var Prob2_CHDIED; where roa=1; run; proc means data=outpp_edu1; var Prob2_CHDIED; where sa=1; run;     Thank you! Stephanie ... View more

I am running a probit regression and want to report marginal effects and the standard errors of the marginal effects. Here is my code for the regression and marginal effects:   proc qlim data=sa.Ipumsi_sa_reg;     model CHDIED=SADC ROA AGE / DISCRETE(d=probit);     where OSEAS=0 and 20 le AGE le 35;     weight WTPER;     output out=outme marginal; run; quit; proc means data=outme n mean;       var Meff_P2_SADC Meff_P2_ROA Meff_P2_AGE;       title 'Average of the Individual Marginal Effects';    run;    quit;   How can I calculate the standard errors of the marginal effects?   Thank you. ... View more

I want to only keep the group of observations (by 'id') if the value of 'gestage' is not missing or </= 20 or >/=43. Delete any group as long as one observation has a missing value in 'gestage' or an implausible value in 'gestage'. Thanks for helping! data test; input id gestage; datalines; 1 36 1 38 2 . 2 41 3 30 3 . 3 32 4 19 4 41 5 35 5 36 5 36 ; run; ... View more

I'm creating a data set of mothers who have given birth at least 2 times during the study period, with summary statistics to be reported for births at time 1 and time 2. I created a variable birth_num that indicates the chronological number of each birth for the same woman during the study period. For example: MOTHER_ID    BIRTH_YR      birth_num 5                     2009              1 7                     2004              1 10                   2005              1 10                   2008              2 10                   2010              3 12                   2006              1 12                   2008              2 14                   2004              1 14                   2006              2 14                   2008              3 I'm finding that as I apply the exclusion criteria to get the sample (only keep records with valid values for first birth, infant health outcomes, prenatal care, mother's race/ethnicity and immigrant status, and women who have given birth at least 2 times), I'm losing either the birth record at time 1 or time 2 for some MOTHER_IDs. For example, if the record for MOTHER_ID=10 where birth_num=2 is excluded, I'm unable to include this record in the summary statistics for births at time 2. However, the record for MOTHER_ID=10 where birth_num=1 is retained and included in summary statistics for births at time 1. Data after exclusion criteria applied: MOTHER_ID    BIRTH_YR      birth_num 10                   2005              1 10                   2010              3 12                   2006              1 12                   2008              2  14                   2004              1 14                   2006              2 14                   2008              3 How can I only retain birth records that meet selection criteria for births at time 1 and time 2? These are the records I am interested in keeping for summary statistics: MOTHER_ID    BIRTH_YR      birth_num 12                   2006              1 12                   2008              2  14                   2004              1 14                   2006              2 Thank you, Stephanie ... View more

Hi, I'm doing some data analysis for a project and my advisor would like me to run simple linear regressions with robust standard errors on all the models. I was planning to use the /acov option in proc reg to calculate the robust standard errors. A few of the models include interaction of variables. I know that I can use proc glm for interacted variables, but I'm not able to generate the robust standard errors with this procedure. Any thoughts on how to get both robust standard errors and include interacted variables in a simple linear regression procedure? I appreciate any help! Stephanie ... View more