1

Rを使用してベクトルのヒストグラムを作成しようとしています。これをプロットすると、x軸は5から9になりますが、ヒストグラムの一部は5より前（最小数は4.414002）で、一部は9より後です。同じことが私のy軸にも当てはまり、0から5000になりますが、タレスの「ステーブ」はこの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に休憩を入れると、問題が解決します。

4

1 に答える 1

0

より低い数値にブレークを配置する (私の例では br = 6)

hist(logIntensityVector, main='Intensity per feature histogram', xlab='logged intensity', ylab='frequencies of features', br=6)

問題を解決しました。

于 2012-02-22T17:12:43.643 に答える