pkmatador

nathan_modularity_issue001.R

Mar 29th, 2019
151
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 2.60 KB | None | 0 0
  1. rdapathgenerator <- function(years, dirpath = 'C:/Users/Nathan/Desktop/vax-p vax-d'){
  2.     errMsg000 = "var 'year' has to be integer!"
  3.     if(is.numeric(years)){
  4.         if(all(years == as.integer(years))){
  5.             years = as.integer(years)
  6.         } else{
  7.             stop(errMsg000) #numeric but not integer
  8.         }
  9.     } else stop(errMsg000) #not even numeric
  10.     ### if necessary conditions are met:
  11.     result = list()
  12.     for(year in years){
  13.         rdafname = paste('va', year, sep='')
  14.         rdafextension = 'RData'
  15.         value = list(
  16.             fullpath=file.path(dirpath, paste(rdafname,'.', rdafextension, sep=''), fsep = .Platform$file.sep),
  17.             filename=rdafname#, year=year, extension=rdafextension
  18.         )
  19.         #print(newvalue)
  20.         result[[length(result)+1]] = value
  21.     }
  22.     return(result)
  23. }
  24.  
  25. years_seq = seq(2000,2014)
  26. for(rdafile in rdapathgenerator(years=years_seq)){
  27.     load(rdafile$fullpath)
  28.     data_va00<-va
  29.     data_va00[-c(415),]<-0
  30.     data_vach_ELET0<-data_va00
  31.     save(data_vach_ELET0, file=paste('data_vach_ELET_', rdafile$filename, '.RData', sep=''))
  32.     load(rdafile$fullpath)
  33.     data_va00<-va
  34.     data_va00[-c(247),]<-0
  35.     data_vabr_ELET00<-data_va00
  36.     save(data_vabr_ELET00, file=paste('data_vabr_ELET_', rdafile$filename, '.RData', sep=''))
  37.     load(rdafile$fullpath)
  38.     data_va00<-va
  39.     data_va00[-c((297),(2369)),]<-0
  40.     data_vausa_ELET00<-data_va00
  41.     save(data_vausa_ELET00, file=paste('data_vausa_ELET_', rdafile$filename, '.RData', sep=''))
  42.     load(rdafile$fullpath)
  43.     data_va00<-va
  44.     data_va00[-c(1361),]<-0
  45.     data_vajap_ELET00<-data_va00
  46.     save(data_vajap_ELET00, file=paste('data_vajap_ELET_', rdafile$filename, '.RData', sep=''))
  47.     load(rdafile$fullpath)
  48.     data_va00<-va
  49.     data_va00[-c(577),]<-0
  50.     data_vager_ELET00<-data_va00
  51.     save(data_vager_ELET00, file=paste('data_vager_ELET_', rdafile$filename, '.RData', sep=''))
  52.     load(rdafile$fullpath)
  53.     data_va00<-va
  54.     data_va00[-c(seq(1137,1232,56),(1417),(2313)),]<-0
  55.     data_va_OT_AS_elet00<-data_va00
  56.     save(data_va_OT_AS_elet00, file=paste('data_vaot_as_ELET_', rdafile$filename, '.RData', sep=''))
  57.     load(rdafile$fullpath)
  58.     data_va00<-va
  59.     data_va00[-c(seq(73,224,56),(353),seq(465,560,56),seq(633,1120,56),seq(1249,1344,56),seq(1473,1624,56),seq(1697,2016,56),seq(2089,2240,56)),]<-0
  60.     data_va_OT_EU_elet00<-data_va00
  61.     save(data_va_OT_EU_elet00, file=paste('data_vaot_eu_ELET_', rdafile$filename, '.RData', sep=''))
  62.     load(rdafile$fullpath)
  63.     data_va00<-va
  64.     data_va00[c((17),(1641),(2089),(2257),(2425)),]<-0
  65.     data_vaROW_elet00<-data_va00
  66.     save(data_vaROW_elet00, file=paste('data_varow_ELET_', rdafile$filename, '.RData', sep=''))
  67. }
Advertisement
Add Comment
Please, Sign In to add comment