BioinfoUNIBA / REDItools

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

Traceback From `get_DE_events.py` dtype='object')] are in the [columns] #23

Open yr542 opened 8 months ago

yr542 commented 8 months ago

I am trying to use get_DE_events.py however, I see to run across an issue using it. I have the following created sif:

Sample,Group,Type
Control-01,GROUPA,Control
Control-02,GROUPA,Control
Control-03,GROUPA,Control
EXPERIMENTAL-01,GROUPB,EXPERIMENTAL
EXPERIMENTAL-02,GROUPB,EXPERIMENTAL
EXPERIMENTAL-03,GROUPB,EXPERIMENTAL 

I get the following error (I removed the paths as I am working on an HPC and they require that I do not show paths):

Traceback (most recent call last):
 File /path/to/get_DE_events.py", line 447, in <module>
 call_differential_editing_sites(samples_informations_file) 
 File /path/to/get_DE_events.py", line 132, in call_differential_editing_sites
 edit_level_table = editing_table[all_people]
 File “/path/to/conda/for/EnvironWithPython2.7/pandas/frame.py", line 2934, in __getitem__
 raise_missing=True)
 File “/path/to/conda/for/EnvrionWithPython2.7/pandas/indexing.py", line 1354, in _convert_to_indexer
 return self._get_listlike_indexer(obj, axis, **kwargs)[1]
 File “/path/to/conda/for/EnvironWithPython2.7/pandas/indexing.py", line 1161, in _get_listlike_indexer
 raise_missing=raise_missing)
 File “/path/to/conda/for/EnvironWithPython2.7/pandas/indexing.py", line 1246, in _validate_read_indexer
 key=key, axis=self.obj._get_axis_name(axis)))
KeyError: u"None of [Index([u'Control-01',\n  u'Control-02',\n  u'Control-03',\n  u'EXPERIMENTAL-01',\n  u'EXPERIMENTAL-02',\n  u'EXPERIMENTAL-03'],
\n  dtype='object')] are in the [columns]"

My script is:

python /path/to/get_DE_events.py \
-c 10 \
-cpval 2 \
 -input_file /path/to/sif/file/Control_vs_EXPERIMENTAL.sif \
-f 0.1 \
-mtsA 50.0 \
-mtsB 50.0 \
-sig yes \
-siglevel 0.05 \
-linear

Interestingly a user on Reddit had virtually the same issue found here it seems to be difficult to solve. How would I get around this error.