R/gencauchy.variog.R

In ProbForecastGOP: Probabilistic weather forecast using the GOP method

"gencauchy.variog" <- function(v,variog,d,w){
   num <- variog - (v[1]+v[2]*(1-(1+(d/v[3])^v[4])^(-v[5]/v[4])))
   den <- (v[1]+v[2]*(1-(1+(d/v[3])^v[4])^(-v[5]/v[4])))
   num.den <- num/den
   num.den2 <- (num.den)^2
   Q <- sum(w*num.den2)
   return(Q)
}