2016-07-31 169 views
1

我正在嘗試使用矩陣佈局組合10個圖塊。我幾乎在那裏,但有一個問題,我需要修復。請看看我到目前爲止:R - 在矩陣佈局中組合10個圖塊

# Reproduce dataset 
obs <- c(8.12903225806452, 4.65483870967742, 4.00967741935484, 2.82903225806452, 
1.02258064516129, 4.94838709677419, 2.03548387096774, 4.56129032258065, 
6.66451612903226, 10.4129032258065, 3.7741935483871, 9.93870967741936) 

uncorrected <- structure(list(CanESM2 = c(3.71099209785461, 6.68828105926514, 
4.2675461769104, 4.29092979431152, 6.47999143600464, 6.30958032608032, 
3.28256869316101, 1.06898427009583, 4.01592206954956, 7.88787031173706, 
10.8663206100464, 9.14756298065186), GFDL.ESM2M = c(4.56643295288086, 
8.36136913299561, 4.6934700012207, 4.62750291824341, 6.34791326522827, 
0.916408598423004, 0.111764870584011, 0.570569038391113, 3.05491662025452, 
3.29908037185669, 5.98263263702393, 7.52161026000977), MRI.CGCM3 = c(3.94400858879089, 
2.48861312866211, 0.663281202316284, 0.422377318143845, 1.49548053741455, 
1.8869035243988, 5.35201597213745, 3.94722390174866, 5.0041446685791, 
6.11909627914429, 6.37435054779053, 5.9401068687439), NorESM1.M = c(8.83542060852051, 
8.36209106445312, 3.22844076156616, 2.86131143569946, 2.05858564376831, 
2.39689040184021, 1.06175291538239, 1.03921580314636, 2.49208211898804, 
6.83605766296387, 7.95251750946045, 4.8852391242981), inmcm4 = c(5.88600730895996, 
8.03514385223389, 3.86212635040283, 1.51728367805481, 3.46821308135986, 
1.30993270874023, 0.497318655252457, 0.589299857616425, 1.63358855247498, 
3.61927103996277, 4.99840354919434, 10.5596942901611)), .Names = c("CanESM2", 
"GFDL.ESM2M", "MRI.CGCM3", "NorESM1.M", "inmcm4"), class = "data.frame", row.names = c(NA, 
-12L)) 

corrected <- structure(list(CanESM2 = c(4.69295692443848, 6.2362003326416, 
4.48602199554443, 4.21368360519409, 6.02497625350952, 5.94498348236084, 
3.92529225349426, 1.07059073448181, 3.92306709289551, 8.00215721130371, 
10.9688482284546, 9.95723056793213), GFDL.ESM2M = c(6.74847173690796, 
10.7781000137329, 6.89084196090698, 7.5780816078186, 8.96726703643799, 
3.57111501693726, 1.20403492450714, 2.76221370697021, 4.98387718200684, 
5.27927255630493, 8.87845706939697, 9.95864582061768), MRI.CGCM3 = c(5.45173072814941, 
3.37272930145264, 2.20000958442688, 1.12101686000824, 2.67246603965759, 
3.10278224945068, 6.02392053604126, 4.70185708999634, 5.86025190353394, 
6.41299057006836, 7.1956262588501, 7.59840393066406), NorESM1.M = c(9.49540519714355, 
9.88559532165527, 4.86537599563599, 4.56189870834351, 3.85974431037903, 
4.07706260681152, 2.55634951591492, 2.63215637207031, 4.21092081069946, 
7.66422271728516, 8.75070095062256, 5.95497798919678), inmcm4 = c(8.21576976776123, 
11.1972188949585, 6.08181619644165, 3.58428621292114, 5.543297290802, 
3.46666693687439, 1.88991332054138, 1.94812619686127, 3.71879911422729, 
5.70409107208252, 6.85269260406494, 11.8774194717407)), .Names = c("CanESM2", 
"GFDL.ESM2M", "MRI.CGCM3", "NorESM1.M", "inmcm4"), class = "data.frame", row.names = c(NA, 
-12L)) 

# Create layout 
m <- cbind(c(1,5,9), c(2,6,10), c(3,7,9), c(4,8,10)) 
layout(m) 

# Plot 
for (i in 1:length(names(uncorrected))) { 
    # Plot uncorrected vs obs 
    qqplot(as.vector(obs), uncorrected[,i], xlim=c(0,14), ylim=c(0,14), 
     xlab=expression("Obs. precip." ~ (mm ~ day^{-1})), 
     ylab=expression("Mod. precip." ~ (mm ~ day^{-1}))) 
    abline(0, 1) 
    # Then plot corrected vs obs 
    qqplot(as.vector(obs), corrected[,i], xlim=c(0,14), ylim=c(0,14), 
     xlab=expression("Obs. precip." ~ (mm ~ day^{-1})), 
     ylab=expression("Mod. precip." ~ (mm ~ day^{-1}))) 
    abline(0, 1) 
} 

這幾乎是我所需要的。但是,最後一個子圖被拉伸到剩餘矩陣行的長度。相反,我希望最後一行中的兩個圖形在中心對齊,與其他圖形的長度和寬度相同。

這可能嗎?

回答

4

問題是與您的佈局矩陣m

#  [,1] [,2] [,3] [,4] 
#[1,] 1 2 3 4 
#[2,] 5 6 7 8 
#[3,] 9 10 9 10 

在最後一排,你有9, 10, 9, 10而非9, 9, 10, 10

使用

m <- cbind(c(1,5,9), c(2,6,9), c(3,7,10), c(4,8,10)) 

代替。

enter image description here

我想有兩個地塊的中心對齊的最後一排,在相同的長度和寬度爲別人。

通過保留相同的副圖大小,我們必須分配12個插槽,但僅使用其中的10個。 layout圖開始放置來自索引1的圖,然後按順序放置其餘圖。所以,我們可以這樣做:

m <- cbind(c(1,5,11), c(2,6,9), c(3,7,10), c(4,8,12)) 

#  [,1] [,2] [,3] [,4] 
#[1,] 1 2 3 4 
#[2,] 5 6 7 8 
#[3,] 11 9 10 12 

enter image description here