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)。
この関数を byexponentialFunc(10,0.5,10000)と呼び、m の一連の値を返します。しかし、出力
を使用したときに0.95の確率を持つこの分布の分位数を見つけたいので、quantile(m,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