我对image
包中的R函数biwavelet
有疑问。此函数用于将相关值绘制为z
函数中的image
。我正在绘制完全相关的数据集(R平方= 1作为我的z值)。绘制图像时,彩色框中存在一些变化。我对此感到困惑,因为我z
矩阵中的所有值都是1.任何想法为什么会发生这种情况,以及如何解决这个问题?
image(x$t, yvals, t(zvals), zlim = zlim, xlim = xlim, ylim = rev(ylim),
xlab = "Time", ylab = "Period", yaxt = "n", xaxt = "n", col = fill.colors)
zvals
是所有1的矩阵
更新:根据以下要求
x$t
[1] 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000
[29] 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012
yvals
[1] 1.046901 1.130235 1.213568 1.296901 1.380235 1.463568 1.546901 1.630235 1.713568 1.796901 1.880235 1.963568 2.046901 2.130235 2.213568
[16] 2.296901 2.380235 2.463568 2.546901 2.630235 2.713568 2.796901 2.880235 2.963568 3.046901 3.130235 3.213568 3.296901 3.380235 3.463568
[31] 3.546901 3.630235 3.713568 3.796901
str(zvals)
num [1:34, 1:40] 1 1 1 1 1 1 1 1 1 1 ...
以上是实际数据。下面是一个示例,但是按预期绘制,而我的数据不是:
a <- 1980:2010
y <- sort( rnorm(30,2,.5) )
z <- matrix(nrow=30,ncol=31,1)
image(a, y, t(z), zlim = range(z), xlim = range(a), ylim = rev(range(y)), col = heat.colors(n=20))
答案 0 :(得分:0)
结果显示zlim的biwavelet默认值导致问题在于颜色缩放沿着非常小的渐变 - (例如0.99-1)改变,这在图中施加了不同的颜色。我假设z值显示的精度小于R存储的精度,这是因为值之间的微小差异?解决方案是为zlim(例如0-1)施加更宽的范围,然后迫使颜色梯度方案适当地打印(单色)。