转录组分析时,计算完FPKM后进行相关性计算并绘图,出现报错:
Error in seq.int(0, 1, length.out = n) :
'length.out' must be a non-negative number
该报错发生在选择颜色的函数colorRampPalette()上,发现是给的值出现问题,而值来自于相关性计算后的矩阵,发现矩阵结果为:
原因是用cor()计算相关性时,矩阵中存在NA,cor()计算前增加na.omit() 即可,修改为:
d<-round(cor(na.omit(fpkm)),4)
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!