1

ファイルからデータをプロットする R スクリプトがあります。スクリプトは現在、ハードコードされた値を ylim に使用しています。プロットされているデータに基づいて、ylim の正しい (つまり、適切な) 値を動的に決定したいと考えています。

を使用して x 軸の値を制限していますxlim。私は、プロット関数が(によって選択されたx軸の値に基づいて)y値を計算できると考えていたでしょうxlim-ylim引数を提供する必要はありませんが、引数xlimなしでplot()を呼び出したときylim、次のエラーが発生しました。

Error in plot.window(...) : need finite 'ylim' values
Calls: plot -> plot.default -> localWindow -> plot.window
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf

私の質問は、 ?ylimに制限を指定した場合に、に指定する値を動的に決定するにはどうすればよいかということxlimです。ylim理想的には、次のように制限を指定したいと思います。

ylim_lower <- (y value for xlim_lower) - [some fixed % distance]
ylim_upper <- (y value for xlim_upper) + [some fixed % distance]

どうすればこれを行うことができますか?

[[編集]]

dat <- read.csv(somefile)

n <- dim(dat)[1]
yvals1 <- rep(0,n)
yvals2 <- rep(0,n)

for(i in 1:n){
  yvals1[i] <- foobar1(dat$X[i])
  yvals2[i] <- foobar2(dat$X[i])
}   

# Note: yvals1 and yvals2 MAY contain NAs

# below is the plot command that barfs:
plot(dat$X, yvals1, typ="l", col="green", xlim=c(lowest_val_cuttoff, highest_val_cuttoff), ylim= c(.2, .60), main=c(the_title, "Title goes here"), xlab="x axis label", ylab="y axis label")
lines(dat$X, yvals2, col="red")
4

1 に答える 1

1

プロットする前に、データセット (または少なくともプロット範囲を見つけるための呼び出し) を有限で欠損値のない値に制限してみてください。

> plot(1:2,1:2, ylim=c(NA,3))
Error in plot.window(...) : need finite 'ylim' values
> plot(1:2,1:2, ylim=c(-Inf,3))
Error in plot.window(...) : need finite 'ylim' values

?is.na ?is.finite

> min(c(NA,3,5))
[1] NA
> min(c(NA,3,5),na.rm=TRUE)
[1] 3
> min(c(-Inf,3,5),na.rm=TRUE)
[1] -Inf
> 
> y <- c(-Inf,NA,3,4,5)
> range(y)
[1] NA NA
> range(y, na.rm=TRUE)
[1] -Inf    5
> range(y[!is.na(y) & is.finite(y)])
[1] 3 5
> 
于 2012-08-01T22:21:36.263 に答える