「統計学関連なんでもあり」の過去ログ--- 040

No.04241　度数分布表について　　【空星】　2007/08/24(Fri) 23:44

マハラノビス距離の正常データの距離の度数をヒストグラムにしたものと，異常データの距離の度数をヒストグラムにしたものを，一つのグラフにしたいのですが，Rによる統計処理の度数分布表のソースを変えればできるのでしょうか？Rについてはまだ初心者で初歩的な質問ですが，アドバイスいただければ幸いです。宜しくお願いします。

No.04243　Re: 度数分布表について　　【青木繁伸】　2007/08/25(Sat) 07:14

barplot を使うのが簡単かも。以下のような例を

set.seed(12345)
x1 <- rnorm(100)*10+50
x2 <- rnorm(100)*12+70
brks <- seq(20, 110, by=10) # 共通の分割点ベクトル
c1 <- hist(x1, breaks=brks, right=FALSE, plot=FALSE)$counts # 度数1
c2 <- hist(x2, breaks=brks, right=FALSE, plot=FALSE)$counts # 度数2
barplot(rbind(c1, c2), beside=TRUE, space=c(0, 0.2), # 棒を並べて描く
        names.arg=brks[-length(c1)], axisnames=TRUE, axis.lty=1) # 横軸の目盛りラベル等

No.04256　Re: 度数分布表について　　【空星】　2007/08/27(Mon) 15:57

うまくいくことができました。ありがとうございました。