我试图用R绘制一个向量的直方图,当我把它画成x轴时,x轴从5到9,但部分直方图在5之前(最低值是4.414002),部分是在9之后。对我的y轴也是一样,从0到5000,但是talles 'stave‘略高于这5000。
为什么它不在向量中的最低值开始呢?请注意,只是条形不够长,图形大小合适(所以mar或oma设置没有问题,我试过更改这些设置,但没有帮助)
下面是我的代码
import rpy2.robjects as R
import R_functions as R_funct
csvData = (R.r['read.csv'](file='/homes/ndeklein/test.csv', head=True, sep='\t'))
hist = R.r.hist
R.r.png('/homes/ndeklein/test_intensity_hist.png', width=300, height=300)
intensityVector = csvData[0]
logIntensityVector = R.r['log10'](intensityVector)
hist(logIntensityVector, main='Intensity per feature histogram', xlab='logged intensity', ylab='frequencies of features', br=20)
R.r['dev.off']()编辑:
我发现了问题所在,在R代码中,它看起来如下所示:
vector = c(5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 5.67750749154535, 6.7399637450174, 6.7399637450174, 6.7399637450174, 6.7399637450174, 6.7399637450174, 6.7399637450174, 6.7399637450174, 6.7399637450174)
hist(vector, breaks=20)但是由于范围在4-9之间(好吧,在本例中较少),有20个中断太大了。6的中断解决了这个问题。
发布于 2012-02-22 17:12:43
对一个较低的数字(在我的示例中为br=6)上插入断点
hist(logIntensityVector, main='Intensity per feature histogram', xlab='logged intensity', ylab='frequencies of features', br=6)解决了问题。
https://stackoverflow.com/questions/9397033
复制相似问题