プロットされたデータが第1(I)象限ではなく、デカルト平面の第4象限( IV )に表示されるように、y軸が反転されたRでプロットを作成する方法を知りたいです。
参考までに、私が作成しようとしているプロットは、次のようになります(ソース)。
y軸の番号を逆にすることに関して、オンラインでいくつかの質問を見つけましたが、これらはすべて、データを第1象限にプロットしています。上に示したようなプロットを作成する方法を誰かが提案できますか?
@timriffeと@joranのコメントに従って、うまくいった答えを提供するだけです...
この回答のマイナー ログ ティックの関数を使用します。
minor.ticks.axis <- function(ax,n,t.ratio=0.5,mn,mx,...){
lims <- par("usr")
if(ax %in%c(1,3)) lims <- lims[1:2] else lims[3:4]
major.ticks <- pretty(lims,n=5)
if(missing(mn)) mn <- min(major.ticks)
if(missing(mx)) mx <- max(major.ticks)
major.ticks <- major.ticks[major.ticks >= mn & major.ticks <= mx]
labels <- sapply(major.ticks,function(i)
as.expression(bquote(10^ .(i)))
)
axis(ax,at=major.ticks,labels=labels,...)
n <- n+2
minors <- log10(pretty(10^major.ticks[1:2],n))-major.ticks[1]
minors <- minors[-c(1,n)]
minor.ticks = c(outer(minors,major.ticks,`+`))
minor.ticks <- minor.ticks[minor.ticks > mn & minor.ticks < mx]
axis(ax,at=minor.ticks,tcl=par("tcl")*t.ratio,labels=FALSE)
}
再現可能なサンプル データを作成します。
x <- 1:8
y <- 10^(sort(runif(8, 1, 10), decreasing = TRUE))
軸なしでプロットする:
plot(x, log10(y), # function to plot
xlab="", # suppress x labels
type = 'l', # specify line graph
xlim = c(min(x), (max(x)*1.3)), # extend axis limits to give space for text annotation
ylim = c(0, max(log10(y))), # ditto
axes = FALSE) # suppress both axes
ファンシーな対数軸を追加し、目盛りラベルを右に向けます (@joran に感謝します!):
minor.ticks.axis(2, 9, mn=0, mx=10, las=1)
上に x 軸を追加します。
axis(3)
x 軸ラベルを追加します (ヒントをありがとう、@WojciechSobala)
mtext("x", side = 3, line = 2)
そして、行末に注釈を追加します
text(max(x), min(log10(y)), "Example", pos = 1)
結果は次のとおりです。