GCbias discrepancy

[Glfrey]

Hi there,

I seem to be encountering a problem with polyester where an error message reports a discrepancy between the length of the gc bias and the sum of numreps "length(extras$gcbias) == sum(num_reps) is not TRUE", however I can't seem to correct this error as R reports they are indeed equal "length(biases) == sum(num_reps)[1] TRUE".

My commands are:

fastapath = ('/path/to/fastafile') <- has 10 sequences within
numtx = count_transcripts(fastapath) <- [1] 10
transcripts = readDNAStringSet(fastapath) <- "A DNAStringSet instance of length 10"
readmat= matrix(readspertx, nrow=numtx, ncol=10)
set.seed(137)
biases = sample(0:7, 10, replace=TRUE)
readmat_biased = add_gc_bias(readmat, as.list(biases), transcripts)

confirming they're equal:
sum(num_reps) - [1] 10
length(biases) - [1] 10

Then I run

simulate_experiment('/path/to/fastafile', reads_per_transcript=readspertx, num_reps=c(5,5), fold_changes = fold_changes, distr = 'empirical',outdir='bias_gc_sim', gcbias = 'biases')

Which returns the error message "length(extras$gcbias) == sum(num_reps) is not TRUE".

I'd be very grateful into any insight into the issues.

Gill