R: brauche endliche 'ylim' Werte in Funktion

Question

Ich möchte die Daten in data.frame xy für jede Gruppe (definiert durch ID) plotten. Wenn ein Jahr vor 1946 in einer Gruppe ist, sollte plot 2 ausgeführt werden. Wenn die Jahre zwischen 1946 und 2014 liegen, sollte plot1 ausgeführt werden.

Mein Problem: Dies funktioniert gut ohne NA-Werte, aber da ich Datenlücken habe, verlasse ich mich auf NAs, um diese Datenlücken zu definieren. Deshalb erhalte ich einen Fehler: error in plot.window(need finite 'ylim' values). Ich habe versucht, finite=T in plot1 auf der y-Achse zu setzen, aber dies ergibt einen subscript out of bounds Fehler. Gibt es eine Möglichkeit, das zu lösen und die Grafiken korrekt zu zeichnen?

Im Folgenden ist mein Code: Es ist lang, aber der Großteil des Codes besteht aus plot() Optionen, auf die ich verlasse.

# read in sample data and split it up by group (defined by ID) 
xy <- data.frame(NAME=c("NAME2","NAME2","NAME2","NAME2","NAME2","NAME3","NAME3","NAME3","NAME3","NAME5","NAME5","NAME5","NAME5"), ID=c(48,48,48,48,48,32,32,32,32,67,67,67,67),YEAR=c(1981,1983,1984,1988,1989,1984,1984,1988,1988,1899,1933,1948,1958),VALUE=c(0,205,-570,0,-310,-3680,-3680,NA,-3680,0,NA,13,-98)) 
ind <- split(x = xy,f = xy[,'ID']) 

# Plot Scenario 1: if only years between 1946 and 2014 are present for each group do this: 
    plot1 <- function(x) { 
    fname <- paste0(x[1, 'ID'], '.png') 
    png(fname, width=1679, height=1165, res=150) 
    par(mar=c(6,8,6,5)) 
    plot(x = c(1946, 2014), 
     y = range(x$VALUE), 
     type='n', 
     main=x[1, 'NAME'], 
     xlab="Time [Years]", 
     ylab="Value") 
    axis(2, at = seq(-100000, 100000, 500), cex.axis=1, labels=FALSE, tcl=-0.3) 
    points(ind[[i]][,c('YEAR','VALUE')], type="l", lwd=2) 
    points(ind[[i]][,c('YEAR','VALUE')], type="p", lwd=1, cex=1, pch=21, bg='white') 
    abline(h=0) 
    dev.off() 
} 

# Plot Scenario 2 if years under 1946 are present do this: 
plot2 <- function(x) { 
    fname <- paste0(x[1, 'ID'], '.png') 
    png(fname, width=1679, height=1165, res=150)  
    par(mar=c(6,8,6,5)) 
    plot(x[,c('YEAR','VALUE')], 
     type='n', 
     main=x[1, 'NAME'], 
    xlab="Time [Years]", 
    ylab="Value [mm]") 
axis(2, at = seq(-100000, 100000, 500), cex.axis=1, labels=FALSE, tcl=-0.3) 
points(ind[[i]][,c('YEAR','VALUE')], type="l", lwd=2) 
points(ind[[i]][,c('YEAR','VALUE')], type="p", lwd=1, cex=1, pch=21, bg='white') 
abline(h=0) 
dev.off() 
} 

# Execute functions 
    lapply(ind, function(x) ifelse(any(x$YEAR < 1946 & x$YEAR < 2014), plot2(x), plot1(x))) 

Answer 1

In plot1 ändern y = range(x$VALUE)-y = range(x$VALUE, na.rm=TRUE)NA fehlende Werte zu entfernen.

Sie haben ein anderes Problem, nämlich, dass Sie in beiden Funktionen auf ind[[i]] verweisen, was bedeutet, dass an einem bestimmten Punkt dieser Code Teil einer Schleife war. Meine Vermutung ist, alle Verweise auf ind[[i]] sollten x sein.