mc2d를 사용하여 나무를 건축하는 데 오류가 있습니다.

Question

나는 각 팔의 끝 부분에있는 값과 함께 각 팔을 추락 할 가능성이있는 팔 세트로 결국 바뀔 나무를 만들고 있습니다.

이것은 내게 마지막 나무 팔다. (나는 뒤를 만들고있다.) mcstoc (...) 나는 그것을 실행하는 나에게 다음과 같은 오류를 제공합니다 - 라인이 <을 probs1 그러나

node1 <- mcprobtree(probs1, list("1"=mcdata(0), "2"=mcdata(1), type="V") 

가 :

는

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)) 

다음 단계는 다음과 같이 트리의 분기를 구축하는 것입니다

Error in function (n, values, prob = NULL) : 
Prob and values should be of same length or have the same number of columns. 

나는이 방정식의 확률값의 = C를 (probs_d_if_fail, 테스트) 부분으로 mcstoc 객체를 사용하여 저를 좋아하지 않는 느낌을 얻을 수 있지만,이 문제를 해결하는 방법을 잘 모르겠습니다.

도움을 주시면 감사하겠습니다. 디모데

Answer 1

rempiricalD

는 prob 인자로 다변량 mcnode이 필요합니다. 이것은 까다 롭습니다.

시도 :

require(mc2d) 
ndvar(10000) 

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 
probs <- mcdata(c(probs_d_if_fail,test),type="V",nvariates=2) 
probs1 <- mcstoc(rempiricalD, values=1:2, type="V", prob= probs) 

node1 <- mcprobtree(probs1, list("1"=mcdata(0), "2"=mcdata(1)), type="V") 
node1 

그러나 사용자의 요구에 맞는 경우 다음 코드는 간단 할 수 있습니다 rbern 확률 (1-probs_d_if_fail)와 1을 제공 할 것입니다

require(mc2d) 
ndvar(10000) 

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(rbeta,type="V",shape1=alpha_dfail, 
shape2=beta_dfail) 

node1 <- mcstoc(rbern, type="V",prob = 1-probs_d_if_fail) 
node1 

하는 것으로 당신으로, 내가 잘 이해할 필요가있다.

RP