Question: Polyester - creating simulated reads based on raw count matrix
gravatar for sseo2
6 months ago by
sseo20 wrote:

I am trying to generate some simulated counts based on a pilot data we have to determine if increasing the number of samples will give us more power in our experiments.

I have a matrix of raw counts containing gene names on rows and sample names in columns and from this I want to use create_read_numbers function from polyester to generate additional simulated samples. However this comes up with an error and I can't understand what is wrong with it. Any ideas?

The code I am running is as below

WTveh_poly <- Counts[,c("WTVEH1IP","WTVEH2IP","WTVEH3IP")]

WTveh_poly <- subset(WTveh_poly, rowMeans(WTveh_poly) > 1)
params = get_params(WTveh_poly)
custom_readmat = create_read_numbers(mu=params$mu, fit=params$fit,
                                     p0=params$p0, m=dim(WTveh_poly)[1], n=dim(WTveh_poly)[2], 

then the output is

Generating data from baseline model.
Error in, size, replace, prob) : 
  cannot take a sample larger than the population when 'replace = FALSE'

Here is my session info

R version 3.4.1 (2017-06-30)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 14.04.5 LTS
 [1] polyester_1.12.0           ballgown_2.8.4             sva_3.24.4                
 [4] BiocParallel_1.10.1        genefilter_1.58.1          mgcv_1.8-18               
 [7] nlme_3.1-131               BiocInstaller_1.26.0       RColorBrewer_1.1-2        
[10] pheatmap_1.0.8             stringr_1.2.0              dplyr_0.5.0               
[13] DESeq2_1.16.1              SummarizedExperiment_1.6.1 DelayedArray_0.2.2        
[16] matrixStats_0.52.2         Biobase_2.36.2             GenomicRanges_1.28.2      
[19] GenomeInfoDb_1.12.0        IRanges_2.10.1             S4Vectors_0.14.1          
[22] BiocGenerics_0.22.0        readr_1.1.1                biomaRt_2.32.0      


Thanks in advances 

ADD COMMENTlink written 6 months ago by sseo20
Please log in to add an answer.


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