{smcl}
{* *! version 1.0.0 12sep2014 N.Orsini}{...}
{hline}

{title:Title}

{p2colset 5 30 35 2}{...}
   Postestimation tool for laplacereg
{p2colreset}{...}

{title:Introduction}

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

{synoptset 11}{...}
{p2coldent :Command}Description{p_end}
{synoptline}
{synopt :{helpb laplacereg postestimation##laplace_surv:laplace_surv}}tabulate and plot the predicted survivor function{p_end}
{synoptline}
{p2colreset}{...}


{marker laplace_surv}
{title:Syntax for laplace_surv} 

{p 8 40 2}
{cmd:laplace_surv}
[{cmd:,} {it:laplace_surv_options}]
{p_end}

{synoptset 35}{...}
{marker laplace_surv_options}{...}
{synopthdr :laplace_surv_options}
{synoptline}
{p2col:{cmd:at#(}{varname}{cmd:=}{it:#} [{varname}{cmd:=}{it:# ...}]{cmd:)}}
value of the specified covariates and mean of unspecified covariates (with {bf:#} = 1,2,...,10){p_end}

{synopt :{opt r:ange(# #)}}range of analysis time{p_end} 
{synopt:{opt gen:erate(varlist)}}saves variables required to reproduce the plot{p_end}
{synopt: {cmdab:f:ormat(%}{it:fmt}{cmd:)}}specifies the display format{p_end}
{synopt: {opt scatter}}specifies the scatter plot{p_end}
{synopt: {opt line}}specifies the line plot{p_end}
{synopt: {opt lowess}}specifies the lowess plot{p_end}
{synopt:{opt lo:ptions(string)}}specifies {helpb lowess} options{p_end}
{synopt: {it:{help twoway_options:twoway_options}}}controls various aspects of the plot{p_end}
{synopt: {opt show}}displays the values of the survivor function{p_end}
{synopt: {opt lpred}}displays the equation of the linear predictor{p_end}


{synoptline}
 

{title:Description for laplace_surv}

{pstd} Tabulate and plot the predicted survival function after fitting quantiles using Laplace regression model.{p_end}

{title:Options for laplace_surv}

{phang}
{cmd:at#(}{varname}{cmd:=}{it:# ...}{cmd:)} requests that the
covariates specified by {it:varname} be set to {it:#}. By default,
{cmd: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 {opt at#()} and at the mean of all unlisted covariates.

{pmore}
{cmd:at1(}{it:varname}{cmd:=}{it:# ...}{cmd:)},
{cmd:at2(}{it:varname}{cmd:=}{it:# ...}{cmd:)}, ...,
{cmd:at10(}{it:varname}{cmd:=}{it:# ...}{cmd:)} 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.  {opt at1()} specifies the values of the covariates
for the first curve, {opt at2()} specifies the values of the covariates for
the second curve, and so on.

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

{phang} {cmd:generate(}{it:varlist}{cmd:)} saves the variables required to
reproduce the plot. This option is useful if one wants to customize the plot using 
{helpb graph twoway}.

{phang}
{cmd:format(%}{it:fmt}{cmd:)} specifies the display format for presenting
numbers.  {cmd:format(%3.2f)} is the default; see help {helpb format}. 

{phang}
{opt scatter} plots the estimates using {helpb twoway scatter}.

{phang}
{opt line} plots the estimates using {helpb twoway line}.

{phang}
{opt line} plots the estimates using {helpb twoway lowess}.

{phang}{opt loptions(string)} specifies {helpb lowess} options
(for instance, {cmd:bwidth()} or {cmd:mean}) when smoothing the regression
coefficients versus the set of specified quantiles.

{phang}
{opt show} displays the values of the survival functions being plotted.

{phang}
{opt lpred} displays the equation of the linear predictor associated with the plot.


{title:Examples}

{stata "webuse catheter, clear"} 

{pstd}Overall survival{p_end}

{stata "quietly laplacereg time , fail(infect) q(1(1)70)"}
{stata "laplace_surv, line"}

{pstd}Survival curve for men and women evaluated at mean age{p_end}

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

{title:Author}

{pstd}Nicola Orsini{p_end}
{pstd}Unit of Biostatistics{p_end}
{pstd}Institute of Environmental Medicine, Karolinska Institutet{p_end}
{pstd}Stockholm, Sweden{p_end}
{pstd}nicola.orsini@ki.se{p_end}