找到约 15 条结果

问题 利用cicro绘图时遇到错误

... cwd: /work/my_gene_family/09.circos       command: /share/work/biosoft/circos/latest/bin/circos -conf config.txt       -outputdir ./ -outputfile HSP90   You have asked to draw [1190] ideograms, but the maximum is currently set at   [200]. To increase this number change max_ideograms...

文章 dQTL.seq方法研究基因定位

...片。Figure=False: no figure output; Figure=True: the output of figures from different methods. The default is True. Figure-resolution:图片分辨率。Figure-resolution=Low: the figure with low resolution; Figure-resolution=High: the figure with high resolution. The default is High. Figure-fi...

问题 利用getGEO函数下载GSE数据时提示“Error in open.connection(x, "rb") : cannot open the connection”下载GDS,GPL和GSM都是正常的。

我就觉得奇怪了,下载GDS.GPL和GSM的时候都是正常的,就是这个GSE下载不了,就这样提示,谢谢!

问题 在鉴定长链非编码RNA分析中,合并转录本之后,如何批量从gffcmp.combine.gtf文件中提取class_code为i, u,和x的转录本?谢谢

文章 淡黄的长裙,蓬松的头发?拿什么拯救你,我的项目

...量。 ![image.png](https://img-blog.csdnimg.cn/20200330142121131.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMTM3NDkz,size_16,color_FFFFFF,t_70) 没听之前我也是不当回事的看过一个视频之后满脑子嗡嗡的一...

文章 PICRUSt2扩增子功能预测与差异比较分析

...: Douglas, G.M., Maffei, V.J., Zaneveld, J.R. et al. for prediction of metagenome functions. Nat Biotechnol 38, 685–688 (2020). https://doi.org/10.1038/s41587-020-0548-6 Chen Yang, Jiahao Mai, Xuan Cao, Aaron Burberry, Fabio Cominelli, Liangliang Zhang, ggpicrust2: an R package for PICRUSt...

文章 Bowtie2使用方法与参数介绍

...2进行段序列比对的步骤一般如下: 1.对参考序列构建index bowtie2-build genome.fasta index 2.序列比对 bowtie2  -p 8 -x index -1reads1.fq -2 reads2.fq -S out.sam 必须参数: -x <bt2-idx>由bowtie2-build所生成的索引文件的前缀。 -1 <m1>双末...

问题 使用plink做pca分析时报错

...析时报错。 Error: Invalid chromosome code 'chr27' on line 30616892 of .vcf file.(This is disallowed for humans.  Check if the problem is with your data, or ifyou forgot to define a different chromosome set with e.g. --chr-set.) 使用的命令如下: plink --vcf  $workdir/2022_clean.vcf....

问题 (ERR): "--dta" does not exist Exiting now ...为什么hisat2有这个选项却报错

...D: hisat2 -p 1 --rg-id=${i} --rg SM:${i} --rg LB:${i} --rg PL:ILLUMINA \ -x $REF_INDEX --dta --rna-strandness RF \ -1 $workdir/data/${i}_1.fastq.gz \ -2 $workdir/data/${i}_2.fastq.gz \ -S ${i}.sam 2>${i}.summary" hisat2 -p 1 --rg-id=${i} --rg SM:${i} --rg LB:${i} --rg PL:ILLUMINA \ -x $RE...

文章 多款软件进行vcf合并--gatk、vcftools、bcftools

... 两个命令执行结果完全一致,结果如下。 ##fileformat=VCFv4.2 ##FILTER=<ID=q10,Description="Quality below 10"> ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth"> ##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality"> ##FORMAT=<ID=GT,Num...

文章 OthroFinder-物种同源基因分析工作原理

...有根树(Rooted species tree)。基于STAG(Species Tree Inference from All Genes)算法从无根基因树上构建无根物种树,再使用STRIDE(Species Tree Root Inference from Gene Duplication Events)算法构建有根物种树;     1.4 (e)通过有根物种树的情况,...

问题 基因家族分析 kaks分析

...parameters: [OK] Checking Epal2nal.pl: [NA] Please make sure Epal2nal.pl is executable (+X) and put into the global directoryy. l total 118M 运行ParaAT批量比对脚本后,显示以上内容。

问题 这是提取SSR,标记的脚本,{ my $redmotif = "([ACGT]{$j})\\1{".($motiflen/$j-1)."};"my $search = "(([acgt]{$motiflen})\\2{$minreps,})";这句代码代表啥意思,特别是\\1{".($motiflen/$j-1)."\\2{$minreps,}这两句代表啥意思

这是提取SSR,标记的脚本,{         my $redmotif = "([ACGT]{$j})\\1{".($motiflen/$j-1)."};"my $search = "(([acgt]{$motiflen})\\2{$minreps,})";这句代码代表啥意思,特别是\\1{".($motiflen/$j-1)."\\2{$minreps,}这两句代表啥意思。

文章 Pacbio三代全长转录组(Iso-Seq ) 纠错软件-HECIL

HECIL(Hybrid Error Correction of Long Reads using Iterative Learning):  这是一款基于二代数据比对,采用迭代学习的方法对pacbio 等长读长序列进行纠错的软件。该软件的第一作者居然是来自IBM。 1. 软件的安装 安装比较简单,从github 上下载...

文章 wgd syn 报错 KeyError: 6.0

...rk/biosoft/python/Python3/bin/wgd", line 33, in <module>     sys.exit(load_entry_point('wgd==1.2', 'console_scripts', 'wgd')())   File "/share/work/biosoft/python/Python-v3.9.16/lib/python3.9/site-packages/click/core.py", line 1130, in __call__     return self.main(*args, **kwargs)  ...