nathan_modularity_issue002

while(t<2015){
    rdapathgenerator <- function(years, prefix, dirpath = 'C:/Users/Nathan/Desktop/vax-p vax-d'){
      errMsg000 = "var 'year' has to be integer!"
      if(is.numeric(years)){
        if(all(years == as.integer(years))){
          years = as.integer(years)
        } else{
          stop(errMsg000) #numeric but not integer
        }
      } else stop(errMsg000) #not even numeric
      ### if necessary conditions are met:
      result = list()
      for(year in years){
        rdafname = paste(prefix, year, sep='')
        rdafextension = 'RData'
        value = list(
          fullpath=file.path(dirpath, paste(rdafname,'.', rdafextension, sep=''), fsep = .Platform$file.sep),
          filename=rdafname, year=year, extension=rdafextension
        )
        #print(newvalue)
        result[[length(result)+1]] = value
      }
      return(result)
    }

    years_seq= c(t)
    for(rdafile in rdapathgenerator(years=years_seq, prefix='dd')){
      load(rdafile$fullpath)
      data_dd<-dd
      data_dd[c((1:392),(449:2464)),]<-0
      data_ddch<-as.matrix(data_dd)
      load(rdafile$fullpath)
      data_dd<-dd
      data_dd[-c((281:336),(2353:2408)),]<-0
      data_ddusa<-as.matrix(data_dd)
      #save(data_ddusa_, file=paste('data_usa_', rdafile$filename, '.RData', sep=''))
      load(rdafile$fullpath)
      data_dd<-dd
      data_dd<-as.matrix(data_dd)
    }

    for(rdafile in rdapathgenerator(years=years_seq, prefix='L')){
      load(rdafile$fullpath)
    }

    for(rdafile in rdapathgenerator(years=years_seq, prefix='a')){
      load(rdafile$fullpath)
    }
    A<- data_a

    for(rdafile in rdapathgenerator(years=years_seq, prefix='va')){
      load(rdafile$fullpath)
      load(rdafile$fullpath)
      data_va<-va
      data_va[-c(225:228),]<-0
      data_vabr_ins<-data_va
      load(rdafile$fullpath)
      data_va<-va
      data_va[-c((229:246),(251),(263)),]<-0
      data_vabr_ind<-data_va
      load(rdafile$fullpath)
      data_va<-va
      data_va[-c((247:250),(252:262),(264:280)),]<-0
      data_vabr_ser<-data_va
    }

    for(rdafile in rdapathgenerator(years=years_seq, prefix='va')){
      va_year <- rdafile$year
      k1<-L%*%data_dd
      k2<-L%*%data_ddch
      k3<-L%*%data_ddusa
      vax_p_brser<- data_vabr_ser%*%k1
    }
}