找到约 15 条结果

文章 使用Python jcvi 时遇到的报错 ,如:FileNotFoundError,如何解决

...rc 或者.bash_profile , 可以先重新启动终端。           2.如果是之前安装的程序,可以对其进行升级,升级之后,旧版的没有删掉的话,新版的安装路             径以及所依赖的库,(尤其像Python,需要依赖很多库),要...

文章 寿命表法计算中位生存期

...vial ~ 1) # # n events median 0.95LCL 0.95UCL ## 307 60 2888 2052 NA  从结果可知,中位生存时间为2888天。生存曲线如下: # 绘制生存曲线 plot(kmfit,conf.int = F,mark.time = F) #绘制一条中位生存期对应红线 abline(h=0.5,lty=2,col="r...

文章 删除linux文件行尾^M符号

...式下输入     :set fileformat=unix     保存退出     :wq 2. 用dos2unix工具:前提Linux安装了dos2unix     dos2unix [filename] 3. 使用sed 工具:需要注意打^M时应输入Ctrl+v+m,长按Ctrl依次点击v和m键,而不是数字6上的^符号。    sed ’s/...

文章 木薯NBS基因家族分析

...读前准备知识点 1. NBS-LRR是最大的抗性基因家族之一 ; 2. 根据N端的Domains不同NBS-LRR分为两个亚族,含TIR的TIR-NBS-LRR(TNL)和含卷曲螺旋(CC)N端结构域的 CC-NBS-LRR(CNL)。 鉴定木薯NBS-LRR基因家族成员 利用从Pfam数据库下载的...

文章 当代博士的真实困境:秃头、延毕、论文难发、母胎单身......

...要求会适当降低,一般也要在竞争激烈的C刊上发表至少2篇小论文,外加10万字毕业论文,通过答辩后方能毕业。 论文被拒十之八九,导师push你凌晨3点改论文也不是没有, 文章这种硬通货,没有,路当然艰难! 困境之三:...

文章 ggplot从数据到图形绘制逻辑

...作图并联合进行展示。分面也称为条件作图或网格作图。2、绘图步骤

问题 利用perl脚本提取基因上游的1500bp(指定位置)序列时,输出文件总是空文件。尝试过很多办法无法解决

输入文件1-基因组文件如下图: 输入文件2-位置信息文件如下图:注,没有空格并且txt文本是excl文件另存制表符分隔 输出文件以及所敲指令-如下图:指令未显示报错,但是输出文件是空文件 所用代码: 所用数据从jgi...

文章 单细胞转录组数据挖掘流程记录-UM葡萄膜黑色素瘤(GSE139829)

... \ --nUMI.min 500 \ --nUMI.max 40000 \ --nGene.min 250 \ --mito.gene.pattern "^MT.*-" \ --percent_mito 50 \ --log10GenesPerUMI 0.5\ -o 01.qc-$i -p GSE139829 --metadata.col.name Sample --metadata.value $i done​ #合并样本 Rscript $scripts/merge_seurat_obj.r -i $(ls 01.qc-*...

问题 叶绿体基因组GetOrganelle组装时出现killed,内存从10G改为30G又改成70G,但是还是不行,程序终止。

1、在get.sh.o文件中,出现:Retrying with more reads。 2、GetOrganelle_out文件夹中只有,没有组装的序列结果。

文章 单细胞转录组数据挖掘流程记录-BRCA乳腺癌(GSE161529)

...i \ --nUMI.min 500 \ --nUMI.max 50000 \ --nGene.min 250 \ --mito.gene.pattern "^MT.*-" \ --percent_mito 25 \ --log10GenesPerUMI 0.7 \ -o 01.qc-$i -p GSE161529 --metadata.col.name Sample type --metadata.value $i ${i%%-*} done Rscript $scripts/merge_seurat_obj.r -i $(ls 01.qc-*/G...

文章 基因家族复制基因查找之手动查找-blast

... 串联重复基因判断条件: 1,两个基因相似度大于75% 2,比对长度大于75%(相对于较长的基因) 3,基因在同一条染色体且距离小于100kb 参考文献如下: https://link.springer.com/article/10.1007/s00425-016-2575-6 blast手动查找方法...

文章 ifconfig不显示ip

...                           TYPE      DEVICE  eno1  1377a295-708b-4cb4-a1cb-04d8d58cc2bc  ethernet  eno1    lo    e9513de2-7f3b-4faa-a7c0-c8544e0bdcfa  loopback  lo      eno2  fe51385f-86b2-470d-a846-27c73ac4eff1  ethernet  --      eno3  f9129c6b-ca95-4d81-b165...

问题 motif绘图时,tbtools图上显示的motif跟meme.txt文件的不一致

...件放进tbtools里,图上显示几乎每个基因都有motif7、1、8、2,但是检查meme.xml文件和meme.html文件里只有部分基因具有这些motif,遇到这种情况怎么办?有没有遇到过相同问题的同学老师,麻烦指点一下~感谢感谢! tbtools生成图里几...

文章 miRTarBase- 一个收录经实验验证的miRNA靶基因数据库

...目前已经更新到7.0 版本。收录的数据也越来越多。 2. 支持多方面的在线查询方式 支持针对miRNA及靶基因,代谢通路,验证方法,相关疾病等多方面的查询 3. 支持数据下载 在线查询,针对少量的miRNA,比较方便,但是如...

问题 基因家族分析时用大豆gff文件报错 killed

...中想要进行“保留蛋白编码基因”出现报错“=> Version of the Bioperl GFF parser selected by AGAT: 3 gff3 reader error level1: No ID attribute found @ for the feature: 1    ena_assembly_gap        biological_region      2584646                             ...