st: RE: gpscore for propensity score of 4 observed groups

[Nick Cox <n.j.cox@durham.ac.uk>]

I don't know the answer, but you should clarify that -gpscore- is

SJ-9-4  st0150_2  The dose-response function adj. for the gen. propensity score
        (help doseresponse, gpscore if installed) . . . . M. Bia and A. Mattei
        Q4/09   SJ 9(4):652
        updated to specify version 10 in order to run correctly
        under Stata 11

SJ-8-4  st0150_1  The dose-response function adj. for the gen. propensity score
        (help doseresponse, gpscore if installed) . . . . M. Bia and A. Mattei
        Q4/08   SJ 8(4):594
        improved handling of predicted probabilities; correction
        of some references to variables named treatment_level_plus
        or similar

SJ-8-3  st0150   The dose-response function adj. for the gen. propensity score
        (help doseresponse, gpscore if installed) . . . . M. Bia and A. Mattei
        Q3/08   SJ 8(3):354--373
        estimates the propensity score with a continuous treatment,
        tests the balancing property of the generalized propensity
        score, and estimates the dose-response function

Why is this request to explain user-written software seemingly so difficult to understand? It is good sense to make clear what you are talking about. 

Nick 
n.j.cox@durham.ac.uk 

oconnor7@u.washington.edu

I have technical question regarding gpscore.

If my treatment variable is type of injury severity, what should I put for "cutpoints", "index", and "nq_gps"? The treatment level has four different levels of injury severity.

What I am seeking is to create a propensity score for the likelihood that a person develops a certain type of injury, in this case TBI, based upon a list of covariates.

This is the syntax I used for a preliminary analysis:

gpscore  age child_gender prev_anxiety prev_depn prev_adhd prev_behavioral prev_other_mh, t(inj_severity) gpscore(mygps) predict(hat_inj_sev) sigma(hat_sd) cutpoints (inj_severity) index(p50) nq_gps(4)

Again, I am unsure of cutpoints, index, and nq_qps given that I have four established groups for the treatment variable.


*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/