GWA_tutorial
GWA_tutorial copied to clipboard
A comprehensive tutorial about GWAS and PRS
Hi, I am running the codes and I got an error in the following code: ### Step 3 ### # Generate bfile with autosomal SNPs only and delete SNPs with...
When generating a plot to assess the type of relationship in the sample, I do get an R script erro. $ Rscript --no-save Relatedness.R Error in legend(1, 1, xjust =...
Hi, I noticed using the command wget ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz from the 2_Main_script_MDS.txt downloads a corrupted file. Hence creating plink files from the vcf resulted in file read failure error. But when...
I am reproducing the whole tutorial. At , ##Permutation: The reduce computational time we only perform this test on a subset of the SNPs from chromosome 22. `awk '{ if...
Hi Marees, Thanks for the pipeline. I have around 5k samples in pihat_min0.2_genome and 390 in zoom_pihat_genome. If I sort them, I could identify the related individuals but don't know...
Thank you for your tutorial in advance. I have some confusion about the "inversion.txt" in the 1_QC_GWAS folder. I know it is for excluding high inversion regions. However, where is...
In part 2, when we convert from vcf to Plink format, I get an error:  Any advice to get around this? I end up with temporary files because of...
Hi, I am new to GWAS, I have 100+ exome SNV annotated files, i want to create same bim, bed and fam files. Can you/anyone pls help me with how...
I want to practice my genetic data following your tutorial, but I do not know how to create the mydata-based inversion.txt file. Could you please tell me in detail. Thank...
Hi Andries, any plan to have an update to GWA_tutorial? Shicheng