error trying to fit a curve (method 2) with coancestry and euclidean distances

[MarinaChi]

Hi!

When I try to fit the IBD trend to a curve (step 2, method 2) using euclidean or coancestry distances, I get this error message:

Error in nls(y ~ cbind(1, 1 - exp(-exp(lrc) * x)), data = xy, start = list(lrc = log(-coef(lm(log(abs(y - :
step factor 0.000488281 reduced below 'minFactor' of 0.000976562

Does it mean that a curve can't be properly fitted to my data? or is it a bug?

Thank you!

Marina

[qt37t247]

Dear Marina,

Sorry for the late response.

In many cases, the non-linear model may not fit your data. In such cases, please use linear model for the IBD by specifying "IBD_method = 1" in the function rdm_IBD().

Best regards,

Tang Qian

[MarinaChi]

Thank you so much!!!
Best
Marina