1

私はこのようなデータセットを持っています

 VegType    87MIN   87MAX   87Q25   87Q50   87Q75   96MIN   96MAX   96Q25   96Q50     96Q75 00MIN   00MAX   00Q25   00Q50   00Q75
 1          0.02    0.32    0.11    0.12    0.13    0.02    0.26    0.08    0.09    0.10    0.02    0.28    0.10    0.11    0.12
 2          0.02    0.45    0.12    0.13    0.13    0.02    0.20    0.09    0.10    0.11    0.02    0.26    0.11    0.12    0.12
 3          0.02    0.29    0.13    0.14    0.14    0.02    0.27    0.11    0.11    0.12    0.02    0.26    0.12    0.13    0.13
 4          0.02    0.41    0.13    0.13    0.14    0.02    0.58    0.10    0.11    0.12    0.02    0.34    0.12    0.13    0.13
 5          0.02    0.42    0.12    0.13    0.14    0.02    0.46    0.10    0.11    0.11    0.02    0.28    0.12    0.12    0.13
 6          0.02    0.32    0.13    0.14    0.14    0.02    0.52    0.12    0.12    0.13    0.02    0.29    0.13    0.14    0.14
 7          0.02    0.55    0.12    0.13    0.14    0.02    0.24    0.10    0.11    0.11    0.02    0.37    0.12    0.12    0.13
 8          0.02    0.55    0.12    0.13    0.14    0.02    0.19    0.10    0.11    0.12    0.02    0.22    0.11    0.12    0.13

実際には、26個の変数と5年があります（列名の87、96、00は年です）。理想的な世界では、変数ごとに1つずつ、26のプロットがあり、各プロットに5つのボックス、つまり1年に1つのボックスが含まれる格子状のグラフが必要です。ラティスは事前定義された統計を受け入れないため、これをラティスで行うことはできないことを理解しています。事前定義された統計を使用してRでこれを行うためのかなり苦痛のない方法はありますか？私はbxp、1つのプロットで1年間すべての変数をプロットする単純な箱ひげ図に使用しました。

Yr01 = read.csv('dat.csv',header=T)
dat01=t(Yr01[,c("01Min","01Q25","01Mean","01Q75","01Max")])
bxp(list(stats=dat01, n=rep(26, ncol(dat01))),ylim=c(0.07,0.2))

でもそこから必要なものに行く方法がわかりません。

ありがとう。

4

1 に答える 1

5

これは、少なくともを使用してggplot2実行できますがreshape、データをかなり処理する必要があります。そして、分位数が実際に意味をなすデータが本当に必要です!! 分位値がすべて台無しになっています！たとえば、Var1持っている01Max = 0.26と01Q75 = .67!!

まず、有効なデータを再作成します。

n  <- c("01Min", "01Max", "01Med", "01Q25", "01Q75", "02Min", 
                            "02Max", "02Med", "02Q25", "02Q75")
v1 <- c(0.03,  0.76,  0.41,  0.13,  0.67,  0.10,  0.43,  0.27,  0.2,   0.33)
v2 <- c(0.03,  0.28,  0.14,  0.08,  0.20,  0.02,  0.77,  0.13,  0.06, 0.44)

df <- data.frame(v1=v1, v2=v2)
df <- as.data.frame(t(df))
names(df) <- n
df <- cbind(var=c("v1","v2"), df)
> df

#    var 01Min 01Max 01Med 01Q25 01Q75 02Min 02Max 02Med 02Q25 02Q75
# v1  v1  0.03  0.76  0.41  0.13  0.67  0.10  0.43  0.27  0.20  0.33
# v2  v2  0.03  0.28  0.14  0.08  0.20  0.02  0.77  0.13  0.06  0.44

次に、データの形状を変更します。

require(reshape2)
df.m <- melt(df, id="var")
# look for a bunch of numbers from the start of the string and capture it
# in the first variable: () captures the pattern. And replace it with the 
# captured pattern with the variable "\\1"
df.m$year <- gsub("^([0-9]+)(.*$)", "\\1", df.m$variable)

# the same but instead refer to the captured pattern in the second 
# paranthesis using "\\2"
df.m$quan <- gsub("^([0-9]+)(.*)$", "\\2", df.m$variable)
df.f <- dcast(df.m, var+year ~ quan, value.var="value")

この形式にするには：

> df.f

#   var year  Max  Med  Min  Q25  Q75
# 1  v1   01 0.76 0.41 0.03 0.13 0.67
# 2  v1   02 0.43 0.27 0.10 0.20 0.33
# 3  v2   01 0.28 0.14 0.03 0.08 0.20
# 4  v2   02 0.77 0.13 0.02 0.06 0.44

column namesこれで、対応する次のように、対応するパラメーターに分位値を直接提供することでプロットできます。

require(ggplot2)
require(scales)
p <- ggplot(df.f, aes(x=var, ymin=`Min`, lower=`Q25`, middle=`Med`, 
                           upper=`Q75`, ymax=`Max`)) 
p <- p + geom_boxplot(aes(fill=year), stat="identity") 
p

ggplot1

# if you want facetting:
p + facet_wrap( ~ var, scales="free")

ggplot2

これで、このコードを使用して、次のように、yearsそれぞれのすべてをvar個別のプロットにプロットするタスクを実行できます。lapplysubsetting

lapply(levels(df.f$var), function(x) {
    p <- ggplot(df.f[df.f$var == x, ], 
            aes(x=var, ymin=`Min`, lower=`Q25`, 
                middle=`Med`, upper=`Q75`, ymax=`Max`))
    p <- p + geom_boxplot(aes(fill=year), stat="identity")
    p
    ggsave(paste0(x, ".pdf"), last_plot())
})

編集：あなたのデータは、いくつかの面であなたが提供した以前のデータとは異なります。したがって、新しいデータのコードのバージョンは次のとおりです。

# change var to VegType everywhere
require(reshape2)
df.m <- melt(df, id="VegType")

df.m$year <- gsub("^X([0-9]+)(.*$)", "\\1", df.m$variable) # pattern has a X
df.m$quan <- gsub("^X([0-9]+)(.*)$", "\\2", df.m$variable) # pattern has a X
df.f <- dcast(df.m, VegType+year ~ quan, value.var="value")
df.f$VegType <- factor(df.f$VegType) # convert integer to factor

require(ggplot2)
require(scales)
p <- ggplot(df.f, aes(x=VegType, ymin=`MIN`, lower=`Q25`, middle=`Q50`, 
                           upper=`Q75`, ymax=`MAX`)) 
p <- p + geom_boxplot(aes(fill=year), stat="identity") 
p

以前と同じコードを使用して、個別のプロットとしてファセット/書き込みを行うことができます。

于 2013-02-07T12:44:28.050 に答える