我正在使用biwavelet
包来进行小波一致性分析。当我想设置自己的x ticklabel时,我发现axis
无效。以下给出了可重复的示例。感谢。
require(biwavelet)
t1 <- cbind(1:100, rnorm(100))
t2 <- cbind(1:100, rnorm(100))
wtc.t1t2 <- wtc(t1,t2,nrands = 10)
plot(wtc.t1t2, plot.cb = TRUE, plot.phase = TRUE,xaxt='n')
axis(1,at = seq(10,100,10),labels = seq(1,10,1))
答案 0 :(得分:2)
破坏你情节的是plot.cb = TRUE
。
在plot.biwavelet
的源代码中,作者注意到以下plot.cb
选项:
##添加颜色条:这必须在一切之后发生,否则会发生混乱 随之而来!
这就是问题所在 - 您在axis()
之后调用了plot.cb
并且随后发生了混乱。 但,您可以使用image.plot
包中的fields
手动添加颜色条,运行plot
而不plot.cb
然后添加了axis()
。
pacman::p_load(biwavelet,fields)
t1 <- cbind(1:100, rnorm(100))
t2 <- cbind(1:100, rnorm(100))
wtc.t1t2 <- wtc(t1,t2,nrands = 10)
plot(wtc.t1t2, plot.phase = TRUE,xaxt='n')
axis(1,at = seq(10,100,10),labels = seq(1,20,2))
image.plot( zlim=c(0,25), legend.only=TRUE)
您可以按照自己喜欢的方式自定义刻度线和颜色条!