Question: LOD Score NaN Error
1
gravatar for Abdul Rafay Khan
13 days ago by
Karachi, PK
Abdul Rafay Khan1000 wrote:
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
R paramlink score lod • 63 views
ADD COMMENTlink modified 13 days ago by zx87547.3k • written 13 days ago by Abdul Rafay Khan1000
1

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 REPLYlink written 13 days ago by zx87547.3k

but my data seems correct

ADD REPLYlink written 6 days ago by Abdul Rafay Khan1000
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1043 users visited in the last hour