R问题的蒙特卡罗模拟

R问题的蒙特卡罗模拟,r,simulation,montecarlo,R,Simulation,Montecarlo,如何解决此练习: 假设U来自均匀0,1分布,Z来自标准正态分布。设X=pZ+1-pU。当X的方差为0.4时,估计p 也许我应该安装“polynom”软件包 p等于0,62,因为我已经分析计算过了,但我不知道如何使用随机数生成器用R来解它,我不会为你解这个练习,但会给你解决它的线索 你需要像实习医生那样做: 模拟U[0,1] 模拟N0,1 使用p的值将它们放在一起 对于给定的p,您可以定义一个函数: 第234集 computeX这里是一个例子 fobj <- function(p) abs

如何解决此练习:

假设U来自均匀0,1分布,Z来自标准正态分布。设X=pZ+1-pU。当X的方差为0.4时,估计p

也许我应该安装“polynom”软件包


p等于0,62,因为我已经分析计算过了,但我不知道如何使用随机数生成器用R来解它,我不会为你解这个练习,但会给你解决它的线索

你需要像实习医生那样做:

模拟U[0,1] 模拟N0,1 使用p的值将它们放在一起

对于给定的p,您可以定义一个函数:

第234集 computeX这里是一个例子

fobj <- function(p)  abs(var(p*rnorm(1e5) + (1-p)*runif(1e5))-0.4)
pval <- optimise(fobj,c(0,1))$minimum

这不是标准问题。需要澄清的是,p等于0,62,因为我已经用分析的方法计算过了,但我不知道如何用随机数生成器用R来求解
> pval
[1] 0.6223968