# NOT RUN {
hp <- mhp(M=diag(2),B=array(c(diag(3),diag(3)),c(3,3,2)),
names=letters[1:3],levels=c("oak","ash"))
M(hp)
B(hp)[1,1,1] <- 30 # try a negative value and see what happens
names(hp)
names(hp) <- c("Alice","Zachy","Annabel")
levels(hp) <- c("squid","snail")
summary(hp)
# }
Run the code above in your browser using DataLab