Skip to content

rgt-TDF AttributeError: 'Statistics' object has no attribute 'tpx_ndef' #292

Description

@SophieEhres

Hello, I was trying to run rgt-TDF and got the following error:

TDF - Regulatory Analysis Toolbox (RGT). Version: 1.0.2

(rgt_env) SAU08:triplex ehresms$ rgt-TDF promotertest -r FENDRR_human_transcripts.fa -de DEG_FENDRR_Names.txt -organism hg38 -rn FENDRR_names -o /Users/ehresms/computational/FENDRR_triplex_promotertest -showdbs

*************** Promoter Test ****************
*** Input RNA sequence: /Users/ehresms/SauvageauLab Dropbox/SauvageauLab/Members/EHRESMANN_Sophie/Projects/FENDRR/triplex/FENDRR_human_transcripts.fa
*** Output directory: computational/FENDRR_triplex_promotertest/FENDRR_names
Step 1: Calculate the triplex forming sites on RNA and DNA.
Dump to file: tdf.dump.hg38.DEG_FENDRR_Names
Running time: 0:04:54
Step 2: Calculate the frequency of DNA binding sites within the promoters.
DEBUG: file_tpx_nde is /Users/ehresms/computational/FENDRR_triplex_promotertest/FENDRR_names/nontarget_promoters.tpx
Counting frequency of promoters on DBD...
620 Binding target promoters
DEBUG: Loaded nontarget_promoters.tpx
34816 Binding non-target promoters
Counting frequency of binding sites on DBD...
620 Binding sites on de promoters
Traceback (most recent call last):
File "/opt/anaconda3/envs/rgt_env/bin/rgt-TDF", line 11, in
sys.exit(main())
File "/opt/anaconda3/envs/rgt_env/lib/python3.7/site-packages/rgt/tdf/Main.py", line 356, in main
file_tpx_de=tpx_de, file_tpx_nde=tpx_nde)
File "/opt/anaconda3/envs/rgt_env/lib/python3.7/site-packages/rgt/tdf/Statistics.py", line 138, in count_frequency_promoters
numdbs_ndef = len(self.tpx_ndef.get_dbs(rm_duplicate=True))
AttributeError: 'Statistics' object has no attribute 'tpx_ndef'

#########

The error in Statistics.py:

  if self.pars.showdbs:
            self.frequency["hits"] = {"de": OrderedDict(), "nde": OrderedDict()}
            numdbs_def = len(self.tpx_def.get_dbs(rm_duplicate=True))
            numdbs_ndef = len(self._tpx_ndef_.get_dbs(rm_duplicate=True))
            for rbs in self.rbss:
                # DE
                l1 = len(self.tpx_def.merged_dict[rbs])
                self.frequency["hits"]["de"][rbs] = [l1, numdbs_def - l1]
                # non-DE
                l2 = len(self._tpx_ndef_.merged_dict[rbs])
                self.frequency["hits"]["nde"][rbs] = [l2, numdbs_ndef - l2]

tpx_ndef should be tpx_nde

thanks!

Sophie

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions