Question: GWAS: access Imputation probabilities: snpStats
gravatar for GENOMIC_region
16 months ago by
GENOMIC_region0 wrote:

Hi there,

I'm interested to work with .gen file generated from IMPUTE2 imputation tool. In particular I need probabilities per SNP. I use snpStats' function:

test_impute<-read.impute("genotyped_snps_chr22.impute", rownames = NULL, nsnp = NULL, snpcol = 2)

But the output I get is snpMatrix, which I fail to unpack and access. Is there a way I can see first 10 rows and column from this matrix? 

attributes(test_impute) This prints all SNP ids on screen

None of these help. 

Is there a better way that could allow me to access these probabilities for my downstream analysis? 

I'm looking at gtx package's read.snpdata.impute function - but this is awfully slow. I've 13K snps in the input file.


snpStats_1.28.0 Matrix_1.2-14   gtx_0.0.8       survival_2.42-3 R version 3.4.2 (2017-09-28)
Platform: x86_64-pc-linux-gnu (64-bit)



snpstats gwas imputation • 217 views
ADD COMMENTlink modified 16 months ago • written 16 months ago by GENOMIC_region0

I haven't really worked with the snpStats package or SnpMatrix but you might try 
as(test_impute, "character") 
as(test_impute, "numeric")
ADD REPLYlink modified 16 months ago • written 16 months ago by shepherl ♦♦ 1.7k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 322 users visited in the last hour