- currently only show effect scan, or plotting phenotype vs genotype (grabbed from `maxmarg()`) - should also show how to use `pull_genoprobpos()` and `fit1()`. - how to make the effect plot? Use `ciplot()` from the [broman package](https://github.com/kbroman/broman)?