sarkka branch: getting the dispersion matrix right
simpler code: Q computation following Sarkka (removed abasset confusing matrix)
merge sarkka
sorry...
better sanitation for logit_ab + threshold for n_failure is now a function of number of time series
merged master