Open cdeboever3 opened 5 years ago
I started doing this in the Process Plink Results
notebook. I think Julia was processing these plink files a bit in her compare_cutoff_pvals
notebook as well. I'm trying to reconcile the number of diseases and alleles left after filtering based on disease/allele frequency right now before making a filtered version of the file.
I think the folder /oak/stanford/groups/mrivas/users/jolivier/repos/hla-assoc/scripts/output/process_output/
should have files of the form _processed.txt
, which should have the odds ratio, p value, and log odds ratio from both R and PLINK (PLINK would be the first columns, not labeled with (R))
On Feb 3, 2019, 7:54 AM +0530, Christopher DeBoever notifications@github.com, wrote:
I started doing this in the Process Plink Results notebook. I think Julia was processing these plink files a bit in her compare_cutoff_pvals notebook as well. I'm trying to reconcile the number of diseases and alleles left after filtering based on disease/allele frequency right now before making a filtered version of the file. — You are receiving this because you were assigned. Reply to this email directly, view it on GitHub, or mute the thread.
Combine the
.hybrid
files in/oak/stanford/groups/mrivas/users/jolivier/repos/hla-assoc/data/PLINK_results
into one file. We only need rows where theTEST
column has valueADD
and we only need the following columns:You'll need to a column with the phenotype code. Make the file a tsv and gzip it before adding to repo.