Как создать комбинированный коробочный график из кадра данных? - PullRequest
Новая
       11

Как создать комбинированный коробочный график из кадра данных?

2 голосов
/ 15 февраля 2011

Я хочу сделать что-то невероятно простое: я хочу создать один блокпост для полного кадра данных. Тем не менее, поиск «комбинированного коробочного графика» и связанных с ним терминов не привел к появлению каких-либо предложений. Если я упустил очевидный способ, дайте мне знать.

У меня есть следующие данные: 

> 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, но дает один "Общий" боксплот?

Ответы [ 2 ]

5 голосов
/ 15 февраля 2011

Попробуйте что-то вроде этого 

boxplot(unlist(theData))
2 голосов
/ 15 февраля 2011

Jura,

Как насчет использования функции 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)
...