找到约 15 条结果

文章 linux之sed用法详解

...、选取等特定工作,下面先了解一下sed的用法 sed命令行格式为:           sed [-nefri] ‘command’ 输入文本 常用选项: -n∶使用silent模式。在一般 sed 的用法中,所有来自 STDIN的信息一般都会被列出到屏幕上。但如果加上 -n ...

文章 IF=10.5 | 转录组代谢组学分析揭示了玉米对盐碱胁迫的响应机制

...H146中MDA含量和相对电解质渗漏率显著高于KN3894,表明H146细胞膜损伤更严重。离子平衡:KN3894中Na⁺积累较少,K⁺和Ca²⁺含量较高,Na⁺/K⁺比值显著低于H146。氧化胁迫:NBT和DAB染色显示H146中O₂⁻和H₂O₂积累更多;KN3894中SOD、...

问题 在计算KaKS时,发现分别有两个基因出现BUG,not squal,然后换成phy,不能转换成axt和转换后在计算时依然出现not equal的现象求解,结果文件肯定也是空的

两对基因,第一对出现phy格式不能识别转换成axt,第二对出现phy转换后依然出现计算时not equal。

文章 小鼠单细胞转录组方法参考 method

Mouse breeding Tissue processing and library construction for scRNA-seq(咨询测序公司)   scRNA-seq data pre-processing Preprocessing of raw data was initially performed using Cell Ranger software pipeline (version 9.0.1)(Zheng et al. 2017). This step included reads alignment ...

文章 R语言使用ggraph绘图导出pdf报错——“invalid font type”

...t) dev.off()#记得关闭画布 cairo主要用于SVG、PDF 和 PostScript 格式图形文件的存储 使用方法: svg(filename = if(onefile) "Rplots.svg" else "Rplot%03d.svg", width = 7, height = 7, pointsize = 12, onefile = FALSE, family = "sans", bg = "white", antialias = c("def...

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

...令仅返回“cab”.但是,我想要所有最大长度的字符串.我该怎么做?a <- c("110","101","abc","cab")要测量字符串的“长度”,您必须使用像nchar这样的东西.如果您希望所有具有最大字符数的元素都可以使用nchar(a)== max(nchar(a))进行过滤.以...

文章 转录组测序分析水稻旱作的重要基因

...合转录因子 (GO: 0001071) 以及氧化还原酶活性 (GO: 0016491),细胞外区域(GO:0005576)和胞质膜结合囊泡(GO:0016023)出现了明显富集。 KEGG富集分析显示,旱作条件下,根部组织的DETs明显富集于代谢以及二级代谢合成等相关途径,嫩...

文章 科研汪专属花式表情包,不懂不是科研人!

...思路就是:换个成分呗,换个比例呗,换个浓度呗,换种细胞呗,换个表征方法呗...... 然后你以为科研是每天高大上研究什么不得了的东西改变全世界, 实际上,就是每天烧炉子、烘箱、换液、跑胶...... 三年过去了,好像啥...

问题 eggNOG 保留蛋白序列时提示同一 ID 被多个一级特征使用,而子特征的 Parent 属性没有得到相应更新,导致注释关系不明确的问题,该怎么解决

老师我在进行如下步骤时,会有warning提示,gff和fasta都是非模式物种ncbi数据库染色体水平的组装,ga_protein_coding.gff3文件正常生成了,请问warning对结果有影响吗,应该如何解决?谢谢! agat_sp_filter_feature_by_attribute_value.pl --gff $gf...

文章 IF=4.2 | 杏黄李关键耐寒代谢途径的转录组代谢组联合分析

...宝贵的候选基因资源,并为未来通过标记辅助选择或基因工程提高杏黄李的抗寒能力奠定了基础。

问题 基因组组装文章发表问题

...隆群,DNA量不够,只测了30G的hifi,没办法补数据了应该怎么回复。(3)survey结果是0.97G,组装结果是1.15G。杂合1.19。请问怎么回答好一点

问题 GEO数据库下载的表达矩阵转换成genesymbol时报错。perl脚本。请问perl脚本的问题出在哪?需要怎么修改。最后有perl及需要处理的数据的百度云链接(链接:https://pan.baidu.com/s/1fKim1pGNsU4DDJXcxchplg 提取码:5ls1)

use strict; use warnings; print STDERR "gene symbol column number: "; my $geneSymbolCol=<STDIN>; chomp($geneSymbolCol); $geneSymbolCol--; my $expFile="probeMatrix.txt"; my $gplFile="ann.txt"; my $expFileWF="geneMatrix.txt"; my %hash=(); my @sampleName=(); open(EXP,"$expFile") or d...

问题 GEO数据库下载的表达矩阵转换成genesymbol时报错。perl脚本。请问perl脚本的问题出在哪?需要怎么修改。最后有perl及需要处理的数据的百度云链接(链接:https://pan.baidu.com/s/1fKim1pGNsU4DDJXcxchplg 提取码:5ls1)

use strict;use warnings;print STDERR "gene symbol column number: ";my $geneSymbolCol=<STDIN>;chomp($geneSymbolCol);$geneSymbolCol--;my $expFile="probeMatrix.txt";my $gplFile="ann.txt";my $expFileWF="geneMatrix.txt";my %hash=();my @sampleName=();open(EXP,"$expFile") or die $!; open(PL,"GEOimmun...

文章 perl脚本批量生成反向互补序列

... input.fa是输入文件,out.fa是反向互补后的序列 输入文件格式入下所示: >bta-26a-2GGCUGUGGCUGGAUUCAAGUAAUCCAGGAUAGGCUGUUUCCAUCUGUGAGGCCUAUUCUUGAUUACUUGUUUCUGGAGGCAGCU>bta-18bCUUGUGUUAAGGUGCAUCUAGUGCAGUUAGUGAAGCAGCUCAGAAUCUACUGCCCUAAAUGCUCCUUCUGGCACA>bta-29aAUGACUGA...

问题 黄老师,能不能在虚拟机上安装一个linux版本的paircoil2,我们用虚拟机好像自己安装不了?然后因为电脑版本较低也没法用cocker

我把文件解压缩了,但是没法运行,如下图,所以不知怎么运行 安装说明中提示要把配置文件.paircoil2放到home dictory 里面,这个文件在哪里,怎么放不知道