R/mng.thin.prob.R

Defines functions thin.prob

Documented in thin.prob

thin.prob <- function( region
                      , skidding.distance.100m
                      , AgeTo5
                      , vuprha.m3.ha
                      , slope.per
                      , SI.m
                      , SI.spp
                      
                       ){

  ## Intercept, V, ....
  AgeTo5[AgeTo5 > -1] <- -1
  pars <- c(-6.3897, 0.5699, -0.038, -0.0693, 0.8227, 1.1354, -21.71,
            -0.0955, -0.0012)
  return(1 / (1 + exp(- (
      pars[1]  + 
      pars[2]  * log(vuprha.m3.ha) +
      pars[3]  * slope.per +
      pars[4]  * skidding.distance.100m +
      pars[5]  * (SI.spp == 2) +
      pars[6]  * (region != "Vest") +
      pars[7]  * 1/SI.m +
      pars[8]  * AgeTo5 +
      pars[9]  * AgeTo5^2
    )))

    )
}

Try the sitree package in your browser

Any scripts or data that you put into this service are public.

sitree documentation built on Oct. 25, 2018, 5:04 p.m.