# NOT RUN {
TotalDateRange <- seq(as.Date("2000-01-01"), as.Date("2002-01-01"), "years")
snapshots <- 1 #This uses a 1 year transition matrices
interval <- 1 #This gives a 1 year transition matrix
lstCnt <-rep(list(list()), length(TotalDateRange)-1)
lstInit <-rep(list(list()), length(TotalDateRange)-1)
kk <- 1
k <- 1
#Create transition and inital countsused as inputs to the function
for (l in 1:(length(TotalDateRange)-1)){
istartDate = POSIXTomatlab(as.POSIXlt(as.Date(TotalDateRange[l],format = "%Y-%m-%d")))
iendDate = POSIXTomatlab(as.POSIXlt(as.Date(TotalDateRange[l+1],format = "%Y-%m-%d")))
DateRange <- as.Date(matlabToPOSIX(cfdates(istartDate,iendDate,snapshots)))
for(i in 1:(length(DateRange)-1)){
sDate <- DateRange[i] # i.e "3/31/1990"
eDate <- DateRange[i+1] # i.e "6/30/1990"
t<-TransitionProb(data,sDate, eDate, 'duration', snapshots, interval)
lstCnt[[k]][[kk]] <- t$sampleTotals$totalsMat
lstInit[[k]][[kk]] <- t$sampleTotals$totalsVec
if(kk>=snapshots){
kk <- 1
k <- k+1
} else {
kk <- kk+1
}
}
}
transtype <- "genAverage"
gYear <- utils::head(TotalDateRange,-1)
nStates <- nrow(as.data.frame(lstInit[[1]][1]))
AverageGeneratorMatrices<-duration.TTC(lstCnt,lstInit,gYear,snapshots,nStates,transtype,
TotalDateRange)
# }
Run the code above in your browser using DataLab