半对数图轴标记和功能不出现在R中

时间:2016-04-16 08:58:53

标签: r

我试图用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)

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))