HBW Destination choice prep file. Do not use, it takes forever.

HBWDCPrep.R

#
# HBWDCPrep.R
#
# Prepares Destination choice files for HBW DC

library(foreign)
source("http://raw.githubusercontent.com/osPlanning/omx/dev/api/r/omx.R")

A1201 = read.dbf("C:\\Modelrun\\Model80\\Base\\A1201.DBF")

timeSkimPK = readMatrixOMX("Data\\HSKIM_PK1.omx",'TIME')
XSkimPK = readMatrixOMX("Data\\HSKIM_PK1.omx",'RXING')

HBWPKCHOICE = data.frame(HHPTID=numeric(),I=integer(),J=integer(),Choice=logical(),Time=numeric(),RXing=integer(),Emp=integer())
print(Sys.time())
for(r in seq(1,nrow(Trips.HBWPK))){ # nrow(Trips.HBWPK))){
  i = Trips.HBW[r,"OTAZ"]
  if(i == 0) next
  j = Trips.HBW[r,"D_TAZ"]
  hhptid = Trips.HBW[r,"HHPERTRPID"]
  for(jj in seq(1,j-1)){
    if(A1201[which(A1201$TAZ == jj),"ATYPE"] == 0) next
    emp = sum(A1201[which(A1201$TAZ == jj),c("AGR","MIN","MAN1","MAN2","WHO","RET","TRA","FIRE","EDU","SOC","ART","PUB")])
    HBWPKCHOICE = rbind(HBWPKCHOICE,cbind(HHPTID=hhptid,I=i,J=jj,Choice=FALSE,Time=timeSkimPK[i,jj],RXing=XSkimPK[i,jj],Emp=emp))
  }
  emp = sum(A1201[which(A1201$TAZ == j),c("AGR","MIN","MAN1","MAN2","WHO","RET","TRA","FIRE","EDU","SOC","ART","PUB")])
  HBWPKCHOICE = rbind(HBWPKCHOICE,cbind(HHPTID=hhptid,I=i,J=j,Choice=TRUE,Time=timeSkimPK[i,j],RXing=XSkimPK[i,j],Emp=emp))
  for(jj in seq(j+1,nrow(A1201))){
    if(A1201[which(A1201$TAZ == jj),"ATYPE"] == 0) next
    emp = sum(A1201[which(A1201$TAZ == jj),c("AGR","MIN","MAN1","MAN2","WHO","RET","TRA","FIRE","EDU","SOC","ART","PUB")])
    HBWPKCHOICE = rbind(HBWPKCHOICE,cbind(HHPTID=hhptid,I=i,J=jj,Choice=FALSE,Time=timeSkimPK[i,jj],RXing=XSkimPK[i,jj],Emp=emp))
  }
  print(paste("Rec: ",r," at ", Sys.time(),sep=""))
}
print(Sys.time())