hfNames <- array(rep("Exponential", 36), dim = c(6,6))
hfNames[3,4:5] <- rep("impossible",2)
hfNames[1:4,6] <- rep("sr.fun", 4)
hfNames[2, 4:5] <- rep("treat.fun", 2)
hfNames[col(hfNames)<=row(hfNames)]<-"NULL"
M <- makeM(hfNames)
Run the code above in your browser using DataLab