weiming-hu / ranenextra

This is a R package that contains helpful functions for the RAnEn package.

Hi there 👋

My name is Weiming Hu. I'm a researcher in ⚡ Renewable Energy ⚡ and ☁️ Weather Forecasting ☁️ .

I'm also an enthusiast for free software and the open-source movement.

• 🚧 I'm working on 🚧

  • Developing R and Python packages for the scientific community
  • Becoming a Bangle.JS contributor.

• 💬 Ask me about 💬

  • R, Python, C++, Supercomputers

• 📫 How to reach me 📫

  • My personal website

ranenextra's Issues

verifyBrier Error

Hi! In using the verifyBrier function I received the following error:

Error in verifyBrier(anenX, obs, threshold, ensemble.func, ### baseline = ozoneForecasts) :
all(ret.flt$check - ret.flt$bs < 1e-06) is not TRUE

The error occurred when using a forecast file that contained "NA" for the first 12 FLTs of data for the parameter of ozone concentrations, while parameters for temp, wind speed, and wind direction included data. When running the same script but with a slightly modified forecast file which contains zeroes for the first 12 FLTs of data for the parameter of ozone concentrations, instead of "NA", the function works as expected.

I can follow the script for verifyBrier, and in running the individual lines of verifyBrier, the verification::brier command completes in line 136. However, it seems as if the ret.flt$check values are all NA (for all monitors and all days).

I can email the input files and scripts that I am using, upon request. Thanks!

Mark

subsetCoordinates POI check

The subsetCoordinates function currently only checks to see if the input Points of Interest are a vector, whereas the prior version included a check to see if the POI were a data frame. The inclusion of a check for POI as a data frame is useful, as observations may be at varying locations other than in a grid format.

subset.coords <- subsetCoordinates(

• fore.coords$X, fore.coords$Y, obs.coords)
  Error in subsetCoordinates(fore.coords$X, fore.coords$Y, obs.coords) :
  The input poi can should be a named vector

is.vector(obs.coords)
[1] FALSE
is.data.frame(obs.coords)
[1] TRUE

Thanks! - Mark

subsetStations dealing with missing StationNames

demo-4_SEA2019 and plotAnalogSample problem

Hi. I am a PhD student and I am working with statistical tools for weather forecast purposes.

I am learning AnEn with your R codes, but in demo-4 something goes wrong.

fcst.id <- which(forecasts$ParameterNames == observations$ParameterNames[config$observation_id])
plotAnalogSample(AnEn, config, i.station = 1, i.test.day = 1,
                 obs.id = config$observation_id, fcst.id = fcst.id, 
                 parameter.names = forecasts$ParameterNames,
                 flts.fraction = 3600, legend.texts = c("Obs","Fcst"))

Error in if (length(parameter.names) != dim(forecasts)[1]) { :
the argument has zero length

I didn't change anything. This problem is also in Binder.

How can I fix it?

Thanks