3

この質問はかなり漠然としていますが、R で biwavelet パッケージを使用して成功した人はいますか? 次のコードがあります。

require(biwavelet)
t1 <- cbind(DecTime,Temp)    
## continuous wavelet transform
wt1 <- wt(t1)
plot(wt1)

そしてそれはエラーを返します：

image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)), のエラー: 'x' と 'y' の値は有限であり、欠落していない必要があります

私のデータには欠損値が含まれておらず、すべて有限であるため、このエラーを理解できません。

例を再現しようとすると：

require(biwavelet)
Date = seq(from=as.POSIXct("2011-01-01 00:00"), 
           to=as.POSIXct("2011-12-31 23:00"), length=8760)
DecTime = julian(Date, Date[1])

data=cbind(as.numeric(DecTime), rnorm(8760))
## Continuous wavelet transform
wt.t1=wt(data)
plot(wt.t1)

面倒くさいので、このパッケージでデータが失敗する理由がわかりません。上記の例のデータは、私のデータ、つまり「マトリックス」と同じクラスです。ヘルプやアドバイスをいただければ幸いです。

変更されました：

require(biwavelet)
Date = seq(from=as.POSIXct("2011-01-01 00:00"), 
           to=as.POSIXct("2011-12-31 23:00"), length=8760)
DecTime = julian(Date, Date[1])
D <- c(4.0267, 4.0211, 4.0005,4.0042,4.0042,4.0191)
data=cbind(as.numeric(DecTime[1:6]),as.numeric(D))
## Continuous wavelet transform
wt.t1=wt(data)
plot(wt.t1)


> data
           [,1]   [,2]
[1,] 0.00000000 4.0267
[2,] 0.04166667 4.0211
[3,] 0.08333333 4.0005
[4,] 0.12500000 4.0042
[5,] 0.16666667 4.0042
[6,] 0.20833333 4.0191
> class(data)
[1] "matrix"
> class(data[,1])
[1] "numeric"
> class(data[,2])
[1] "numeric"
> wt.t1=wt(data)
> plot(wt.t1)
Error in image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)),  : 
  invalid z limits
4

2 に答える 2

1

他の投稿への返信を見逃した方のために、wt.R 関数のバグ (誤った括弧) を発見したことをお知らせします。このバグはバージョン 0.12 の biwavelet パッケージで修正されているため、コードは動作するようになります。

エラーを見つけてくれてありがとう！

于 2012-10-07T07:19:28.030 に答える
0

これは、少し光を当てることを望んでいる部分的な答えです。動作する例からの最初のN行を取得するdataと、 のときに失敗しN < 2762ます。次の場合に機能しますN >= 2762(おそらく、このしきい値は ? の実際のランダム値に依存しdataます):

> wt.t1=wt(data[1:2762,])
> plot(wt.t1, plot.cb=T)
> wt.t1=wt(data[1:2761,])
> plot(wt.t1, plot.cb=T)
Error in image.default(x$t, yvals, t(zvals), zlim = zlims, ylim = rev(range(yvals)),  : 
  invalid z limits

重要な点は、すべてのwt.t1$power値が 1 未満であることです。これにより、zlims変数の計算方法が混乱します。とはいえ、何が何だかさっぱりわかりませんwt.t1$power。

于 2012-06-25T18:31:31.023 に答える