removedSnpFemaleChrXhweControl | R Documentation |
Hardy Weinberg Equilibrium test for SNPs on the chromosome X in female controls.
removedSnpFemaleChrXhweControl( plink, inputPrefix, pval = 1e-06, outputPvalFile, outputSNPfile, outputPrefix )
plink |
an executable program in either the current working directory or somewhere in the command path. |
inputPrefix |
the prefix of the input PLINK binary files. |
pval |
the p-value cutoff for controlling HWE test in female control subjects. Only chromosome X SNPs are considered. The default value is 0.000001. |
outputPvalFile |
the output pure text file that stores chromosome X SNPs and their sorted HWE p-values. |
outputSNPfile |
the output pure text file that stores the removed SNPs, one per line. |
outputPrefix |
the prefix of the output PLINK binary files. |
The output PLINK binary files after HWE test on chromosomal X in female controls.
Junfang Chen
removedSnpHWEauto
## In the current working directory bedFile <- system.file("extdata", "genoUpdatedData.bed", package="Gimpute") bimFile <- system.file("extdata", "genoUpdatedData.bim", package="Gimpute") famFile <- system.file("extdata", "genoUpdatedData.fam", package="Gimpute") system(paste0("scp ", bedFile, bimFile, famFile, " .")) inputPrefix <- "genoUpdatedData" outputPvalFile <- "2_12_snpHwePvalfemaleXct.txt" outputSNPfile <- "2_12_snpRemovedHweFemaleXct.txt" outputPrefix <- "2_12_removedSnpHweFemaleX" ## Not run: Requires an executable program PLINK, e.g. ## plink <- "/home/tools/plink" ## removedSnpFemaleChrXhweControl(plink, inputPrefix, pval=0.000001, ## outputPvalFile, outputSNPfile, ## outputPrefix)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.