在linux中用同一个版本的R 同时安装 Seurat 2 和 Seurat 3

Seurat  作为单细胞分析中的重量级R包，有多好用用，用过的人都知道。Seurat 分析流程基本涵盖了单细胞分析中的所有常见分析方法，包括filtering，tSNE，UMAP降维及画图等。还有一个重量级功能就是矫正不同实验之间的批次效应。然而Seurat 2和Seurat 3的矫正方法完全不一样，得到的结果也不一致。Seurat 2是基于CCA (典型相关性)的，可以矫正肿瘤，外周血及癌旁组织间由于实验带来的批次效应，也能很好的矫正用不同的单细胞实验平台进行试验带来的批次效应。虽然速度慢，效果还是不错的。而Seurat 3 则是基于样本间具有相似表达谱的细胞群来进行矫正，对于同一种性质的实验，由于不同单细胞技术造成的实验批次效应，seurat 3 能够很好的矫正。从官网给的pancers矫正结果就可以看到其矫正能力多么强大。然而正式因为如此强大的矫正能力，对于肿瘤和外周血样本的矫正却过了头，导致不该分在一起的细胞具有了相似的基因表达谱。本人也是做了好几个课题，发现都存在这样的问题，因此果断放弃Seurat 3的矫正方法，继续用Seurat 2的。但是Seurat 3的 findmarker 这个功能可以一次计算10万以上的细胞不报错，而Seurat 2就不行，折衷的方案是同时安装 Seurat 2和 Seurat 3的包，在内存里切换数据，而不用写到本地后再用Seurat 3读取后升级。尤其是对于动辄10几万个细胞来说，保存数据到本地这个操作要花费至少30min, 读取也要30min.

下面我就告诉大家不用读写到本地就可以在Seurat 2 和 Seurat 3之间完美切换，。

其实方法很简单，将Seurat 2和 Seurat 3 安装在不同的 library 里面就行了。

我已经安装好了，以我自己进行的自由切换为例:

> R.version

_

platform x86_64-conda_cos6-linux-gnu

arch x86_64

os linux-gnu

system x86_64, linux-gnu

status

major 3

minor 6.1

year 2019

month 07

day 05

svn rev 76782

language R

version.string R version 3.6.1 (2019-07-05)

nickname Action of the Toes

我用的是最新的R版本 3.6.1很好用。

默认的library 是conda 自带的

> .libPaths()

[1] "/data/home/heshuai/anaconda3/lib/R/library"

默认的Seurat是最新版的 Seurat 3

> library(Seurat)

Registered S3 method overwritten by 'R.oo':

method from

throw.default R.methodsS3

> packageVersion("Seurat")

[1] ‘3.0.2’

我在另一个library 里安装了 Seurat 2

/data/home/heshuai/R/x86_64-conda_cos6-linux-gnu-library

在两者间自由切换

1. 首先将 Seurat 2 所在的library 加载进来

> .libPaths("/data/home/heshuai/R/x86_64-conda_cos6-linux-gnu-library")

> .libPaths()

[1] "/data/home/heshuai/R/x86_64-conda_cos6-linux-gnu-library" "/data/home/heshuai/anaconda3/lib/R/library"

>

2. detach Seurat 3 后加载 Seurat 2, 因为这个时候Seurat 2 所在的library 已经在Seurat 3 之前了，系统会默认先加载Seurat 2

> detach("package:Seurat", unload = T)

> library(Seurat)

Loading required package: ggplot2

RStudio Community is a great place to get help: https://community.rstudio.com/c/tidyverse.

Loading required package: cowplot

********************************************************

Note: As of version 1.0.0, cowplot does not change the

default ggplot2 theme anymore. To recover the previous

behavior, execute:

theme_set(theme_cowplot())

********************************************************

Loading required package: Matrix

> packageVersion("Seurat")

[1] ‘2.3.4’

>

现在Seurat 3已经成功的切换成Seurat 2了. 想要加载Seurat 3的时候，将默认library 换到Seurat 2的前面即可。

是不是 so easy !

Linux,Unix各种版本的操作系统在线安装软件命令

摘自:http://blog.csdn.net/zjg555543/article/details/8278266 linux和unix,各个版本的操作系统都有自己的软件安装方式,最方便的莫过于在线安 ...

Linux mint 18版本开启SSH服务

