Revised reduce_markers() so that it can handle the case of many markers, by working with them in smaller batches.
fit1() now returns both fitted values and residuals.
fit1() can be run with genotype probabilities omitted, in which case an intercept column of 1's is used (Issue #151).
Updated mouse gene database with 2020-09-07 data from MGI.
Implemented Issue #184, to make calc_het() multi-core.
Made the vdiffr package optional: only test the plots locally, and only if vdiffr is installed.
calc_sdp() can now take a plain vector (Issue #142).
Added a lodcolumn argument to maxlod() (Issue #137).