Singular results of biascorrection()

[Fisico2k]

I am trying to use biascorrection() with isimip3 algorithm. Here the code:

pr.proj.BC <- biasCorrection(pr.arcis.ORIG.hist, pr.hist, pr.proj, precipitation=TRUE, wet.threshold=1,method="isimip3")

Now: pr.arcis.ORIG.hist is a grid 2x2km for period 1971-2005. It is an interpolation of station data. Variable: daily cumulative precipitation pr.hist is a grid 12x12km for period 1971-2005. It is the output of a CORDEX model in historical period. Variable: daily cumulative precipitation pr.proj is a grid 12x12km for period 2006-2015. It is the output of a CORDEX model in forecast period. Variable: daily cumulative precipitation.

The issue I realize is that in some grid points the resulting grid pr.proj.BC have negative values for pr variable. How is it possible? Is there an error in the R script. Is the mistake mine?

[Alireza-Amani]

@Fisico2k , Hey there,

I think it would be very helpful if you could provide a minimal re-producible code that also includes small part of your data. Im not sure what is the proper way to include data in a github issue, maybe you can upload it somewhere?

Note that im just a user of this nice library, not a developer. But It d be great if I can be of any use for your problem.