recording covariates in separate matrix to be bound later

changed sd_avg fxn to record 0 if a species mass was not recorded.
Minimize data loss, shouldn't alter results. Need to denote
appropriately in methods since its not a TRUE zero.