This function plots the estimated results from the 3 parameter vitality model. It is used within the function vitality.k.
plotting.k(r.final,s.final,k.final,mlv,time,sfract,x1,x2,Ni,
pplot,tlab,lplot,cplot,Iplot,gfit,rc.data)r estimate
s estimate
k estimate
TODO mlv
time vector giving the ages at which calculate the estimated value
observed survival fraction
Time 1
Time 2
Initial population
Boolean. Plot cumulative survival fraction?
Boolean. Plot incremental survival?
Boolean. Plot likelihood functions? Provides likelihood function plotting. Defaults to FALSE. Note: these plots are not "likelihood profiles" in that while one parameter is varied, the others are held fixed, rather than re-optimized. (must also have pplot=T.)
Boolean. Plot likelihood contour plot? Provides a likelihood contour plot for a range of r and s values (can be slow so default is FALSE). Must also have lplot=T and pplot=T to get contour plots.
Character, label for time axis. Defaults to "days".
Provides a Pearson C type test for goodness of fit. Default is FALSE. Must provide the initial study population to compute goodness of fit.
Boolean, right-censored data?
See vitality.k for further description of function arguments.