Skip to content

Instantly share code, notes, and snippets.

View irudnyts's full-sized avatar
🏠
Working from home

Iegor Rudnytskyi irudnyts

🏠
Working from home
View GitHub Profile
while (any(teamR$name == teamR$target)) {
teamR <- teamR %>%
mutate(target = sample(name))
}
library(tidyverse)
teamR <- tribble(
~name, ~email,
"Bobby", "xxx@epfl.ch",
"Enrico", "xxx@epfl.ch",
"Harry", "xxx@epfl.ch",
"Giulia", "xxx@epfl.ch",
"Daniel", "xxx@epfl.ch",
"Iegor", "xxx@epfl.ch"
@irudnyts
irudnyts / set-up.R
Last active December 18, 2020 10:42
set-up
# install.packages("blastula")
library(blastula)
create_smtp_creds_key(
id = "my_gmail",
user = "iegor.rudnytskyi@gmail.com",
provider = "gmail"
)
@irudnyts
irudnyts / tr2.R
Last active November 16, 2020 14:40
glimpse(fatalities)
library(rvest)
library(tidyverse)
library(janitor)
library(lsr)
library(scales)
library(patchwork)
import::from(magrittr, extract)
wiki_page <- str_c(
@irudnyts
irudnyts / remove_keys.R
Created December 6, 2019 11:11
Remove `blastula` keys from a keyring
library(dplyr)
library(stringr)
keys <- keyring::key_list() %>%
filter(str_detect(string = service, pattern = "^blastula-v1"))
apply(
keys,
1,
function(key) keyring::key_delete(service = key[1], username = key[2])
library(tidyverse)
probs <- readRDS(file = "/Users/irudnyts/Downloads/probs.rds")
(p1 <- ggplot(probs) +
geom_ribbon(
aes(x = time_horizon,
ymin = lower_bound,
ymax = upper_bound,
fill = factor(initial_capital)),
# Step 1
library(ruin)
library(actuar)
library(ggplot2)
probs <- data.frame(
initial_capital = rep(0:10, each = 50),
time_horizon = rep(1:50, times = 11),
lower_bound = NA,
# Step 1
library(ruin)
library(actuar)
library(ggplot2)
probs <- data.frame(
initial_capital = rep(0:10, each = 50),
time_horizon = rep(1:50, times = 11),
lower_bound = NA,
library(ruin)
library(actuar)
model0 <- CramerLundbergExtended(
initial_capital = 0,
premium_rate = 1,
claim_poisson_arrival_rate = 1,
claim_size_mixing_parameter = 0.01,
claim_size_light_tail_distribution = rexp,
claim_size_light_tail_parameters = list(rate = 1),