Add tests for benchmarking family against old reference implementation. Thus, we need a reference implementation stored in a data file, and a reference data set (or better a function for simulating data of arbitrary dimension).

• test log-density
• test scores
• test loglik after backfitting