- genomeID=read.csv("genomeID.csv",row.names = 1)
- gene_to_ko=read.csv("gene_to_ko.csv")
- genomeIDbackup=genomeID
- gene_to_kobackup=gene_to_ko
- for(genomeX in rownames(genomeID)){
- for (geneY in list(colnames(genomeID))) {
- if(genomeID[genomeX,geneY] %in% gene_to_ko$genomeX){
- genomeID[genomeX,geneY]<-1
- }else{
- genomeID[genomeX,geneY]<-0
-
- }}}
- ##
- Error in if (genomeID[genomeX, geneY] %in% gene_to_ko$genomeX) { :
- the condition has length > 1


雷达卡


京公网安备 11010802022788号







