run error

16 views
Skip to first unread message

er qw

unread,
May 12, 2025, 4:41:28 PMMay 12
to picrust-users
Rscript /home/tools/picrust2-2.6.2/picrust2/Rscripts/castor_nsti.R qiime_realsc/arc.tre /tmp/tmpgk8m80hl/known_tips.txt /tmp/tmpgk8m80hl/nsti_out.txt

Rscript /home/tools/picrust2-2.6.2/picrust2/Rscripts/castor_hsp.R qiime_realsc/arc.tre /tmp/tmpgtrf8_3n/subset_tab_0 mp 0.5 FALSE FALSE /tmp/tmp1rehla6h/predicted_counts.txt /tmp/tmp1rehla6h/predicted_ci.txt 100


Finished getting marker and NSTI predictions for both domains: qiime_realsc/bac_marker_predicted_and_nsti.tsv.gz, qiime_realsc/arc_marker_predicted_and_nsti.tsv.gz
Now finding the best one for each sequence.
Picking the best domains for all sequences.
4902 sequences were kept for bac
1 sequences were kept for arc
Traceback (most recent call last):
  File "/home/tools/anaconda3/envs/picrust2-sc/bin/picrust2_pipeline.py", line 312, in <module>
    main()
  File "/home/tools/anaconda3/envs/picrust2-sc/bin/picrust2_pipeline.py", line 272, in main
    func_outfiles, pathway_outfiles = full_pipeline_split(study_fasta=args.study_fasta,
  File "/home/tools/picrust2-2.6.2/picrust2/pipeline.py", line 306, in full_pipeline_split
    prune_tree(list(nsti_dom1_df.index.values)+ref_d1_seqs, out_tree_ref1, out_tree_ref1_red)
  File "/home/tools/picrust2-2.6.2/picrust2/util.py", line 496, in prune_tree
    tree.prune(names)
  File "/home/tools/anaconda3/envs/picrust2-sc/lib/python3.10/site-packages/ete3/coretype/tree.py", line 531, in prune
    to_keep = set(_translate_nodes(self, *nodes))
  File "/home/tools/anaconda3/envs/picrust2-sc/lib/python3.10/site-packages/ete3/coretype/tree.py", line 2607, in _translate_nodes
    raise TreeError("Invalid target node: "+str(n))
ete3.coretype.tree.TreeError: 'Invalid target node: 180227'

this is some ERROR information, my command is 
picrust2_pipeline.py -s rep-seqs_single/dna-sequences.fasta -i table_single/otu_table_single.tsv -o qiime_realsc -p 180 --verbose

but when i run 
picrust2_pipeline_oldIMG.py -s rep-seqs_single/dna-sequences.fasta -i table_single/otu_table_single.tsv -o qiime_realsc -p 180 --verbose
it is OK, no ERROR, that's why?
otu_table_single.tsv
dna-sequences.fasta
Reply all
Reply to author
Forward
0 new messages