title: multiple GLMs in a for loop
body: <p>I am looking to estimate parameters for a large set (274) of correlated response variables which follow a NB dist. The goal is to use the parameters for a generalized linear model for each of the variables with a two level categorical predictor. I am familiar with the MASS package in R, but there doesn't seem to be a way to do a multivariate analyses. Is there a way this can be achieved (for instance a for loop) in R.</p>
  
  <p>Thanks in advance for any help/suggestions</p>