R/mode.lands.t.R

mode.lands.t <-
function (nu, zeta) 
{
    if (nu < 2) 
        stop("'nu' must be greater than or equal to 2")
    fcn <- function(theta, nu, zeta) {
        -(nu - 1) * tan(theta) + (1 + nu) * zeta * cos(theta)
    }
    theta.max <- uniroot(fcn, interval = c(-pi/2, pi/2), nu = nu, 
        zeta = zeta)$root
    sqrt(nu) * tan(theta.max)
}

Try the EnvStats package in your browser

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

EnvStats documentation built on Aug. 22, 2023, 5:09 p.m.