我试图用x中的对数轴创建一个图,从10 ^ -8到10 ^ -3,只标记其他每个刻度,y轴从0到100。 /> 我无法标记每个其他标记,并且该功能不会出现在图表上 有谁知道我哪里出错了?
c4=1e-6
c5=2
K3=1e-6
R1=0.02
R2=0.05
R3=0.2
g <- function(x,c5,K3,R,c4)
{
z = (c5*K3*x)/(R*c4-c5*x)
return(z)
}
z1 <- g(x, c5, K3, R1, c4)
z2 <- g(x, c5,K3,R2,c4)
z3 <- g(x, c5,K3,R3,c4)
plot(x, z1, log="x", lty=3, xaxt="n", yaxt="n", xlim=c(10^-8, 10^-3), ylim=c(0, 100),
main="Fig 2B",
xlab="[Ca2+] in cytosol (M)",
ylab="[InsP3] (nM))")
axis(side=2, at= seq(0,100, by=10), labels= c(0, FALSE, 20, FALSE, 40, FALSE, 60, FALSE, 80, FALSE, 100), srt=90)
xticks <- seq(-7, -3, by=2)
xlabels <- sapply(xticks, function(i) as.expression(bquote(10^ .(i))))
axis(side=1, at=c(10^-7, 10^-5, 10^-3), labels=xlabels)
lines(x, z2, lty=2)
lines(x, z3, lty=1)
答案 0 :(得分:0)
以下是此问题的两种解决方案:
1)将每个其他标签设置为NA(当你尝试时不是FALSE)
plot(1:10, axes = FALSE)
axis(1, at = 1:10, labels = c(NA, 2, NA, 4, NA, 6, NA, 8, NA, 10))
2)添加轴两次。一旦没有标签就可以得到刻度线,一次带有你想要的标签
plot(1:10, axes = FALSE)
axis(1, at = 1:10, labels = FALSE)
axis(1, at = seq(2, 10, 2))