myrna generates all p-values = 1.000000e+00

[Adam Witney]

I am trying to get myrna working on our bacterial RNA-seq data.

I first modified a reftools shell script to access EnsemblBacteria for the Streptococcus genome and this seems to work ok.

After running myrna all the directories and files seem to be created but all genes have a p-value of 1.000000e+00 in the pvals.txt file.

The dataset has no replicates so only has one fastq for each of 2 conditions, but I wouldn't expect all p-values to be exactly 1. Do i need more replicates or is something not running correctly within the process?

thanks

adam

[Ben Langmead]

Hi Adam,

Can you send me the log of your run? (Myrna prints out a message at the end indicating what the output-dump file is named.) Also, can you send me the count_table.txt file from the run?

Thanks,
Ben

I am trying to get myrna working on our bacterial RNA-seq data.

I first modified a reftools shell script to access EnsemblBacteria for the Streptococcus genome and this seems to work ok.

After running myrna all the directories and files seem to be created but all genes have a p-value of 1.000000e+00 in the pvals.txt file.

The dataset has no replicates so only has one fastq for each of 2 conditions, but I wouldn't expect all p-values to be exactly 1. Do i need more replicates or is something not running correctly within the process?

thanks

adam

[Rachelly]

Hi,

I'm facing a similar problem, but I get all the p-values as 0.000000e+00, which is also unlikely..
Any conclusions on this issue?

Thanks,
Rachelly.

[Adam Witney]

In my case the p-value's of 1.000000e+00 were being generated as there were no replicates in this dataset.

I don't know enough about the system to answer your question, but maybe Ben can help you.

Adam

[csmatyi]

Same problem for us over here! All of our p-values come out as 1.