生物信息-软件安装

软件安装——conda 下载

从终端进wsl:

wsl

WSL2下安装conda
1.wget: unable to resolve host address解决方法

vim /etc/resolv.conf
nameserver 8.8.8.8 #google域名服务器
nameserver 8.8.4.4 #google域名服务器

2.Windows中WSL2(子系统)设置默认root用户登入修改conf文件出现报错:

E212: Can't open file for writing Press ENTER or type command to continue

原因:

1.当前用户的权限不足

2.此文件可能正被其他程序或用户使用。
一般错误原因都是前者,解决方案是在使用vi命令打开文件时,前面加上sudo来临时提供管理员权限,比如使用命令“sudo vi hosts”打开编辑文件。
或者直接设置root登录
执行命令:

ubuntu2204.exe config --default-user root

值得一提的是,退出root用户命令为:

ubuntu2204.exe config --default-user 用户名

之后重启WSL2
3.使用wget命令下载miniconda.sh
http://repo.anaconda.com/miniconda/Miniconda3-py38_23.1.0-1-Linux-x86_64.sh

wget http://repo.anaconda.com/miniconda/Miniconda3-py38_23.1.0-1-Linux-x86_64.sh

之后我们使用bash命令来运行安装包,并将conda添加到环境变量

vi .bashrc
在.bashrc文件中添加 export PATH=/home/wuhouo/miniconda3/bin:$PATH

如果在安装的时候选择了Init操作,则不需要上面的命令

软件安装——配置channels

1.编写.condarc文件

channels:
 - bioconda
 - conda-forge
 - r
 - defaults

2.搜索fastqc

conda search fastqc

3.安装

conda install fastqc

软件安装:安装blast,hisat2,trinity

软件安装——R Rstudio

来自CRAN:BiocManager,devtools,tidyverse,gridExtra,pheatmap,ggpubr,ggplot2,ggsci

options(repos=" ")#换镜像channels
install.packages("tidyverse")

来自bioconductor:DESeq2,edgeR,GEOquery,ClusterProfiler,limma,org.Hs.eg.db,PCAtools,EnhancedVolcano

