vcontact2:病毒聚类(失败)

ops/2024/10/18 21:22:10/

Bitbucket

安装 

mamba create --name vContact2 biopython=1.78
mamba install -c bioconda vcontact2=0.11.3vim ~/envs/vContact2/lib/python3.9/site-packages/vcontact2/exports/summaries.py
把 np.warnings.filterwarnings('ignore')
改成
import warnings
warnings.filterwarnings('ignore')## 下载聚类软件,移动到conda/bin路径
wget -c http://www.paccanarolab.org/static_content/clusterone/cluster_one-1.0.jar## 下载MCL
mamba install -c bioconda mcl

使用

#使用自带函数构建[gene-to-genome mapping file]
vcontact2_gene2genome -p (faa文件) -o (输出文件) -s Prodigal-FAA#正式运行
运行失败报错
----------------------------------Loading data----------------------------------
INFO:vcontact2: Read 422588 entries (dropped 3534 singletons) from fetal_vc2/vConTACT_profiles.csv
Traceback (most recent call last):File "/home/zhongpei/mambaforge/envs/vContact2/bin/vcontact2", line 834, in <module>main(options)File "/home/zhongpei/mambaforge/envs/vContact2/bin/vcontact2", line 602, in mainmatrix, singletons = vcontact2.pcprofiles.build_pc_matrices(profiles, contigs_csv_df, pcs_csv_df)File "/home/zhongpei/mambaforge/envs/vContact2/lib/python3.9/site-packages/vcontact2/pcprofiles.py", line 358, in build_pc_matricesmatrix = sparse.coo_matrix(([1]*len(profiles), (zip(*profiles.values))), shape=(len(contigs), len(pcs)),File "/home/zhongpei/mambaforge/envs/vContact2/lib/python3.9/site-packages/scipy/sparse/_coo.py", line 99, in __init__self._check()File "/home/zhongpei/mambaforge/envs/vContact2/lib/python3.9/site-packages/scipy/sparse/_coo.py", line 188, in _checkraise ValueError('mismatching number of index arrays for shape; '
ValueError: mismatching number of index arrays for shape; got 0, expected 2
算了不用它了


http://www.ppmy.cn/ops/20897.html

相关文章

Redis分布式锁 - 基于Jedis和LUA的分布式锁

先基于单机模式&#xff0c;基于Jedis手工造轮子实现自己的分布式锁。 首先看两个命令&#xff1a; Redis 分布式锁机制&#xff0c;主要借助 setnx 和 expire 两个命令完成。 setnx命令: setnx 是 set if not exists 的简写。将 key 的值设为 value &#xff0c;当且仅当…

Llama-7b-Chinese本地推理

Llama-7b-Chinese 本地推理 基础环境信息&#xff08;wsl2安装Ubuntu22.04 miniconda&#xff09; 使用miniconda搭建环境 (base) :~$ conda create --name Llama-7b-Chinese python3.10 Channels:- defaults Platform: linux-64 Collecting package metadata (repodata.js…

模拟LinkedList实现的双向链表

1. 前言 前文我们用java语言实现了无哨兵的单向链表.稍作修改即可实现有哨兵的单向链表.有哨兵的单向链表相较与无哨兵的而言&#xff0c;其对链表的头结点的增删操作更为方便.而在此我们实现了带有头节点和尾节点的双向链表(该头节点和尾节点都不存储有效的数据). 2. 带有头…

手机照片传到电脑的最快方法有哪些?手把手教程来啦!

当我们需要将手机中的照片导入到电脑时&#xff0c;往往会因为速度慢、操作复杂而感到困扰。那么&#xff0c;手机照片传到电脑的最快方法有哪些呢&#xff1f;本文将为您详细介绍2个超快的方法&#xff0c;并附上详细教程&#xff0c;针对每种方法提供详细的操作步骤&#xff…

学习操作系统路线

操作系统 简介 本课程为计算机专业学生量身定制&#xff0c;补足计算机操作系统相关知识&#xff0c;查漏补缺&#xff0c;也可用于考研复习。内容包括&#xff1a;操作统概述、进程管理、内存管理、文件管理、输入/输出管理等章节。内容力求精炼、重点突出、条理清晰、深入浅…

GNU Radio之Frequency Mod底层C++实现

文章目录 前言一、频率调制原理二、Frequency Mod 模块三、底层 C 代码实现 前言 频率调制&#xff08;Frequency Modulation, FM&#xff09;是一种重要的调制技术&#xff0c;广泛应用于无线广播和通信&#xff0c;本文对 GNU Radio 中的 Frequency Mod 模块进行深入剖析。 …

亚远景科技-什么是R.A.S.I.C角色职权矩阵

什么是 R.A.S.I.C角色职权矩阵 在流程定义过程中&#xff0c;亚远景科技推荐使用RASIC 矩阵。RASIC 矩阵是一个非常有用的管理方法。可以明确流程定义中的角色和其相关责任。 "RASIC" 是"Responsible" 、"Accountable" 、"Supportive"…

【论文阅读】EgoPCA: A New Framework for Egocentric Hand-Object Interaction

论文主要贡献 提出一种新的框架&#xff1a;Ego-HOI recognition by Probing, Curation and Adaption (EgoPCA)。构建了全面的预训练集&#xff0c;平衡的测试集&#xff0c;以及一个包含了微调策略的baseline。 在Ego-HOI达到了SOTA&#xff0c;并且建立了有效的机制方法。 …