regi0008
commented
5 years ago Based on example of 2018_ECOMS-UDG example in climate4R notebook, under their Verifications (Correlation Analysis), there are errors:
Correlation Analysis
object cor contains correlation coefficients for each point.
need to convert list into array
summary(NCEP_Iberia_tas) summary(CFS_Iberia_tas)
obs <- NCEP_Iberia_tas[["Data"]]
obs_array <- array(as.numeric(unlist(obs)))
dim(obs_array) = 86640
fcst <- CFS_Iberia_tas[["Data"]]
fcst_array <- array(as.numeric(unlist(fcst)))
dim(fst_array) = 2599200
cor <- veriApply(verifun = "EnsCorr", fcst = fcst, obs = obs, tdim = length(dim(fcst))-1, ensdim = length(dim(fcst)))
above error:
Error in veriApply(verifun = "EnsCorr", fcst = fcst, obs = obs, tdim = length(dim(fcst)) - :
odims[-otdim] == dim(fcst)[-c(ensdim, tdim)] are not all TRUE
easyVeri2grid is a bridging function that performs transformation
from the raw verification matrices produced by library(easyVerification)
to the grid structure with all required metadata for data collection and representation by climate4R bundle
cor.grid <- easyVeri2grid(easyVeri.mat = cor, obs.grid = NCEP_Iberia_tas[["Data"]], verifun = "EnsCorr")
above error:
$ operator is invalid for atomic vectors
plot correlation map
spatialPlot(climatology(cor.grid), backdrop.theme = "countries", at = seq(-1,1,0.05), colorkey = list(space = "bottom"), main = "Ensemble correlation", sub = "CFSv2 24 member - DJF Mean Temperature (1982-2010")
above error:
Error in getDim(grid) : object 'cor.grid' not found
--------------------------------------
Notes: 1) Am still figuring out if there is anything wrong with the veriApply() function used, and whether the inputs "obs" and "fcst" is correct (as an array/vector)
rykang11
To do by Tuesday, 20th August 2019 (everyone)