我使用biwavelet软件包进行小波相干分析。当我想设置我自己的x标签时,我发现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))发布于 2016-07-09 04:31:19
破坏你的阴谋的是plot.cb = TRUE。
在plot.biwavelet的源代码中,作者注意到关于plot.cb选项的以下内容:
添加颜色栏:这必须发生后,一切,否则混乱接踵而来!
这就是问题所在--在plot.cb和混乱之后,你调用了plot.cb。但是,您可以使用image.plot从fields包手动添加颜色条,在没有plot.cb的情况下运行plot,然后添加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)

您可以自定义的滴答和颜色条,以您喜欢这样!
https://stackoverflow.com/questions/38278211
复制相似问题