Changelog -> NEWS.md
Fixed “uninitialized variable” warning in C code.
Small bug fix re || vs | in function for gamma probabilities
Use markdown in documentation
Small bug fix in simStahl.c re over-running max_nxo.
Small changes to avoid Note about “R_registerRoutines”
Remove calls to missing() and use NULL argument defaults instead.
Added some error messages for fitStahl and stahlLoglik, if input is not a list.
Modified find.breaks to handle an F2 intercross.
Revised stahlLoglik to handle the case of an intercross.
stahlLoglik returns just the vector of log likelihoods (rather than a data.frame with the parameters as well). The parameter values are saved as attributes of the result.
Removed use of R_zeroin() by incorporating code within the package as the function Rxoi_zeroin().
Reduced computation time for examples in several help files.
Added NAMESPACE.
Compressed bssbsb data.
Fixed a bug in find.breaks (changing call to locate.xo to locateXO).
Fixed links in help files.
Added functions stahlLoglik and fitStahl, for calculating the log likelihood under the Stahl model and for obtaining MLEs and the log LR for testing p=0.
Revised convertxoloc() so that the output is a data frame.
In est.recrate, modified the default positions at which the recombination rate is calculated (that is, if the “pos” argument is not specified) so that calculations will be done at the marker positions and at a grid with 4 positions per Mbp.
Added a function recrate2scanone for converting est.recrate results for multiple chromosomes to the form of the output of scanone (from R/qtl), but I’ll leave it undocumented for now.
Added a function joint.given.two for calculating the joint density of the crossover locations, given that there are two.
Added a function est.coi, for estimating the three-point coincidence along a chromosome.
Added function est.intensity and est.kfunc to estimate the intensity and K_inhom functions for an inhomogeneous point process.
Revised most of the functions for calculating densities and such so that you may specify (via the argument “x”) the exact values at which these functions are to be calculated. If x is unspecified, we calculate things at n equally spaced points between 0 and L.