R/calcul.xi.R

Defines functions calcul.xi

#' @import ggplot2
#' @import rstan
#' @import methods
#' @import stats
#' @import graphics
#' @import grDevices
#' @import stats4
calcul.xi <-
function(wm, meanbeta, a = NULL, model){
    if(model == "logistic"){
        x <- -(a +log(1/wm -1))/exp(meanbeta)
    }else{
        if(model == "power"| model == "power_log"){
            x <- wm^exp(meanbeta)
        }else{
            stop("Error: Please insert a valid model")	
        }
    }
    return(x)
}

Try the dfped package in your browser

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

dfped documentation built on May 2, 2019, 8:36 a.m.