How to simulate a count matrix using Polyester
1
0
Entering edit mode
@fagholizadeh89-12402
Last seen 5.0 years ago

Hi,

I'm gonna simulate some RNA-seq read counts data from a real data set for publication purposes. My problem is I can't download big Fasta files. So I want to obtain the count matrix using galaxy and then simulate data using that count matrix. There are some R packages such as SimSeq and Polyester which simulate data. SimSeq can simulate data by getting as input a count matrix which is what I'm looking for. But the problem is I want to simulate a time course study and SimSeq can't accommodate dependence over time. Technically, it looks at each time point as an independent treatment. Polyester can handle time course studies but can I give it a count matrix as input instead of big Fasta and GTF files? I'm a master student in biostatistics and I'm very new to simulation. So I don't know how to write the codes on my own.

Thanks a lot.

simulation count matrix RNA-seq Polyester • 1.3k views
ADD COMMENT
0
Entering edit mode
Alyssa Frazee ▴ 210
@alyssa-frazee-6710
Last seen 3.4 years ago
San Francisco, CA, USA

Yes, you can use a count matrix as input for polyester! Please see the "Simulating Reads" section of the vignette (subsection = "simulate_experiment_countmat example", which is available in the package vignette or at https://github.com/alyssafrazee/polyester#simulate_experiment_countmat-example. ?simulate_experiment_countmat will give you the details of how to use the function. Good luck!

ADD COMMENT
0
Entering edit mode

Thank you very much. I try to apply it. But from the manual, it says that the output is the fasta files ("create FASTA files containing RNA-seq reads simulated from provided transcripts, with optional differential expression between two groups (designated via read count matrix)"). does it mean that I have to obtain the count matrix using those fasta files? Or can I have the count matrix directly by running some other codes in the package?

ADD REPLY

Login before adding your answer.

Traffic: 609 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6