I'm wondering if someone could please help.
I have a population of individuals, some are in families and some are not in families. I have a binary phenotype. I want to do an association analysis between phenotype and genotype, controlling for population stratification and family structures in the population (and other co-variates that I have in the --covar
file).
The command I'm using is:
plink --bfile myfile --logistic --ci 0.95 --out TestPlink --noweb --covar CovariateFile --covar-name SEX,AGE --hide-covar --model-rec
First, I'm checking for population stratification using the command:
plink --bfile MyFile --ibs-test --out TestPlink --noweb
This is currently running.
Second, I want to account for the family structure of my population. I know for sure that my population is made up of families (i.e. it was a community based sample of maybe 300 families + individuals not in families). I've tried to read thru the doc, and also other papers that use PLINK, but I'm just really confused on what the command is that says "look in my population for evidence of family clustering, and account for this family clustering when doing the association analysis". For example, do I use permutation "controlling for identified substructure or familial relationships by permuting only within cluster" , multidimensional scaling plots or the DFAM procedure in PLINK (implements the sib-TDT and also allows for unrelated individuals to be included (via a clustered-analysis using the Cochran-Mantel-Haesnzel)).
Would the correct command for "look in my population for evidence of family clustering, and account for the fact that there are families in my population when doing the association analysis" be something like this:
plink --bfile myfile --logistic --ci 0.95 --out TestPlink --noweb --covar CovariateFile --covar-name SEX,AGE --hide-covar --model-rec --dfam
And then would the command to also control for population stratification, should there be any be:
plink --bfile myfile --logistic --ci 0.95 --out TestPlink --noweb --covar CovariateFile --covar-name SEX,AGE --hide-covar --model-rec --dfam --cluster # (--mds-plot 4 ???)
Many thanks