无法在MASS包中绘制内核密度 - 在min(x)中:min没有非缺失参数;返回Inf

时间:2012-07-20 11:34:56

标签: r plot kernel-density

在尝试绘制数据时,我收到以下错误消息:

1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf

这个可重复的代码有效:

# some dummy data
library(mvtnorm)
dat <- rmvnorm(10000, mean=c(x=4,y=4), sigma=matrix(c(1,0.5,0.5,1), ncol=2))
dat <- as.data.frame(dat)

# 2d density plot
library(MASS)
kdexy <- kde2d(dat$x,dat$y, n=50)
image(kdexy, col=grey(seq(1,0.2,length=10)))

但是使用我的真实数据不会:

kdexy <- kde2d(temp$V1, temp$V2, n=50)
image(kdexy, col=grey(seq(1,0.2,length=10)))

然而两个数据集的结构是相同的(虚拟数据[dat]和真实数据[temp]):

> str(dat$x)
 num [1:80000] 0.669 0.609 -0.633 0.565 0.559 ...
> str(temp$V1)
 num [1:823180] 0 0 0 0 0.0146 ...
> summary(dat$x)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
-4.2270 -0.1902  0.4841  0.4900  1.1600  5.3570 
> summary(temp$V1)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
-1.0000  0.0000  0.0000 -0.0289  0.0000  0.9844 
> range(dat$x)
[1] -4.227400  5.357184
> range(temp$V1)
[1] -1.000000  0.984375
> str(temp$V2)
 num [1:823180] 1 1 15.5 15.5 18.5 ...
> summary(temp$V2)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
   1.00    1.00   18.45   23.55   35.96  116.10 
> range(temp$V2)
[1]   1.0000 116.0829

它们都存储在数据框中,我所知道的唯一区别是长度,temp$V1-11有界。

kdexy <- kde2d()的输出在两个数据集之间有所不同。在示例数据中,Z部分填充了非常小的数字;在真实数据集中,每个点都填充'NaN'。

1 个答案:

答案 0 :(得分:1)

我的数据主要集中在0(total n=850,000n==0 ~ 650,000)。使用x<0n=150,000)的数据子集,图表可以正常工作。