1

我有一个从文件中绘制数据的 R 脚本。该脚本当前对 ylim 使用硬编码值。我想根据绘制的数据动态确定 ylim 的正确（即合理）值。

我正在使用 限制 x 轴值xlim。我原以为 plot 函数将能够计算出 y 值（基于 选择的 x 轴值xlim）——而我不必提供 ylim 参数，但是，当我调用 plot() 时，xlim没有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 回答