Mark-recapture models for closed populations
|
|
- Sabina Morgan
- 5 years ago
- Views:
Transcription
1 Mark-recapture models for closed populations A standard technique for estimating the size of a wildlife population uses multiple sampling occasions. The samples by design are spaced close enough in time so that the change in population size during the study is negligible. A population with no changes due to births, deaths, or migration is termed a closed population. During each sample occasion, the individuals in the population face some chance of being caught in the sample. Those animals caught are marked with an individually identifiable mark and released back into the population. The identification of individuals in the samples allows modeling of capture histories which leads to improved estimates of population size. In a typical small mammal study, live traps are baited and laid out in the study area in a large grid. The traps are visited each morning and the animals caught are recorded, marked if necessary, and released. Such a study might last several days or weeks and can provide a reliable estimate of population size during that time if care is taken during design, conduct, and analysis of the study. The data for, say, a study with $ sampling occasions is summarized in the form of frequency counts for capture histories. An animal's capture history can be represented as a string of!'s and "'s, for instance, "!" signifies that a particular animal was captured on the first day, not captured on the second day, and captured on the third day. Each animal in the population is thereby categorized into
2 one of ) possible capture histories, with a resulting frequency count of animals displaying each capture history: Capture day frequency " # $ count probability Capture " " " ] " 1 " history " "! ] # 1 # "! " ] $ 1 $ "!! ] % 1 %! " " ] & 1 &! "! ] ' 1 '!! " ] ( 1 (!!! ] ) 1 ) The sampling scheme induces a multinomial distribution for the frequency counts: ], ],..., ] µ multinomial( >, 1, 1,..., 1 ). " # ) " # ) Here > is the total population size, and 1 4 is the probability that an animal will have capture history 4. The last count, ] ), is the number of animals that did not appear in the sample and is unobserved. The capture history probabilities depend on the probabilities of events which cause animals to end up or not in the samples. The various closed population mark-
3 recapture models represent different constructions for the 1 4 's. For instance, suppose every animal has an equal probability of being captured on any given day, but also that the common capture probability varies in value from day to day. Let ) " œ capture probability on the first day, )# œ capture probability on the second day, ) $ œ capture probability on the third day. The different possible capture histories have the following probabilities: 1 " œ ))) " # $ 1 # œ )")# ") $ 1 $ œ )" ") # ) $ 1 % œ )" " )# ") $ 1 & œ ") " ) #) $ 1 ' œ " )" ) # ") $ 1 ( œ " )" ") # ) $ 1 œ " ) " ) ") ) " # $ Parameter estimation In general a closed population capture recapture study will have < capture occasions, producing 5 œ # < possible capture histories and corresponding categories in the multinomial likelihood. The category probabilities 1", 1#,..., 15 are modeled as functions 1 " œ1" ), 1 # œ1# ),..., 15 œ15 ) of underlying parameters ) œ ) ", )#,..., ) 6, representing different structures for capture probabilities. The essential structural problem in parameter estimation is that the number > of observations in the
4 multinomial distribution (and along with it, the last count ] 5 ) is unknown. The data are the observed counts C", C#,..., C5". The likelihood is a multinomial probability: > C C C C, C,..., C " # 5 " # 5 P ), > œ 1 1 â1. Here " # 5 C5 œ > C" C# âc5", 1 œ"1 1 â1, 5" " # 5" > >x C C C C xc xâc x,,..., œ. " # 5 " # 5 The integer ML trick can be used to find the relationship between the integer ML estimate of > and the other model parameters. Set P ), > > P ), >" "œ >815 "œ!, where 8 œ C" C# âc5" (the number of unique animals sampled), to find that s 8 > œ floor. "s 1 5
5 The computational task is to find the ML estimates of the other parameters. One can substitute the expression for >s into the multinoial log-likelihood and numerically maximize the log-likelihood for the remaining parameters. ML estimates for the two-sample mark recapture model with time-dependent capture probabilities )" and )# can be found symbolically. The model has four possible capture histories ((1 1), (1 0), (0 1), (0 0)) and three parameters: ) ", ) #, and >. Thus, the model is a saturated multinomial model, with 1 " œ )")#, 1 œ ) "), # " # 1 œ ") ), $ " # 1 œ " ) "). % " # ML estimates of )", )#, and > are found by equating 1", 1#, and with their saturated multinomial ML estimates: 1 $ C" > œ )) " #, C# > œ )" ")#, C$ > œ " " # ) ).
6 Solving the equations algebraically gives the ML estimates: )s œ " )s œ # C CC " $ ", C CC " # ", s CC " # CC " $ > œ C ". The estimate for > is recognized as the Lincoln-Peterson estimate. The integer ML estimate for > turns out to be floor C C C C ÎC. " # " $ " Wald intervals The large sample results for ML estimates in ordinary reduced-parameter multinomial models must be modified for these mark-recapture models, because the number of trials > is itself an unknown parameter. The relevant results were proved and reported by Sanathanan (1972 Annals of Mathematical Statistics 43: ). According to Sanathanan's results, the Fisher information matrix for the scaled ML estimates given by s) ), s) ),..., s) ), s> > Î > " " # # 6 6 is in the following form:
7 M ), > œ E + w " Here E is the information matrix per observation for the parameters )", )#,..., ) 6 when > is known, with the element in the 3th row and 7th column given by 5 # ` log P ) E œ " `14 ) `14 ). `) `) 1 ) `) `) 4œ" Also, + is an 6 " column vector in which the 3th element is + œ 3 " `1 ) 1 ) `) k. 5 3 Finally, + 6" is the scalar quantity given by + œ 6" "15 ) 1 ). 5 As >Ä, the joint distribution of the scaled ML parameter estimates approaches a multivariate normal distribution with mean vector! and variance-covariance " matrix Z ), > œ M ), > Þ Asymptotic confidence intervals are formed using as standard errors the square roots of the diagonal elements of the matrix given by s>z s), s>.
8 It is straightforward to work out the Fisher information matrix for the two-sample mark-recapture model with timedependent capture probabilities ), ) : " # " " )! " " )" ") " " " M > œ! ), )# " )# ") #. " " " " )" ") # " ) " ) " ) ") " # " # The asymptotics for mark-recapture models are especially slow, and the resulting coverage rates for the Wald intervals tend to be inaccurate. An essential problem is that the speed of convergence depends on the number 8 of unique animals appearing in the samples, which is usually low in comparison to > in wildlife studies. Profile likelihood A more informative approach is to construct a profile likelihood for >. Using the multinomial log-likelihood function, one fixes the value of > and numerically maximizes the function with respect to the other parameters. One repeats the mazimization for a succession of integer values of >, starting say at its minimum possible value of 8 and proceeding upward through the integers. The resulting graph of the maximized log-likelihood at different values of > is the profile log-likelihood. The range of values of > for which the profile log-likelihood is within 1.92 ( œ " # #; Þ!&, 1 df) of its maximum value is a valid asymptotic 95% confidence interval for >. Simulations
9 have suggested that the profile likelihood confidence intervals, as well as parametric bootstrap confidence intervals, for mark-recapture estimates of population size have much better coverage properties than symmetric Wald intervals (Evans et al Journal of Agricultural, Biological and Environmental Statistics). Models for 1", 1#,..., 15 that are in wide use include: (1) time varying capture probabilities (as in the above example), (2) capture probabilities that vary behaviorally, such as enhanced (trap-happy) or diminished (trap-shy) chance of capture after being captured, (3) capture probabilities that are heterogeneous among individuals, (4) various combinations of the above. Mark recapture estimation for closed populations was unified in a wellknown paper by Otis el al. (1978 Wildlife Monograph); those authors produced and freely distributed a computer program known as CAPTURE that made mark recapture models for closed populations widely accessible in wildlife management. # R program to calculate profile log-likelihood for population size # using data from a two-sample capture recapture survey. # Counts are entered into the vector yy here: capture histories (1 1), # (1 0), (0 1), in that order. yy=c(7,80,7); # Data in example are from: Skalski et al Ecology 64: # Initial patameter values are calculated from multinomial ml estimates. theta10=yy[1]/(yy[1]+yy[3]); theta20=yy[1]/(yy[1]+yy[2]); psi10=log(theta10/(1-theta10)); psi20=log(theta20/(1-theta20));
10 nn=sum(yy); tt0=nn/(1-(1-theta10)*(1-theta20)); # Range of values for profile plot. Adjust "tlo" and "thi" to get a # good plot. tlo=110; thi=370; t.vals=tlo:thi; loglike.vals=t.vals; # ML objective function "negloglike.ml" is negative of log-likelihood; # the optimization routine in R, "optim", is a minimization # routine. The three function arguments are: psi = vector of # parameters (transformed to real line), ts = fixed population size for # profile calculation, ys = vector of frequencies. negloglike.ml=function(psi,ts,ys) { theta1=1/(1+exp(-psi[1])); # Constrains 0 < theta1 < 1. theta2=1/(1+exp(-psi[2])); # Constrains 0 < theta2 < 1. k=length(ys)+1; p=rep(0,k); p[1]=theta1*theta2; p[2]=theta1*(1-theta2); p[3]=(1-theta1)*theta2; p[4]=1-(p[1]+p[2]+p[3]); nn=sum(ys); ys=c(ys,ts-nn); ofn=-(lfactorial(ts)-sum(lfactorial(ys))+sum(ys*log(p))); return(ofn); } for (tt in tlo:thi) { ii=tt-tlo+1; MULTML=optim(par=c(psi10,psi20), negloglike.ml,null,method="nelder-mead",ts=tt,ys=yy); loglike.vals[ii]=-multml$val; } ci.height=max(loglike.vals)-3.843/2; # = 95th chisquare(1) percentile ci.cut=rep(ci.height,length(t.vals)); plot(t.vals,loglike.vals,type="l",xlab="population size",
11 ylab="profile log-likelihood"); points(t.vals,ci.cut,type="l",lty=2); # Calculate Wald confidence interval. that=tt0; theta1=theta10; theta2=theta20; I.mat=matrix(0,3,3); I.mat[1,1]=1/(theta1*(1-theta1)); I.mat[2,2]=1/(theta2*(1-theta2)); I.mat[3,3]=(1-(1-theta1)*(1-theta2))/((1-theta1)*(1-theta2)); I.mat[1,3]=1/(1-theta1); I.mat[3,1]=I.mat[1,3]; I.mat[2,3]=1/(1-theta2); I.mat[3,2]=I.mat[2,3]; Vhat=solve(I.mat); Vhat; me=1.96*sqrt(that*vhat[3,3]); c((that-me),that,(that+me));
Maximum Likelihood Estimation
Maximum Likelihood Estimation The likelihood and log-likelihood functions are the basis for deriving estimators for parameters, given data. While the shapes of these two functions are different, they have
More informationIntroduction to Population Modeling
Introduction to Population Modeling In addition to estimating the size of a population, it is often beneficial to estimate how the population size changes over time. Ecologists often uses models to create
More informationMeasuring Financial Risk using Extreme Value Theory: evidence from Pakistan
Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan Dr. Abdul Qayyum and Faisal Nawaz Abstract The purpose of the paper is to show some methods of extreme value theory through analysis
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationMaximum Likelihood Estimation
Maximum Likelihood Estimation EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #6 EPSY 905: Maximum Likelihood In This Lecture The basics of maximum likelihood estimation Ø The engine that
More informationEcon 8602, Fall 2017 Homework 2
Econ 8602, Fall 2017 Homework 2 Due Tues Oct 3. Question 1 Consider the following model of entry. There are two firms. There are two entry scenarios in each period. With probability only one firm is able
More informationThe method of Maximum Likelihood.
Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed
More informationLecture 3: Factor models in modern portfolio choice
Lecture 3: Factor models in modern portfolio choice Prof. Massimo Guidolin Portfolio Management Spring 2016 Overview The inputs of portfolio problems Using the single index model Multi-index models Portfolio
More informationData that can be any numerical value are called continuous. These are usually things that are measured, such as height, length, time, speed, etc.
Chapter 8 Measures of Center Data that can be any numerical value are called continuous. These are usually things that are measured, such as height, length, time, speed, etc. Data that can only be integer
More informationThe normal distribution is a theoretical model derived mathematically and not empirically.
Sociology 541 The Normal Distribution Probability and An Introduction to Inferential Statistics Normal Approximation The normal distribution is a theoretical model derived mathematically and not empirically.
More informationGamma Distribution Fitting
Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics
More informationChapter 8 Estimation
Chapter 8 Estimation There are two important forms of statistical inference: estimation (Confidence Intervals) Hypothesis Testing Statistical Inference drawing conclusions about populations based on samples
More informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationLog-linear Modeling Under Generalized Inverse Sampling Scheme
Log-linear Modeling Under Generalized Inverse Sampling Scheme Soumi Lahiri (1) and Sunil Dhar (2) (1) Department of Mathematical Sciences New Jersey Institute of Technology University Heights, Newark,
More informationDepartment of Agricultural Economics. PhD Qualifier Examination. August 2010
Department of Agricultural Economics PhD Qualifier Examination August 200 Instructions: The exam consists of six questions. You must answer all questions. If you need an assumption to complete a question,
More informationAgricultural and Applied Economics 637 Applied Econometrics II
Agricultural and Applied Economics 637 Applied Econometrics II Assignment I Using Search Algorithms to Determine Optimal Parameter Values in Nonlinear Regression Models (Due: February 3, 2015) (Note: Make
More informationSTRESS-STRENGTH RELIABILITY ESTIMATION
CHAPTER 5 STRESS-STRENGTH RELIABILITY ESTIMATION 5. Introduction There are appliances (every physical component possess an inherent strength) which survive due to their strength. These appliances receive
More informationMVE051/MSG Lecture 7
MVE051/MSG810 2017 Lecture 7 Petter Mostad Chalmers November 20, 2017 The purpose of collecting and analyzing data Purpose: To build and select models for parts of the real world (which can be used for
More informationLogit Models for Binary Data
Chapter 3 Logit Models for Binary Data We now turn our attention to regression models for dichotomous data, including logistic regression and probit analysis These models are appropriate when the response
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More information9. Logit and Probit Models For Dichotomous Data
Sociology 740 John Fox Lecture Notes 9. Logit and Probit Models For Dichotomous Data Copyright 2014 by John Fox Logit and Probit Models for Dichotomous Responses 1 1. Goals: I To show how models similar
More informationA New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
More informationThe Delta Method. j =.
The Delta Method Often one has one or more MLEs ( 3 and their estimated, conditional sampling variancecovariance matrix. However, there is interest in some function of these estimates. The question is,
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationProbability Distributions: Discrete
Probability Distributions: Discrete Introduction to Data Science Algorithms Jordan Boyd-Graber and Michael Paul SEPTEMBER 27, 2016 Introduction to Data Science Algorithms Boyd-Graber and Paul Probability
More informationStatistical Intervals (One sample) (Chs )
7 Statistical Intervals (One sample) (Chs 8.1-8.3) Confidence Intervals The CLT tells us that as the sample size n increases, the sample mean X is close to normally distributed with expected value µ and
More informationTABLE OF CONTENTS - VOLUME 2
TABLE OF CONTENTS - VOLUME 2 CREDIBILITY SECTION 1 - LIMITED FLUCTUATION CREDIBILITY PROBLEM SET 1 SECTION 2 - BAYESIAN ESTIMATION, DISCRETE PRIOR PROBLEM SET 2 SECTION 3 - BAYESIAN CREDIBILITY, DISCRETE
More informationNCSS Statistical Software. Reference Intervals
Chapter 586 Introduction A reference interval contains the middle 95% of measurements of a substance from a healthy population. It is a type of prediction interval. This procedure calculates one-, and
More informationRobust Critical Values for the Jarque-bera Test for Normality
Robust Critical Values for the Jarque-bera Test for Normality PANAGIOTIS MANTALOS Jönköping International Business School Jönköping University JIBS Working Papers No. 00-8 ROBUST CRITICAL VALUES FOR THE
More informationAnalysis of Microdata
Rainer Winkelmann Stefan Boes Analysis of Microdata Second Edition 4u Springer 1 Introduction 1 1.1 What Are Microdata? 1 1.2 Types of Microdata 4 1.2.1 Qualitative Data 4 1.2.2 Quantitative Data 6 1.3
More informationAn analysis of momentum and contrarian strategies using an optimal orthogonal portfolio approach
An analysis of momentum and contrarian strategies using an optimal orthogonal portfolio approach Hossein Asgharian and Björn Hansson Department of Economics, Lund University Box 7082 S-22007 Lund, Sweden
More informationu (x) < 0. and if you believe in diminishing return of the wealth, then you would require
Chapter 8 Markowitz Portfolio Theory 8.7 Investor Utility Functions People are always asked the question: would more money make you happier? The answer is usually yes. The next question is how much more
More informationAn Intertemporal Capital Asset Pricing Model
I. Assumptions Finance 400 A. Penati - G. Pennacchi Notes on An Intertemporal Capital Asset Pricing Model These notes are based on the article Robert C. Merton (1973) An Intertemporal Capital Asset Pricing
More informationThis homework assignment uses the material on pages ( A moving average ).
Module 2: Time series concepts HW Homework assignment: equally weighted moving average This homework assignment uses the material on pages 14-15 ( A moving average ). 2 Let Y t = 1/5 ( t + t-1 + t-2 +
More informationAnnual risk measures and related statistics
Annual risk measures and related statistics Arno E. Weber, CIPM Applied paper No. 2017-01 August 2017 Annual risk measures and related statistics Arno E. Weber, CIPM 1,2 Applied paper No. 2017-01 August
More informationRisk Measuring of Chosen Stocks of the Prague Stock Exchange
Risk Measuring of Chosen Stocks of the Prague Stock Exchange Ing. Mgr. Radim Gottwald, Department of Finance, Faculty of Business and Economics, Mendelu University in Brno, radim.gottwald@mendelu.cz Abstract
More informationIOP 201-Q (Industrial Psychological Research) Tutorial 5
IOP 201-Q (Industrial Psychological Research) Tutorial 5 TRUE/FALSE [1 point each] Indicate whether the sentence or statement is true or false. 1. To establish a cause-and-effect relation between two variables,
More informationQuantile Regression. By Luyang Fu, Ph. D., FCAS, State Auto Insurance Company Cheng-sheng Peter Wu, FCAS, ASA, MAAA, Deloitte Consulting
Quantile Regression By Luyang Fu, Ph. D., FCAS, State Auto Insurance Company Cheng-sheng Peter Wu, FCAS, ASA, MAAA, Deloitte Consulting Agenda Overview of Predictive Modeling for P&C Applications Quantile
More informationThe Two Sample T-test with One Variance Unknown
The Two Sample T-test with One Variance Unknown Arnab Maity Department of Statistics, Texas A&M University, College Station TX 77843-343, U.S.A. amaity@stat.tamu.edu Michael Sherman Department of Statistics,
More informationIntro to GLM Day 2: GLM and Maximum Likelihood
Intro to GLM Day 2: GLM and Maximum Likelihood Federico Vegetti Central European University ECPR Summer School in Methods and Techniques 1 / 32 Generalized Linear Modeling 3 steps of GLM 1. Specify the
More informationLecture 21: Logit Models for Multinomial Responses Continued
Lecture 21: Logit Models for Multinomial Responses Continued Dipankar Bandyopadhyay, Ph.D. BMTRY 711: Analysis of Categorical Data Spring 2011 Division of Biostatistics and Epidemiology Medical University
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationPortfolio Construction Research by
Portfolio Construction Research by Real World Case Studies in Portfolio Construction Using Robust Optimization By Anthony Renshaw, PhD Director, Applied Research July 2008 Copyright, Axioma, Inc. 2008
More informationMulti-armed bandits in dynamic pricing
Multi-armed bandits in dynamic pricing Arnoud den Boer University of Twente, Centrum Wiskunde & Informatica Amsterdam Lancaster, January 11, 2016 Dynamic pricing A firm sells a product, with abundant inventory,
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationROM Simulation with Exact Means, Covariances, and Multivariate Skewness
ROM Simulation with Exact Means, Covariances, and Multivariate Skewness Michael Hanke 1 Spiridon Penev 2 Wolfgang Schief 2 Alex Weissensteiner 3 1 Institute for Finance, University of Liechtenstein 2 School
More informationREINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS
REINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS By Siqi Chen, Madeleine Min Jing Leong, Yuan Yuan University of Illinois at Urbana-Champaign 1. Introduction Reinsurance contract is an
More informationThe application of linear programming to management accounting
The application of linear programming to management accounting After studying this chapter, you should be able to: formulate the linear programming model and calculate marginal rates of substitution and
More informationCourse objective. Modélisation Financière et Applications UE 111. Application series #2 Diversification and Efficient Frontier
Course objective Modélisation Financière et Applications UE 111 Application series #2 Diversification and Efficient Frontier Juan Raposo and Fabrice Riva Université Paris Dauphine The previous session
More informationMean-Variance Portfolio Choice in Excel
Mean-Variance Portfolio Choice in Excel Prof. Manuela Pedio 20550 Quantitative Methods for Finance August 2018 Let s suppose you can only invest in two assets: a (US) stock index (here represented by the
More informationMAKING SENSE OF DATA Essentials series
MAKING SENSE OF DATA Essentials series THE NORMAL DISTRIBUTION Copyright by City of Bradford MDC Prerequisites Descriptive statistics Charts and graphs The normal distribution Surveys and sampling Correlation
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
More informationConover Test of Variances (Simulation)
Chapter 561 Conover Test of Variances (Simulation) Introduction This procedure analyzes the power and significance level of the Conover homogeneity test. This test is used to test whether two or more population
More informationVERSION 7.2 Mplus LANGUAGE ADDENDUM
VERSION 7.2 Mplus LANGUAGE ADDENDUM This addendum describes changes introduced in Version 7.2. They include corrections to minor problems that have been found since the release of Version 7.11 in June
More informationThe following content is provided under a Creative Commons license. Your support
MITOCW Recitation 6 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To make
More informationAn Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process
Computational Statistics 17 (March 2002), 17 28. An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Gordon K. Smyth and Heather M. Podlich Department
More informationEstimation of Volatility of Cross Sectional Data: a Kalman filter approach
Estimation of Volatility of Cross Sectional Data: a Kalman filter approach Cristina Sommacampagna University of Verona Italy Gordon Sick University of Calgary Canada This version: 4 April, 2004 Abstract
More information23.1 Probability Distributions
3.1 Probability Distributions Essential Question: What is a probability distribution for a discrete random variable, and how can it be displayed? Explore Using Simulation to Obtain an Empirical Probability
More information8. International Financial Allocation
8. International Financial Allocation An Example and Definitions... 1 Expected eturn, Variance, and Standard Deviation.... S&P 500 Example... The S&P 500 and Treasury bill Portfolio... 8.S. 10-Year Note
More informationExact Inference (9/30/13) 2 A brief review of Forward-Backward and EM for HMMs
STA561: Probabilistic machine learning Exact Inference (9/30/13) Lecturer: Barbara Engelhardt Scribes: Jiawei Liang, He Jiang, Brittany Cohen 1 Validation for Clustering If we have two centroids, η 1 and
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationChapter 6 Part 3 October 21, Bootstrapping
Chapter 6 Part 3 October 21, 2008 Bootstrapping From the internet: The bootstrap involves repeated re-estimation of a parameter using random samples with replacement from the original data. Because the
More informationRisk and Return and Portfolio Theory
Risk and Return and Portfolio Theory Intro: Last week we learned how to calculate cash flows, now we want to learn how to discount these cash flows. This will take the next several weeks. We know discount
More information1 The Solow Growth Model
1 The Solow Growth Model The Solow growth model is constructed around 3 building blocks: 1. The aggregate production function: = ( ()) which it is assumed to satisfy a series of technical conditions: (a)
More informationExam 2 Spring 2015 Statistics for Applications 4/9/2015
18.443 Exam 2 Spring 2015 Statistics for Applications 4/9/2015 1. True or False (and state why). (a). The significance level of a statistical test is not equal to the probability that the null hypothesis
More informationLecture 9. Probability Distributions. Outline. Outline
Outline Lecture 9 Probability Distributions 6-1 Introduction 6- Probability Distributions 6-3 Mean, Variance, and Expectation 6-4 The Binomial Distribution Outline 7- Properties of the Normal Distribution
More informationExpected utility theory; Expected Utility Theory; risk aversion and utility functions
; Expected Utility Theory; risk aversion and utility functions Prof. Massimo Guidolin Portfolio Management Spring 2016 Outline and objectives Utility functions The expected utility theorem and the axioms
More informationBayesian Multinomial Model for Ordinal Data
Bayesian Multinomial Model for Ordinal Data Overview This example illustrates how to fit a Bayesian multinomial model by using the built-in mutinomial density function (MULTINOM) in the MCMC procedure
More informationBEST LINEAR UNBIASED ESTIMATORS FOR THE MULTIPLE LINEAR REGRESSION MODEL USING RANKED SET SAMPLING WITH A CONCOMITANT VARIABLE
Hacettepe Journal of Mathematics and Statistics Volume 36 (1) (007), 65 73 BEST LINEAR UNBIASED ESTIMATORS FOR THE MULTIPLE LINEAR REGRESSION MODEL USING RANKED SET SAMPLING WITH A CONCOMITANT VARIABLE
More informationChapter 10 Inventory Theory
Chapter 10 Inventory Theory 10.1. (a) Find the smallest n such that g(n) 0. g(1) = 3 g(2) =2 n = 2 (b) Find the smallest n such that g(n) 0. g(1) = 1 25 1 64 g(2) = 1 4 1 25 g(3) =1 1 4 g(4) = 1 16 1
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationLecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall Financial mathematics
Lecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall 2014 Reduce the risk, one asset Let us warm up by doing an exercise. We consider an investment with σ 1 =
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationEstimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function
Australian Journal of Basic Applied Sciences, 5(7): 92-98, 2011 ISSN 1991-8178 Estimating the Parameters of Closed Skew-Normal Distribution Under LINEX Loss Function 1 N. Abbasi, 1 N. Saffari, 2 M. Salehi
More informationTutorial 4 - Pigouvian Taxes and Pollution Permits II. Corrections
Johannes Emmerling Natural resources and environmental economics, TSE Tutorial 4 - Pigouvian Taxes and Pollution Permits II Corrections Q 1: Write the environmental agency problem as a constrained minimization
More informationLecture 9. Probability Distributions
Lecture 9 Probability Distributions Outline 6-1 Introduction 6-2 Probability Distributions 6-3 Mean, Variance, and Expectation 6-4 The Binomial Distribution Outline 7-2 Properties of the Normal Distribution
More informationAlgebraic Problems in Graphical Modeling
Algebraic Problems in Graphical Modeling Mathias Drton Department of Statistics University of Chicago Outline 1 What (roughly) are graphical models? a.k.a. Markov random fields, Bayesian networks,... 2
More informationTwo-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models
Economic Review (Otaru University of Commerce), Vo.59, No.4, 4-48, March, 009 Two-term Edgeworth expansions of the distributions of fit indexes under fixed alternatives in covariance structure models Haruhiko
More informationvalue BE.104 Spring Biostatistics: Distribution and the Mean J. L. Sherley
BE.104 Spring Biostatistics: Distribution and the Mean J. L. Sherley Outline: 1) Review of Variation & Error 2) Binomial Distributions 3) The Normal Distribution 4) Defining the Mean of a population Goals:
More informationThe Baumol-Tobin and the Tobin Mean-Variance Models of the Demand
Appendix 1 to chapter 19 A p p e n d i x t o c h a p t e r An Overview of the Financial System 1 The Baumol-Tobin and the Tobin Mean-Variance Models of the Demand for Money The Baumol-Tobin Model of Transactions
More informationBlack-Litterman Model
Institute of Financial and Actuarial Mathematics at Vienna University of Technology Seminar paper Black-Litterman Model by: Tetyana Polovenko Supervisor: Associate Prof. Dipl.-Ing. Dr.techn. Stefan Gerhold
More informationSubject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018
` Subject CS1 Actuarial Statistics 1 Core Principles Syllabus for the 2019 exams 1 June 2018 Copyright in this Core Reading is the property of the Institute and Faculty of Actuaries who are the sole distributors.
More information1 Residual life for gamma and Weibull distributions
Supplement to Tail Estimation for Window Censored Processes Residual life for gamma and Weibull distributions. Gamma distribution Let Γ(k, x = x yk e y dy be the upper incomplete gamma function, and let
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationThe 2 nd Order Polynomial Next Bar Forecast System Working Paper August 2004 Copyright 2004 Dennis Meyers
The 2 nd Order Polynomial Next Bar Forecast System Working Paper August 2004 Copyright 2004 Dennis Meyers In a previous paper we examined a trading system, called The Next Bar Forecast System. That system
More informationRISK ADJUSTMENT FOR LOSS RESERVING BY A COST OF CAPITAL TECHNIQUE
RISK ADJUSTMENT FOR LOSS RESERVING BY A COST OF CAPITAL TECHNIQUE B. POSTHUMA 1, E.A. CATOR, V. LOUS, AND E.W. VAN ZWET Abstract. Primarily, Solvency II concerns the amount of capital that EU insurance
More informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationthe display, exploration and transformation of the data are demonstrated and biases typically encountered are highlighted.
1 Insurance data Generalized linear modeling is a methodology for modeling relationships between variables. It generalizes the classical normal linear model, by relaxing some of its restrictive assumptions,
More informationChapter 14 : Statistical Inference 1. Note : Here the 4-th and 5-th editions of the text have different chapters, but the material is the same.
Chapter 14 : Statistical Inference 1 Chapter 14 : Introduction to Statistical Inference Note : Here the 4-th and 5-th editions of the text have different chapters, but the material is the same. Data x
More informationBackground. opportunities. the transformation. probability. at the lower. data come
The T Chart in Minitab Statisti cal Software Background The T chart is a control chart used to monitor the amount of time between adverse events, where time is measured on a continuous scale. The T chart
More informationTHE DISTRIBUTION OF LOAN PORTFOLIO VALUE * Oldrich Alfons Vasicek
HE DISRIBUION OF LOAN PORFOLIO VALUE * Oldrich Alfons Vasicek he amount of capital necessary to support a portfolio of debt securities depends on the probability distribution of the portfolio loss. Consider
More informationChapter 8. Markowitz Portfolio Theory. 8.1 Expected Returns and Covariance
Chapter 8 Markowitz Portfolio Theory 8.1 Expected Returns and Covariance The main question in portfolio theory is the following: Given an initial capital V (0), and opportunities (buy or sell) in N securities
More information1.1 Some Apparently Simple Questions 0:2. q =p :
Chapter 1 Introduction 1.1 Some Apparently Simple Questions Consider the constant elasticity demand function 0:2 q =p : This is a function because for each price p there is an unique quantity demanded
More informationPARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS
PARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS Melfi Alrasheedi School of Business, King Faisal University, Saudi
More informationStatistics 431 Spring 2007 P. Shaman. Preliminaries
Statistics 4 Spring 007 P. Shaman The Binomial Distribution Preliminaries A binomial experiment is defined by the following conditions: A sequence of n trials is conducted, with each trial having two possible
More informationModelling, Estimation and Hedging of Longevity Risk
IA BE Summer School 2016, K. Antonio, UvA 1 / 50 Modelling, Estimation and Hedging of Longevity Risk Katrien Antonio KU Leuven and University of Amsterdam IA BE Summer School 2016, Leuven Module II: Fitting
More informationPhD Qualifier Examination
PhD Qualifier Examination Department of Agricultural Economics May 29, 2015 Instructions This exam consists of six questions. You must answer all questions. If you need an assumption to complete a question,
More informationMissing Data. EM Algorithm and Multiple Imputation. Aaron Molstad, Dootika Vats, Li Zhong. University of Minnesota School of Statistics
Missing Data EM Algorithm and Multiple Imputation Aaron Molstad, Dootika Vats, Li Zhong University of Minnesota School of Statistics December 4, 2013 Overview 1 EM Algorithm 2 Multiple Imputation Incomplete
More information1 Economical Applications
WEEK 4 Reading [SB], 3.6, pp. 58-69 1 Economical Applications 1.1 Production Function A production function y f(q) assigns to amount q of input the corresponding output y. Usually f is - increasing, that
More information