2016-05-18 12 views
0

私は、毎時のデータセットに対して作業しているコードがあり、適切なウェーブレットプロットを表示したいのですが、好きなように調整しています。私は、x軸、y軸を変更し、カラーバーを設定したいと思います。ただし、カラーバーがアクティブ化されていない場合にのみ、x軸を正常に変更しました。そうであれば、私はx軸を変更できません。さらに、私はy軸を変更して、自動的に生成される値の間にさらに多くの値を持たせることができます。あなたの助けウェーブレットプロット:x軸、y軸、およびカラープロットの変更

library(biwavelet) # used for wavelets 

    n <- 141696 
    d <- data.frame(1:n, round(runif(n, 38, 100),2)) 

    # X-Axis for plotting 
    TIME1 <- as.POSIXlt("2000-01-01 00:00:00 PST", format = '%Y-%m-%d %H:%M:%S') 
    TIME2 <- as.POSIXlt("2016-02-29 23:00:00 PST", format = '%Y-%m-%d %H:%M:%S') 
    LABELS <- seq(from = TIME1, to = TIME2, by = "3 months") 
    xAxis <- seq(from = TIME1, to = TIME2, by = "hour") 
    Location <-NA 
    for (i in 1:length(LABELS)) { Location[i] <- which(LABELS[i] == xAxis) } 
    LABELS <- format(LABELS, "%b %Y") 

    # Wavelet 
    WAV <- wt(d) 

これは正しいx軸を持っていますが、私は、引数としてプロットでplot.cb = TRUEを入れていなかったので、カラーバーを表示しないために、事前にありがとうございます。

# PLOT (Has no legend but correct x-axis) 
    par(oma=c(0, 0, 0, 1), mar=c(5, 4, 4, 5) + 0.1) 
    plot(WAV, type="power.corr.norm", main="Bias-corrected wavelet power ", ylab="Period(hourly)", xlab="Time", lwd.sig=1, xaxt='n') 
    axis(side = 1, at = Location, labels = LABELS, tick = TRUE, las = 2) 

No Color Bar, correct x-axis

これは、カラーバーが表示されますが、適切なx軸はラベルがありません。

# PLOT (Has legend but no x-axis) 
    par(oma=c(0, 0, 0, 1), mar=c(5, 4, 4, 5) + 0.1) 
    plot(WAV, type="power.corr.norm", main="Bias-corrected wavelet power ", ylab="Period(hourly)", xlab="Time", lwd.sig=1, xaxt='n', plot.cb=TRUE) 
    axis(side = 1, at = Location, labels = LABELS, tick = TRUE, las = 2) 

No x-axis, but color bar present

答えて

0

plot.biwaveletを入力して、機能自体に取得スクリプトエディタにコピーし、機能を編集し、それをこのようmyplotなどの新しい名前を付け、myplot代わりのplot.biwaveletを実行し、それを編集します。 plot.biwaveletファンクション内の任意のものを変更できます。たとえば、xlimラベルの数を増やすには、locs <- pretty(range(xlim), n = 10)を実行します。現在、n=5

関連する問題