commit 2f1a3be417f5d9e8482642c2026f1147470156b1
Author: Filipe Rodrigues <filipejacintorodrigues1@gmail.com>
Date:   Sun Jun 4 01:07:39 2023 +0100

    Initial commit.

diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..166edb4
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1,2 @@
+# Output files
+output.*
diff --git a/.vscode/settings.json b/.vscode/settings.json
new file mode 100644
index 0000000..0d78ddc
--- /dev/null
+++ b/.vscode/settings.json
@@ -0,0 +1,3 @@
+{
+  "cSpell.words": ["ddesemp", "ggplot", "ggsave", "readxl"]
+}
diff --git a/1/econ.xlsx b/1/econ.xlsx
new file mode 100644
index 0000000..2929c7e
Binary files /dev/null and b/1/econ.xlsx differ
diff --git a/1/script.R b/1/script.R
new file mode 100755
index 0000000..07a3191
--- /dev/null
+++ b/1/script.R
@@ -0,0 +1,38 @@
+#!/bin/env Rscript
+
+# Libraries
+library("readxl")
+library("ggplot2")
+
+
+# Read the data
+df <- read_xlsx("econ.xlsx")
+
+
+# Filter for years >= 1971
+df <- df[df$tempo >= "1971/01/01", ]
+
+
+# Create a new data frame with adjusted `ddesemp` and `tpp` columns
+ddesemp_mean <- mean(df$ddesemp)
+ddesemp_sd <- sd(df$ddesemp)
+
+tpp_mean <- mean(df$tpp)
+tpp_sd <- sd(df$tpp)
+
+df_adjusted <- data.frame(
+    tempo = lapply(df[, c("tempo")], function(tempo) tempo),
+    ddesemp = lapply(df[, c("ddesemp")], function(ddesemp) {
+        (ddesemp - ddesemp_mean) / ddesemp_sd
+    }),
+    tpp = lapply(df[, c("tpp")], function(tpp) {
+        (tpp - tpp_mean) / tpp_sd
+    })
+)
+
+#  Then plot them and save the output
+plot <- ggplot(df_adjusted, aes(x = tempo)) +
+    geom_line(aes(y = ddesemp, color = "ddesemp")) +
+    geom_line(aes(y = tpp, color = "tpp"))
+
+ggsave(plot = plot, "output.png")