提交: 1b8f8131ba0dc385fea00ebe92177b6143f274e7 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 1b8f8131ba0dc385fea00ebe92177b6143f274e7 (tree)
• 时间: 2022-12-12 22:59:07
• 作者: Lorenzo Isella <lorenzo.isella@gmai...>
• Commiter: Lorenzo Isella

Log Message

Code seld-contained.

更改概述

• modified: R-codes/arima_simple2_parallel.R

差异

diff -r 753a0d2ce999 -r 1b8f8131ba0d R-codes/arima_simple2_parallel.R
--- a/R-codes/arima_simple2_parallel.R	Mon Dec 12 14:58:47 2022 +0100
+++ b/R-codes/arima_simple2_parallel.R	Mon Dec 12 14:59:07 2022 +0100
@@ -14,7 +14,6 @@
 tidymodels_prefer()
 
 
-source("/home/lorenzo/myprojects-hg/R-codes/stat_lib.R")
 
 
 extract_ms_data <- function( ms, log_trans){
@@ -47,6 +46,43 @@
 }
 
 
+su <- function(x, decreasing=F, ...){
+
+    ## res <- x %>% unique %>% sort(decreasing=decreasing, ...)
+
+    res <- x  |>
+        unique() |>
+        sort(decreasing=decreasing, ...)
+
+    
+    return(res)
+
+    
+}
+
+
+save_excel <- function(output, fileName, sheetName="data", na_yes = TRUE,...){
+tryCatch({
+
+    wb <- loadWorkbook(fileName)
+    addWorksheet(wb = wb, sheet = sheetName)
+  writeData(wb = wb, sheet = sheetName, x = output, colNames = T, rowNames = F,
+       keepNA=na_yes    ,...)
+   saveWorkbook(wb = wb, file = fileName, overwrite = T)
+}, 
+error = function(err){
+
+  wb <- createWorkbook(fileName)
+  addWorksheet(wb = wb, sheet = sheetName)
+    writeData(wb = wb, sheet = sheetName, x = output, colNames = T, rowNames = F, keepNA=na_yes ,  ... )
+  saveWorkbook(wb = wb, file = fileName, overwrite = T)  
+})
+}
+
+
+
+
+
 ###################################################################
 ###################################################################
 ###################################################################
@@ -65,6 +101,24 @@
 
 choose_log <- 1
 
+iso_map_eu27 <- tibble(iso3=c("AUT", "BEL", "BGR", "CYP", "CZE", "DEU", "DNK", "ESP", "EST", "FIN",
+"FRA", 
+"GRC", "HRV", "HUN", "IRL", "ITA", "LTU", "LUX", "LVA", "MLT",
+"NLD", "POL", "PRT",
+"ROM", "SVK", "SVN", "SWE"),
+iso2 = c("AT", "BE",  "BG", "CY", "CZ","DE",
+              "DK",   "ES" ,"EE",  "FI",  "FR",  "EL", "HR", "HU", "IE",
+              "IT" , "LT","LU", "LV",  "MT", "NL", "PL",
+              "PT", "RO", "SK", "SI", "SE"),
+
+country=c("Austria", "Belgium", "Bulgaria", "Cyprus", "Czechia", "Germany",
+  "Denmark", "Spain", "Estonia", "Finland", "France", "Greece",
+  "Croatia", "Hungary", "Ireland","Italy", "Lithuania", "Luxembourg",
+  "Latvia", "Malta", "Netherlands", "Poland", "Portugal",
+  "Romania", "Slovakia", "Slovenia", "Sweden")
+)
+
+
 
 present_month <- month(today())