私は10列のcsvファイルを持っており、それらの7つをテキストファイルの時系列と一緒に読み込んでプロットしたいと思っています。次のコードでは、私はそれらを描画することができますが、いくつかの問題に直面しています:rとパネルプロット
1]プロット間の空白は各プロットを非常に小さくしています。ほとんどのスペースの使い方は? 2]凡例が7番目のプロットにオーバーレイされています。共通のxおよびy軸を作り、各パネルプロット
の個々のタイトルを定義する方法]3(パネルプロット領域8および9を使用していないので)8と9のパネルのプロットエリアにそれを配置する方法
library(qmap)
obs <- read.table("D:/test.txt")
mod <- read.csv("D:/pr.csv", header = FALSE)
par(mfrow=c(3,3), oma = c(4, 1, 1, 1))
for(i in 4:10){
modp <- mod[,i]
qm1.fit <- fitQmap(obs,modp,
method="PTF",
transfun="expasympt",
cost="RSS",wett.day=TRUE)
qm1 <- doQmap(modprecip,qm1.fit)
# repeated same procedure for qm2, qm3, qm4, and qm5
# QQ-plot of all mappings
obs_v <- as.vector(as.matrix(obs))
mod_v <- as.vector(as.matrix(modp))
qqplot(mod_v, obs_v, col = 1, xlab = 'Pm [mm/day]', ylab = 'Po [mm/day]')
qqline(qm1, col = 2, lwd = 2)
qqline(qm2, col = 3, lwd = 2)
qqline(qm3, col = 4, lwd = 2)
qqline(qm5, col = 5, lwd = 2)
}
legend("bottom", inset=c(-0.2,0), legend=c("A","B"), pch=c(1,3), title="Group")
すべてのヘルプは、プロットのこの種のは非常に簡単にggplot2
を使用して生成することができます
結果のデータや数字はありません。私たちは何をすると思われますか? –
@RomanLuštrikは図を追加 – Ibe
これは['ggplot2 :: facet_wrap'](http://docs.ggplot2.org/current/facet_wrap.html)を使って簡単に解決できます。 –