Reference Distributions

library(CPTtools)

Model setup

skill1 <- c("H","M","L")
troph <- c("gold","silver","none")

sDist <- calcDNTable(list(),skill1,numeric(),0,log(.5))
round(sDist,3)

##          H     M     L
## [1,] 0.194 0.611 0.194

r <- .8
se <- sqrt(1-r^2)
rT <- calcDNTable(list(s1=skill1),skill1,log(r),0,log(se))
round(rT,3)

##          H     M     L
## [1,] 0.716 0.261 0.022
## [2,] 0.236 0.527 0.236
## [3,] 0.022 0.261 0.716

Psx <- calcDPCTable(list(s1=skill1),troph,log(1),
                    list(gold=1,silver=-1))
round(Psx,3)

##       gold silver  none
## [1,] 0.478  0.505 0.018
## [2,] 0.134  0.732 0.134
## [3,] 0.018  0.505 0.478

simTable <- function(n,Ps,Pss,Psy) {
  
}