(reproducible exampleがなければ、それはこれで多くを行うのは難しいです。)
簡単な解決策はで1つの散布図を作り、その後、作られたプロットを受け入れるために他の変数をループをPDFファイルを開くことであろう時間。条件を設定したい因子の2つの異なるレベルを取る観測値を示すために、異なる記号と色を使用することができます。
pdf(file=<something>, <other settings>)
for(i in 1:13){ # or perhaps: for(i in c(1,2,5,10,11,12,13)){
plot(AF[,i], AF[,14], pch=<something>, col=<something>)
}
dev.off()
アップデート:私は、一例としてiris
データセットを使用して上記の私の答えを説明することができます。
data(iris)
add.legend <- function(...){ # taken from: https://stackoverflow.com/a/21784009/1217536
opar <- par(fig=c(0, 1, 0, 1), oma=c(0, 0, 0, 0),
mar=c(0, 0, 0, 0), new=TRUE)
on.exit(par(opar))
plot(0, 0, type='n', bty='n', xaxt='n', yaxt='n')
legend(...)
}
pdf(file="Scatterplots of variables against petal width.pdf")
for(i in 1:3){
plot(iris[,i], iris[,4], pch=as.numeric(iris$Species), col=as.numeric(iris$Species),
xlab=names(iris)[i], ylab="Petal Width")
add.legend("top", legend=levels(iris$Species), horiz=TRUE, pch=1:3, col=1:3)
}
dev.off()
新しいページ上の各プロットでPDFファイルを取得します。 @ Hack-Rと@Alexの答えは問題がありません(より良い言葉がないので、あまりにも厳しくないというわけではありません)。プロットは3つで絞ればOKですが、13のときは非常に有益な情報にはなりません。
「AF」というオブジェクトはどこにありますか?再現可能な例が必要です。 http://stackoverflow.com/help/mcveそれはまた、あなたが何を求めているのか不明瞭です。 –
@ Hack-Rに追加する:あなたはオブジェクトを検査して、プロットの完全な行列を与える必要があります(私はあなたの最後の文から理解しています、その行列のスライスが一番右の列であることを理解してください)これをさらに処理します.... – Dilettant
以前使用していたAFデータセットは、再現性のために 'iris'データセットに変更されました。お問い合わせ –