I am currently programming a cokriging between PM10 and PM2.5, but we would like to add external drifts to this classic cokriging.
Is there any possibility to add different external drift to the different variables (meaning PM10 external drift to PM10 observations and PM2.5 external drift to PM2.5 observations)?
Here is the code I use:
- Code: Select all
tabRValC=data.frame(x1=tabvario$x,x2=tabvario$y,z1=tabvario$Mesure,z2=tabvario$Mesure_Co,f1=tabvario$Modele,f2=tabvario$Modele_Co)
dataVC.db = db.create(tabRValC,flag.grid=F,ndim=2,nvar=2,autoname=F)
## Variables polluants
dataVC.db=db.locate(dataVC.db,"z1","z",1)
dataVC.db=db.locate(dataVC.db,"z2","z",2)
## Derives
dataVC.db=db.locate(dataVC.db,"f1","f",1)
dataVC.db=db.locate(dataVC.db,"f2","f",2)
I have no error issues but I am not sure of what the function is doing with the two variables (z1 and z2) and the two external drifts (f1 and f2).
Thank you for your help!
Charline