Recent Enhancements to the MULTIFAN-CL Software

Similar documents
4.8 Improper Integrals

LAPLACE TRANSFORMS. 1. Basic transforms

Estimating the population parameter, r, q and K based on surplus production model. Wang, Chien-Hsiung

Bipartite Matching. Matching. Bipartite Matching. Maxflow Formulation

e t dt e t dt = lim e t dt T (1 e T ) = 1

Optimality of Myopic Policy for a Class of Monotone Affine Restless Multi-Armed Bandit

Positive and negative solutions of a boundary value problem for a

S Radio transmission and network access Exercise 1-2

A new model for limit order book dynamics

f t f a f x dx By Lin McMullin f x dx= f b f a. 2

graph of unit step function t

3. Renewal Limit Theorems

Minimum Squared Error

Minimum Squared Error

( ) ( ) ( ) ( ) ( ) ( y )

1.0 Electrical Systems

GEOMETRIC EFFECTS CONTRIBUTING TO ANTICIPATION OF THE BEVEL EDGE IN SPREADING RESISTANCE PROFILING

SOME USEFUL MATHEMATICS

Chapter 2: Evaluative Feedback

Flow Networks Alon Efrat Slides courtesy of Charles Leiserson with small changes by Carola Wenk. Flow networks. Flow networks CS 445

Randomized Perfect Bipartite Matching

CSC 373: Algorithm Design and Analysis Lecture 9

Sph3u Practice Unit Test: Kinematics (Solutions) LoRusso

Transformations. Ordered set of numbers: (1,2,3,4) Example: (x,y,z) coordinates of pt in space. Vectors

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.

u(t) Figure 1. Open loop control system

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

Solutions to assignment 3

Motion. Part 2: Constant Acceleration. Acceleration. October Lab Physics. Ms. Levine 1. Acceleration. Acceleration. Units for Acceleration.

A Kalman filtering simulation

A Point Estimate of Natural Mortality Coefficient of Southern Minke Whales Using JARPA data

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

Contraction Mapping Principle Approach to Differential Equations

USP. Surplus-Production Models

Introduction to Congestion Games

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

The solution is often represented as a vector: 2xI + 4X2 + 2X3 + 4X4 + 2X5 = 4 2xI + 4X2 + 3X3 + 3X4 + 3X5 = 4. 3xI + 6X2 + 6X3 + 3X4 + 6X5 = 6.

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER

0 for t < 0 1 for t > 0

Introduction to SLE Lecture Notes

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

MATH 124 AND 125 FINAL EXAM REVIEW PACKET (Revised spring 2008)

Physics 2A HW #3 Solutions

Physic 231 Lecture 4. Mi it ftd l t. Main points of today s lecture: Example: addition of velocities Trajectories of objects in 2 = =

A Numerical Model of Commodity Moisture Loss and Temperature Distribution During Refrigerated Storage

Thermal neutron self-shielding factor in foils: a universal curve

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems

EECE 301 Signals & Systems Prof. Mark Fowler

EXERCISE - 01 CHECK YOUR GRASP

Stat13 Homework 7. Suggested Solutions

A 1.3 m 2.5 m 2.8 m. x = m m = 8400 m. y = 4900 m 3200 m = 1700 m

CHAPTER 7: SECOND-ORDER CIRCUITS

An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples.

Mathematics 805 Final Examination Answers

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

RESPONSE UNDER A GENERAL PERIODIC FORCE. When the external force F(t) is periodic with periodτ = 2π

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

Properties of Logarithms. Solving Exponential and Logarithmic Equations. Properties of Logarithms. Properties of Logarithms. ( x)

Price Discrimination

September 20 Homework Solutions

Algorithmic Discrete Mathematics 6. Exercise Sheet

Temperature Rise of the Earth

Laplace Examples, Inverse, Rational Form

ANSWERS TO EVEN NUMBERED EXERCISES IN CHAPTER 2

can be viewed as a generalized product, and one for which the product of f and g. That is, does

ARTIFICIAL INTELLIGENCE. Markov decision processes

Curvature. Institute of Lifelong Learning, University of Delhi pg. 1

Graduate Algorithms CS F-18 Flow Networks

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

Maximum Flow. Flow Graph

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

REAL ANALYSIS I HOMEWORK 3. Chapter 1

Let. x y. denote a bivariate time series with zero mean.

Assignment 16. Malaria does not affect the red cell count in the lizards.

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Investment-specific Technology Shocks, Neutral Technology Shocks and the Dunlop-Tarshis Observation: Theory and Evidence

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t)

What is maximum Likelihood? History Features of ML method Tools used Advantages Disadvantages Evolutionary models

2D Motion WS. A horizontally launched projectile s initial vertical velocity is zero. Solve the following problems with this information.

How to Prove the Riemann Hypothesis Author: Fayez Fok Al Adeh.

SLOW INCREASING FUNCTIONS AND THEIR APPLICATIONS TO SOME PROBLEMS IN NUMBER THEORY

