LOD Score NaN Error
0
1
Entering edit mode
4.9 years ago
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 • 1.2k views
ADD COMMENT
1
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

but my data seems correct

ADD REPLY

Login before adding your answer.

Traffic: 1771 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6