Fit MARSS models as usual, but use MARSS_tmb() instead of MARSS(). Not R and Q can only be diagonal, unconstrained or fixed. equalvarcov is not available nor are other custom Q matrices available with the EM algorithm in MARSS(). But TMB is very fast.