Posts Tagged ‘ PLINK ’

Prune GWAS data in R

March 29, 2011
By
Prune GWAS data in R

Hansong Wang, our biostats professor here at the Hawaii Cancer Center, generously gave me some R code that goes through a SNP annotation file (i.e. a mapfile) and selects SNPs that are at least a certain specified distance apart. You might want to do t...

Read more »

Prune GWAS data in R

March 29, 2011
By
Prune GWAS data in R

Hansong Wang, our biostats professor here at the Hawaii Cancer Center, generously gave me some R code that goes through a SNP annotation file (i.e. a mapfile) and selects SNPs that are at least a certain specified distance apart. You might want to do t...

Read more »

Efficient Mixed-Model Association eXpedited (EMMAX) to Simutaneously Account for Relatedness and Stratification in Genome-Wide Association Studies

June 9, 2010
By

A few months ago I covered an algorithm called EMMA (Efficient Mixed-Model Association) implemented in R for simultaneously correct for both population stratification and relatedness in an association study. This method/software is very useful because ...

Read more »

Happy Birthday GGD! The 10 Most Popular Posts Since GGD’s Launch

February 23, 2010
By

The first post on Getting Genetics Done was one year ago today. To celebrate, here are the top 10 most viewed posts since GGD launched last year. Incidentally, nine of the ten are tutorials on how to do something in R. Thanks to all the readers and all...

Read more »

LocusZoom: Plot regional association results from GWAS

February 10, 2010
By

Update Friday, May 14, 2010: See this newer post on LocusZoom. If you caught Cristen Willer's seminar here a few weeks ago you saw several beautiful figures in the style of a manhattan plot, but zoomed in around a region of interest, with several ot...

Read more »

GWAS Manhattan plots and QQ plots using ggplot2 in R

January 20, 2010
By

Will posted earlier this week about how to produce manhattan plots of GWAS results using Stata, so I thought I'd share how I do this in R using ggplot2. First, if you've never used ggplot2, you'll need to add it to your R installation by typing: install.packages("ggplot2") Once you've done that, copy and paste this command to download the functions I wrote...

Read more »

Visualizing sample relatedness in a GWAS using PLINK and R

October 9, 2009
By

Strict quality control procedures are extremely important for any genome-wide association study.  One of the first steps you should take when running QC on your GWAS is to look for related samples in your dataset.  This does two things for you.  First, you can get an idea of how many related samples you have in your dataset, and second,...

Read more »