The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: Using PEER to identify hidden confounders in RNA-seq data (not working)
gravatar for rodrigo.duarte88
23 days ago by
rodrigo.duarte880 wrote:

Hi all,

I am trying to use PEER in R (on conda) to identify hidden confounders in the raw counts of my RNA-seq experiment. I am following their tutorial ( ), but without much luck, and there's not much help available unfortunately. Anyone used this tool before and could advise?

I have installed R-PEER on Conda, but am having issues with following with the PEER pipeline, to identify 10 hidden confounders in my expression data. I was wondering if anyone could advise? I couldn't find much help online.

> $ source activate r-peer (r-peer) 
> $ R  
> expr <- read.table("~/raw_counts.tsv", header=TRUE)
> dim(expr) 
>[1] 14211   538      # 14211 row-genes, 538 col-samples 
> expr <- t(expr)    # transpose 
> model = PEER()
> PEER_setPhenoMean(model,as.matrix(expr))
> [1] NA                      # in the tutorial, it says this should be "NULL" 
> PEER_setNk(model,10)  # so PEER identifies 10 hidden confounders 
>PEER_getNk(model) [1] 10
> PEER_update(model)
>         iteration 0/1000
>         iteration 1/1000 
>Converged (bound) after 1 iterations NULL     ###### why isn't it iterating?
> sessionInfo() R version 3.4.1 (2017-06-30) Platform:
> x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 18.04.1 LTS
> Matrix products: default BLAS:
> /home/rodrigo/miniconda2/envs/r-peer/lib/R/lib/ LAPACK:
> /home/rodrigo/miniconda2/envs/r-peer/lib/R/lib/
> locale:  [1] LC_CTYPE=C.UTF-8       LC_NUMERIC=C          
> LC_MESSAGES=C.UTF-8  [7] LC_PAPER=C.UTF-8       LC_NAME=C             
> attached base packages: [1] stats     graphics  grDevices utils    
> datasets  methods   base
> other attached packages: [1] peer_1.0
> loaded via a namespace (and not attached): [1] compiler_3.4.1
ADD COMMENTlink modified 23 days ago by Martin Morgan ♦♦ 22k • written 23 days ago by rodrigo.duarte880

PEER isn't a Bioconductor package so you'll have more luck directing questions to the maintainer or other resource identified in packageDescription("PEER")

ADD REPLYlink written 23 days ago by Martin Morgan ♦♦ 22k

Oops, sorry about that! To be fair, I think am going to try "sva", looks like there's more documentation for that one. Thanks, Martin!

ADD REPLYlink written 23 days ago by rodrigo.duarte880
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: 194 users visited in the last hour