找到约 15 条结果

文章 蛋白质组和转录组相关性图绘制

...id <- function(x){ return(substr(x, start = 0, stop = 16)) } dep$protein_id <- as.character(llply(dep$protein, extr_id)) gene_ids <- deg$`#ID` protein_id <- dep$protein_id # 共有ID common_id <- intersect(gene_ids,protein_id) # 提取数据列表 dep_common = dep[dep$protein_id %...

文章 运行命令makeblastdb建库时报错,生成pdb-lock 空文件。

运行命令makeblastdb -in all.pep.fa -dbtype prot -title all.pep.fa 时卡住报错,生成all.pep.fa.pdb-lock 空文件。 解决:makeblastdb 工具版本问题,从2.13.0版本改成2.6.0版本,就可以正常运行

文章 蛋白质亚细胞定位method

Prediction of secondary structure and subcellular localization of apple BES1 proteins The secondary structure of apple BES1 Proteins was analyzed with NPS@: SOPMA secondary structure prediction (https://www.predictprotein.org/signin). The subcellular localization prediction of apple BES1 proteins wa...

文章 进化树构建方法及原理

...Comparison of the accuracies of several phylogenetic methods using protein and DNA sequences. Mol Biol Evol. 2005 Mar;22(3):792-802. doi: 10.1093/molbev/msi066. Epub 2004 Dec 8. Erratum in: Mol Biol Evol. 2005 Apr;22(4):1160. PMID: 15590907.

文章 perl读写gz压缩文件

...要获取文件句柄,只是稍有变化。 输入文件句柄: open IN,"gzip -dc infile.gz|"    输出文件句柄: open OUT, "| gzip > outfile.gz" or die $!; 也可以用包:PerlIO::gzip use PerlIO::gzip;   #Perl  IO open FQ,"<:gzip",file1.gz||die"$!"; open FN,">...

文章 Graphpad Prism 绘制的直方图美化

...设置界面。 1.Data set 选择直方图中的柱子。 2.Bars and boxes 可以对柱子的填充色,网格,网格颜色,边框粗细,边框颜色进行设置。 3.Error bars 可以对误差线颜色,形状,粗细,显示等进行设置。 4.修改图注,图中...

文章 TCGA数据库中临床样品编号详解(Barcode)

...等 TSS Tissue source site 样品来自哪个组织机构:01 代表International Genomics Consortium, 更多见:TSS Participant Study participant 样品唯一编号(可以理解为一个病人唯一编号) Sample Sample type 样品来自人体组织类型,如:01代表Primary...

文章 blast 建库报错:blast options error: *.fa does not exist

使用blast通过如下指令建库的时候遇到报错: makeblastdb -in *.fa -dbtype prot -title *.fa 解决方法,在库文件前面加上文件的绝对路径 建库成功。

文章 R 相关性系数计算 cor() 报错

...,计算完FPKM后进行相关性计算并绘图,出现报错: Error in seq.int(0, 1, length.out = n) : 'length.out' must be a non-negative number 该报错发生在选择颜色的函数colorRampPalette()上,发现是给的值出现问题,而值来自于相关性计算后的矩阵,...

问题 juicer 报错: Failure during chimera handling of

... 老师好,在运行juicer时报错: awk: calling undefined function and input record number 1116, file /home/data/t090519/JF/20.Jui/juicer/work/splits/JF2.fastq.gz.sam source line number 298***! Failure during chimera handling of /home/data/t090519/JF/20.Jui/juicer/work/splits/JF2.fastq.gz...

文章 samtools之index

必须对bam文件进行默认情况下的排序后,才能进行index。否则会报错。 Usage: samtools index <in.bam> [out.index] Options:    -b    Generate BAI-format index for BAM files[default]    samtools index -b xx.bam  建立索引后将产生后缀为.bai的文件...

文章 evolview 3.0最新更新美化进化树

... Hu, Wei-Hua Chen, Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees, Nucleic Acids Research,(2019) (https://academic.oup.com/nar/advance-article/doi/10.1093/nar/gkz357/5494715) 1.新建上传自己的tree文件: Upload tree control iconProject name tex...

文章 R批量读取文件

...; 把文件的名字(前缀)赋给每个读取到的数据 > for (i in 1:length(ABC)){+ assign(make.names(gsub(".txt","",tem[i])),data.frame(ABC[[i]]))+ }> A V1 V21 A 1> B V1 V21 B 2> C V1 V21 C 3

问题 老师,请问sample type出现这种情况怎么解决,麻烦老师解惑一下,谢谢啦。

...#####################################################> # A. 普通R包: install.packages('caret')> # B. bioconductor 包参考: https://bioconductor.org/install/> # > # BiocManager::install("TCGAbiolinks")> # BiocManager::install("SummarizedExperiment")> ##########################...

文章 安装R包ggpubr 报错

/usr/lib64/R/library/cpp11/include/cpp11/protect.hpp: In member function ‘constexpr R cpp11::protect::function<R(A ...)>::operator()(A ...) const [with R = SEXPREC*; A = {SEXPREC*, SEXPREC*}]’:/usr/lib64/R/library/cpp11/include/cpp11/protect.hpp:186:5: warning: control reaches end of non-v...