找到约 15 条结果

文章 NCBI下载SRR并转换为fastq文件

...ited at GenBank under the accession PRJNA356989.  想下载该篇文献某一样品的fastq原始数据,看了一些下载教程,方法各异,看完思路还是不清晰,以下记录我的操作步骤,仅供参考。其红色圆圈为步骤记录,蓝色为补充或其他选择...

文章 ConsensusClusterPlus 包对基因表达数据进行一致性聚类

...Consensus Clustering 的基本原理假设:从原数据集不同的子类提取出的样本构成一个新的数据集,并且从同一个子类有不同的样本被提取出来,那么在新数据集上聚类分析之后的结果,无论是聚类的数目还是类内样本都应该和原...

文章 群体遗传进化PCA分析

#方法1: ## plink分析PCA plink --vcf  $workdir/00.filter/clean.sorted.vcf.gz --pca 10 --out  plink_pca   \     --allow-extra-chr --set-missing-var-ids @:#    --vcf-half-call missing #绘图 pca_plink_plot.r -i plink_pca.eigenvec -f $GROUP -g group --name plink_pca #方...

文章 java 安装以及环境变量配置

...右击我的电脑->属性->高级->环境变量->户变量单击[新建(N)]添加以下环境变量 1.配置JAVA_HOME 新建JAVA_HOME 变量名        JAVA_HOME 变量值        C:\Program Files\Java\jdk1.8.0_152  你的jdk安装地址 2.配置PATH         ...

问题 老师,我在学习基因家族视频,perl脚本get_gene_weizhi预测基因顺式作元件时出现了问题,老师请看

老师,是这样的,我的物种是川桑(Morus notabilis),没有收录在在Ensembl数据库,所以我是在川桑数据库下载的genome,gene,cds及protein。如下图 然后在NCBI上下载了川桑的GFF注释文件,如下图 老师,下图是我根据视频学习找到...

问题 老师,我在基因家族鉴定时,新生成的hmm文件做hmmsearch时,报TC bit thresholds unavailable on model domain的错误。最开始以为是e值太高1e-5,但我后面改为1e-10还是报这个错误。请问老师是什么原因呢?

问题 老师,我hmmsearch 没搜索出结果,然后拟南芥的query.fa和我要分析的物种cds进行blast分析,比对出来的基因只有几个,我老师说我要做的基因有200多个,完全对不上,不知道该怎么办了,请教下老师。

文章 IF=5|干眼症微生物多样性分析案例解析

...身的基因组多100倍,有1000到1150种细菌定植在肠道内,其约160种是人体肠道的优势菌群。已有研究表明人体微生物与很多人类疾病息息相关。如肠道微生物分泌的神经递质,可调节肠道神经系统的神经信号,导致帕金森病、...

问题 HISAT2比对时报错|sam文件为空白,summary文件为(ERR): hisat2-align died with signal 9 (KILL)

HISAT2比对时报错|sam文件为空白,summary文件为(ERR): hisat2-align died with signal 9 (KILL),后续检查了内存什么的下图请问为何执行不了,是运行空间不够还是什么?

问题 不同胁迫处理以及不同类型的RNAseq数据denovo组装

...RNAseq数据做分析,之前不同时间处理的样品,call表达量的基因ID不一样。现在想统一下。 由于不同批次测序类型也不一样,有双端测序也有单端测序结果。样品是各种胁迫处理,比如说高温低温干旱处理还有对照 这些样...

问题 请问相同的文件,只是把名字改了一下,然后docker自动把重复的序列删除后,为什么出来的motif不一样啊,然后我又把名字改成不一样的,让他不能过滤掉重复序列,之后生成的就和之前的一样了,是什么回事

问题 请问为什么GWAStassel mlm模型画的manhattan图什么都没有呢 也没有报错 能生成一个什么都没有的图片 只有坐标和题

 ../../script/gwas_manhattan_plot.r -i mlm_pvalue.txt -F msu7.0_bowtie2.fa.fai -T Grains_per_panicle -n grains_per_panicle_manhattan -c 1.3e-6 --> Q&A for bioinformatics, please visit the website: www.omicsclass.com --> Recommended to learn R language: --> https://study.163.com/cour...

文章 使DNA序列做主成分分析(PCA)——R语言adegenet包

...据。fasta2genlight()函数的主要作就是从比对好的fasta文件提取snp数据,存入到genlight对象。 三、从DNA序列提取多态性可以通过 genind 对象或 genlight 对象来有效处理。对于较大的(接近全基因组)数据集,可以使 fasta2genlight...

文章 基因家族分析后续试验怎么做?

...下面就通过一篇2018年9月10日发表在Frontiers in plant science(科院分区: 2区;IF=3.678)题为Genome-Wide Identification and Expression Profiling of the TCP Family Genes in Spike and Grain Development of Wheat (Triticum aestivum L.)的文章,一起来学习下基因家族分...

问题 请问我跑一个perl脚本,怎么限定它在服务器的线程呢?因为没有限定线程服务器好响啊,之前有些脚本会提前-T指定线程,可是这个脚本是我自己写的,突然不知道怎么限定线程了,谷歌没找到答案,求救于此