- #新建clonal.subclonal,存储56个重要基因的克隆、亚克隆、野生型数目
- clonal.subclonal <- matrix(NA, nrow = length(gene), ncol = 3, dimnames= list(gene, c("Clonal", "Subclonal", "WT")))
- for(i in gene)
- {
- clonal.subclonal[i,] <- table(clonal.table[,i])
- }
- Warning message:
- In strsplit(code, "\n", fixed = TRUE) :
- input string 1 is invalid in this locale
- Error in clonal.subclonal[i, ] <- table(clonal.table[, i]) :
- number of items to replace is not a multiple of replacement length
也就是说运行到一半就结束了,我观察了数据发现,可能是SAMRCD1这个基因没有亚克隆的状态,但是这个问题应该怎么解决