Hi
Im using the kinship package to draw a pedigree. On my data set this works fine
but when i add indivudals to the pedigree i keep getting an error i hope
someone can help me!
This is the code im using:
Data<-read.table("Tree.txt", header=T, sep=",")
attach(Data)
ped<-pedigree(id, dadid, momid, sex, aff)
par(xpd=T)
plot.pedigree(ped)
This is my data looks like without the last 3 individuals it works perfect,when
i add them i get the following error:
Error in if (min(pos2) < 0) { : missing value where TRUE/FALSE needed
famid,id,dadid,momid,sex,aff
1,8860,9972,8856,2,0
1,8855,9972,8856,2,0
1,8859,9976,8860,2,0
1,8854,9971,8855,2,0
1,8863,9975,8859,2,0
1,8858,9975,8859,2,0
1,8865,9975,8859,2,0
1,9969,9970,8854,1,0
1,8864,9980,8863,2,0
1,8862,9980,8863,2,0
1,23834,9981,8865,2,0
1,9968,9969,8853,1,0
1,21141,9974,8858,1,0
1,21142,9974,8858,2,0
1,23172,265,21142,2,0
1,25409,265,21142,1,0
1,23171,265,21142,2,0
1,265,NA,NA,1,0
1,9981,NA,NA,1,0
1,8853,NA,NA,2,0
1,9974,NA,NA,1,0
1,9980,NA,NA,1,0
1,9972,NA,NA,1,0
1,8856,NA,NA,2,0
1,9976,NA,NA,1,0
1,9971,NA,NA,1,0
1,9975,NA,NA,1,0
1,9970,NA,NA,1,0
1,9979,NA,NA,1,0
1,20644,9979,8862,2,0
1,20670,9979,8862,1,0
I tried putting nuclear families closer together in de pedigree but this
doesn't seem to help. All individuals have both parents in the pedigree or have
2 missing parents.
I hope someone can help me with this!
Iris
____________________________________________________________________________________
Never miss a thing. Make Yahoo your home page.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