Average & instantaneous velocity and acceleration Motion with constant acceleration

FM Applications of Integration 1.Centroid of Area

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Probability, Estimators, and Stationarity

Abdul-Aziz Yakubu. Department of Mathematics Howard University Washington, D.C

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER

Chapter 2. First Order Scalar Equations

PHYSICS 1210 Exam 1 University of Wyoming 14 February points

Solutions for Nonlinear Partial Differential Equations By Tan-Cot Method

Time series Decomposition method

Math 2214 Solution Test 1 B Spring 2016

Physics 240: Worksheet 16 Name

Math 10B: Mock Mid II. April 13, 2016

The order of reaction is defined as the number of atoms or molecules whose concentration change during the chemical reaction.

1. Introduction. 1 b b

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Performance Comparison of LCMV-based Space-time 2D Array and Ambiguity Problem

Transcription:

SCTB15 Working Pper MWG-2 Recen Enhncemen o he MULTIFAN-CL Sofwre John Hmpon 1 nd Dvid Fournier 2 1 Ocenic Fiherie Progrmme Secreri of he Pcific Communiy Noume, New Cledoni 2 Oer Reerch Ld. PO Box 2040 Sidney, BC, Cnd V81 3

1 Inroducion The lengh-bed, ge-rucured ock emen known MULTIFAN-CL i currenly ued for un nd oher pelgic ock emen by he SPC Ocenic Fiherie Progrmme (kipjck, yellowfin, bigeye nd Souh Pcific lbcore un in he weern nd cenr l Pcific), he Nionl Mrine Fiherie Service (NMFS) Honolulu Lborory (blue hrk, blue mrlin nd wordfih), he NMFS L Joll Lborory (Norh Pcific lbcore) nd he Jpn Nionl Reerch Iniue of Fr Se Fiherie (wh). Technicl documenion of mo of he model feure i provided in Hmpon nd Fournier (2001). Some ddiionl informion on recen developmen o he ofwre were provided he MWG eion SCTB14 (Hmpon nd Fournier 2001b). Thi pper briefly ouline developmen o he ofwre ince SCTB14. 2 MULTIFAN-CL developmen ince SCTB14 2.1 Reference poin nlyi A reference poin nlyi h been dded o he yield nlyi ecion of he model. The reference poin h hve been choen re he ol equilibrium biom ( B ), he dul dul equilibrium biom ( B ) nd he ggrege fihing morliy ( F ). To ummrie, he erie of ep involved i follow: 1. Eime populion model prmeer, including he prmeer of Beveron nd Hol ockrecruimen relionhip (SRR). Aume h he eimed prmeer of he SRR (α nd β) cn be ued o decribe he relionhip beween equilibrium recruimen (R ) nd equilibrium pwning biom ( B ): αb R = (1) β + B 2. Compue be ge-pecific fihing morliy vecor, F, vriou muliple of which re umed o be minined ino he fuure; normlly, he be i compued he verge over ome recen period of ime. 3. For vriou muliplier (x) of F compue he equilibrium ol biom ( B ) nd equilibrium pwning biom. Le = R( x) (2) where 1 exp ( M + xf ) = wm, (3) = 1 ' = 1 M i he nurl morliy re ge, w i men weigh ge nd m i he proporion pwning ge. Wih ll oher prmeer fixed, vlue of φ (x) cn be deermined for ny vlue of x. From ubiuion of (2) ino (1), we cn hen obin = α β (4) which pecifie he equilibrium pwning biom ocied wih ny fihing morliy muliplier x. R ( x ) my hen be deermined from rerrngemen of erm in (2). Similrly, he ol 2

ol equilibrium biom ocied wih he fihing morliy muliplier x my hen be deermined by: = R ( x) (5) where 1 exp ( M + xf ) = w. (6) = 1 ' = 1 4. Compue he equilibrium yield, Y ( x ) funcion of he fihing morliy muliplier x: Y ( x) = R ( x) = 1 xf xf + M 1 xf M ( 1 e ) exp ( M + xf ) ' = 1 w (7) Le z be he vlue of x h mximize Y ( x ). The i hen given by Y ( z). 5. The reference poin of inere re: B dul B ( z ) B F = ol = = B ( z ) Y ( z ) B ( z) 6. Compre he cul eimed biom nd fihing morliy level ime wih hee reference ol ol dul dul poin. Thi i done by compuing he rio B B, B B, F F nd heir 95% confidence inervl nd compring hem wih 1.0. Vlue of F F ignificnly greer hn 1.0 would indice overfihing, while vlue of would indice n overfihed e. 2.2 Incorporion of weigh-frequency d B B nd/or ol ol B B of le hn 1.0 A SCTB14, weigh-frequency d from he longline fihery unloding in Gum were preened by Kikkw nd Cuhing (2001). Thee d ppered o be exremely informive wih repec o ge rucure nd growh. A erch of OFP dbe nd conulion wih vriou nionl gencie reveled h here re ubnil weigh-frequency d vilble for longline fiherie for yellowfin nd bigeye un. Thee d ppered o be imilrly informive he Gum d nd in mny ce he mple ize repreened very lrge proporion of he ol cch (pproching 100% in he ce of he Aurlin e-co longline fihery). I w herefore decided o incorpore weigh-frequency d ino MULTIFAN -CL o ke dvnge of he vilbiliy of hee d. Thi required he following chnge/ddiion o he ofwre: 1. The inpu d file ( frq file) rucure w chnged o llow he pecificion of he rucure of he weigh-frequency d nd o include weigh-frequency d record, if vilble, for ech fihery record. 2. The MULTIFAN-CL model i ge-rucured he compuionl level bu he probbiliy diribuion of lengh for ech ge cl i lo compued o llow he model o be fi o he lengh-frequency d. We oped o pecify he probbiliy diribuion of weigh--ge deerminiic funcion of he lengh--ge. The ddiion of weigh-frequency d o he model herefore did no involve he eimion of ddiionl prmeer. The inegrion of he weigh dul dul 3

