This is a test version of Biostars. For the public version, visit https://www.biostars.org.
LOD Score NaN Error
library(paramlink) 

y = read.table("markers.txt", header=TRUE) 
x = linkdat(y, model=1)

lod(x)

Result:

LOD:    NaN

Data:

FAMID   ID  FID MID SEX AFF M_A1 M_A2
1   1   0   0   1   1   2   3
1   2   0   0   2   1   1   2
1   3   1   2   2   2   2   2
1   4   1   2   1   2   2   3
1   5   1   2   2   1   1   2
1   6   1   2   2   1   1   2
1   7   1   2   2   1   1   2
1   8   1   2   1   1   1   2
lod score r paramlink

From the manual:

If a marker has incompatible values (i.e. if a child of homozygous 1/1 parents has a 2 allele), the corresponding output entries are NaN.

0 answers

No answers yet.

Log in to answer this question.