linux mint 18版本默认是没有安装ssh server的 需要手动安装 安装ssh server: 此命令需要联网,会自动下载安装 安装之后看是否开始了ssh, 看到ssh-agent 和s ...

linux中用shell获取昨天、明天或多天前的日期

linux中用shell获取昨天.明天或多天前的日期 时间 -- :: BlogJava-专家区 原文 http://www.blogjava.net/xzclog/archive/2015/12/0 ...

linux中用shell获取时间,日期

linux中用shell获取昨天.明天或多天前的日期:在Linux中对man date -d 参数说的比较模糊,以下举例进一步说明:# -d, --date=STRING display time d ...

linux部署不同版本mysql

测试环境部署过程中经常会遇到同一个服务器上部署两个不同版本的mysql数据库,在部署过程中也会有各种各样的问题,现将部署多版本mysql的方法总结如下: 1.下载mysql版本 http://down ...

linux查看内核版本、系统版本、系统位数(32or64)

linux查看内核版本.系统版本.系统位数(32or64) 2011-05-01 22:05:12 标签:linux 内核版本 休闲 系统版本 系统位数 1. 查看内核版本命令: 1) [root@ ...

Linux下指定版本编译安装LAMP

说明: 操作系统:CentOS 6.5 64位 需求: 编译安装LAMP运行环境 各软件版本如下: MySQL:mysql-5.1.73 Apache:httpd-2.2.31 PHP:php-5.2 ...

Linux 查看系统版本及位数

1. 查看内核版本命令: 1) [root@www ~]# cat /proc/version    Linux version 2.6.9-22.ELsmp (bhcompile@crowe.dev ...

Linux 查看系统版本和内核

查看系统内核版本 [root@11e71db4a00e /]# cat /proc/version Linux version -.el7.x86_64 (builder@kbuilder.dev.c ...

随机推荐

不停止MySQL服务增加从库的两种方式

不停止MySQL服务增加从库的两种方式 转载自:http://lizhenliang.blog.51cto.com/7876557/1669829 现在生产环境MySQL数据库是一主一从,由于业务量访 ...

你好，C++(5)如何输出数据到屏幕、从屏幕输入数据与读写文件？

2.2  基本输入/输出流 听过HelloWorld.exe的自我介绍之后,大家已经知道了一个C++程序的任务就是描述数据和处理数据.这两大任务的对象都是数据,可现在的问题是,数据不可能无中生有地产生 ...

BZOJ 1651: [Usaco2006 Feb]Stall Reservations 专用牛棚( 线段树 )

线段树.. -------------------------------------------------------------------------------------- #includ ...

洛谷 P1083 借教室

传送门:Probem 1083 https://www.cnblogs.com/violet-acmer/p/9721160.html 一.暴力简述 首先我们不难看出,这道题--并不是一道多难的题,因 ...

Android长度单位详解(dp、sp、px、in、pt、mm、dip)

Android中定义的dimension单位有以下这些:px(Pixels ,像素):对应屏幕上的实际像素点.in(Inches ,英寸):屏幕物理长度单位.mm(Millimeters ,毫米):屏 ...

django配置超级管理员

Admin站点管理: 概述: 内容发布:负责添加,修改,删除内容的 公告访问 配置Admin应用: 在settings.py文件中的INSTALLED_APPS中添加'django.contrib.a ...

oozie 重新提交作业

在oozie的运行过程当中可能会出现错误,比如数据库连接不上,或者作业执行报错导致流程进入suspend或者killed状态,这个时候我们就要分析了,如果确实是数据或者是网络有问题,我们比如把问题解决 ...

Ajax棵

ajax 1.什么是ajax?(异步请求,局部刷新) ajax是一个改善用户体验的技术,实质上是利用浏览器端ajax对象()向服务器发送异步(ajax对象在向服务器发送请求的时候,用户可以继续其他操作 ...

《opencv学习》 之 OTSU算法实现二值化

主要讲解OTSU算法实现图像二值化:    1.统计灰度级图像中每个像素值的个数. 2.计算第一步个数占整个图像的比例. 3.计算每个阈值[0-255]条件下,背景和前景所包含像素值总个数和总概率(就 ...

2017-2018-1 20155220 《信息安全系统设计基础》课下实践——实现mypwd

学习pwd命令 输入pwd命令 于是man 1 pwd查看pwd详细 然后查看pwd实现需要的系统调用man -k; grep 在这发现了一个功能相同的内核函数getcwd 到这步就很简单了,先查看这 ...