2012-03-29 10 views
2

は、私は、データファイルからヒストグラムをこのように作成しました:R、ヒストグラムドメインを絞り込む方法は?私はR言語で取り組んできた

< - (「_ datafile.txt」)read.tableを。

b < - 表(a)。

c < - as.numeric(名前(b));

hist(c、100);

これは、添付の画像で、結果である:

Histogram http://www.imagehost.it/dm/FILK/istogrammaBis.png

私の問題は、私は、ヒストグラムの1rstバーを除外したいということです。

あなたはそれを行うためのRコマンドを知っていますか?

ありがとう

答えて

3

histデフォルトの区切り位置として(nclass.Sturgesによって決定される数の等しい間隔をおい)「スタージェス」ブレークを使用しますが、過剰に乗っているのです。第ブレーク値は

配列であろう(MIN(X)、MAX(X)=差分によって(範囲(X))/ 10)[2]を

hist(c[ c >= seq(min(c), max(c), by= diff(range(c))/10)[2], 9 ]) 

私はカテゴリーを同じに保つためにビンの数を減らしました。デフォルトのブレーク数の計算式を見たい場合は、単に入力します。

nclass.Sturges 

それは等間隔の観測数の少ないhistの行動によって困惑の人を見つけることが珍しくありません。結果は不規則に見えるためです。試してみてください:

hist(0:8) 
+0

ありがとう、しかし、あなたの中にxは何を表していますか?それは私のCですか? –

+1

右。修正されます。 ( "c"は変数名の非常に悪い選択ですが) –

+0

そうです、それはちょうどabc ... Thanx :-) –

0

をCRAN上HistogramToolsパッケージには、ヒストグラムからビンを削除するか、ヒストグラムのドメインを絞り込むことが簡単なヘルパー関数を多数提供しています。たとえば、SubsetHistogramおよびTrimHistogramの機能を参照してください。

例えば、我々は、この最初のヒストグラムがあるとします。私たちは、このようなTrimHistogramと空のバケツをオフにトリミングすることができ

hist.1 <- hist(runif(100,min=2,max=4), breaks=seq(0,6,by=.2), plot=FALSE) 

enter image description here

を:

hist.trimmed <- TrimHistogram(hist.1) 

enter image description here

私たちは、毛皮ができるTHERは、このようなSubsetHistogram機能付き(オリジナルポスター用語でドメインを狭める)、ヒストグラムの左または右のビンを除外するために、最大または最小のブレークポイントを指定します。

plot(SubsetHistogram(hist.trimmed, 2.2)) 

enter image description here

関連する問題