#Chapter 7
#7.3
#multivariate predictor
x=weight_behavior[,c(2:3,5:14)]
pred=weight_behavior[,4]
y=weight_behavior[,15]
data.mb.b <-data.org(x,y,mediator=5:12,jointm=list(n=1,j1=c(5,7,9)),
                     pred=pred,predref="OTHER", alpha=0.4,alpha2=0.4)
summary(data.mb.b)

#multivariate responses
x=weight_behavior[,c(2:3,5:14)]
pred=weight_behavior[,4]
y=weight_behavior[,c(1,15)]
data.mb.mb<-data.org(x,y,mediator=5:12,jointm=list(n=1,j1=c(5,7,9)),
                     pred=pred,predref="OTHER", alpha=0.4,alpha2=0.4)
summary(data.mb.mb)