- Fixed a regression error involving the external program MINX (part of MERLIN). The new version of
`checkMerlin()`

by default checks that both`merlin`

and`minx`

are available on the system, and that they both come from the latest MERLIN version.

`pedprobr`

now depends on R 4.1 and`pedtools`

1.1.0.The deprecated argument

`loop_breakers`

has been completely removed (renamed to`loopBreakers`

)

`likelihood.ped()`

gains a new argument`lump`

, which activates allele lumping. This is TRUE by default, and should remain so in most cases.Using theta correction in pedigrees with inbred founders now gives an error.

The internal

`peelingProcess()`

now checks for unbroken loops in the pedigree. (Of interest for developers mainly.)

- Fixed bug affecting likelihood computations with theta correction.

When calling MERLIN, the window-centric and unnecessary

`.exe`

extensions has been removed.Many examples have been rewritten in simpler code, taking advantage of

`pedtools::addMarker`

and the pipe`|>`

.

- Removing/disabling mutation models is now easier, with
`setMutationModel(..., model = NULL)`

.

- Fixed bug in
`setMutationModel()`

affecting lists of multiple pedigrees.

`likelihood()`

has been refactored, moving the treatment of (two) linked markers to a separate function,`likelihood2()`

.

Theta-correction is implemented in likelihood calculations, through the argument

`theta`

of`likelihood()`

.`merlin()`

and`likelihoodMerlin()`

have been overhauled.A new function

`checkMerlin()`

checks if MERLIN is installed and available.`merlin()`

now has an argument`linkageMap`

facilitating analysis of linked markers.`likelihoodMerlin()`

gains arguments`rho`

and`logbase`

.New function

`lumpAlleles()`

.

The

`likelihood()`

function now handles multiple marker inputs. For example, the call`likelihood(x, 1:2)`

results in a vector of length 2 with the likelihoods of the first two markers attached to`x`

.Recombination parameter

`theta`

is renamed to`rho`

everywhere, to align with other ped suite packages.`theta`

still works, though.

Many input checks have been added, giving more useful error messages.

The data structure used in likelihood calculations is simplified and more efficient.

Likelihood computation of two linked markers (and therefore also

`twoMarkerDistribution()`

) is much faster now.

New general MERLIN wrapper,

`merlin()`

.Both

`likelihoodMerlin()`

and the new`merlin()`

now accepts ped lists.

- Fixed a bug affecting lumped mutation models.

- Initial CRAN release.