Compute correlation for initializing the active set
cor_init(m0, mb)
m0
modFUN(formula, data)
mb
modFUN(mandatory, data)