注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

wangyufeng的博客

祝愿BB 健康开心快乐每一天

 
 
 

日志

 
 

R/fGWAS: Functional GWAS package for R (Version 1.0)  

2013-04-13 10:10:43|  分类: GWAS |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

Description:

The R/fGWAS Package (Functional Genome-wide Association Studies) is developed as a new method for genome-wide association studies based on a single SNP analysis [1]. The proceduce called by Bayesian Lasso model includes two steps, which is characterized by the production of a preconditioned reponse variable by a supervised principle component analysis and the formulation of Bayesian lasso for selecting a subset of significant SNPs. Based on this theoretical model, the package integrates two steps into one function to analyse the genotypical data and phenotypical data. The analysis flow can be started by the data loading, and then the data goes through the main function of the Bayesian Lasso model. Finally the results will be generated. By the result summarization, the information about significant SNP positions can be reported.

For more details, please refer to the document.

[1] Jiahan Li, Kiranmoy Das, Guifang fu, Runze Li and Rongling Wu. The Bayesian Lasso for Genome-wide Associations Studies. Bioinformatics (2010). First published online: December 14, 2010. doi: 10.1093/bioinformatics/btq688

Download:


The packages for Windows, Linux and Mac OS are available. The document also can be downloaded by the following link.


Installation:


The R/fGWAS package depends on 3 packages, including mvtnorm, glmnet and bigmemory. Before R/fGWAS is installed, these packages should be installed in advance. After you download this package file, please type the following command or click the menu item "Install packages from local zip files".

Windows OS:

>install.packages("x:/fullpath/fgwas_1.0-1.zip", repos=NULL)

Linux/Mac OS:

>install.packages("/fullpath/ fgwas_1.0-1.tar", repos=NULL)

Before the package is used in R, It is necessary to loading package by the following command:

> library(fgwas)

After it is loaded, all functions within R/fGWAS will be readily available to the user.

Sample Script:


The following source shows how to call the main function in R.

dat1 <- LS.load_plink("bmi_chr1_t1000.tped", "bmi_chr1_t1000.tfam", "bmi_phenos.csv", 1 );
res1 <- LS.run(dat1);

dat2 <- LS.load_plink("bmi_chr2_t1000.tped", "bmi_chr2_t1000.tfam", "bmi_phenos.csv", 2 );
res2 <- LS.run(dat2);

ret <- LS.merge_result(res1, res2);
summary(ret);


Sample Figures:


R/fGWAS: Functional GWAS package for R (Version 1.0) - 喜欢吃桃子 - wangyufeng的博客

Last updated:


Version 1.0 01/19/2011
  评论这张
 
阅读(1343)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017