Raw curves & Outcome estimates

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:plyr':
## 
##     arrange, count, desc, failwith, id, mutate, rename, summarise,
##     summarize

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

## Type 'citation("pROC")' for a citation.

## 
## Attaching package: 'pROC'

## The following objects are masked from 'package:stats':
## 
##     cov, smooth, var

## Package: CFC, Version: 1.1.2

## Cause-specific Framework for Competing-risk survival analysis

## School of Public Heath, Imperial College London &

## Davidson Kempner Capital Management, NY

## 
## Attaching package: 'purrr'

## The following object is masked from 'package:plyr':
## 
##     compact

## Loading required package: viridisLite

## NOTE: Either Arial Narrow or Roboto Condensed fonts are required to use these themes.

##       Please use hrbrthemes::import_roboto_condensed() to install Roboto Condensed and

##       if Arial Narrow is not on your system, please see https://bit.ly/arialnarrow

## Loading required package: prodlim

## riskRegression version 2022.03.22

## 
## Attaching package: 'riskRegression'

## The following objects are masked from 'package:pec':
## 
##     ipcw, selectCox

## The following object is masked from 'package:InformationValue':
## 
##     plotROC

## Loading required package: Hmisc

## Loading required package: lattice

## Loading required package: Formula

## 
## Attaching package: 'Hmisc'

## The following objects are masked from 'package:DescTools':
## 
##     %nin%, Label, Mean, Quantile

## The following objects are masked from 'package:dplyr':
## 
##     src, summarize

## The following objects are masked from 'package:plyr':
## 
##     is.discrete, summarize

## The following objects are masked from 'package:base':
## 
##     format.pval, units

## Loading required package: SparseM

## 
## Attaching package: 'SparseM'

## The following object is masked from 'package:base':
## 
##     backsolve

## Loading required package: plotrix

## Loading required package: gtools

## 
## Attaching package: 'cutpointr'

## The following objects are masked from 'package:InformationValue':
## 
##     npv, precision, sensitivity, specificity

## The following objects are masked from 'package:pROC':
## 
##     auc, roc

## 
## Attaching package: 'matrixStats'

## The following object is masked from 'package:dplyr':
## 
##     count

## The following object is masked from 'package:plyr':
## 
##     count

## 
## Attaching package: 'mice'

## The following object is masked from 'package:stats':
## 
##     filter

## The following objects are masked from 'package:base':
## 
##     cbind, rbind

## 
## Attaching package: 'ROCit'

## The following object is masked from 'package:gtools':
## 
##     logit

## Loading required package: doSNOW

## Loading required package: foreach

## 
## Attaching package: 'foreach'

## The following objects are masked from 'package:purrr':
## 
##     accumulate, when

## The following object is masked from 'package:DescTools':
## 
##     %:%

## Loading required package: iterators

## Loading required package: snow

## Loading required package: mexhaz

## 
## Attaching package: 'mexhaz'

## The following object is masked from 'package:stats':
## 
##     vcov

## Loading required package: Matrix

## 
## Attaching package: 'Matrix'

## The following objects are masked from 'package:tidyr':
## 
##     expand, pack, unpack

## Loading required package: gridExtra

## 
## Attaching package: 'gridExtra'

## The following object is masked from 'package:dplyr':
## 
##     combine

Transition from CDR 0.5 AD to CDR >=1

Transition from CDR 1 AD to CDR >=2

Transition from CDR 2 AD to 3

main.text0="CDR 3 AD mortality"
x.text="years from CDR 3 AD"
dat0=cdr3.death

main.text=paste("All ages: ",main.text0,sep="")
dat1=dat0
mort.plot()

age.table=table(dat0$age.group10)
age.groups=names(age.table)
t.max=10

par(mfrow=c(1,1))
for(i in 1:length(age.groups)){
  main.text=paste("Age ",age.groups[i],": ",main.text0,sep="")
  dat1=dat0[dat0$age.group10==age.groups[i],]
  mort.plot()}