Could I please ask two questions, both related to SNPs.
I have .bim/.bed file for a set of SNPs. I want to remove SNPs from my data set that have a >2% Mendelian error rate.
./plink --bfile IN_FILE --me 0.02 0.02 --noweb --make-bed --out X.
But the output tells me 0 SNPs were removed. I am working with 500K,100K and 50K arrays and I would have expected at least one SNP to be removed. I messed around with increasing and decreasing the -me 0.02, 0.02 parameters, and nothing is ever removed.
Could someone tell me the correct command to remove SNPs from a data set that have >2% Mendelian error rate?
I have a set of SNPs (it's an Affymetrix 100K). How do I tell whether a SNP is on the plus or minus strand?
For example, I have a set of SNPs. The information I have for each SNP is:
This is an example of a SNP (let's call it SNPx):
And I have a set of genes and their location in the genome, for example:
EntrezID Chr # GeneStart GeneEnd Strand
1 12 744000 9067900 minus
4 12 744001 130887675 plus
111 3 123282296 123449077 minus
142 1 226360691 226408100 minus
185 3 148697871 148743003 plus
This is entrez ID, chromosome number, their start and end position in the genome, and whether they are on the plus/minus strand.
I need to find out whether SNPx lies within any of these genes.
I am confused as to how to tell whether the SNP is on the plus or minus strand?
(1) , in this case, would SNPx belong to Entrez gene ID 1, or 4? They both are on the same chromosome in about the same position, just one is on the plus strand and one is on the minus strand?
(2) Do I need to account for ss to rs orientation, and rs to genome orientation? How should I do this, if so?