2

信じられないほどシンプルなことをしたい。完全なデータフレーム用に1つの箱ひげ図を作成したい。それでも、「結合された箱ひげ図」および関連する用語を検索しても、提案は見つかりませんでした。明らかな方法を見落とした場合は、お知らせください。

私は次のデータを持っています:

> theData
     X20.7    X21.7    X22.7    X23.7    X24.7    X25.7    X26.7    X27.7    X28.7    X29.7    X30.7    X31.7    X32.7    X33.7    X34.7    X35.7
1 99.64920 99.49319 99.49319 99.49319 99.49319 99.49319 99.80837 99.29348 99.29348 99.29348 99.29348 99.29348 99.29348 99.46376 99.46376 99.51554
2 98.76469 98.60867 98.60867 98.60867 98.60867 98.60867 99.41553 98.40896 98.40896 98.40896 98.40896 98.40896 98.40896 98.74975 98.74975 98.54527
3 98.37824 98.22222 98.22222 98.22222 98.22222 98.22222 98.70900 98.13767 98.13767 98.13767 98.13767 98.13767 98.13767 98.47846 98.47846 98.01791
4 98.11356 97.95754 97.95754 97.95754 97.95754 97.95754 97.82447 97.93003 97.93003 97.93003 97.93003 97.93003 97.93003 98.27083 98.27083 97.81027
5 97.80027 97.64424 97.64424 97.64424 97.64424 97.48632 97.43801 97.40158 97.40158 97.40158 97.40158 97.40158 97.40158 97.74239 97.74239 97.28181
6 97.47825 97.32222 97.32222 97.32222 97.43795 97.12131 97.17333 97.03658 97.10158 97.10158 97.10158 97.10158 97.10158 97.44239 97.44239 96.98180
> dput(theData)
structure(list(X20.7 = c(99.6492, 98.7646913866934, 98.3782376564915, 
98.1135635544627, 97.8002672890352, 97.4782549804011), X21.7 = c(99.4931928571429, 
98.6086741582754, 98.2222160140822, 97.9575388921788, 97.6442390541023, 
97.3222230681959), X22.7 = c(99.4931928571429, 98.6086741582754, 
98.2222160140822, 97.9575388921788, 97.6442390541023, 97.3222230681959
), X23.7 = c(99.4931928571429, 98.6086741582754, 98.2222160140822, 
97.9575388921788, 97.6442390541023, 97.3222230681959), X24.7 = c(99.4931928571429, 
98.6086741582754, 98.2222160140822, 97.9575388921788, 97.6442390541023, 
97.437947563131), X25.7 = c(99.4931928571429, 98.6086741582754, 
98.2222160140822, 97.9575388921788, 97.4863155584865, 97.121313307238
), X26.7 = c(99.8083714285714, 99.415530164398, 98.7090041774867, 
97.8244717838903, 97.4380076185552, 97.173326388931), X27.7 = c(99.2934828571429, 
98.4089615689001, 98.1376722694449, 97.9300324124538, 97.401583100132, 
97.03657716757), X28.7 = c(99.2934828571429, 98.4089615689001, 
98.1376722694449, 97.9300324124538, 97.401583100132, 97.1015782240536
), X29.7 = c(99.2934828571429, 98.4089615689001, 98.1376722694449, 
97.9300324124538, 97.401583100132, 97.1015782240536), X30.7 = c(99.2934828571429, 
98.4089615689001, 98.1376722694449, 97.9300324124538, 97.401583100132, 
97.1015782240536), X31.7 = c(99.2934828571429, 98.4089615689001, 
98.1376722694449, 97.9300324124538, 97.401583100132, 97.1015782240536
), X32.7 = c(99.2934828571429, 98.4089615689001, 98.1376722694449, 
97.9300324124538, 97.401583100132, 97.1015782240536), X33.7 = c(99.4637585714286, 
98.7497473555799, 98.478463763926, 98.2708282766442, 97.7423900760775, 
97.4423915096353), X34.7 = c(99.4637585714286, 98.7497473555799, 
98.478463763926, 98.2708282766442, 97.7423900760775, 97.4423915096353
), X35.7 = c(99.5155421428571, 98.5452656069643, 98.0179127183643, 
97.81026932055, 97.2818110000344, 96.9818010094329)), .Names = c("X20.7", 
"X21.7", "X22.7", "X23.7", "X24.7", "X25.7", "X26.7", "X27.7", 
"X28.7", "X29.7", "X30.7", "X31.7", "X32.7", "X33.7", "X34.7", 
"X35.7"), row.names = c(NA, 6L), class = "data.frame")

このすべてのデータを1つの箱ひげ図にまとめたいのですが、箱ひげ図をプロットしようとすると(つまりboxplot(theData))、Rは列名に基づいて自動的にグループを作成します。

また、完全なデータフレームをベクトルに入れようとしましたが、私の(完全な)データセットにもNA値が含まれているため、これは成功しませんでした。これまでのところ、箱ひげ図にプロットできるように、データフレームのベクトルを作成するための次の関数があります。

for(i in 1:ncol(allTheData)) {
        tmpData <- allTheData[,i]
        for(j in 1:length(tmpData)){
            if(!is.na(j)){
                tmpVector <- c(tmpVector, j)
            }
        }
    }

しかし、私はこの問題を複雑にしすぎていると思います。そのようなループ構造がRのパフォーマンスに役立つかどうかは疑問です。

では、完全なデータフレームの1つの箱ひげ図で構成される箱ひげ図を作成するにはどうすればよいですか?それで、X20.7からX35.7で構成される箱ひげ図は得られませんが、1つの「全体的な」箱ひげ図が得られますか?

4

2 に答える 2

5

このようなものを試してください

boxplot(unlist(theData))
于 2011-02-15T14:09:52.603 に答える
2

ジュラ、

meltの関数を使用しreshapeてデータを「長い」形式に変換し、それを使用boxplotするのはどうですか。データが次の名前のオブジェクトにあると仮定しますdf

> library(reshape)
> df.m <- melt(df)
Using  as id variables
> head(df.m)
  variable    value
1    X20.7 99.64920
2    X20.7 98.76469
3    X20.7 98.37824
4    X20.7 98.11356
5    X20.7 97.80027
6    X20.7 97.47825
> 
> boxplot(df.m$value)
于 2011-02-15T14:08:06.187 に答える