From 2fecc9f2bc8ca6d6ad2476e6847168e906e240e6 Mon Sep 17 00:00:00 2001 From: Gaddis Date: Tue, 10 Mar 2020 11:44:46 -0400 Subject: [PATCH] [#3] Updated QC workflow README --- genotype_array_qc/README.md | 258 +++++++++++++++++++++++++++++++++++- 1 file changed, 257 insertions(+), 1 deletion(-) diff --git a/genotype_array_qc/README.md b/genotype_array_qc/README.md index 5005752..69ea5ee 100644 --- a/genotype_array_qc/README.md +++ b/genotype_array_qc/README.md @@ -107,6 +107,26 @@ Parameters:
3. Remove duplicate IDs (based on call rate) +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
@@ -135,7 +155,6 @@ Input Files: | FILE | DESCRIPTION | | --- | --- | -| `[PREFIX_LIST]` | List of prefixes of files to be merged | Output Files: @@ -197,66 +216,303 @@ Parameters:
7. Format phenotype data to standard format +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
8. Structure workflow (separate supporting workflow) +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
9. Partition data by ancestry +Sample command: +``` shell +plink \ + --bfile [INPUT_BED_BIM_FAM_PREFIX] \ \ + --keep [KEEP_LIST] \ + --make-bed \ + --out [OUTPUT_BED_BIM_FAM_PREFIX] +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | +| `[INPUT_BED_BIM_FAM_PREFIX].bed` | PLINK format bed file for input genotypes | +| `[INPUT_BED_BIM_FAM_PREFIX].bim` | PLINK format bim file for input genotypes | +| `[INPUT_BED_BIM_FAM_PREFIX].fam` | PLINK format fam file for input genotypes | +| `[KEEP_LIST]` | List of subjects to keep | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | +| `[OUTPUT_BED_BIM_FAM_PREFIX].bed` | PLINK format bed file for output genotypes | +| `[OUTPUT_BED_BIM_FAM_PREFIX].bim` | PLINK format bim file for output genotypes | +| `[OUTPUT_BED_BIM_FAM_PREFIX].fam` | PLINK format fam file for output genotypes | +| `[OUTPUT_BED_BIM_FAM_PREFIX].log` | PLINK log file | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- | +| `--bfile [INPUT_BED_BIM_FAM_PREFIX]` | Prefix for input genotypes in PLINK bed/bim/fam format | +| `--keep [KEEP_LIST]` | List of subjects to keep | +| `--make-bed` | Flag indicating to generate genotypes in PLINK bed/bim/fam format | +| `--out [OUTPUT_BED_BIM_FAM_PREFIX]` | Prefix for output genotypes in PLINK bed/bim/fam format |
10. Call rate filter +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
11. HWE filter +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
12. Subject call rate filter (based on autosomes) +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
13. Relatedness workflow (separate supporting workflow) +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
14. Remove samples based on relatedness +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
15. Sex check and sample removal +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
16. Excessive homozygosity filtering +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |
17. Set het haploids to missing +Sample command: +``` shell +``` + +Input Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Output Files: + +| FILE | DESCRIPTION | +| --- | --- | + + +Parameters: + +| PARAMETER | DESCRIPTION | +| --- | --- |