# Consider the knowledge space and the parameters used in Brancaccio,
# de Chiusole, Stefanutti (2023) in Example 1
states<-matrix(c( 0,0,0,0,0,
0,0,0,0,1,
0,0,1,0,1,
0,0,0,1,1,
0,0,1,1,1,
1,0,1,0,1,
0,1,0,1,1,
1,0,1,1,1,
0,1,1,1,1,
1,1,0,1,1,
1,1,1,1,1), byrow=TRUE, ncol=5)
beta <-c(.004,.03,.02,.01,.007)
eta <-c(5e-06, 5e-05, 4e-05,.007,.08)
likelihood <-c(0,0,0,0,0,1/3,0,1/3,0,0,1/3)
Q_pool <- c(2,4,5)
half_split(likelihood,states,beta,eta,Q_pool)
Run the code above in your browser using DataLab