от а = б * с
Мне нужно найти среднее значение, P10, P50 и P90 "a", происходящие из дистрибутивов "b" и "c", я пытался использовать для этого Монте-Карло, но получаю ошибку.
Я знаю, что это немного глупо, но вы можете мне помочь с этим
Пожалуйста, найдите код R, который я использовал ниже
library(MonteCarlo)
mult<-function(a,b) {sample<-rnorm(b,c)
a<-sample*sample
mean<-mean(a)
return(list("mean"=mean))}
b_grid<-c(50,100,250,500)
c_grid<-seq(1,2,0.2)
param_list=list("b"=b_grid, "c"=c_grid)
MC_result<-MonteCarlo(func=mult, nrep=1000, param_list=param_list)