search

magnusdv / pedprobr

Pedigree probabilities in R
GNU General Public License v2.0
4 stars 0 forks source link

`likelihoodMerlin` problem #13

Closed thoree closed 4 years ago

thoree commented 4 years ago

I don't understand the error (I updated from GitHub this morning):

library(pedprobr)
#> Loading required package: pedtools
x = nuclearPed(4)
# Two identical markers
m1 = marker(x, geno = c("a/a", "a/b", "a/b", "a/b", "a/b", "a/b"))
# Attach markers
m2 = m1
x = setMarkers(x, list(m1, m2))
likelihoodMerlin(x, markers = 1:2, rho = 0, verbose = F)
#> Error: Replacement vector larger than the number of markers: 1, 1

Created on 2020-08-17 by the reprex package (v0.3.0)

magnusdv commented 4 years ago

My mistake, this is due to a mistake in pedtools which I fixed but forgot to commit. Should be ok now.