Problems with "'x' and 'y' must be given in ascending order" after using subsetGrid

[rbalmaceda]

Hi,
Im using subsetGrid function to select some stations. Then I want to calculate bias, or correlations using valueMeasure from package climate4R.value but I get this error:

[2020-11-26 01:32:58] Computing member 1 out of 11
Error in mat2Dto3Darray(mat, xy$x, xy$y) :
Coordinates 'x' and 'y' must be given in ascending order

Im not very clear what the problem is,maybe you can guide me. The error only appears when I use subsetGrid first.
Thanks a lot :)

[jorgebanomedina]

Hi Rocio, Could you please share to us some data so we can reproduce the error? Also the chunk of code you are running please. Thank you! Jorge El jue, 26 nov 2020 a las 1:55, Rocio B (<notifications@github.com>) escribió:

…

Hi, Im using subsetGrid function to select some stations. Then I want to calculate bias, or correlations using valueMeasure from package climate4R.value but I get this error: [2020-11-26 01:32:58] Computing member 1 out of 11 Error in mat2Dto3Darray(mat, xy$x, xy$y) : Coordinates 'x' and 'y' must be given in ascending order Im not very clear what the problem is,maybe you can guide me. The error only appears when I use subsetGrid first. Thanks a lot :) — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#35>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AE4DRYJEQNV2EIOSUFDESU3SRWROPANCNFSM4UDCILZA> .

[rbalmaceda]

Thanks Jorge for your quick response. I realized that I was a bit confused about the error. I think I have two different problems. I attach some data with 2 objects that are BindGrids after using downscale.CV: 1 analog and 1 GLM
In the case of analogs, I already have them in a specific region, and the error come up after ScaleGrid .
In the case of the GLM that I have them in a larger domain, I first calculate the anomalies but then I cannot select the stations. I was able to solve it by changing the attributes of the data matrix.
I hope it is understood, and thank you very much!

This is my code:
`#########analogs in Region 1
ANOM_ERA_analogos.cv_R1=scaleGrid(ERA_analogos.cv_Tn_R1_tt,time.frame = "daily")
ANOM_cor_ERA.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE)

############################################ GLM
estaciones_grupo1=c( "10178", "10332", "10221", "10156", "10204", "10190", "10159", "10166",
"10334", "10134", "10453", "10170","10148", "10210", "10369", "10113",
"10456", "10206", "10177", "10133", "10154", "10192", "10451", "10311",
"10490", "10446", "10400", "10216", "10900", "10358", "10902", "10421",
"10502", "10907", "10114", "10145", "10492", "10912")
tt=attributes(ERA_GLM_cv_tt$Data)
anom_modelos_ERA_GLM.cv=scaleGrid(grid=ERA_GLM_cv_tt,base=ERA_GLM_cv_tt,time.frame = "daily")
attributes(anom_modelos_ERA_GLM.cv$Data)=tt
ANOM_ERA_GLM.cv_R1=subsetGrid(anom_modelos_ERA_GLM.cv,station.id = estaciones_grupo1)
ANOM_cor_ERA_GLM.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE)#, measure.cor(deseason = 1, method = pearson))
`

example.zip

[jorgebanomedina]

Hi Rocio, Maybe I am missing the error or not understanding it correctly. If I run the following chunk of code I do not find any errors (note that I am commenting the 'attributes' lines to show that you do not need to touch them manually). #########analogs in Region 1 ANOM_ERA_analogos.cv_R1=scaleGrid(ERA_analogos.cv_Tn_R1_tt,time.frame = "daily") ANOM_cor_ERA.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE) ############################################ GLM estaciones_grupo1=c( "10178", "10332", "10221", "10156", "10204", "10190", "10159", "10166", "10334", "10134", "10453", "10170","10148", "10210", "10369", "10113", "10456", "10206", "10177", "10133", "10154", "10192", "10451", "10311", "10490", "10446", "10400", "10216", "10900", "10358", "10902", "10421", "10502", "10907", "10114", "10145", "10492", "10912") # tt=attributes(ERA_GLM_cv_tt$Data) anom_modelos_ERA_GLM.cv=scaleGrid(grid=ERA_GLM_cv_tt,base=ERA_GLM_cv_tt,time.frame = "daily") # attributes(anom_modelos_ERA_GLM.cv$Data)=tt ANOM_ERA_GLM.cv_R1=subsetGrid(anom_modelos_ERA_GLM.cv,station.id = estaciones_grupo1) ANOM_cor_ERA_GLM.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE)#, measure.cor(deseason = 1, method = pearson)) Am I missing something? Jorge El vie, 27 nov 2020 a las 23:56, Rocio B (<notifications@github.com>) escribió:

…

Thanks Jorge for your quick response. I realized that I was a bit confused about the error. I think I have two different problems. I attach some data with 2 objects that are BindGrids after using downscale.CV: 1 analog and 1 GLM In the case of analogs, I already have them in a specific region, and the error come up after ScaleGrid . In the case of the GLM that I have them in a larger domain, I first calculate the anomalies but then I cannot select the stations. I was able to solve it by changing the attributes of the data matrix. I hope it is understood, and thank you very much! This is my code: `#########analogs in Region 1 ANOM_ERA_analogos.cv_R1=scaleGrid(ERA_analogos.cv_Tn_R1_tt,time.frame = "daily") ANOM_cor_ERA.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_analogos.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE) ############################################ GLM estaciones_grupo1=c( "10178", "10332", "10221", "10156", "10204", "10190", "10159", "10166", "10334", "10134", "10453", "10170","10148", "10210", "10369", "10113", "10456", "10206", "10177", "10133", "10154", "10192", "10451", "10311", "10490", "10446", "10400", "10216", "10900", "10358", "10902", "10421", "10502", "10907", "10114", "10145", "10492", "10912") tt=attributes(ERA_GLM_cv_tt$Data) anom_modelos_ERA_GLM.cv=scaleGrid(grid=ERA_GLM_cv_tt,base=ERA_GLM_cv_tt,time.frame = "daily") attributes(anom_modelos_ERA_GLM.cv$Data)=tt ANOM_ERA_GLM.cv_R1=subsetGrid(anom_modelos_ERA_GLM.cv,station.id = estaciones_grupo1) ANOM_cor_ERA_GLM.cv_R1=valueMeasure(subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices =1),subsetDimension(ANOM_ERA_GLM.cv_R1,dimension = "member", indices = c(2:3)), measure.code = "ts.rp",return.NApercentage=FALSE)#, measure.cor(deseason = 1, method = pearson)) ` example.zip <https://github.com/SantanderMetGroup/transformeR/files/5609762/example.zip> — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#35 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AE4DRYNVHY34P2GG2P5UZNTSSAVCXANCNFSM4UDCILZA> .

[rbalmaceda]

Thank u Jorge, I was able to solve the problem by reinstalling transformeR to version 2.0.1 after heeding this warning:

WARNING: Your current version of transformeR (v2.0.2) is ahead of the master branch version (2.0.1)
Development version may have an unexpected behaviour