qcdhRoutine: QCDH routine 2
In kindlychung/CollapsABEL: An R package for quantitative CDH analysis

This routine generates arguments to be passed to PLINK and task name automatically, then pass these to routine 1.

1 2	qcdhRoutine(wDir = ".", pheno, pheno_name, covar = NULL, covar_name, plinkArgs, initGwas = FALSE, initGwasArgs, pFilter, nMaxShift)

`wDir`	Working directory. All paths within this function will be relative to this.
`pheno`	Phenotype file (usually also containing covariates)
`pheno_name`	Name of phenotype in the phenotype file
`covar`	Covariates file, default to the same as `pheno`
`covar_name`	Name of covariates in the covariates file
`plinkArgs`	Arguments to be passed to PLINK (in QCDH analysis)
`initGwas`	Whether to perform an initial GWAS
`initGwasArgs`	Arguments to pass to PLINK for the initial GWAS
`pFilter`	Filter out SNPs with p value higher than this (in the initial GWAS)
`nMaxShift`	Maximum shift number (for genotype collapsing)