老师，您好。我再用vcftools计算Fst和Pi时，设置的窗口步长是一致的，但是跑出来的结果有一些区域出现了跳过的情况，导致两个文件的窗口不一致，如图所示，pi值的范围直接从15w到了50我，怎么解决这一问题呢？
...有500多个，可是文献里鉴定出来的是14个左右，想问一下筛选的cdd条件怎么确定呢？使用STKc_MAPK作为筛选条件只有3个，STKc_MAP作为筛选有7个，以STKc作为筛选有163个，我查阅了文献，他们的筛选方法是保证有STK，丝氨酸和组氨酸...
运行水稻谷子均可 但是玉米就结果只显示一行，是不是需要修改代码？谢谢老师~
piranha输出的bed文件会把每一个bin给出一个p值，大部分bin的p值都小于0.01，要怎么筛选出显著的peak去找motif？如果一个peak跨两个bin，给出的p值是bin的而不是peak的，这种要怎么处理？
按着步骤下来，这是什么原因呢？是我哪里弄错了吗？ 该怎么解决呢？老师们
limma 筛选差异基因的时候，上调基因80多个，很少，但是下调基因多大900 多个，这样正常吗？请教大家一下，谢谢 # 以logFC>1，Qvalue<0.05为阈值进行过滤up <- subset(DEG1, adj.P.Val < 0.05 & DEG1$logFC >= 1)down <- subset(DEG1, ad...
...因家族发表的文章数量也是逐年递增，小编到NCBI的pubmed搜索了一下基因家族分析相关的文章，历史文章总共有11948篇，非常多！2020年有1000多篇，今年2021刚刚过去一个季度就有300多篇，按这个速度，今年发表1000多篇想必也是妥...
...带来一些涉及到行业关键词的基金统计结果。 小编一共搜索了14个关键词（个别词只是细微差别），同时统计了17年和18年两年的搜索结果。 为了看起来更直观一些，小编特意做了一个柱状图，如下： 从图中可以看出： RNA...
#根据序列ID特点筛选对应ID序列$ zcat hairpin.fa.gz | seqkit grep -r -p ^hsa >hsa-let-7a-1 MI0000060 Homo sapiens let-7a-1 stem-loop UGGGAUGAGGUAGUAGGUUGUAUAGUUUUAGGGUCACACCCACCACUGGGAGAUAACUAU ACAAUCUACUGUCUUUCCUA >hsa-let-7a-2 MI0000061 Homo sapiens let-7a-2 stem-loop AGGUUGAGGUAGU...
做WGCNA分析，数据量比较大，有45783个基因，但样本数只有12个。运行后得到的cluster dendrogram图中的module colors的颜色只能显示出蓝色，但我table(moduleColors)后得到的颜色模块还是很多的。想问下我这个操作是哪里出现问题了呢？万...
...到变量选择中去。 说简单点：在回归分析中因素筛选主要用到逐步回归stepwise、向前、向后等等方法，这些方法比较传统，而对于共线性问题比较严重的数据，或者变量个数大于观测值个数例如基因测序数据，基因个数远...
st 与 pi 联合 筛选受选择区域 fst_pi_select_sweep.r --fst Fst.wild.cultivated.windowed.weir.fst \ --pi1 pi.wild.windowed.pi --pi2 pi.cultivated.windowed.pi --zscore --log2 \ -A wild -B cultivated -c 0.05 -n fst-pi.wild-vs-cultivated -f pdf fst_pi_select_sweep.r --fst Fst...