r - 通过 R 中的随机数生成查找分布

Question

X与参数呈指数分布lambda=0.5。
我想找Pr(1/mean(X)< K)=0.95 因为我不知道的分布1/mean(X)我通过下面的代码生成一个分布。

m<-c()
exponentialFunc<-function(n,lambda,nsim){
    for(i in seq(nsim)){
    x<-rexp(n,lambda)
    y<-mean(x)
    m<-c(m,1/y)
    }
    m

}

现在m是分布1/mean(x)。
现在我调用了这个函数exponentialFunc(10,0.5,10000)，它给出了一系列 m 的值。
但是当我想找到这个分布的分位数时，当我使用quantile(m,0.95)输出时它的概率为 0.95

> quantile(m,0.95) 95% NA
为什么会产生NA？我的代码有什么问题吗？

score 2 · Accepted Answer

似乎您没有将函数的输出保存在名为的对象中m：

> m <- exponentialFunc(10,0.5,10000)
> quantile(m,0.95)
     95% 
0.939102

您正在调用您的函数并且没有将其保存在任何地方，看来您正在做的事情（基本上）是这样的：

> m<-c()
> quantile(m,0.95)
95% 
 NA

1 回答 1