R/clean.station.R

In metno/esd: Climate analysis and empirical-statistical downscaling (ESD) package for monthly and daily data

#' Remove stations with a lot of missing data
#'
#' @param x a \code{station} object
#' @param miss fraction of data that are allowed to be missing,
#' e.g., miss=.1 means that stations with more than 10\% missing data
#' will be excluded
#' @param verbose a boolean; if TRUE print information about progress
#'
#' @export
clean.station <- function(x,miss=.1,verbose=TRUE) {
  if(verbose) print("clean.station")
  stopifnot(inherits(x,"station"))
  loc <- loc(x)
  if (verbose) par(mfrow=c(2,1))
  if (verbose) image(coredata(annual(x)))
  if (verbose) title('Original data')
  # compute number of finite values
  n <- apply(coredata(x),2,FUN=nv)

  if (verbose) summary(n)
  keep <- n > (dim(x)[1]*(1-miss))
  x <- subset(x,is=keep)
  if (verbose) print("Removed stations ...")
  if (verbose) print(loc[!keep])
  if (verbose) image(coredata(annual(x)))
  if (verbose) title('Filtered data')
  return(x)
}