BioinfoUNIBA / REDItools

REDItools are python scripts to investigate RNA editing at genomic scale.
MIT License
57 stars 35 forks source link

sample_status_file_creator.py #22

Closed sixingchen718 closed 9 months ago

sixingchen718 commented 9 months ago

Hi, Thank you develop this cool tool. I've completed running REDItoolDnaRna.py on my RNA-seq results, and I want to proceed with differential analysis. However, I'm stuck. Some protocols suggest creating a sample information file and then running the get_DE_events.py program. I tried this but wasn't successful. Others mention using sample_status_file_creator.py to generate the required sample information file (.sif), but I haven't seen this script mentioned in your guidance. Could you please guide me through the correct process? Thank you very much. Best, Sixing

gianfilippo commented 8 months ago

Hi, same issue here. How did you solve it ? Thanks

yr542 commented 8 months ago

@gianfilippo Did you resolve your issue with the sif file. I have been have a lot of issues with get_DE_events.py and I figured out how the sample informations file should look but I cannot seem to get the get_DE_events.py so much so I have had to convert to python3 and then am still debugging (whether in python 2 or 3 it seems to have quite a few issues. If you haven't figured out the sample informations file format (SIF) it seems that CTRL or DIS are required in a Status column and the sample names in the Sample column. The scaffold of what I used is below:

Sample,Status
SRRExperimental-01,DIS
SRRExperimental-02,DIS
SRRExperimental-03,DIS
SRRControl-01,CTRL
SRRControl,CTRL
SRRControl,CTRL

Note(s):

The notes are just what I have observed however, I just cannot seem to get the get_DE_events.py to work.

gianfilippo commented 8 months ago

Hi, no, sorry. Still stuck

yr542 commented 8 months ago

I have been working on/editing the get_DE_events script as it seems to have quite a few requirements and might have a couple of bugs.

Requirements Seem To Be:

Scaffold:

SRRsampleName/editing/outTable_12345678

This is all I have so far but I do not know if it helps @gianfilippo as sadly I still cannot seem to get it to work.

gianfilippo commented 8 months ago

thanks!!! I will give it a try as soon as I go back to this project. Will keep you posted.