R/weibull.surface.R

In fzwaeustc/pcrfn: Peng's Collection of R Functions at Nielsen

weibull.surface <-
function(x, y, a1_intercept=0.025, a1_slope=0.0001,
                            b1_intercept=0.576, b1_slope=0.0021,
                            a2_intercept=0.011, a2_slope=0.0023,
                            b2_intercept=0.134, b2_slope=0.00005){
  (
  exp(1)^(-(a1_intercept+a1_slope*y)*(x-1)^(b1_intercept+b1_slope*y))+
    exp(1)^(-(a2_intercept+a2_slope*y)*(y+1-x)^(b2_intercept+b2_slope*y))-(
    exp(1)^(-(a2_intercept+a2_slope*y)*y^(b2_intercept+b2_slope*y))+
    (exp(1)^(-(a1_intercept+a1_slope*y)*y^(b1_intercept+b1_slope*y))-
    exp(1)^(-(a2_intercept+a2_slope*y)*y^(b2_intercept+b2_slope*y)))*(x-1)/y)
  )
}