分析的文件样本数很少（个位数）。样本数多的文件，分析就不会报错。所以这个应该怎么改呢（我前一段时间做都不会报错，难道是数据需要更新？）
群体遗传课程57课PSMC那里，我看视频是一个.bam文件就能做一次分析，出一张图了。那我的数据有200多个样本400个bam文件怎么办呢？？需要合并吗？怎么合并啊？ 下面这个是例子代码：
为什么我从数据库下载的相关基因的氨基酸序列，格式也是fasta,进行非模式作物注释时，提交上去总是会出现这样的错误
如果我做水稻重测序的话，table_annovar.pl注释的时候，-protocol refGene，这个我需要自己下水稻的基因组数据库吗？还是有官方的？如果需要自己下载的话，是下载参考基因组的吗？
首先是准备材料： 1.R包 2.绘图数据 install.packages('pheatmap')#安装R包 library(pheatmap)#加载R包 画图 画个基础图形 Rscript heatmap.r -i fpkm.txt -s -n heatmap绘制图形，并保存文件。 画的图拍你可以根据自己的喜好进行调整，包括横, ...
宏基因组binning后的mags，也可以用kegg注释出功能是嘛？是拿到mags的碱基序列，直接扔进官网去注释吗，还是也要基因预测，然后蛋白质序列再扔进官网去注释？或者用本地数据库？这方面的课程会有教学吗？
在mfuzz分析完后，通过ClusterGVis可视化，但是出现了图片中的报错，开始以为是自己数据格式有问题，我替换了课程资料提供的示例文件，也出现同样的报错内容，请问怎么解决？
TCGA 数据库中的基因编号采用的Esembl 的编号，但是有些分析软件，需要输入的基因编号是 gene symbol ，这就需要将Esemble 的ID 转换成gene symbol 。 今天介绍采用clusterProfiler 进行转换： # 加载相关软件包 > library(clusterProfiler) > ...
我们课题组测了基因组但是还未曾发表，但同时让公司以基因组的测序结果为参照测了有参转录组，那现在手机的这个转录组数据能不能发sci？？？跪求翻牌～～
...命令行参数input_file <- args[1]output_prefix <- args[2]# 读取数据文件, 不转换列名venn_dat <- read.delim(input_file, check.names = FALSE)# 从文件中创建一个列表venn_list <- list(venn_dat[, 1], venn_dat[, 2], venn_dat[, 3], venn_dat[, 4], venn_dat[, 5], venn_dat[, ...
...变积累系（MA）的基因组变异，并基于野莲资源的重测序数据重建了野莲群体历史，探讨了野莲迁徙的历史路径及亚洲野莲的受选择分析。 本文的研究数据来自于两种材料，一个是课题组种植的突变积累系，同一个L0分芽后无性...
这个是您示例Arabidopsis_thaliana.TAIR10.41.gff3文件最后的几段文件，但是我下载的物种数据库下面没有这些，我想要补齐，想问问老师，最下面的这些是什么意思。谢谢！
从公共数据库中的下载的二代基因组测序文件（双端测序：.fastq_1，.fastq_2文件），下载后多个样本后，如何从这些测序文件中提取特定某个基因序列信息进行后续分析特异性分析。
各位老师，我想问下我用得到的某一物种的特异蛋白的CDS序列去NCBI-blast，鉴定是否是这个物种的特异蛋白序列。选择reference genomic sequences数据库时为什么出现这种错误，比对时需要设置一些参数吗，谢谢老师
请问我的数据是一个实验设计有5个处理，每个处理做了3个重复样测序，想做一下每个处理的单独的otu网络，只有3个重复，800多个otu,在用conet构建网络，在过滤那步跑了几天也没跑完，是不是样品量太少的时候不能用conet啊？