Version with sample-specific catchability variables

This version incorporates a new design matrix (Q_ik) that contains variables affecting catchability but not underlying density. The model then controls for these variables and predicts abundance for the reference value (Q_ik = 0). This is necessary to account for difference in catchability within season, between passes, etc.