Dear freinds
I want to do LD analysis on my data (Illumina BovinSNP50 BeadChip) by plink
I have two problem :
1- the pedigree of my data are very crowded and I don't know if I put zero for parental and maternal column in ped file, i can trust to my result or not. because, as you know in livestock population we have a many relationship and it affected the LD extent and pattern
2-in Ped file, what can I put in phenotype column? I know if I put zero I can use --no-pheno in plink. but i have EBV for many traits such as, EBV for milk production, protein percentage and ....
Is there any one who can help me?
1 answer
Hi,
1- I remember that once I had to hack a few lines in the assoc.cpp file of PLINK, so I had a look in the code (Plink::correlation2SNP() funtions, called from Plink::calcLDStatistics()) and found the following lines to keep only founders in the analysis:
// Iterate over every individual but only consider founders
... {
// Allow for 1 or both SNPs to be non-autosomal
for (int i=0; i < n; i++) {
Individual * person = sample[i];
// Only consider founders
if ( ! person->founder ) continue;
-> yes, I suppose that parental and maternal informations really matter !
2- on the contrary i suppose you do not have to care about the pheotype column in the LD analysis, just put zero for every sample. The --no-pheno option is needed only if you skip the column.
Log in to answer this question.