имитировать lapply / цикл по списку параметров

Question

Я хотел бы смоделировать частоту и серьезность по списку параметров.Вот случай для первого элемента в списке:

data <- data.frame(
lamda = c(5, 2, 3),
meanlog = c(9, 10, 11), 
sdlog = c(2, 2.1, 2.2))

freq <- rpois(s, data$lamda[1])

freqsev <- lapply(freq, function(k) rlnorm(k, data$meanlog[1], sdlog = data$sdlog[1]))

freq
freqsev

Как настроить цикл или оператор lapply для перебора всех элементов в данных?(не только первый).

Спасибо.

Answer 1

Мы можем использовать map (из пакета purrr, часть пакета tidyverse) следующим образом для создания столбцов списка.Содержимое теперь хранится в столбцах freq и freqsev.

library(tidyverse)

set.seed(123)

s <- 2

data2 <- data %>%
  mutate(freq = map(lamda, ~rpois(s, .x)),
         freqsev = map(freq, ~map(.x, function(k) rlnorm(k, meanlog, sdlog))))

data2$freq
# [[1]]
# [1] 4 7
# 
# [[2]]
# [1] 2 4
# 
# [[3]]
# [1] 6 0

data2$freqsev
# [[1]]
# [[1]][[1]]
# [1]    9330.247   28897.323 2605520.369   20370.283
# 
# [[1]][[2]]
# [1]    645.4047   5206.2183  22461.1778  93729.0634  46892.3129 144595.7492  10110.8606
# 
# 
# [[2]]
# [[2]][[1]]
# [1]   2665.955 938950.074
# 
# [[2]][[2]]
# [1]  21931.9763    354.2858 280122.6952   3147.6681
# 
# 
# [[3]]
# [[3]][[1]]
# [1]   957.5257 13936.3063  6265.3530  1886.0077  5927.8540  1464.5081
# 
# [[3]][[2]]
# numeric(0)

Обновление

Вот способ замены значений, больших чем равных 500.

data3 <- data2 %>%
  mutate(capat500 = map(freqsev, ~map(.x, function(y) ifelse(y >= 500, 500, y))))