找到约 15 条结果

文章 boxplot绘图中的离群值Outliers 如何去除

...boxplot完成: 得到离群值的索引: OutVals = boxplot(x)$out which(x %in% OutVals) 如果想删除离群值: x[! x %in% OutVals] 不想绘图,只想得到结果可以这样: OutVals = boxplot(x, plot=FALSE)$out 生物信息入门到精通必修基础课:linux系统...

问题 tophat2建立索引报错:Error: gtf_to_fasta returned an error.

...--------------------------------------------[2019-04-28 05:49:42] Checking for Bowtie                  Bowtie version:        2.2.7.0[2019-04-28 05:49:44] Checking for Bowtie index files (genome)..        Found both Bowtie1 and Bowtie2 indexes.[2019-04-28 05:49:44] Checking for refe...

问题 做有参转录组分析,在进行基因表达定量的时候,出现如下警告,请问老师这是什么原因呀

...3900000 SAM alignment record pairs processed.Warning: Mate records missing for 270457 records; first such record: <SAM_Alignment object: Paired-end read 'A00224:87:H5WHYDSXX:1:2474:9118:3098' aligned to utg9_pilon_pilon_pilon:[3553256,3553424)/+>.24000000 SAM alignment record pairs processed.2...

文章 Linux批量修改文件名

...量文件的时候该如何进行批量修改呢? 方法一:mv配合for循环方式进行修改 $ for name in `ls *.html`;do echo $name ${name%.html}.jpg;done00.html 00.jpg01.html 01.jpg02.html 02.jpg03.html 03.jpg04.html 04.jpg05.html 05.jpg06.html 06.jpg07.html 07.jpg08.html 08.jpg09.html ...

文章 扩增子序列直接物种注释代码

...es = re.findall(r'(Query=.*?)(?=Query=|$)', text, re.DOTALL) results = [] for match in matches: query_id = re.search(r'Query= (\S+)', match).group(1) species_matches = re.findall(r'>([^ ]+) ([^>]+?)\nLength', match, re.DOTALL) identities = re.findall(r'Identities = (.*?)\,', match...

文章 vt软件安装

...Each undeclared identifier is reported only oncecram/cram_io.c:684: error: for each function it appears in.)cram/cram_io.c:684: warning: implicit declaration of function ‘lzma_easy_buffer_encode’cram/cram_io.c:684: error: ‘LZMA_CHECK_CRC32’ undeclared (first use in this function)cram/cram_io...

文章 ggtern 安装报错

** inst ** byte-compile and prepare package for lazy loading Error in get(x, envir = ns, inherits = FALSE) :    object 'message_wrap' not found Error: unable to load R code in package ‘ggtern’ Execution halted ERROR: lazy loading failed for package ‘ggtern’ * removing ‘/share/bios...

问题 bioconductor一直安装不对

...standard repositories, see 'help("repositories", package = "BiocManager")' for details. Replacement repositories: CRAN: https://cran.rstudio.com/ BioCsoft: https://bioconductor.org/packages/3.16/bioc Bioconductor version 3.16 (BiocManager 1.30.20), R 4.2.2 (2022-10-31 ucrt) Installing packag...

文章 HGMD人类变异数据库ANNOVAR注释使用转换格式

官方使用说明:http://resources.qiagenbioinformatics.com/tech-notes/HGMD_Tech_Note.pdf 但是我们需要转换一下,使ANNOVAR能够支持vcf批量注释,转换perl脚本如下: perl ./prepare_hgmd.pl HGMD_PRO_2016.4_hg19.vcf hg19_HGMD_PRO_2016.4.txt my ($dbfile) = $ARGV[0];my ($o...

文章 使用peer根据表达量数据推断隐藏因子

...r of PEER factors was selected as function of sample size (N): 15 factors for N < 15030 factors for 150 ≤ N < 25045 factors for 250 ≤ N < 35060 factors for N ≥ 350 根据文章的提示,使用的是peer这个工具: peer有两个版本,python和R,我们这里使用R包封...

问题 Linux安装anaconda后,加载镜像包,提示失败

...d not found: default_channels: manager@bl8vbox[manager]   - https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/main zsh: no such file or directory: https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/main manager@bl8vbox[manager]   - https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/free zsh...

文章 TCGA数据下载—TCGAbiolinks包参数详解

...the TCGA Data Portal hosted by the TCGA Data Coordinating Center (DCC), in which uses as references GRCh37 (hg19) and GRCh36 (hg18). GDC harmonized database: data available was harmonized against GRCh38 (hg38) using GDC Bioinformatics Pipelines which provides methods to the standardization of biosp...

文章 multi_cox.r 基因表达量做多因素cox分析

...                        set event column name must 0 or 1 code format                        [required]  -v variate [variate ...], --variate variate [variate ...]                        variate for cox analysis [required]  -P predict.time [predict.time ...], --p...

文章 GOChord_plot.r GOplot包绘制GO弦图

... GO enrichment analysis data,the data with columns for 'category','ID','term',adjusted p-value('adj_pval') and 'genes'[required] -d DEG_DATA, --deg_data DEG_DATA input the result of differential expression analysis ...

问题 如何用荧光定量数据绘制热图?

... by qRT-PCR. There were three biological and three technicalreplicates performed for each sample. The relative expression levelswere evaluated by 2-ΔΔCT. GraphPad software was used for the statistical analyses. 作者是用的2-ΔΔCT的方法做的数据分析,可是结果怎么会是负数...