Skip to contents

Calculate the negative log-likelihood of the expected profit model (EPM) with log-normal catch function. For more information on the EPM lognormal model see section 8.4.5 in the FishSET user manual. https://docs.google.com/document/d/1dzXsVt5iWcAQooDDXRJ3XyMoqnSmpZOqirU_f_PnQUM/edit#heading=h.ps7td88zo4ge

Usage

epm_lognormal(starts3, dat, otherdat, alts, project, expname, mod.name)

Arguments

starts3

Starting parameter values as a numeric vector. The order of parameters in the vector is:
c([catch-function params], [travel-dist params], [stdev], [common scale param]),
where the length of catch-function parameters is the # of alternatives * # of catch-function variables, length of travel-distance parameters is the # of travel-distance variables, length of standard deviation defaults to 1 but alternative- specific standard deviation values can be specified (length = # of alternatives), and the common scale parameter is a single value.

dat

Data matrix, see output from shift_sort_x, alternatives with distance.

otherdat

List that contains other data used in the model, see section 8.4.5 in the FishSET user manual for more details (link in the description above): (1) 'griddat': catch-function variables that interact with alternative-specific catch-function parameters and do not vary across alternatives (e.g., vessel gross tonnage). (2) 'intdat': travel-distance variables that interact with travel-distance parameters and the distance matrix and do not vary across alternatives. (3) 'prices': price in terms of $/landings units. This is typically a vector with prices for each observation, but can be a single value representing price for the entire dataset.

alts

Number of alternative choices in model

project

Name of project

expname

Expected catch table (optional)

mod.name

Name of model run for model result output table

Value

ld: negative log likelihood

Details

This function is called in discretefish_subroutine when running an EPM model with a log-normal catch function.