我有一个 5x5 矩阵示例,其中包含以下值:
dat <- matrix(seq(1,13,0.5), nrow=5, byrow=TRUE)
dat[seq(2,25,2)] <- NA
1 | NA | 2 | NA | 3
NA | 4 | NA | 5 | NA
6 | NA | 7 | NA | 8
NA | 9 | NA | 10 | NA
11 | NA | 12 | NA | 13
由于缺少值,我一生都无法使用 persp3d() 等方法获得 3D 曲面图。有没有一种方法可以让 R 只对值进行插值并仍然绘制它们?
虽然以前可能有人问过这个问题,但我找不到这种情况的简洁示例。尝试一下,这将给出一个结果,然后您可以将其传递给persp
, image
etc...
#install.packages("akima")
library(akima)
nas <- !is.na(dat)
interp(
row(dat)[nas], #row index - 'x' values
col(dat)[nas], #col index - 'y' values
dat[nas], #height data - 'z' values
xo=seq(1,nrow(dat)), #'x' values for output
yo=seq(1,ncol(dat)) #'y' values for output
)
#$x
#[1] 1 2 3 4 5
#
#$y
#[1] 1 2 3 4 5
#
#$z
# [,1] [,2] [,3] [,4] [,5]
#[1,] 1.0 1.5 2.0 2.5 3.0
#[2,] 3.5 4.0 4.5 5.0 5.5
#[3,] 6.0 6.5 7.0 7.5 8.0
#[4,] 8.5 9.0 9.5 10.0 10.5
#[5,] 11.0 11.5 12.0 12.5 13.0
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)