data(daPe)
ma <- glm(grade ~ gpa + tuce + psi, x = TRUE,
data = daPe, family = binomial(link = "probit"))
summary(ma)
(ea <- maBina(w = ma, x.mean = TRUE, rev.dum = TRUE))
(ta <- maTrend(q = ea, nam.c = "gpa", nam.d= "psi"))
plot(ta)
Run the code above in your browser using DataLab