# NOT RUN {
tmat = matrix(c(0,0,1,0,0,0,0,1,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0),
5,5, byrow= TRUE ) #for "X" (5 nodes) shaped tree
states = c("P","N") #"P" represent cases(or positive) and "N" represent controls(or negative)
hmmA = initHMM(states,list(c("L","R")), tmat) #one feature with two discrete levels "L" and "R"
hmmB = initHMM(states, list(c("X","Y")),tmat, c(0.5,0.5), matrix(c(0.7,0.3,0.3,0.7),2,2))
# }
Run the code above in your browser using DataLab