タグ付けされた質問 「beta-distribution」

間隔で定義された一変量分布の2パラメーターファミリー [0,1]

2
UMVUE
ましょう密度からのランダムサンプルである(X1,X2,…,Xn)(X1,X2,…,Xn)(X_1,X_2,\ldots,X_n)fθ(x)=θxθ−110&lt;x&lt;1,θ&gt;0fθ(x)=θxθ−110&lt;x&lt;1,θ&gt;0f_{\theta}(x)=\theta x^{\theta-1}\mathbf1_{00 のUMVUEを見つけようとしています。θ1+θθ1+θ\frac{\theta}{1+\theta} の結合密度は(X1,…,Xn)(X1,…,Xn)(X_1,\ldots,X_n) fθ(x1,⋯,xn)=θn(∏i=1nxi)θ−110&lt;x1,…,xn&lt;1=exp[(θ−1)∑i=1nlnxi+nlnθ+ln(10&lt;x1,…,xn&lt;1)],θ&gt;0fθ(x1,⋯,xn)=θn(∏i=1nxi)θ−110&lt;x1,…,xn&lt;1=exp⁡[(θ−1)∑i=1nln⁡xi+nln⁡θ+ln⁡(10&lt;x1,…,xn&lt;1)],θ&gt;0\begin{align} f_{\theta}(x_1,\cdots,x_n)&=\theta^n\left(\prod_{i=1}^n x_i\right)^{\theta-1}\mathbf1_{00 \end{align} 母集団pdfは1パラメータ指数ファミリに属しているため、これは完全な十分な統計がfθfθf_{\theta}θθ\thetaT(X1,…,Xn)=∑i=1nlnXiT(X1,…,Xn)=∑i=1nln⁡XiT(X_1,\ldots,X_n)=\sum_{i=1}^n\ln X_i 以降、最初の考えで、の私に与えるUMVUEによってレーマン・シェッフェの定理。この条件付き期待値が直接見つかるか、条件付き分布を見つける必要があるかどうかはわかりません 。E(X1)=θ1+θE(X1)=θ1+θE(X_1)=\frac{\theta}{1+\theta}E(X1∣T)E(X1∣T)E(X_1\mid T)θ1+θθ1+θ\frac{\theta}{1+\theta}X1∣∑ni=1lnXiX1∣∑i=1nln⁡XiX_1\mid \sum_{i=1}^n\ln X_i 一方、私は次のアプローチを検討しました: 我々は、つまり。Xi∼i.i.dBeta(θ,1)⟹−2θlnXi∼i.i.dχ22Xi∼i.i.dBeta(θ,1)⟹−2θln⁡Xi∼i.i.dχ22X_i\stackrel{\text{i.i.d}}{\sim}\text{Beta}(\theta,1)\implies -2\theta\ln X_i\stackrel{\text{i.i.d}}{\sim}\chi^2_2−2θT∼χ22n−2θT∼χ2n2-2\theta\, T\sim\chi^2_{2n} したがって、カイ2乗pdfを使用して計算された次の未加工モーメントはrrr−2θT−2θT-2\theta\,TE(−2θT)r=2rΓ(n+r)Γ(n),n+r&gt;0E(−2θT)r=2rΓ(n+r)Γ(n),n+r&gt;0E(-2\theta\,T)^r=2^r\frac{\Gamma\left(n+r\right)}{\Gamma\left(n\right)}\qquad ,\,n+r>0 したがって、異なる整数の選択に対して、異なる整数の累乗の不偏推定量(およびUMVUE)が得られるようです。たとえば、およびとのUMVUEを直接指定してください。rrrθθ\thetaE(−Tn)=1θE(−Tn)=1θE\left(-\frac{T}{n}\right)=\frac{1}{\theta}E(1−nT)=θE(1−nT)=θE\left(\frac{1-n}{T}\right)=\theta1θ1θ\frac{1}{\theta}θθ\theta さて、我々が持っている。θ&gt;1θ&gt;1\theta>1θ1+θ=(1+1θ)−1=1−1θ+1θ2−1θ3+⋯θ1+θ=(1+1θ)−1=1−1θ+1θ2−1θ3+⋯\frac{\theta}{1+\theta}=\left(1+\frac{1}{\theta}\right)^{-1}=1-\frac{1}{\theta}+\frac{1}{\theta^2}-\frac{1}{\theta^3}+\cdots などのUMVUEを確実に取得できます。したがって、これらのUMVUEを組み合わせると、の必要なUMVUEを取得できます。この方法は有効ですか、それとも最初の方法から続行しますか?UMVUEが存在する場合、UMVUEは一意であるため、どちらも同じ答えを返すはずです。1θ,1θ2,1θ31θ,1θ2,1θ3\frac{1}{\theta},\frac{1}{\theta^2},\frac{1}{\theta^3}θ1+θθ1+θ\frac{\theta}{1+\theta} 明確にするために、E(1+Tn+T2n(n+1)+T3n(n+1)(n+2)+⋯)=1−1θ+1θ2−1θ3+⋯E(1+Tn+T2n(n+1)+T3n(n+1)(n+2)+⋯)=1−1θ+1θ2−1θ3+⋯E\left(1+\frac{T}{n}+\frac{T^2}{n(n+1)}+\frac{T^3}{n(n+1)(n+2)}+\cdots\right)=1-\frac{1}{\theta}+\frac{1}{\theta^2}-\frac{1}{\theta^3}+\cdots つまり、E(∑r=0∞Trn(n+1)...(n+r−1))=θ1+θE(∑r=0∞Trn(n+1)...(n+r−1))=θ1+θE\left(\sum_{r=0}^\infty \frac{T^r}{n(n+1)...(n+r-1)}\right)=\frac{\theta}{1+\theta} 場合、必要なUMVUEがである可能性はありますか?∑r=0∞Trn(n+1)...(n+r−1)∑r=0∞Trn(n+1)...(n+r−1)\displaystyle\sum_{r=0}^\infty \frac{T^r}{n(n+1)...(n+r-1)}θ&gt;1θ&gt;1\theta>1 用、私はなるだろう、及びUMVUEが異なることになるので。0&lt;θ&lt;10&lt;θ&lt;10<\theta<1g(θ )= θ (1 + θ + θ2+ ⋯ )g(θ)=θ(1+θ+θ2+⋯)g(\theta)=\theta(1+\theta+\theta^2+\cdots) 最初のアプローチの条件付き期待値を直接見つけることができなかったと確信しており、、私は先に進みました条件付き分布を検索します。そのため、の結合密度が必要でした。E(X1| Σ LNバツ私= t )= E(X1| Π X私= et)E(X1∣∑ln⁡Xi=t)=E(X1∣∏Xi=et)E(X_1\mid \sum\ln …

2
しきい値処理されたベータ分布を効率的にサンプリング
次の分布から効率的にサンプリングするにはどうすればよいですか? x∼B(α,β), x&gt;kx∼B(α,β), x&gt;k x \sim B(\alpha, \beta),\space x > k が大きすぎない場合、リジェクションサンプリングが最善のアプローチである可能性がありますが、kが大きい場合の処理​​方法がわかりません。おそらく、適用できる漸近近似がありますか?kkkkkk

1
相互に排他的でないカテゴリを分類できる深層学習モデル
例:仕事の説明に「英国のJavaシニアエンジニア」という文があります。 私は2つのカテゴリとして、それを予測することは、深い学習モデルを使用したい:English とIT jobs。従来の分類モデルを使用する場合softmax、最後のレイヤーで機能を持つ1つのラベルのみを予測できます。したがって、2つのモデルのニューラルネットワークを使用して、両方のカテゴリで「はい」/「いいえ」を予測できますが、さらに多くのカテゴリがあると、コストがかかりすぎます。では、2つ以上のカテゴリを同時に予測するためのディープラーニングまたは機械学習モデルはありますか? 「編集」:従来のアプローチによる3つのラベルでは、[1,0,0]によってエンコードされますが、私の場合、[1,1,0]または[1,1,1]によってエンコードされます 例:3つのラベルがあり、文がこれらすべてのラベルに収まる場合。したがって、softmax関数からの出力が[0.45、0.35、0.2]である場合、3つのラベルまたは2つのラベルに分類する必要がありますか、それとも1つにすることができますか?それを行うときの主な問題は、1、2、または3つのラベルに分類するための適切なしきい値は何ですか?
9 machine-learning  deep-learning  natural-language  tensorflow  sampling  distance  non-independent  application  regression  machine-learning  logistic  mixed-model  control-group  crossover  r  multivariate-analysis  ecology  procrustes-analysis  vegan  regression  hypothesis-testing  interpretation  chi-squared  bootstrap  r  bioinformatics  bayesian  exponential  beta-distribution  bernoulli-distribution  conjugate-prior  distributions  bayesian  prior  beta-distribution  covariance  naive-bayes  smoothing  laplace-smoothing  distributions  data-visualization  regression  probit  penalized  estimation  unbiased-estimator  fisher-information  unbalanced-classes  bayesian  model-selection  aic  multiple-regression  cross-validation  regression-coefficients  nonlinear-regression  standardization  naive-bayes  trend  machine-learning  clustering  unsupervised-learning  wilcoxon-mann-whitney  z-score  econometrics  generalized-moments  method-of-moments  machine-learning  conv-neural-network  image-processing  ocr  machine-learning  neural-networks  conv-neural-network  tensorflow  r  logistic  scoring-rules  probability  self-study  pdf  cdf  classification  svm  resampling  forecasting  rms  volatility-forecasting  diebold-mariano  neural-networks  prediction-interval  uncertainty 

4
コックスハザードモデルの生存曲線を解釈するにはどうすればよいですか?
コックス比例ハザードモデルから生存曲線をどのように解釈しますか? このおもちゃの例ではage、kidneyデータの変数にcox比例ハザードモデルがあり、生存曲線を生成するとします。 library(survival) fit &lt;- coxph(Surv(time, status)~age, data=kidney) plot(conf.int="none", survfit(fit)) grid() たとえば、時間、どのステートメントが正しいですか?または両方が間違っていますか?200200200 ステートメント1:被験者は20%残ります(たとえば、人がいる場合、200日目までに、およそ200人が残っているはずです)。 100010001000200200200200200200 ステートメント2:特定の人に対して、彼/彼女は200日目に生存する可能性がます。20%20%20\%200200200 βTxβTx\beta^Tx


1
が独立したベータの場合もベータであることを示します
数年前に私たちの大学の学期試験で出てきた問題を解決しようとしています。 場合独立している密度を有するランダム変数とをそれぞれその表示以下の。X1,X2X1,X2X_1,X_2ββ\betaβ(n1,n2)β(n1,n2)\beta(n_1,n_2)β(n1+12,n2)β(n1+12,n2)\beta(n_1+\dfrac{1}{2},n_2)X1X2−−−−−√X1X2\sqrt{X_1X_2}β(2n1,2n2)β(2n1,2n2)\beta(2n_1,2n_2) ヤコビアン法を使用して、の密度が次のようになることを確認しました: Y=X1X2−−−−−√Y=X1X2Y=\sqrt{X_1X_2}fY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫1y1x2(1−x2)n2−1(1−y2x2)n2−1dxfY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫y11x2(1−x2)n2−1(1−y2x2)n2−1dxf_Y(y)=\dfrac{4y^{2n_1}}{B(n_1,n_2)B(n_1+\dfrac{1}{2},n_2)}\int_y^1\dfrac{1}{x^2}(1-x^2)^{n_2-1}(1-\dfrac{y^2}{x^2})^{n_2-1}dx この時点で私は実際に迷っています。さて、メインの論文で、ヒントが提供されていました。ヒントを使ってみましたが、希望の表現が得られませんでした。ヒントは次のとおり逐語的です。 ヒント:と与えられた密度の観点からの密度の式を導き出し、で変数の変更を使用してみます。Y=X1X2−−−−−√Y=X1X2Y=\sqrt{X_1X_2}X1X1X_1X2X2X_2z=y2xz=y2xz=\dfrac{y^2}{x} したがって、この時点で、この変数の変更を考慮して、このヒントを利用しようとします。したがって、簡略化後、(を書き込む)fY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫yy2z2y4(1−y4z2)n2−1(1−y2.z2y4)n2−1y2z2dzfY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫y2yz2y4(1−y4z2)n2−1(1−y2.z2y4)n2−1y2z2dzf_Y(y)=\dfrac{4y^{2n_1}}{B(n_1,n_2)B(n_1+\dfrac{1}{2},n_2)}\int_{y^2}^y\dfrac{z^2}{y^4}(1-\dfrac{y^4}{z^2})^{n_2-1}(1-y^2.\dfrac{z^2}{y^4})^{n_2-1}\dfrac{y^2}{z^2}dzxxxzzzfY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫yy21y2(1−y4x2)n2−1(1−x2y2)n2−1dxfY(y)=4y2n1B(n1,n2)B(n1+12,n2)∫y2y1y2(1−y4x2)n2−1(1−x2y2)n2−1dxf_Y(y)=\dfrac{4y^{2n_1}}{B(n_1,n_2)B(n_1+\dfrac{1}{2},n_2)}\int_{y^2}^y\dfrac{1}{y^2}(1-\dfrac{y^4}{x^2})^{n_2-1}(1-\dfrac{x^2}{y^2})^{n_2-1}dx どうすればいいのか分かりません。ヒントを適切に解釈しているかどうかさえわかりません。とにかく、残りのヒントを次に示します。 変数の変更を使用することで、平均化することで、必要な密度を2つの方法で表すことができます。今への統合の範囲を分割し、、書き込みおよび。z=y2xz=y2xz=\dfrac{y^2}{x}fY(y)=constant.y2n1−1∫1y2(1−y2x)n2−1(1−x)n2−1(1+yx)1x−−√dxfY(y)=constant.y2n1−1∫y21(1−y2x)n2−1(1−x)n2−1(1+yx)1xdxf_Y(y)=constant.y^{2n_1-1}\int_{y^2}^1(1-\dfrac{y^2}{x})^{n_2-1}(1-x)^{n_2-1}(1+\dfrac{y}{x})\dfrac{1}{\sqrt{x}}dx(y2,y)(y2,y)(y^2,y)(y,1)(y,1)(y,1)(1−y2x)(1−x)=(1−y)2−(yx−−√−x−−√)2(1−y2x)(1−x)=(1−y)2−(yx−x)2(1-\dfrac{y^2}{x})(1-x)=(1-y)^2-(\dfrac{y}{\sqrt{x}}-\sqrt{x})^2u=yx−−√−x−−√u=yx−xu=\dfrac{y}{\sqrt{x}}-\sqrt{x} まあ、正直なところ、私はこれらのヒントの使い方を理解できません。助けていただければ幸いです。前もって感謝します。

2
混合モデルのパラメトリック、セミパラメトリック、ノンパラメトリックブートストラップ
以下の移植片は、この記事から引用したものです。私はブートストラップの初心者であり、R bootパッケージを使用した線形混合モデルのパラメトリック、セミパラメトリック、ノンパラメトリックのブートストラップブートストラップを実装しようとしています。 Rコード これが私のRコードです: library(SASmixed) library(lme4) library(boot) fm1Cult &lt;- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=Cultivation) fixef(fm1Cult) boot.fn &lt;- function(data, indices){ data &lt;- data[indices, ] mod &lt;- lmer(drywt ~ Inoc + Cult + (1|Block) + (1|Cult), data=data) fixef(mod) } set.seed(12345) Out &lt;- boot(data=Cultivation, statistic=boot.fn, R=99) Out ご質問 …
9 r  mixed-model  bootstrap  central-limit-theorem  stable-distribution  time-series  hypothesis-testing  markov-process  r  correlation  categorical-data  association-measure  meta-analysis  r  anova  confidence-interval  lm  r  bayesian  multilevel-analysis  logit  regression  logistic  least-squares  eda  regression  notation  distributions  random-variable  expected-value  distributions  markov-process  hidden-markov-model  r  variance  group-differences  microarray  r  descriptive-statistics  machine-learning  references  r  regression  r  categorical-data  random-forest  data-transformation  data-visualization  interactive-visualization  binomial  beta-distribution  time-series  forecasting  logistic  arima  beta-regression  r  time-series  seasonality  large-data  unevenly-spaced-time-series  correlation  statistical-significance  normalization  population  group-differences  demography 

1
ベータ分布の注文統計
ましょう IIDから引くこと。最小および最大次数統計はそれぞれどのように分散されますか? B e t a (kバツ1、… 、 xんx1,…,xnx_1,\dots,x_nB e t a ( K2、k − p − 12)Beta(k2,k−p−12)Beta\left(\frac{k}2,\frac{k-p-1}{2}\right) できれば参考にしていただければ幸いです。一般的に、私は次数統計の導出に慣れていません。 編集:ベータ分布が均一分布の統計として解釈できることを考えると、ベータ分布の最小値または最大値は別のベータ分布に従って分布していると思います。kkk Edit_2:たまたま気になっている、もう少し正確な設定を追加しました。結局、私は最小値と最大値のテールバウンドを探しているので、どんな形でもこれらにつながるので満足です。私は最終的に漸近的なケースにも興味がありますが、いわば私の次の関心事です。

1
以下のための閉じた形のために、
我々が知っている場合p∼Beta(α,β)p∼Beta(α,β)p \sim Beta(\alpha, \beta)、次いで E[lnp]=ψ(α)−ψ(α+β)E[ln⁡p]=ψ(α)−ψ(α+β) \mathbb{E}[\ln p] = \psi(\alpha) - \psi(\alpha + \beta) ここでψ(.)ψ(.)\psi(.)ディガンマ関数です。\ mathbb {E} [\ ln(1-p)]の簡単な形式はあり E[ln(1−p)]E[ln⁡(1−p)] \mathbb{E}[\ln (1-p)]ますか?

2
HMM Baum-Welchでの後方変数のスケーリング
スケーリングされたBaum-Welchアルゴリズムを実装しようとしていますが、スケーリング後の後方変数が1の値を超えているという問題に遭遇しました。これは正常ですか?結局のところ、確率は1を超えてはなりません。 私はフォワード変数から取得したスケール係数を使用しています: ct=1/∑s∈Sαt(s)ct=1/∑s∈Sαt(s) c_t = 1 / \sum_{s\in S}\alpha_t(s)\\ ここで、c_tは時間tのスケーリング係数、alphaは前方変数、sはhmmの状態です。 後方アルゴリズムについては、以下のJavaで実装しました。 public double[][] backwardAlgo(){ int time = eSequence.size(); double beta[][] = new double[2][time]; // Intialize beta for current time for(int i = 0; i &lt; 2; i++){ beta[i][time-1] = scaler[time-1]; } // Use recursive method to calculate beta double tempBeta …

1
Jeffreysの以前のベータ配布
可能性がベータ分布の形をしていて、そのパラメーターにジェフリーズの事前分布を使用したい場合、事前分布の形式は何ですか? 一部のディストリビューションでは、計算が非常に簡単です。たとえば、2項の場合、2次導関数の期待値は明確に与えます。しかし、可能性自体にすでにベータフォームがある場合、それを導き出そうとして迷った。誰かが私を助けてくれますか? Beta(0.5,0.5)Beta⁡(0.5,0.5)\operatorname{Beta}(0.5, 0.5)

1
コサインカーネルはベータ分布の場合と理解できますか?
Wand and Jones(1995)が述べたように、ほとんどの標準カーネルは、 K(x;p)={22p+1B(p+1,p+1)}−1(1−x2)p1{|x|&lt;1}K(x;p)={22p+1B(p+1,p+1)}−1(1−x2)p1{|x|&lt;1} K(x;p) = \{ 2^{2p+1} \; \mathrm{B}(p+1,p+1) \}^{-1} \; (1-x^2)^p \;\boldsymbol{1}_{\{|x|<1\}} ここで、B(⋅,⋅)B(⋅,⋅)\mathrm{B}(\cdot,\cdot)はベータ関数です。pの異なる値はppp、長方形(p=0p=0p=0)、エパネチニコフ(p=1p=1p=1)、バイウェイト(p=2p=2p=2)、およびトライウェイト(p=3p=3p=3)カーネルにつながります。 余弦カーネル(R density関数で理解できる) 12(1+cos(πx))1{|x|&lt;1}12(1+cos⁡(πx))1{|x|&lt;1} \frac{1}{2} (1 + \cos(\pi x)) \;\boldsymbol{1}_{\{|x|<1\}} この家族の一員としても考えられますか?もしそうなら、それのためのpの適切な値は何pppですか?いくつかのシミュレーションを実行した後、≈2.35≈2.35\approx 2.35はかなり近いと思いますが、(方法)シミュレーションなしで適切なものを見つけるにはどうすればよいですか?そうでない場合、ベータ分布を使用して概算できますか? ワンド、MPおよびジョーンズ、MC(1995)。 カーネル平滑化。 チャップマンとホール、ロンドン。

3
分布
XXXがベータ分布Beta (1,K−1)(1,K−1)(1,K-1)持ち、YYY が2K2K2K度のカイ2乗に従うと仮定します。さらに、XXXとYYYは独立していると仮定します。 製品の分布はどのようなものですかZ=XYZ=XYZ=XY。 私の試みを更新: fZ=∫y=+∞y=−∞1|y|fY(y)fX(zy)dy=∫+∞01B(1,K−1)2KΓ(K)1yyK−1e−y/2(1−z/y)K−2dy=1B(1,K−1)2KΓ(K)∫+∞0e−y/2(y−z)K−2dy=1B(1,K−1)2KΓ(K)[−2K−1e−z/2Γ(K−1,y−z2)]∞0=2K−1B(1,K−1)2KΓ(K)e−z/2Γ(K−1,−z/2)fZ=∫y=−∞y=+∞1|y|fY(y)fX(zy)dy=∫0+∞1B(1,K−1)2KΓ(K)1yyK−1e−y/2(1−z/y)K−2dy=1B(1,K−1)2KΓ(K)∫0+∞e−y/2(y−z)K−2dy=1B(1,K−1)2KΓ(K)[−2K−1e−z/2Γ(K−1,y−z2)]0∞=2K−1B(1,K−1)2KΓ(K)e−z/2Γ(K−1,−z/2)\begin{align} f_Z &= \int_{y=-\infty}^{y=+\infty}\frac{1}{|y|}f_Y(y) f_X \left (\frac{z}{y} \right ) dy \\ &= \int_{0}^{+\infty} \frac{1}{B(1,K-1)2^K \Gamma(K)} \frac{1}{y} y^{K-1} e^{-y/2} (1-z/y)^{K-2} dy \\ &= \frac{1}{B(1,K-1)2^K \Gamma(K)}\int_{0}^{+\infty} e^{-y/2} (y-z)^{K-2} dy \\ &=\frac{1}{B(1,K-1)2^K \Gamma(K)} [-2^{K-1}e^{-z/2}\Gamma(K-1,\frac{y-z}{2})]_0^\infty \\ &= \frac{2^{K-1}}{B(1,K-1)2^K \Gamma(K)} e^{-z/2} \Gamma(K-1,-z/2) \end{align} それが正しいか?はいの場合、どのようにこの分布を呼び出しますか?

1
ベータ分布での異常値の検出
私は値の大きなサンプルを持っていると言う[0,1][0,1][0,1]。基礎となるBeta(α,β)Beta(α,β)\text{Beta}(\alpha, \beta)分布を推定したいと思います。サンプルの大部分は、この想定されるBeta(α,β)Beta(α,β)\text{Beta}(\alpha, \beta)分布からのものですが、残りは、αα\alphaおよび推定で無視したい外れ値ですββ\beta。 これについて進める良い方法は何ですか? なる標準:Inliers={x∈[Q1−1.5IQR,Q3+1.5IQR]}Inliers={x∈[Q1−1.5IQR,Q3+1.5IQR]}\text{Inliers} = \left\{x \in [Q1 - 1.5\, \text{IQR}, Q3 + 1.5 \,\text{IQR}] \right\}箱ひげ図で使用される式は、悪い近似ですか? これを解決するためのより原則的な方法は何でしょうか?この種の問題でうまく機能する、αα\alphaと特定の事前分布はありますββ\betaか?

1
ベータ分布の積
私はトリガーの効率を調べています。つまり、イベントのうちイベントで起動するデバイスがあるということです。結局、ランダムに与えられたイベントで発砲する確率である効率推定に興味があります。を超える均一な事前分布でベイジアンアプローチを使用すると、確率分布をベータ分布としてモデル化できます。kkknnnϵϵ\epsilon[0,1][0,1][0,1]ϵϵ\epsilonβ(ϵ;k+1,n−k+1)β(ϵ;k+1,n−k+1)\beta(\epsilon; k+1, n-k+1) ここで問題が発生します。「ブートストラップ」を使用して効率を計算します。つまり、最終的なトリガー効率は2つのトリガー効率の積であり、どちらもベータ分布としてモデル化できます。 大きな値に対して2つのベータPDFのこの積を計算するにはどうすればよいですかk1,2k1,2k_{1,2}及び効率?閉じた形の製品はありますか(AFAIKはありません)?現時点ではこれを数値的に行っていますが、かなり時間がかかります。n1,2n1,2n_{1,2} この質問には、大きな引数値のベータ分布の積分を評価する方法の答えがありますが、これはここでは役に立ちません。 私の質問が明確で完全に愚かではないことを願っています...

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.