mirror of
https://github.com/Zenithsiz/ist-pe-proj.git
synced 2026-02-03 22:23:56 +00:00
39 lines
873 B
R
Executable File
39 lines
873 B
R
Executable File
#!/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")
|