ログインしてさらにmixiを楽しもう

コメントを投稿して情報交換!
更新通知を受け取って、最新情報をゲット!

Rコミュのヒストグラム

  • mixiチェック
  • このエントリーをはてなブックマークに追加
こんにちは。
ちょっとお聞きします。
2つのグループのヒストグラムを1つのOutputに描きたいのですがRで可能ですか?
例えば、Rにて2組(XとY)のRandom Numberを作り出します。

X <- rnorm(100,0,1)
Y <- rnorm(100,10,1)

このXとYをもとに2つのヒストグラムを1つのOutputに描けたらいいなって思うんですけど、どなたかそれが可能かどうか知ってる方いらっしゃいますか?
またそれが可能ならどのようにRでするかも教えていただけたらうれしいです。
よろしくお願いします。

コメント(5)

はじめまして

ざっと見たところそういう引数は見当たらなかったので,ちょっと強引ですが,barplotで無理矢理描いてしまおう,と

tmp1<-rnorm(100,0,1)
tmp2<-rnorm(100,5,1)
barplot(rbind(hist(tmp1,plot=F,breaks=-3,9)$counts,rbind(hist(tmp2,plot=F,breaks=-3,9)$counts),beside=TRUE,space=0,names=as.character(-3:8))

こんな感じでそれっぽいものは出来ますが,如何なものでしょうか?(x軸のラベルはとりあえず貼ってみた感じなので,ちょっと修正が必要)
>丸太さん

お返事ありがとうございます。
試してみたんですが、エラーメッセージがでました。
invalid number of breaksだそうです。
何か原因はわかりますか?
ちなみに、"breaks"というのがどのような役割を果たすのか教えていただけたら幸いです。
>Hiroさん

ああ,申し訳ないちと間違っておりました…
別のPCのコマンドをタイプするときになんか間違っていた風です,申し訳ない

barplot(rbind(hist(tmp1,plot=F,breaks=-3:9)$counts,hist(tmp2,plot=F,breaks=-3:9)$counts),beside=TRUE)

Breaksはヒストグラムの区間幅を決める引数です
データをみて適宜指定するという感じで(histならある程度自動でやってくれますが・・・,今回は手動ということで)
>丸太さん

できました。
とてもうまくできてますね。
今までこのようなヒストグラムを作りたくていろいろ試してみたんですが、全然うまくいなかったんです。
俺が作ったのはいつも2つのヒストグラムがくっついてしまって、2回ピークがある1つのヒストグラムになってしまってたんです。
今回のはとてもいいアイデアですね。
ありがとうございました。
>Hiroさん
上手くいって何よりです
どうも,hist自体区間幅を決めてカウントしてそれをbarplotに渡してグラフを書かせる,という挙動をしているっぽいので,カウントさえあれば三つでも四つでもでも書けそうな感じでした〜

高水準な関数だと,複数のプロットを出すっていうのが無いんですかね・・・?三次元のグラフを同じとこに重ねがきしたりとか出来るんでしょうか・・・?

ログインすると、みんなのコメントがもっと見れるよ

mixiユーザー
ログインしてコメントしよう!

R 更新情報

Rのメンバーはこんなコミュニティにも参加しています

星印の数は、共通して参加しているメンバーが多いほど増えます。

人気コミュニティランキング