2

我想做一些非常简单的事情:我想为一个完整的数据框创建一个箱线图。然而,搜索“组合箱线图”和相关术语并没有提出任何建议。如果我忽略了一个明显的方式,请告诉我。

我有以下数据:

> 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")

我希望将所有这些数据汇总在一个箱线图中,但是,当我尝试绘制箱线图(即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 的性能。

那么,如何制作一个包含一个完整数据框的箱线图的箱线图?那么,我没有得到一个由 X20.7 到 X35.7 组成的箱线图,而是给出了一个“整体”箱线图?

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 回答