diribuion w done uing he me inervl pecified in he weigh-frequency d. A weighlengh relionhip ( W = L b ) w ued o compue he weigh diribuio n from he lengh diribuion. The weigh-lengh prmeer nd b mu be pecified. 3. The likelihood funcion for he weigh d i idenicl o h for he lengh d. However, we dded ome flexibiliy in pecifying he effecive mple ize for boh he lengh nd weigh d. Previouly, he effecive mple ize (which influence he vrince erm in he likelihood funcion) w he cul mple ize up o mximum of 1,000 divided by 10. Thi w done o cknowledge he likelihood h rel ize mple re unlikely o be ruly rndom umed by he iicl model. However, ome of he weigh mple repreen very high proporion of he ol cch, reuling in n increed likelihood of he mple being rndom (or repreenive). In hee ce, i would be pproprie for he cul mple ize o more ccurely reflec he effecive mple ize. We hve herefore dded more flexibiliy o boh he weigh nd lengh likelihood funcion by llowing he relionhip beween cul nd effecive mple ize o be pecified on fihery by fihery bi. 4. The Jv-bed MULTIFAN-CL viewer w upded o genere oberved nd prediced weigh frequencie. A cpbiliy o red he d file (o diply he oberved lengh nd weigh d only) w lo dded o he viewer. 2.3 Cobb-Dougl cchbiliy model A pr of bioeconomic modelling projec, we were ked o e he hypohee h: () cchbiliy i independen of he level of effor, gin he lernive h he level of effor impc cchbiliy hrough, for exmple, informion hring (poiive impc) or crowding (negive impc) (b) cchbiliy i independen of ock ize, gin he lernive h ock ize impc cchbiliy (he noion of n underlyin g populion h feed urfce chool which would enble pure eine CPUE o remin relively conn he ol populion declined) Thi w done follow: () ln( F ) = ln( ) + ln( q) + α ln( e), where F i fihing morliy, i eleciviy, q i cch biliy, e i fihing effor nd α i he prmeer. Thi effec cn be grouped cro fiherie, in which ce n index of relive grouped effor i ued. Eime of α ignificnly differen o 1 would provide evidence for he lernive hypohei over he null hypohei. (b) ln( F ) = ln( ) + ln( q) + ln( e) + β ln( B), where B i he biom index nd β i he prmeer. Thi effec cn lo be grouped cro fiherie if deired. Eime of β ignificnly differen o 0 would provide evidence for he lernive hypohei over he null hypohei. 2.4 Seonl nd ime -erie movemen Currenly, movemen coefficien re conn over ime. We re currenly ( he ime of wriing) incorporing rucure o llow eonl vriion in movemen for ech e of coefficien. There i ome indicion in he d for Souh Pcific lbcore h uch eonliy could be imporn. Vriou nlye hve uggeed h lrge-cle ENSO-ype ocenogrphic vriion influence he diribuion nd movemen of un. We herefore pln o inroduce rucure ino he model o llow movemen coefficien o vry over ime ccording o n environmenl correle. 4

3 Reference Hmpon, J., nd Fournier, D.A. 2001. A pilly-diggreged, lengh-bed, ge-rucured populion model of yellowfin un (Thunnu lbcre) in he weern nd cenrl Pcific Ocen. Mr. Frehw. Re. 52:937 963. Hmpon, J., nd Fournier, D.A. 2001b. Recen enhncemen o he MULTIFAN-CL ofwre. Working Pper MWG-1, SCTB 14, 9 16 Augu 2001, Noume, New Cledoni. Kikkw, B., nd Cuhing, J. 2001. Vriion in growh nd morliy of bigeye un (Thunnu obeu) in he equoril weern Pcific Ocen. Working Pper BET-2, SCTB 14, 9 16 Augu 2001, Noume, New Cledoni. 5