.-
help for ^stgtcalc^, ^stgtplot^                    (Patrick Royston, Peter Sasieni)
.-

Time-varying regression coefficients in Cox PH models
      (via Grambsch and Therneau residuals)
-----------------------------------------------------


    ^stgtcalc^


    ^stgtplot^ varlist [^, ci le^vel^(^#^)^ ^me^an graph_options ]



Description
-----------

^stgtcalc^ calculates local estimates of the regression coefficients in the Cox
model most recently fitted for using @stcox@.

^stgtplot^ plots the estimates and confidence intervals against the failure
times. varlist is (a subset of) the independent variables in model fitted by
^stcox^.


Options
-------

^ci^ causes a ^level()^% confidence band to be plotted around the estimates of the
    regression coefficients. The pointwise confidence intervals are derived
    from the running line smoother @running@.

^level(^#^)^ sets the width of a pointwise confidence band for the estimated
    b(t). Default: current value of system macro ^$S_level^ (usually 95%).

^mean^ uses the ^mean^ option of @running@ to produce a running mean smoother.

graph_options are any of Stata's ^graph, twoway^ options


Remarks
-------

^stgtcalc^ implements a suggestion of Grambsch and Therneau (1994). For an xvar
with estimated regression coefficient b, ^stgtcalc^ calculates 

	g(T) = b + #deaths * V^^-1 * r(T)

at each failure time T, where V is the estimated dispersion matrix of the
regression coefficients, #deaths is the number of uncensored observations and
r(T) is the vector of Schoenfeld residuals at T. The g(T)'s are smoothed and
the smoothed values are saved. g(T) is a local estimate of beta, so the shape
of the smoothed profile may indicate whether beta depends on time. If so, the
dataset has non-proportional hazards and the model will need to be refined.

Note that the smoothed g(T) may be unreliable in the initial and final phases
of follow-up.  This may be partly due to lack of robustness in the smoother.

The values of g(T) for each xvar are saved with names ^GT^xvar1, ^GT^xvar2, etc.

If more than one xvar is to be plotted by ^stgtplot^, a panel of plots is produced.


Example
-------

	. ^stset survtime cens^
	. ^stcox x1 x2^
	. ^stgtcalc^
	. ^stgtplot x1^
	. ^stgtplot x1 x2, ci mean^


Reference
---------

Grambsch PM, Therneau TM. Proportional hazards tests and diagnostics based
    on weighted residuals. Biometrika 81: 515-526 (1994).


Authors
-------

        Patrick Royston
        Imperial College School of Medicine, London
        FAX:  (011)-44-181-383-8573

        Peter Sasieni
        Imperial Cancer Research Fund, London
        FAX:  (011)-44-171-269-3429


Also see
--------

 Manual:  [R] st stcox
On-line:  help for @stcox@, @stbtcalc@
