Передача векторных и одиночных значений для функции с pmap

Question

Как мои ноги мокрые с purrr.

Я оцениваю различные комбинации параметров для бета-распределения, и я пытаюсь передать вектор значений x и параметры для распределения.

Смотрите код ниже:

library(tidyr) 
library(purrr) 
library(dplyr) 

x <- seq(0.01, 0.99, 0.01) 
alpha <- c(0.1, 0.5, 1:5, 10) 
beta <- c(0.5, 1, 2, 5) 


df <- expand.grid(alpha = alpha, beta = beta, x = x) %>% 
    nest(x) %>% 
    mutate(plotdata = pmap(list(data$x, alpha, beta), dbeta)) %>% 
    unnest() 

Что такое правильный синтаксис для ртар в этом случае?

Answer 1

Давайте рассмотрим первый случай в качестве примера.

x = 0.01 
alpha = 0.1 
beta = 0.5 

dbeta(shape1 = alpha, shape2 = beta, x = x) 
# [1] 5.600381 

Второй пример был бы:

x = 0.01 
alpha = 0.5 
beta = 0.5 

dbeta(shape1 = alpha, shape2 = beta, x = x) 
# [1] 3.199135 

Если это вывод, который вы ожидаете, то мы могли бы обобщать:

x <- seq(0.01, 0.99, 0.01) 
alpha <- c(0.1, 0.5, 1:5, 10) 
beta <- c(0.5, 1, 2, 5) 

df <- expand.grid(alpha = alpha, beta = beta, x = x) 

args_beta <- with(df, list(x = x, shape1 = alpha, shape2 = beta)) 
args_beta %>% 
    pmap(dbeta) %>% 
    unlist() %>% 
    as.data.frame() %>% 
    cbind(df, .) 

Тогда результат будет:

alpha beta x   . 
1 0.1 0.5 0.01 5.600381e+00 
2 0.5 0.5 0.01 3.199135e+00 
3 1.0 0.5 0.01 5.025189e-01 
4 2.0 0.5 0.01 7.537784e-03 
5 3.0 0.5 0.01 9.422230e-05 
6 4.0 0.5 0.01 1.099260e-06