1

データをプロットしようとすると、次のエラー メッセージが表示されます。

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)))

ただし、2 つのデータ セットの構造は同じです (ダミー データ [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は と で制限され-1て1います。

の出力はkdexy <- kde2d()、2 つのデータセット間で異なります。サンプル データでは、Zセクションに非常に小さな数が入力されています。実際のデータセットでは、すべてのポイントが「NaN」で満たされています。

4

1 に答える 1

1

total n=850,000私のデータは 0 ( , )付近に集中していますn==0 ~ 650,000。データのサブセットを使用すると、x<0( n=150,000) でプロットが機能します。

于 2012-07-27T11:35:09.030 に答える