data("weight_behavior")
x=weight_behavior[,2:14]
y=weight_behavior[,15]
temp.b.b.glm<-mma(x,y,pred=2,contmed=c(8:10,12:13),binmed=c(7,11),binref=c(1,1),
catmed=6,catref=1,predref="M",alpha=0.4,alpha2=0.4, jointm=NULL,margin=1, n=2,
seed=1,mart=FALSE,nu=0.001,D=3,distn="bernoulli",
family1=binomial(link = "logit"),n2=2)
plot(temp.b.b.glm)
Run the code above in your browser using DataLab