options(BioC_mirror=" ")
if (!require("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("limma")

来自github
1.使用devtools:

 if (!requireNamespace('devtools', quietly = TRUE))
        install.packages('devtools')

    devtools::install_github('kevinblighe/PCAtools')

2.自行打包安装

install.packages("pkgbuild")
pkgbuild::build("需要打包的文件路径" "打包后存入的文件路径")#获得的文件是.tar.gz文件
install.packages("打包后存入的文件路径",repos=NULL)

升级R包

remove.packages("")
library("rvcheck")
uodate_all()

生信linux基础

目录操作

ls
ll = ls -l
cd

文件查看

less #查看文件(可翻页)
cat #一次性查看所有文件

j - 向下移动
k - 向上移动
q - 退出 less 
- -chop-long-lines或 - -ch  - 开启水平滚动
/ - 搜索
&something - 只显示文件中包含某些内容的行  #为了移除过滤器,只需键入 &<ENTER> 

常用less -s:显示连续几行空行为一行

管道和重定向

做一个六边形战士
本文转载于网络 如有侵权请联系删除

相关文章

  • Apple 老设备(系统)精简记录

    为了打造HomeLab桌面监控和个人消息推送终端,我偶然入手了一台发布于十年前的老设备(iPad2)。为了能够更方便的使用这台设备,需要对这台设备进行一些特殊操作,比如系统降级、系统运行服务精简。本篇内容,将作为监控终端打造的前篇,先聊聊如何对这台老设备进行系统精简。写在前面这台设备的最终效果图类似下面这样。http://mpvideo.qpic.cn/0bc3fyab4aaaveaam63a5zrfalwddyxaahqa.f10003.mp4?dis_k=6b73e44bbd2cfe6f5769ea6866f55a4e&dis_t=1652760133&vid=wxv_2363450576543367169&format_id=10003&support_redirect=0&mmversion=false接近完成的效果图出于方便后续进行系统降级和持久化越狱,来安装或运行自定义程序,我推荐使用搭载了“过时的”32位芯片的苹果设备,比如:iPad2~4、iPadMini1、iPhone3GS~5。这些设备在二手市场中供应充足,因为一般场景的使用

  • 麒麟绝唱--华为Mate X2-买还是不买

    华为仍然活着尽管美国进行了多轮制裁,但在全球消费者和合作伙伴的支持下,华为仍然还在。在2020年,华为在智能手机和可穿戴设备的市场份额将在中国排名第一,而TWS耳机将实现80%的增长。华为拥有超过10亿的智能互联设备,超过5000万的HiLink生态用户以及每月5.8亿的HMS活跃用户。该公司在全国拥有10,000多家授权体验店,拥有2,000多家授权服务中心,客户服务满意度为99%,连续六年在中国市场排名第一。华为X2折叠屏智能手机这款手机将以前的MateXs的外部折叠模式更改为向内折叠。现在,它带有一个新的双螺旋下降铰链。后者允许屏幕以弧形存储,并避免硬折叠。华为MateX2使用8:7.1黄金比例的柔性内屏。它已经覆盖了业界首个磁控管纳米光学膜。这样可以将反射率降低到1.5%以下。两种屏幕均支持90Hz高刷新率和DCI-P3宽色域。其中,8英寸可折叠柔性内屏具有2480×2200的超高分辨率和412ppi的像素密度。此外,内部屏幕还具有180Hz的触摸采样率。手机的外屏尺寸为6.45英寸,分辨率为2700×1160。像素密度达到456ppi,触摸采样率提高到240Hz。外屏的主流屏

  • PHP数组迭代器的使用方法

    近来在开发一个视力筛查电子报告系统的产品,这个产品的作用是自动提取视力筛查过程中得到的屈光检查数据,并结合数据自动生成通俗易懂且专业的电子报告,以方便家长可以通过公众号或H5链接查阅。要实现这个需求,第一步是对验光设备里打印出来的纸质报告做OCR,图片识别接口返回的是二维数组,报告的原图是这样的: OCR接口返回的数据是这样的array(3){ ["words_result"]=> array(36){ [0]=> array(1){ ["words"]=> string(8)"FA-6000A" } [1]=> array(1){ ["words"]=> string(10)"2022-09-16" } [2]=> array(1){ ["words"]=> string(7)"04:00" } [3]=> array(1){ ["words"]=> string(8)&

  • 深度优先遍历--最大的岛屿

    问题描述:给定一个二维矩阵,0表示水,1表示陆地,一个岛屿是指相邻的上下左右的陆地面积,求最大的岛屿a=[[1,1,1,0,0,0], [1,1,1,0,0,0], [1,0,0,0,1,1], [0,1,1,0,1,0], [0,1,1,0,0,0]] area=0 defmaxAreaOfIsland(a): #记录地图的行,列 row=len(a) col=len(a[0]) foriinrange(row): forjinrange(col): ifa[i][j]==1: #存储当前岛屿的面积 cur=1 #深度优先遍历 dfs(i,j,cur,a) returnarea defdfs(i,j,cur,a): #定义全局变量 globalarea #将以已经遍历过的标记 a[i][j]=2 ifi>0anda[i-1][j]==1: cur=dfs(i-1,j,cur+1,a) ifi<len(a)-1anda[i+1][j]==1: cur=dfs(i+1,j,cur+1,a) ifj>0anda[i][j-1]==1: cur=dfs(i,j-1,cur+

  • SFFAI分享 | 崔泽宇:套装搭配推荐在图神经网络上的应用【附PPT与视频资料】

    作者简介崔泽宇,中国科学院自动化研究所博士一年级在读,本科毕业于华北电力大学(保定)。研究方向为多模态数据挖掘与推荐系统。在WWW会议发过一篇论文。背景 搭配问题是电商系统中的一个非常重要的任务,它需要综合利用图像文本,协同信息等等商品及用户关系,对商品进行深入的理解整合,并且建立一个优秀的美感评估机制,才能够获得一个优秀的搭配系统。当下的互联网电商平台中,阿里京东等公司,都设计了专门独立的搭配平台,比如时尚大脑,滴搭系统。这些搭配平台在给予用户以更高质量的商品推荐的同时,也在辅助商家,帮助商铺提供更好的营销手段有很高的商业价值。 问题定义搭配任务主要的定义为一个打分问题。就是输入一套商品输出一个打分,通过这个打分的高低来评价这一套商品搭配的优劣情况;目前很多工作为了简化任务,把这个任务转化成了两两配对的问题,即输入两个商品,评价两者之间的搭配程度,如下图所示: 过去工作总结近年来,服装搭配任务在深度学习的框架下,获得了更好的学习和表达能力。如下图为主流的工作脉络。在深度学习的基础框架下上的从度量学习的角度提出了基本的框架,把商品特征映射到风格空间进行距离度量,再后来加入多专家模型提高

  • python统计前十出现最多的词

    一、描述    这是一道python面试题:    “一个可读文件,有一万行,一行只有一个单词,单词可以重复的,求出这一万行中出现频繁次数最多的前10个单词”二、思路      先读取文件变为列表,再用集合去重得到一个参照的列表,逆排序取前10(最大即最多的的10个元素),再用参照列表中的每个元素从文件中去统计,把参照列表中的元素作为键,统计到的结果为值,放入字典,打印出来。三、代码#!/usr/bin/python #coding:utf-8 all_C = [] with open("words.txt",'r') as f:     for line in f.readlines():         all_C.append(line) #获取无重复元素 all_set=set(sorted(all_C)) #统计为字典 counts={} for key in all_set:      counts[key] = all_C.count(key) #获取前10个元素的个数变为列表 tens = sorted(counts.values

  • 零基础VB教程027期:各种分支语句的课题小结

    Elseif语句语法格式If条件1then执行代码Elseif条件2then执行代码Elseif条件3then执行代码。。。。【else执行代码】Endif为什么要引入elseif语句?所有的if语句都是必须要进行判断的。而elseif是指一旦满足一个条件之后,将不再判断后面的条件。使得elseif语句的执行效率要高于if语句。If语句与elseif语句的重要区别?If语句条件判断机制:有几个if就判断几次Elseif语句判断机制:至少判断1次Selectcase的引入是为了当条件超多的时候使用!Selectcase中case所要表示的东西Selectcase变量case条件表达式/值执行代码case条件表达式/值执行代码caseelse执行代码endselectcase条件表达式/值这个位置该如何写?1、直接写值“aa”“bb”902、使用is关键词,必须要搭配比较符号><>=<=<>3、使用to关键词,注意to的两边的值是包含的4、使用符号,可以将多个case的条件进行合并错误写法:N>10andn<89在vb6中的case写法是错误

  • 【AI-1000问】为什么CNN中的卷积核半径都是奇数?

    点击边框调出视频工具条为什么CNN中的卷积核一般都是奇数*奇数?熟悉CNN应该都知道常见的卷积核都是3*3或者5*5等,也就是奇数*奇数,似乎都没看过偶数的,这是为什么呢?作者/编辑汤兴旺在CNN中,卷积核的大小是3*3或者5*5是最常见的。也就是说我们见到的卷积核几乎都是奇数*奇数的。在LeNet5中两个卷积层的卷积核都是5*5。而在AlexNet中共有5个卷积层,conv1的卷积核大小是11*11;conv2的卷积核大小是5*5;conv3的卷积核大小是3*3;conv4的卷积核大小是3*3;conv5的卷积核大小是3*3。看到了吗?都是奇数!这是为什么呢? 解答1:更容易padding!在卷积时,我们有时候需要卷积前后的尺寸不变。这时候我们就需要用到padding。假设图像的大小,也就是被卷积对象的大小为n*n,卷积核大小为k*k,padding的幅度设为(k-1)/2时,卷积后的输出就为(n-k+2*((k-1)/2))/1+1=n,即卷积输出为n*n,保证了卷积前后尺寸不变。但是如果k是偶数的话,(k-1)/2就不是整数了。 解答2:更容易找到卷积锚点!在CNN中,进行卷积操

  • .Net Core 自定义配置源从配置中心读取配置

    配置,几乎所有的应用程序都离不开它。.NetFramework时代我们使用App.config、Web.config,到了.NetCore的时代我们使用appsettings.json,这些我们再熟悉不过了。然而到了容器化、微服务的时代,这些本地文件配置有的时候就不太合适了。当你把本地部署的服务搬到docker上后,你会发现要修改一个配置文件变的非常麻烦。你不得不通过宿主机进入容器内部来修改文件,也许容器内还不带vi等编辑工具,你连看都不能看,改都不能。更别说当你启动多个容器实例来做分布式应用的时候,一个个去修改容器的配置,这简直要命了。 因为这些原因,所以“配置中心”就诞生了。配置中心是微服务的基础设施,它对配置进行集中的管理并对外暴露接口,当应用程序需要的时候通过接口读取。配置通常为Key/Value模式,然后通过http接口暴露。好了,配置中心不多说了,感觉要偏了,这次是介绍怎么自定义一个配置源从配置中心读取配置。废话不多说直接上代码吧。模拟配置中心我们新建一个asp.netcorewebapi站点来模拟配置中心服务,端口配置到5000,并添加相应的controller来模拟配置

  • 公交车到站预测2----数据后处理

    之前从csv获取了数据,但是我们的目的是用机器学习的方式对其分类。目测使用sklearn的机器学习库,所以要把数据处理成符合要求的格式。 importtime importnumpyasnp fromsklearnimportcluster,datasets复制 首先是我们需要的一些模块,time就是时间处理的模块,这里的作用就是把表示时间的字符串转换成时间戳,这样就可以计算字符串代表的时间的含义了。numpy个人理解就是可以让python实现基本的matlab功能,也就是说是矩阵处理和基本数学处理的模块。sklearn就是机器学习模块了。十分强大。 defclearDic(dic,number):#delettheunfitablerecords dicb=dic.copy() forkindic: iflen(dic[k].keys())!=20: deldicb[k] returndicb#retrunthecleardictionary复制这是将之前处理好的字典格式的记录中的残缺数据删除,只保留下具有完整数据组的记录。defdic2list(dic):#transferthe

  • 【死磕Java并发】-----J.U.C之阻塞队列:SynchronousQueue

    原文出处http://cmsblogs.com/『chenssy』【注】:SynchronousQueue实现算法看的晕乎乎的,写了好久才写完,如果当中有什么错误之处,忘各位指正作为BlockingQueue中的一员,SynchronousQueue与其他BlockingQueue有着不同特性:SynchronousQueue没有容量。与其他BlockingQueue不同,SynchronousQueue是一个不存储元素的BlockingQueue。每一个put操作必须要等待一个take操作,否则不能继续添加元素,反之亦然。因为没有容量,所以对应peek,contains,clear,isEmpty...等方法其实是无效的。例如clear是不执行任何操作的,contains始终返回false,peek始终返回null。SynchronousQueue分为公平和非公平,默认情况下采用非公平性访问策略,当然也可以通过构造函数来设置为公平性访问策略(为true即可)。若使用TransferQueue,则队列中永远会存在一个dummynode(这点后面详细阐述)。SynchronousQueu

  • 从WGS测序得到的VCF文件里面提取位于外显子区域的【直播】我的基因组84

    首先要下载并且得到人类基因组的外显子坐标记录文件这里我用的参考基因组版本仍然是hg19,所以去CCDS数据库里面下载对应版本,并且格式化成BED文件。wgetftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/archive/Hs37.3/CCDS.20110907.txtcatCCDS.20110907.txt|perl-alne'{/[(.*?)]/;nextunless$1;$exons=$1;$exons=~s/\s//g;$exons=~s/-/\t/g;print"$F[0]\t$"foreachsplit/,/,$exons;}'>hg19exon.bed复制制作好的bed格式的人类全部的exon区域坐标文件如下:1801942802433186132186139218655348657151866418866468187115187127518744198745081874654874839187652387668518775158776301877789877867复制从VCF文件里面根据BED文件进行

  • 互联网应用的状态线程库

    http://state-threads.sourceforge.net/docs/st.html StateThreadsforInternetApplications Introduction StateThreadsisanapplicationlibrarywhichprovidesafoundationforwritingfastandhighlyscalableInternetApplicationsonUNIX-likeplatforms.Itcombinesthesimplicityofthemultithreadedprogrammingparadigm,inwhichonethreadsupportseachsimultaneousconnection,withtheperformanceandscalabilityofanevent-drivenstatemachinearchitecture. 状态线程库是一个应用程序库,是编写快速、高可扩展类UNIX平台互联网应用的基础。他利用了多线程编程泛型和具有高性能和可扩展的事件驱动的状态机架构。 1.Definitio

  • 入门到放弃:理清前端技术概念

    什么是JavaScript JavaScript是一种脚本语言,1995年时,由Netscape公司在网景导航者浏览器上实现。浏览器内置JavaScript引擎来解释执行,用来给网页增加动态功能。Netscape希望它看起来像Java,因此取名为JavaScript,实际上毫无共通之处。Java和JavaScript的关系,就好比“老婆”和“老婆饼”的关系。 什么是ECMAScript 由于微软等公司也推出浏览器上运行的脚本语言进行竞争,1996年Netscape公司决定将JavaScript提交给标准化组织ECMA(“EuropeanComputerManufacturersAssociation”,欧洲计算机制造商协会),希望这种语言能够成为国际标准。次年,ECMA基于JavaScript发布ECMA-262标准文件第一版,即ECMAScript1.0。为体现开放性和中立性,并避免侵犯Sun公司的商标权,该标准被称为ECMAScript而非JavaScript。 ECMAScript和JavaScript的关系 ECMAScript和JavaScript的关系是,前者是后者的规格

  • AI_Trade_Model

    importmatplotlib.pyplotaspltfrommatplotlib.animationimportFuncAnimationfromcollectionsimportdeque#坐标系fig,ax=plt.subplots()等价于fig,ax=plt.subplots(11)。"""fig=plt.figure()matpltlib.pyplot.figure(num=None,#设定figure名称。系统默认按数字升序命名的figure_num(透视表输出窗口)e.g.“figure1”。可自行设定figure名称,名称或是INT,或是str类型;figsize=None,#设定figure尺寸。系统默认命令是rcParams["figure.fig.size"]=[6.4,4.8],即figure长宽为6.4*4.8;dpi=None,#设定figure像素密度。系统默命令是rcParams["sigure.dpi"]=100;facecolor=None,#设定figure背景色。系统默认命令是rcParams["figure.facecolor"]='w',

  • Switch.88E6321/88E6320

    物理Switch数据流:   该设备接受以太网帧,或者丢弃它们,或者从Switch的端口中的一个或多个端口发送它们。   ... 物理接口:   每个端口包含某种物理接口,以接收和发送端口到MAC的帧。   一些端口支持许多不同的物理接口选项,而其他端口只支持一个。   如果端口支持许多接口选项,则一次只能使用一个选项。   每个端口支持的物理接口选项覆盖在每个部分的相应数据表的应用示例部分中。   设备特性参照寄存器的控制特性。   Switch设备中的寄存器被组织成称为端口的三组,全局1、全局2和全局3,另外一组用于访问PHY们,称为PHY。   这些组支持32个16位寄存器,每个端口有属于自己的32个端口寄存器。   一个组中的32个特定的寄存器可以用“偏移”来表示。   举例:     端口控制寄存器被引用为端口偏移0x4,它出现在寄存器地址0x4的端口设备地址空间中。      Switch寄存器:   通过MDIO_CPU与MDC_CPU访问寄存器,PHY设备支持IEEE串口管理接口(SMI)或远程管理的以太网帧。   设备支持2种SMI地址使用模型。   1.使

  • 不定参数宏

    引入 在工程代码中,都会有较为规范的日志信息,除了表达函数是否正常执行、变量值为多少等调试信息外,还有功能模块、打印等级(debug打印、error打印)、所在文件、行数以及函数等信息。 有两种方法可实现,如下: 定义自己的打印函数voidprint_log(intfunctionModule,intprintLevel,char*format,...); 复制 可以在函数内加入时间、文件名、行号以及函数名等信息。 定义各个功能模块的日志信息不定参数宏/*FTP为宏定义*/ #defineFTP_DEBUG(level,arg...) DEBUG(FTP,level,##arg) #defineDEBUG(module,level,format,...) print_log(module,level,format,##__VA_ARGS__); 复制 虽然还是通过不定参数函数实现,但是通过给不同功能模块定义不同的打印宏,使得日志信息表达含义更加明确。 是什么 通过把宏参数列表中最后的参数写成省略号(...),使其可以接受数量可变的宏参数。 表现形式 不定参数宏__VA_ARGS_

  • DC平衡双向控制解串器 转接IC GM8914:FPD-LINK III转LVTTL芯片

    1概述GM8914型DC平衡双向控制解串器,其主要功能是实现将2.8Gbps高速串行数据转换成10或12位并行控制信号,并同步输出一路时钟信号;同时低速通道将芯片控制信息调制到高速差分信号上传输给前级驱动器模块的功能。芯片内部集成终端电阻,可通过外部I/0或I2C总线进行配置,支持powerdown模式。芯片core电源VDDn为1.8V,I/O电源VDDIO可支持3.3V和1.8V两种电压。该芯片的主要应用领域是汽车AdvancedDriverAssistanceSystems(ADAS)中ECU(电子控制单元)视频处理器与防碰撞系统前端摄像机、后视镜摄像机和停车系统成像仪模块之间的无缝、独立双向、低时延通讯。2特征a)工作温度范围:-40℃~105℃;b)电源电压VDDn:1.8V;c)电源电压VDDIO:3.3V或1.8V;d)封装形式:QFN48;e)器件等级:工业级。3封装及引脚功能说明本器件采用48引线的方形扁平无引脚封装(QFN48),引脚排序如下所示。 该芯片的各引脚功能描述见表1: ... 4功能描述   功能框图按图3规定。器件主要实

  • swift 使用匿名函数初始化属性

    swift使用匿名函数初始化属性 匿名函数格式: { //代码体 }() 复制 大括号是定义函数体的,小括号()是调用函数体的。 为什么不写一个函数然后调用函数呢?为何用匿名函数呢? 主要还是代码的简洁性. 举个? privatelazyvarguideScrollView:UIScrollView={ letview=UIScrollView.init() view.backgroundColor=UIColor.lightGray view.bounces=false view.isPagingEnabled=true view.showsHorizontalScrollIndicator=false view.delegate=self returnview }() 复制 上述代码简单清晰明了,如果使用函数,要先写一个函数,然后再去调用,麻烦许多,代码也不明朗。 当时,你也可以使用闭包来这样初始化 privatelazyvarguideScrollView={()->UIScrollViewin letview=UIScrollView.init() view.backgr

  • like 内容转义

    如题,当SQL语句中使用Like查询,且期望匹配的结果中含有"\"的,应当把"\"替换为"\\\\"。   比如数据库中text字段有以下三行: idtext1au1224su5b9e1u9a8c2\24\asu5b9e1u9a8c3\24\\asu5b9e1u9a8c 当我们使用下面这三种方式查询时,都会得到id为1,2,3三行,原因后面会讲。 select*fromtest_tablewheretextlike"%\24%";select*fromtest_tablewheretextlike"%\\24%";select*fromtest_tablewheretextlike"%\\\24%"; 只有像下面这样使用四个反斜杠"\\\\"查询时,才会得到期望的包含"\24"的结果(id为2、3的两行)。 select*fromtest_tablewheretextlike"%24\\\\%";进一步,如果期望查询到的结果更准确,比如只得到id为2的第二行,应该像下边这样查询。 select*fromtest_tablewheretextlike"%24\\\\a

  • linux 扩容磁盘空间

    1、fdisk-l  #查看磁盘情况    2、fdisk/dev/sda PS:n是创建新分区,p是设置成主分区,回车、回车、w保存即可   3、fdisk-l #查看磁盘情况  PS:可以看到sda4就是刚才创建的主分区,但是sda4还没挂载无法使用,此时先reboot重启系统   4、mkfs-text3/dev/sda3 #格式化新分区为ext3格式   5、lvs# 查看lv信息     6、pvcreate/dev/sda4# 将物理硬盘分区初始化为物理卷,以便LVM使用,输入指令  PS:这里会询问是否要初始化,输入y回车   7、df-h #查看lvm组名  PS:红框处就是组名   8、vgextendcentos/dev/sda4#增加卷组的容量   9、vgdisplay #查看可以扩展空间大小  PS:红框处就是可扩展的空间,扩展

相关推荐

推荐阅读