私は、最終的に各腕を下る確率が付いた腕のセットになり、各腕の端に値があるツリーを作成しています。
これは私の木の最後の腕です(私はそれを後方に構築しています)。
require(mc2d)
ndvar(1000)
sderror_dfail<-(0.30-0.86)/1.96
alpha_dfail<-0.62*{{0.62*{1-0.62}/{sderror_dfail*sderror_dfail}}-1}
beta_dfail<-alpha_dfail*{{1-0.62}/0.62}
probs_d_if_fail<-mcstoc(rbetagen,type="V",shape1=alpha_dfail, shape2=beta_dfail,
min=0, max=1)
lots1<-mcdata(1)
test<-lots1-probs_d_if_fail
probs1 <- mcstoc(rempiricalD, values=1:2, type="V", prob= c(probs_d_if_fail,test))
次のステップは、次のようにツリーのブランチを構築することです。
node1 <- mcprobtree(probs1, list("1"=mcdata(0), "2"=mcdata(1), type="V")
ただし、行probs1 <-mcstoc(....)を実行すると、次のエラーが発生します。
Error in function (n, values, prob = NULL) :
Prob and values should be of same length or have the same number of columns.
方程式のprob=c(probs_d_if_fail、test)部分としてmcstocオブジェクトを使用するのは気に入らないように感じますが、これを回避する方法がわかりません。
どんな助けでも大歓迎です。
ティモシー