...A 1>
把文件的名字(前缀)赋给每个读取到的数据
> for (i in 1:length(ABC)){+ assign(make.names(gsub(".txt","",tem[i])),data.frame(ABC[[i]]))+ }> A V1 V21 A 1> B V1 V21 B 2> C V1 V21 C 3
运行core_pan_gene_compare.r ,报错Error in contrasts<-(*tmp*, value = contr.funs[1 + isOF[nn]]) :
contrasts can be applied only to factors with 2 or more levels
Calls: aov ... model.matrix.terms -> model.matrix.default -> contrasts<-
Execution halted。我在泛基因组分类的...
[root@30b37d57431e 00:02:31 /work/3.map/dadi]# gtfToGenePred -genePredExt /work/ref/Bamboo.Hic.gtf unknown_refGene.txt
conflicting frame for PH02Gene38416.t1 exon index 0, was 0, trying to assign 2
...with this analysis:Some pairwise distances could not be estimated.For example,between sequences 23 and 11.
这是因为你的序列差异太大,无法构建进化树,应当将序列差异较大的gap给删除掉,再进行进化树分析,如下图所示:
当然手动删除...
...able),或在硬件上增加内存。
Error: package ‘XYZ’ 1.2.3 is required by ‘somePackage’ but will not be installed
R的生态系统包含大量的扩展包,这些包彼此之间可能有依赖关系。有时,一个包需要的依赖包版本与另一个包需要的版本不...
...: will use block size 2938.
pickSoftThreshold: calculating connectivity for given powers...
..working on genes 1 through 2938 of 2938
Error in datk[c(startG:endG), ] <- foreach(t = actualThreads, .combine = rbind) %dopar% :
被替换的项目不是替换值长度的倍数
谢...
各位老师,我基因家族分析的文章被审稿人问道:It is not clear to me why authors search protein from the genome, not proteome. Moreover, the authors should underline that their results report the finding of candidate genes encoding protein and not ...