*! version 1.0.0 12sep2014 N.Orsini

Title

Postestimation tool for laplacereg

Introduction

The following postestimation commands are of special interest after the most
recently fit Laplace regression model using laplacereg:

Command Description

laplace_surv tabulate and plot the predicted survivor function

Syntax for laplace_surv

laplace_surv
[, laplace_surv_options]

laplace_surv_options

at#(varname=# [varname=# ...])
value of the specified covariates and mean of unspecified covariates (with # = 1,2,...,10)

range(# #) range of analysis time

generate(varlist) saves variables required to reproduce the plot

format(%fmt) specifies the display format

scatter specifies the scatter plot

line specifies the line plot

lowess specifies the lowess plot

loptions(string) specifies lowess options

twoway_options controls various aspects of the plot

show displays the values of the survivor function

lpred displays the equation of the linear predictor

Description for laplace_surv

Tabulate and plot the predicted survival function after fitting quantiles using Laplace regression model.

Options for laplace_surv

at#(varname=# ...) requests that the
covariates specified by varname be set to #. By default,
laplace_surv evaluates the function by setting each covariate to its mean
value. This option causes the function to be evaluated at the value of the
covariates listed in at#() and at the mean of all unlisted covariates.

at1(varname=# ...),
at2(varname=# ...), ...,
at10(varname=# ...) specify that multiple curves
(up to 10) be plotted on the same graph. They request that the function
be evaluated at the value of the covariates specified and at the mean of
all unlisted covariates. at1() specifies the values of the covariates
for the first curve, at2() specifies the values of the covariates for
the second curve, and so on.

range(# #) specifies the range of the time axis to be plotted. If this
option is not specified, laplace_surv plots the desired curve on an
interval expanding from 0 to the latest time in the data.

generate(varlist) saves the variables required to
reproduce the plot. This option is useful if one wants to customize the plot using
graph twoway.

format(%fmt) specifies the display format for presenting
numbers. format(%3.2f) is the default; see help format.

scatter plots the estimates using twoway scatter.

line plots the estimates using twoway line.

line plots the estimates using twoway lowess.

loptions(string) specifies lowess options
(for instance, bwidth() or mean) when smoothing the regression
coefficients versus the set of specified quantiles.

show displays the values of the survival functions being plotted.

lpred displays the equation of the linear predictor associated with the plot.

Examples

webuse catheter, clear

Overall survival

quietly laplacereg time , fail(infect) q(1(1)70)
laplace_surv, line

Survival curve for men and women evaluated at mean age

quietly laplacereg time female age, fail(infect) q(1(1)70)
laplace_surv, at1(female=0) at2(female=1) line

Author

Nicola Orsini

Unit of Biostatistics

Institute of Environmental Medicine, Karolinska Institutet

Stockholm, Sweden

nicola.orsini@ki.se