找到约 15 条结果

文章 immune_infiltrates_rnaseq.r RNA-seq 转录组数据免疫侵润分析

...疫侵润分析. In general values should be TPM-normalized ,not log-transformed.optional arguments:  -h, --help            show this help message and exit  -i expset, --expset expset                        input gene expression set matrix from RNA-seq data tsv            ...

文章 WGCNA 中Relationship between module eigengenes图片绘制--plotMEpairs函数

... interpreted as a biologically motivated data reduction scheme that allows for dependency betweenthe resulting components. In contrast, principal component analysis imposes orthogonality between the components.Since modules may represent biological pathways there is no biological reason why modules ...

文章 pennCNV软件linux系统安装

PennCNV安装,官方安装说明:http://penncnv.openbioinformatics.org/en/latest/user-guide/install/ wget https://github.com/WGLab/PennCNV/archive/v1.0.4.tar.gztar xvfz v1.0.4.tar.gzcd kext/make 但是会报错: gcc `perl -MExtUtils::Embed -e ccopts` -fPIC -c -o khmm_wrap.o khmm_wra...

问题 安装docker toolbox 报错:Looks like something went wrong in step ´Looking for vboxmanage.exe´… Press any key to continue…

win10家庭版,CPU虚拟化已经开启。

文章 安装picrust2 保存:pkg_resources.DistributionNotFound: The 'cython' distribution was not found and is required by biom-format

...hub.com/picrust/picrust2.git cd Next, install the required dependencies listed here. Install PICRUSt2 into it's own environment using conda: conda env create -f picrust2-env.yaml source activate picrust2 pip install --no-deps --editable . Finally, run the tests to verify the install was succ...

文章 泛基因组可视化

...aths.txt#转换 GFA odgi view -i subregion518.og -g > subregion518.gfa #for  Bandage view#绘图odgi viz -i subregion518.og -o subregion518.png -x 5000 #-z strand Black is forward, red is reverse.#获取FASTA 文件:odgi paths -i subregion518.og -f > subregion518.paths.fasta 结果图如下...

问题 gatk报错:java.lang.IllegalStateException: Read position too high for BAI bin indexing.

...91) Caused by: htsjdk.samtools.SAMException: Exception creating BAM index for record A00821:774:H5KM3DSX2:3:1628:6551:31281 1/2 150b aligned to Chr1A:536904672-536904820.         at htsjdk.samtools.BAMIndexer.processAlignment(BAMIndexer.java:119)         at htsjdk.samtools.BAMFileWriter.w...

文章 tcga_gene_exp_download.r TCGA数据下载脚本使用说明。

...oject project                        input project ID of TCGA, for example TCGA-STAD,more                        project ID:https://www.omicsclass.com/article/1061                        [required]  -f files.per.chunk, --files.per.chunk files.per.chunk       ...

文章 maf_oncoplot.r突变注释文件分类可视化maftools

...se] -o outdir, --outdir outdir output file directory [default cwd] -H height, --height height the height of pic inches [default 8] -W width, --width width the width of pic inches [default 8] 参数说明: -i   突变...

文章 linux用户管理命令介绍

...y.  If you _do_ define both, MAIL_DIR takes precedence.#   QMAIL_DIR is for Qmail##QMAIL_DIR      MaildirMAIL_DIR        /var/spool/mail  注:创建用户时,要在目录/var/spool/mail中创建一个用户mail文件;#MAIL_FILE      .mail# Password aging controls:##   ...

文章 miRNA与靶基因相关性作图

...rary(ggplot2)library('getopt');spec = matrix(c('help' , 'h', 0, "logical","for help",'input' , 'i', 1, "character","input file ,required",'name' , 'n', 1, "character","photo name"), byrow=TRUE, ncol=5);opt = getopt(spec);print_usage <- function(spec=NULL){cat(getopt(spec, usage=TRUE));q(status=1)...

问题 老师,根据您的课程分析拟南芥WRKY基因家族,进化树构建时出现No common sites found for computing distances,该怎样解决?已把大的gap删除

文章 使用R语言安装包编译时gcc标准报错——“error: for loop initial declarations are only allowed in C99 mode”

...装Biostrings包的时候报错了: BiocManager::install("Biostrings",force=TRUE) 报错信息提示是gcc标准的问题,应该使用C99标准,可以使用withr包指定使用标准: with_makevars(c(PKG_CFLAGS = "-std=c99"), + BiocManager::install("Biostrings",force=TRUE), + assignment =...

文章 提取R中字符向量中的所有最大长度值,提取向量最大值

...=max(nchar(a))][1]"110""101""abc""cab" 得到最长字符,可以使用which.max() > a <- c("110", "101", "abc", "cab")> a <- c("110", "101", "abc", "cabs")> which.max(nchar(a))[1] 4 更多生物信息课程:1. 文章越来越难发?是你没发现新思路,基因家族分...

文章 BSA分析中QTL-seq 分析:95% 99%阈值置信区间是如何计算的

...RIL"){x<-runif(1) if (x<=0.5){count<-1}else{count<-0} }else{for(i in 1:2){x<-runif(1) if (x<=0.5){number<-0.5}else{number<-0} if(number == 0.5){count<- count+0.5}}}return(count)}#######################################################################caluclate of genotyp...