本土化Linux系统,科学网—linux本地化进行lefse分析

注：参考来自网络，如侵权则删。

##对应于上述 A-F 6 个模块，本地版的命令行操作示例如下

#A，设置 LEfSe 的数据格式，详情 format_input.py -h

#-c，指定 class 的行(必须指定)；-s，指定 sub_class 的行(可缺省)；

#-u，指定 subject_id 的行(可缺省)；-o，设置归一化值，默认 -1 即不执行标准化

#注：版本问题，有时 format_input.py 命令找不到，可能为 lefse-format_input.py

#python format_input.py lefse_test.txt A_lefse_test.in -c 1 -s 2 -u 3 -o 1000000

format_input.py lefse_test.txt A_lefse_test.in -c 1 -s 2 -u 3 -o 1000000

#B，LEfSe 分析，详情 run_lefse.py -h

#-l 2.0，设定 LDA 得分的对数值的最低阈值为 2

run_lefse.py A_lefse_test.in B_lefse_test.res -l 2.0

run_lefse.py A_lefse_test.in B_lefse_test.res -l 4.0

#备注：对于 B 的输出，我们可以选择从中删一些不必要(不显著)的数据，

#以增强 C、D、E 作图时的美感

#C，绘制 LEfSe 得分值，详情 plot_res.py -h

#注：版本问题，有时 plot_res.py 命令找不到，可能为 lefse-plot_res.py

#lefse-plot_res.py B_lefse_test.res C_lefse_test.lda.pdf --format pdf --dpi 150 --width 16

plot_res.py B_lefse_test.res C_lefse_test.lda.pdf --format pdf --dpi 150 --width 16

#D，绘制进化分支图，详情 plot_cladogram.py -h

#注：版本问题，有时 plot_cladogram.py 命令找不到，可能为 lefse-plot_cladogram.py

plot_cladogram.py B_lefse_test.res D_lefse_test.cladogram.pdf --format pdf --dpi 150

#E，单张图的展示略，直接使用 F 绘制所有的图

#F，绘制差异特征，详情 plot_features.py -h

#注：版本问题，有时 plot_features.py 命令找不到，可能为 lefse-plot_features.py

mkdir -p F_out

plot_features.py A_lefse_test.in B_lefse_test.res F_out/lefse_test --format pdf --dpi 200

转载本文请联系原作者获取授权，同时请注明本文来自林国鹏科学网博客。

链接地址：http://blog.sciencenet.cn/blog-3448646-1278372.html

上一篇：实验工程菌，无害菌，为什么不能流通在大自然中。

下一篇：conda环境命令

本土化Linux系统,科学网—linux本地化进行lefse分析 - 林国鹏的博文相关推荐

w10系统自带linux系统,科学网—Windows10自带Linux系统 - 陈超的博文
从Build 14393版本(Win10 1607周年更新版)开始,Windows10已经原生支持Linux了,在控制面板的"启用或关闭Windows功能"中即可看到"适 ...
r语言算巢式设计方差分析_科学网—巢式设计之R篇 - 林元震的博文
巢式设计(Nested design),也称为嵌套设计或系统分组设计.把研究对象分成若干组,每组内又分若干亚组,每个亚组又有若干观测值的设计,称为巢式设计.根据因素数的不同,巢式设计可分为二因素(二级 ...
连接linux系统的mysql,Linux系统MySQL开启远程连接
1.远程连接上Linux系统,确保Linux系统已经安装上了MySQL数据库.登陆数据库.mysql -uroot -p(密码). 2.创建用户用来远程连接 GRANT ALL PRIVILEGES ...
imx6ull移植Linux系统第二篇——Linux内核的移植
imx6ull移植Linux系统第二篇--Linux内核的移植花了大概两周的时间,把Linux的移植认真学了一遍,期间踩了不少坑,花费了不少时间去解决各种奇奇怪怪的问题,最终完成了uboot.系统镜 ...
Alpine linux 小内存版,Alpine Linux系统-Alpine Linux系统【轻量级Linux】下载v3.11.0 最新版-西西软件下载...
Alpine Linux系统[轻量级Linux],Alpine Linux其实就是一个清凉版的Linux系统,通过牺牲一些不常用的功能和特性来达到更加轻巧的体积以及更快的响应速度.Alpine Lin ...
Linux系统之部署Linux管理面板1Panel
Linux系统之部署Linux管理面板1Panel 一.1Panel介绍 1.1Panel简介 2.1Panel特点二.本地环境规划 1.本此实践目的 2.本地环境规划三.检查本地环境 1.检查操 ...
linux系统清理内存,Linux实验室：五个清理磁盘空间小妙招
1删除安装包:清空系统APT缓存玩儿电脑最怕的就是卡慢,那么电脑卡慢应该怎么解决呢?对于windows系统来说,你可能有各种免费的杀毒软件.全家桶帮你清空系统空间,那么Linux系统怎么办?今天笔者 ...
linux系统的笔记本,Linux：开源的免费操作系统_笔记本电脑_笔记本评测-中关村在线...
■linux Linux是一套免费使用和自由传播的类Unix操作系统,是一个基于POSIX和UNIX的多用户.多任务.支持多线程和多CPU的操作系统.它能运行主要的UNIX工具软件.应用程序和网络协议 ...
p4vasp的linux版,科学网-p4vasp on Linux-单光存的博文
p4vasp從0.2.0發展到目前的0.3.9除了可以讀POSCAR,CONTCAR以視覺化分子結構之外,更可以xml格式的vasprun.xml的輸出檔,以分析DOS,band structure等 ...
linux去重文件第一列,科学网—Linux实用命令 - 刘洪波的博文
1.查看系统访问日志(实时) tail --lines=100 -f /var/log/secure 2.按列合并文件 paste file1 file2 >combind.txt 3.按行合并 ...

本土化Linux系统,科学网—linux本地化进行lefse分析 - 林国鹏的博文

本土化Linux系统,科学网—linux本地化进行lefse分析 - 林国鹏的博文相关推荐

最新文章

热门文章