APPLICATION OF LESLIE MATRIX MODELS TO WILD TURKEY POPULATIONS
|
|
- Earl Collins
- 6 years ago
- Views:
Transcription
1 APPLICATION OF LESLIE MATRIX MODELS TO WILD TURKEY POPULATIONS by Wen-Ching Li Nrth Carlina State University Department f Statistics Bimathematics Graduate Prgram Raleigh 1994
2 CONTENTS ABSTRACT 1 THE LESLIE MATRIX MODEL 2 Intrductin 2 Mdel structure 2 Parameter estimatin 7 Eigensystem 9 Stable age distributin 11 Applicatin and generalizatin 14 APPLICATION TO EASTERN WILD TURKEY IN lowa Intrductin Objectives Mdel structure Parameter estimatin Result and discussin Future wrk
3 ABSTRACT The Leslie matrix mdel is a ppular and useful methd fr ppulatin prjectin based n age specific survival and fecundity rates. The first part f this prject will describe the therems, prperties, applicatin, and develpment f the Leslie matrix mdel in detail. In the secnd part, specific Leslie mdels are develped fr wild turkeys in Iwa. The Leslie matrix mdel has becme very ppular because f its easy understanding by bilgists and because many theretical prperties can be easily studied using matrix algebra. In the example fr Iwa wild turkey ppulatins, a simple, linear, deterministic mdel based n the Leslie mdel is built. The mdel's simplicity makes it easy t understand and analyze. Hwever, its simplicity als means that it is far frm bilgical reality. There are t few parameters t adequately describe the cmplexity f a real turkey ppulatin. Hwever, it is a gd mdel n which t build mre cmplex and realistic mdels which may allw stchasticity and density dependence f survival and fecundity rates. 1
4 THE LESLIE MATRIX MODEL Intrductin Births and deathsare age-dependent in mst species except very simple rganisms. It is usually necessary t take accunt f a ppulatin's age structure t get a reliable descriptin f its dynamics. The English bilgist P.R.Leslie (1945) intrduced a mdel which uses the agespecific rates f fertility and mrtality f a ppulatin. He develped the mdel based n the use f matrices. The matrix frm makes the mdel flexible and mathematically easy t study. This basic mdel is an increasingly ppular tl fr describing ppulatin dynamics f plants (Usher 1969) and animals (Yu 1990, Cruse 1987). The Leslie matrix mdel is widely used t prject the present state f a ppulatin int the future, either as an attempt t frecast the age distributin, r as a way t evaluate life histry hyptheses by cnsidering different sets f survival and fecundity parameters. Mdel structure In the simplest Leslie matrix mdel (Leslie 1945), nly females are cnsidered. The female ppulatin can be divided int several categries by age r by size. The classificatin styles and their differences will be discussed later in this sectin. If it's gruped by age, the grup intervals are suppsed t have equal length f time. Fr example, it may be 5 years fr human and 2 years fr large whales (Cullen 1985). We assume that the survival and fecundity rates 2
5 f each categry are cnstant ver time and therefre nt dependent n ppulatin density. Define nj(t) = the number f females alive in the grup i at time t, Pj = the prbability that a female f grup i at time t will be alive in grup i + 1 at time t+ 1, (t+ 1 means t+ne unit f time), Pi lies between 0 and 1, Fj = the number f daughters brn per female in grup i frm time t t time t+ 1, let i = 1,2,...,x, then nl(t+l) = F1*n](t) + F 2 *n 2 (t) Fx*nx(t) n 2 (t+l) = n 2 (t+ 1) = P1*n)(t) P 2 *n 2 (t) (1.1) These equatins can be represented in matrix frm as n(t+ 1) = A n(t) (1.2) where n(t) is a clumn vectr with cmpnents f the number f individuals in each recgnized categry at time t. A is a nnnegative, square matrix f rder x with all the elements zer except thse in the first rw and in the subdiagnal immediately belw the principal diagnal. 3
6 F i F 2 F 3 F X - 1 F x P 1 a 0 a a A = a P 2 0 a a a 0 P 3 a a (1.3) a 0 0 P X A is knwn as the "Leslie matrix" based n age classificatins r the "ppulatin prjectin matrix" mdeled n stage classificatins (Lefkvitch 1965, Grenendae et al. 1988, Cruse et al. 1987). The prjectin matrix was intrduced by Lefkvitch (1965). The stages in the prjectin matrix mdel are nt necessarily related t age. Its main assumptin is that all individuals in a given stage are subject t identical mrtality, grwth, and fecundity schedules. These matrices are mre cmplicated in appearance than the Leslie matrices, but have almst the same analytical methd. Fig. 1.1 shws the life cycle graph depending n different decmpsitin f the life cycle. The arrws in the graph indicate the transitins which are pssible fr an individual frm ne time t the next. The Leslie matrix crrespnding with the life cycle graph (A) in Fig 1.1 is just the same as (1.3) and much simpler than the prjectin matrix fr the stage-classified ppulatin, (8) f Fig 1.1, btained as fllws: 4
7 A F u F 12 F 13 F 14. F 1X P 21 P 22 P P 32 P 33 P P X - 1 x-2 P X - 1 x-l P X - 1 x 0 0 P x x-l P xx (1. 4) A.,/, J, F, 1F 2 Fx ~Pl Age class Age class P 2 Px- 1 Age class,- -, ~ x B.,.l J, F ll l -.J P 21 F l2 Fix Stage Stage P 32.. Pxx- 1 Stage ;, I... ~ J 1 2 P r- Px- 1 x IJ' ~ r x Fig 1.1 Life cycle graphs fr tw chices f categries. (A) Age-classified ppulatin. (B) One example f stage-classified ppulatin. 5
8 In the stage-classified mdel, the pssibility f an individual f grup i mving backwards t the previus grup, mving frwards t the next grup, r staying the same grup at next time perid is cnsidered. The chice f categries is a very imprtant step in establishing a ppulatin matrix mdel (Grenendael et al. 1988). Since smetimes it is hard t determine age-specific parameters f the mdel fr sme rganisms with variable develpment, stage categries may be mre apprpriate than age categries (e.g. fr turtles, Cruse 1987). Fr ther rganisms like mammals and birds, their develpment is mre dependent n age than size and therefre age categries are used frequently. If age and size are bth imprtant and significant interactins exist between them, a twdimensinal mdel using bth stages and age classes simultaneusly can be develped(yu 1990). In Yu' s study, develpmental stages are as meaningful as chrnlgical age classes in the crn earwrm (CEW) management. He divided the CEW ppulatin int five stages and split the first stage int fur age grups, the remaining stages int six age grups, except taking the last stage as a single age grup. Thus, his prjectin matrix is a 23*23 matrix taking bth stage and age int accunt. In this prject, I will mainly discuss the prperties f an age-classified matrix mdel based n frm (1.3). The mdel has linear, time invariant functins which prject the current state f the ppulatin int the next state. The female ppulatin alne will be cnsidered. 6
9 Parameter estimatin Survival prbability Pi The values f Pi can be btained frm a life table r experimental data. It is usually assumed that (Leslie 1945, Ple 1974) L i + 1 p.=- ~ L ~ (1.5) where Li = the number alive in the age grup i t i+1 in the statinary r life table age distributin. Pi is assumed t be cnstant ver a unit f time and nt dependent n the ttal number in the ppulatin. Actually, seasnal differences and density dependence f the survival rate d exist, but we assume that they are nt significant enugh t cunt. Fecundity rates F j F; is determined nt nly by the average number f female births per mther per time perid, m j in the life table, but als by infant survival rates (Leslie 1945, Grenendael et al. 1988). Ppulatins can be divided int tw cases f breeding systems, 'birth-flw' ppulatin with cntinuus reprductin and 'birth-pulse' ppulatins with discrete reprductin (Caughley 1977). Ppulatins, e.g. human, f the first case prduce ffspring at a rate almst cnstant thrughut the year and thse f the secnd case prduce ffspring ver a restricted seasn, e.g. 7
10 blue whale. In birth-flw ppulatins, the estimatin f F i invlves integratin f cntinuus fertility and infant mrtality functins ver the time interval (Leslie 1945). Fr birth pulse ppulatins reprductin is during a shrt perid f time per year. It is imprtant t knw the exact breeding perid with respect t the time step in birth-pulse ppulatins. Tw different ways f cunting the ppulatin, census just befre r after reprductin, will influence the estimatin f F i (Cullen 1985). Sme studies cunt the ppulatin just befre reprductin and we illustrate that here. Define M j = average number f female ffspring per female in categry i brn between t and t+i, S = the ffspring survival rate between birth and the time when they are cunted as part f the ppulatin, then F i = S * M i 1 = 1,2,...,x (1,6) F i, the age specific fecundity rate, is assumed t be cnstant and nt dependent n density just as the age specific survival rate, Pi' is. 8
11 Eigensystem One advantage f using Leslie matrices is the flexibility f their mathematical frmulatin. By applying linear algebra techniques t the mdel, several bilgical phenmena can be interpreted. Cnsider the hmgeneus discrete-time system (1.2), n(t+ 1) = A net), where A IS a square matrix (1.3). If A has the frmula, (1. 7) then AI,...,Axare called the eigenvalues f A with I AI I > I A2 I ~ I A3 I ~... ~ I Ax I and el'0..,exare the crrespnding linearly independent eigenvectrs. The eigenvalue Al f greatest magnitude is the dminant eigenvalue. Then the slutins t n(t+ 1) = A net) have the frm (Ducet and Slep 1992), nl(t) = all All + a l2 A2' + n2(t) = a21 All + a 22 A2 1 + (1.8) njt) = axl All + ax2 A a xx A/ Fr large t, althugh the influence f A2.3...,x des nt necessarily disappear except that their abslute values are smaller than 1, the first term with the dminant eigenvalue will grw faster I than the thers. Therefre, nj(t) will becme similar t the single term aila I. Each age class eventually will grw expnentially at the rate f the dminant eigenvalue, AI, per time perid. The system (1.2) is asympttically stable if and nly if the eigenvalues f A all have 9
12 magnitude less than ne (Luenberger 1979), that js, I A; I < I fr every i. The state vectr n(t) will tend t an equilibrium pint (the rigin) fr any initial cnditin if the system is asympttically stable. Once the system state vectr is equal t an equilibrium pint, it will remain equal t that fr all future time. The system is called marginally stable jf n eigenvalues have magnitude greater than ne but ne r mre is exactly equal t 1. eigenvectr is an equilibrium pint. If unity is an eigenvalue f A, then the crrespnding If unity is nt an eigenvalue f A, the rigin is the nly equilibrium pint f the system (Luenberger 1979). Therefre, if unity is the dminant eigenvalue, the state vectr eventually will be a specific multiple f the assciated eigenvectr and wn't change in future. It means that ppulatin size will remain the same fr each age class after it accmplishes that pint. When the dminant eigenvalue AI = 1, we say the ppulatin is statinary. Each eigenvalue defines nt nly a characteristic grwth rate but als a characteristic frequency f scillatin (Luenberger 1979). Fr a discrete-time system, n scillatins are derived frm a psitive eigenvalue, in ther wrds, scillatins are due t a negative r cmplex eigenvalue with perid 2 r with a larger perid respectively (Ducet et al. 1992). Sme eigenvalue-eigenvectr therems fr Leslie matrices are useful in determining the develpment f a ppulatin (Cullen 1985). (1) A Leslie matrix has at least ne psitive real eigenvalue. (2) If there are at least tw cnsecutive age classes that are fertile, a psitive real dminant eigenvalue always exists. 10
13 Stable age distributin Accrding t the utstanding prperties f the dminant eigenvalue, a ppulatin grwing based n the Leslie matrix will have an age cmpsitin which is entirely determined by its Leslie matrix and des nt depend n the initial cmpsitin after a perid f time. Then frm sme year n, next year's ppulatin will be a multiple f this year's. This age cmpsitin is called a stable age distributin. (1. 9) The multiple Al is the dminant eigenvalue f the Leslie matrix A (1.3) and if e l IS the assciated eigenvectr, the stable age distributin can be btained frm el' Let (1.10) and v = VI + v vx' then the vectr S is called the stable age distributin f the ppulatin (Cullen 1985). s= (1.11) 11
14 Tw different Leslie matrices may have the same AI and the same stable age cmpsitin, that is, given a value f AI and a stable age vectr, mre than ne Leslie matrix can be determined (Pielu 1977). Fr example, , 6 6~. 0] x= [ 0.5 and ~. 5] y= X and Y bth have the same dminate eigenvalue, AI = 1.5, and the same stable vectr S =[ ]T. Since they are actually unlike, they differ in ther eigenvalues. The apprach t S depends n hw much larger I AI I is than I A 2 I,..., I Ax I. The larger the difference between them, the faster the ppulatin mves tward stability. Therefre, X and Y have varius appraching rates t the stable age cmpsitin. The cnvergence twards the stable age cmpsitin ccurs nly if AI is "strictly" larger than the ther eigenvalues. Beginning with an unstable r stable age distributin, a ppulatin in an unlimited envirnment will reach a stable age distributin in time, whether increasing, decreasing, r staying at the same density (Ple 1974). Unfrtunately, there is an exceptin t this rule (Ple 1974, Ducet and Slep 1992). When there exists ne f the negative r cmplex eigenvalues has the same abslute value as AI, the cntributin f All will never 12
15 dminate the ther cntributins and the ppulatin will scillate frever. This happens when the entries f the first rw f Leslie matrix are all zer except the rightmst ne. Fr example (Ple 1974), A= The eigenvalues f A are and The dminant eigenvalue is 1, s the ppulatin will neither increase r decrease as a trend. Nte that A 2 = a a - 6 A 3 =I and s frth. The ppulatin will scillate with perid 3 time units and never apprach stability unless it begins with a stable age distributin. 13
16 Applicatin and Generalizatin In the riginal Leslie matrix mdel, the last age class is assumed t be remved frm the ppulatin after a time unit. S the entry with rw x and clumn x is 0 in the matrix (1.3). Actually, individuals in the last age class may nt die after ne time unit and their survival prbability and fecundity rate will just stay the same until they die. It indicates that we allw members f the last grup t live and reprduce fr sme years. Besides, lder individuals are assumed t be essentially the same as thse f the last age grup. Let P x = the prbability that an individual f grup x at time t will alive and stay in the grup x at time t+ 1. Then the Leslie matrix can be written as F 1 F 2 F 3 F X - 1 F x P 1 a a a a A= a P a a P 3 a a (1.12 ) a a a P X - 1 P x and the last equatin f (1.1), n x (t+l) = Px-l nx_l(t), is replaced by which is an aggregatin f all remaining age classes. 14
17 The basic Leslie matrix mdel is a deterministic mdel. The parameters, survival and reprductive rates, are cnstant ver time. In real ppulatin, thse parameters may change with many factrs. In the cases f human ppulatin, the survival rates, especially amng the elderly, will change with new medical advances r better feeding habits, and the fecundity rates will be affected by changing scial attitudes tward marriage and family. The effect f harvest and seasnal differences n survival prbabilities shuld be cnsidered fr animal r plant ppulatins. If the parameter values fr years are randmly generated frm a specific prbability distributin, the mdel will be stchastic. Only females are cnsidered in Leslie's mdel (1945). In bisexual ppulatins, the basic mdel can be mdified t keep track f bth sexes (Pielu 1977, Cullen 1985). The Leslie matrix in its classical frm makes n allwance fr density dependent ppulatin grwth and it is fundamentally linear. Cnsidering the effect f density dependence n fertility and survival rates will turn the mdel int a nnlinear ne (Pielu 1977). The size f the current ppulatin will play an imprtant rle in determining thse parameters and the initial size has a lasting effect n the ppulatin's future chances f reprducing and surviving. 15
18 APPLICATION TO EASTERN WILD TURKEY IN lowa Intrductin The eastern wild turkey (Meleagris gallpav silvestris), ne f five subspecies, inhabits rughly the eastern half f the United States. Turkey hunting has a substantial ecnmic effect in many rural cmmunities. It nt nly brings in revenue because f the actual turkey hunting, but it als enhances the develpment f the related industries f turkey-hunting clthes and equipment. Imprvement f the knwledge f turkey ppulatin dynamics is imprtant t frmulate hunting regulatins and ther turkey management practices. Much research abut the rates f reprductin, mrtality, and survival, and the mvement f wild turkeys has been dne (Dicksn 1992). Hwever, there were few mdels established fr evaluating turkey ppulatins. Fr a specific cmbinatin f survival and reprductive rates, a mdel may tell us whether a wild turkey ppulatin will grw, decline, r remain stable. The survival rates will bviusly als depend n the level f harvest. By changing the harvest level we can see the effect n the ppulatin size and cnsider changes t the hunting regulatins and enhance the management f wild turkeys. A Leslie matrix mdel is develped fr the ppulatin dynamics f eastern wild turkeys in Iwa. The parameters used in the mdel are based n Iwa wild turkey studies (Suchy et al. 16
19 1983, Dicksn 1992), and the ppulatin characteristics will be interpreted by the mdel analysis. Objectives 1. T build a mdel in rder t predict the ppulatin size and age structure f Iwa wild turkeys in the future. 2. T see the effect f harvest level (percent f ttal ppulatin) n the ppulatin age structure and grwth. 3. T find the mst apprpriate wild turkey hunting seasns. Mdel structure An age-classified mdel is built here. There culd be several stages in the age structure f the wild turkey ppulatin. A three-stage mdel is chsen in rder t simplify the mdeling prcedure. The first categry is "pults", aged frm 0 t 1, the secnd categry is "yearlings", aged frm 1 t 2, and the last categry is " adults", aged 2 and lder. Reprductin ccurs frm yearlings nwards. The time unit is ne year. Only females are cnsidered in this mdel. Let Pi = the prbability that a female f grup i at time t will be alive in grup i+ 1 fr i= 1,2 r will stay in grup i fr i =3 at time t+ 1. F j = the number f daughters brn per female in grup i frm time t t time t+ 1. i = 1, 2, r 3. 17
20 The life cycle graph is as fllws, l I r Pults 1 PI Yearlings 2 P 2 Adults 3,. i \ 10- The state variables are the numbers f pults, yearlings, and adults, dented as n l, n 2, and n3 respectively. Then, n 1 (t+l) = F z n z (t)+f 3 n 3 (t) n z ( t +1) = PIn I ( t) (2. 1 ) n 3 (t+l) = P z n Z (t)+p 3 n 3 (t) A matrix frm can be derived frm thse equatins. n(t+l)==an(t) (2.2) I.e. n 1 (t+l) 0 F z F 3 n 1 (t) n z (t+l) PI 0 a n 3 (t+l) 0 Pz P 3 n 3 ( t) * n z ( t) (2.3) 18
21 The Leslie matrix A, a square matrix f rder 3, is defined in this three-age-class mdel as A= (2.4) Estimating parameters annually we verlk sme cmplicatins, fr example, different mrtality rates in different seasns. A number f studies n estimating the survival and reprductive prcesses were well dne fr wild turkeys. It is nt t hard t estimate the parameters, Pi and F j, t use in the mdel. This is the reasn why an age-classified mdel is chsen rather than a stage (size)-classified mdel in this prject. Assumptins The basic assumptins are as fllws: 1. The survival and reprductive rates are cnstant ver time within each age-class and are therefre nt dependent n ppulatin density. 2. The individuals in each grup always g t next grup frm ne time t the next except thse in the last age-class which may stay fr sme years befre they die. 3. The cnditins f each individual in the same stage are similar, that is, all individuals in that stage have the same survival and reprductive rates. In practice survival and fecundity d vary frm year t year due t envirnmental factrs. If the Leslie mdel is still applied using averages there is a degree f apprximatin invlved. In practice survival and fecundity may als be density dependent. 19
22 Parameter estimatin Survival prbability Pi 1983). The survival rates are derived frm a study n Iwa wild turkey ppulatins (Suchy's et al. Pl=0.445 P P They are the average estimates f the data frm 1977 t Fecundity rates F j Frm Dicksn (1992), sme infrmatin abut Iwa turkey ppulatins is imprtant fr estimating the reprductive rates. Estimates f nesting and renesting rates (Table 2.1) indicates that a higher prprtin f adult hens attempt t nest. Clutch size (Table 2.2) is als larger fr adult hens than fr yearlings. Since reprductin starts frm yearling-stage, the value f F I is zer. RecalI that F j shuld be determined nt nly by the number f daughters brn per female per time perid but als by infant survival prbabilities (hatching success). Hatching success = number eggs h~tching = 0.85 clutch s~ze fr successful nests. 20
23 Table 2.1 Nesting rates and nest success fr Iwa wild turkey hens. (Dicksn 1992) First nest Renest nesting nest renesting nest rate success rate success yearling adults Table 2.2 Clutch size fr Iwa wild turkey hens. (Dicksn 1992) clutch SIze First nest Renest yearlings adults Using the values in Table 2.1, Table 2.2, and the hatching success, F 2 and F 3 can be estimated as fllws, 21
24 F 2 estimatin F 2 = number f female pul ts number f females = number f female eggs * number f female pul ts number f females number f female eggs = female eggs f first nesting + female eggs f renesting number f females * number f femal e pul ts number f female eggs (0.42) (0.56) (1) (~) +0 = *0.85.E 3 estimatin Use the similar way as abve. (0.97) (0.38) (1) (2...:..!)+(0.97) (0.62) (0.32) (0.52) (1) (1...:1.) * 0.85 =1.860 Therefre, the riginal Leslie matrix f Iwa wild turkey ppulatin in this prject is btained. A= (2.5) 22
25 Result and discussin 1. Stability I used a FORTRAN prgram t find eigenvalues and eigenvectrs f the Leslie matrix A. The initial values f state variables are n1(o) = the initial number f female pults in the ppulatin = n2(o) = the initial number f female yearlings in the ppulatin = n,(o) = the initial number f female adults in the ppulatin = which are the data derived frm Suchy et al. (1983). Fr the riginal matrix (2.5), eigenvalues are A)=1.15, A2,3=-0.27±0.40i and I Al I =1.15>1 which shw that the ppulatin is unstable. Actually, sme parameters will fluctuate at randm in the real situatin. Fr example, nest success fluctuates widely frm year t year in weather cnditins r rainfall (Dicksn 1992). T study the sensitivity f the mdel t individual parameter values and t manage the turkey ppulatins, I will change the parameter values ne by ne until I get the stable situatin. I emphasize that this is dne t investigate mdel prperties but smetimes the results may nt be bilgically reasnable. Let Al be the dminant eigenvalue. Case1 The adult reprductive rate, F3, is decreased frm 1.86 t The Leslie matrix changed t be A= a
26 Then eigenvalues are Al =0.999, A2=0.05, A3=-0.42, I Ai I < 1. It is asympttically stable. In this case, a big change in F3 is required t btain the stable situatin. Case2 The pult survival rate, PI' is decreased frm t Other parameters remain the same values as in the riginal matrix. Thatis A= ~60] The eigenvalues f this mdel are Al =0.998, A2. 3 =-0.19±0.35i and I Ai I < 1. It is a stable system. Case3 The yearling survival rate, P2' is decreased frm t Other parameters remain the same values as in the riginal matrix. Then Al =0.996, A2=0.02, A 3 = Case4 The adult survival rate, P 3, is decreased frm 0.61 t Other parameters remain the same values as in the riginal matrix. Then A( =0.997, Av =-0.42±0.52i. Case5 It's hard t get a stable state if nly F2 is decreased. I A( I is still bigger than I when F 2 had been decreased t 0.01 which is t small fr the yearlings fertility. Case6 I tried t change bth F 2 and F3 simultaneusly. Then when F 2 is reduced frm 0.88 t 0.40 and F3 is changed frm 1.86 t 1.13, AJ =0.999, A2,3 =0.19±0.4li. This is mre bilgically reasnable than case 1. There are several ways t change the values f F 2 and F3 in this case. In casel,2,3,4,and 6, the abslute values f Ai are less than ne. Thse systems are asympttically stable and because unity is nt an eigenvalue f A, the rigin is the nly equilibrium pint f system. The ppulatin is ging t be extinct after a lng run. 24
27 If unity is the dminant eigenvalue f A, instead f ging t be extinct, the ppulatin will eventually achieve a stable age distributin, which is prprtinal t the crrespnding eigenvectr f unity. This situatin is called the marginal stability. In rder t get the marginal stability, I change the parameter values ne by ne again frm the riginal matrix (Table 2.3). Table 2.3 The eigenvalues f marginally stable ppulatins btained by adjusting parameters separately riginal mdified Al A 2 A 3 value value PI i i P P i 0.52i F Sensitivity Analysis I used the change f the ttal ppulatin size based n the change f the parameter as an indicatin f parameter sensitivity. If Q(p) means the ttal ppulatin size (Q) is a functin f the parameter (p), then fr example, 25
28 5% sensitivity f Q t p = Q(p+5%p) -Q(p-5%p) 2*5%*Q(p) Table 2.4 shws that PI' the pult survival rate, is the mst sensitive parameter. In Suchy et al. (1983), they reasned that the mst sensitive parameters wuld be thse that prduced a statinary ppulatin with the smallest change in parameter value. Recall the data in Table 2.3. Table 2.4 The 5 %, 10%, and 15 % sensitivities f the ttal ppulatin size t each parameter Sensitivities Parameters 5 % 10% 15 % PI P P F F A decrease f 41 % in pult survival rate (PI) is the smallest reductin f all parameters t prduce statinary ppulatins. S PI shuld be the mst sensitive parameter which is shwn in Table 2.4. Changes in survival rates f female pults have a great effect n the ppulatin's net rate f change. The fecundity rate f adults is the secnd mst sensitive parameter. 26
29 3. Stable age distributin Recall that any ppulatin beginning with an unstable r stable age distributin will reach a stable age cmpsitin in time, the exceptin being a matrix with the first rw which nly has ne nnzer entry at the rightmst psitin. Frtunately, the turkey mdel is nt that kind f matrix and it will cme t a stable age distributin sner r later depending n the differences between the dminant eigenvalue and the ther nes. Fr the unstable, riginal matrix, the dminant eigenvalue is and the crrespnding eigenvectr is [ ]T. Then the stable age distributin is S=[ ]T. Fig 2.1 shws the simulatin f 15 years. Althugh the ppulatin keeps grwing, the age cmpsitin reaches a cnstant rati after nine years and wn't change in the future. One asympttically stable system is chsen (case 2) t d the simulatin f frecasting the ppulatin size. A= ] Fig 2.2 shws the result f fifteen-year simulatin. After a shrt perid (within five years) f scillatin, the age distributin reaches stability althugh the ppulatin is decreasing. We are interested in the statinary situatin. When the dminant eigenvalue becmes ne, nt nly the age cmpsitin appraches a stable state but als the ppulatin size wn't change eventually. Then the ppulatin will never vergrw r die ut. Pick the statinary cases I have run befre and d the simulatin fr each ne. Fig 2.3, 2.4, 2.5, 2.6 shw the result f thse marginally stable systems gtten by adjusting PI' P 2, P 3, F 3 respectively. There are scillatins in all f them but it des nt take a lng time. 27
30 a aa ('t) pults yearlings adults Q) N '0 :: a C\I 1a ::::J Q. Q. a aa ~ /... /.. ' /..../.. ',...,~...,,...,.. ' ----:-:.::: :-::.-::..':.:':;/.. ~.-:-:-:'... :- /... /.. ' / / /...' ' / r r ~ 5 year / Fig 2.1 Simulatin f riginal ppulatin mdel 28
31 LO Q) V N "00 c: "ca ::J 0 a C") pults yearlings adults C\I /', / '----- /"" /"" , year Fig 2.2 Simulatin f the asympttically stable ppulatin mdel btained by adjusting P1 29
32 LO 0 Q) 'I::t N en c: 0 pults 1a... ::l yearlings a adults 0 0 a. C"') C\I // /.,..--- T'"" ' year Fig 2.3 Simulatin f the stable ppulatin mdel btained by adjusting P1 30
33 LO Q) 0 N... "en c: 1a ::J a. a. 0 ('I) pults yearlings adults C\I... 'f, / :'...,/ year Fig 2.4 Simulatin f the stable ppulatin mdel btained by adjusting P2 31
34 LO 0 Q) V N.c;; c: 0 pults 1a... yearlings ~ adults 0- ('I) C\I..- '.,,,,.,.....t 1\'. " '\ II \\. /"""-_ j\ / -' \J year Fig 2.5 Simulatin f the stable ppulatin mdel btained by adjusting P3 32
35 LO pults yearlings adults Q) 0 N ~.Ci) c: 15 ~ a. 0 0 a. ('I) C\I,... /"- / "-, ~ ,. /,~. ', ' '.. '/.: /'.., ',:~..J.,...r, year Fig 2.6 Simulatin f the stable ppulatin mdel btained by adjusting F3 33
36 4. The effect f harvest Harvest will influence the parameters in yearlings and adults stages directly and in pults stage indirectly. We assume that nly the survival prbabilities f yearlings and adults are affected by the harvest level which is defined as the percentage f the ttal ppulatin. Only fall hunting seasns are cnsidered because the mdel is fr female ppulatins and it's usually nly males can be hunted in the spring seasns. Data frm the bk named "The Wild Turkey" (Dicksn 1992) suggests pssible reductins in P2 and P 3 due t harvest (Table 2.5). Table 2.5 The survival rates f adult hens and yearling hens fr different fall harvest level Survival rates Fall harvest level female adults female yearlings ( % f Ttal ppulatin) (P 3 ) (P 2 ) I changed the parameters f the riginal matrix based n this pssible effect f harvest and tried t figure ut the trend. Accrding t the different harvest level, P 2 and P 3 are replaced by the values in the table. Then the dminant eigenvalue is decreasing with the increase f harvest level (Table 2.6). 34
37 Table 2.6 The eigenvalues f different turkey systems based n the riginal Leslie matrix accrding t different harvest levels Harvest level Al A ± 0.388i ± 0.365i ± 0.342i ± 0.317i These results are surprising and suggest sme inadequacy f ur mdel because when the 20 percent f ttal ppulatin is remved by hunting ne might expect the ppulatin size shuld be smaller and smaller, that is, the dminant eigenvalue shuld be smaller than 1. The result didn't tell us that. Maybe it is because the mdel is t simple r there is smething wrng with ur parameter estimatin. As it has been mentined befre, the survival and fecundity rates will vary frm time t time depending n many factrs, like weather cnditin. Amng thse parameters, the survival rate f pults, PI' will change much widely than the thers. S PI is chsen t be adjusted t d sme mre analysis. If PI changed frm t 0.330, the result f simulatin will be as fllws (Table 2.7). 35
38 Table 2.7 The eigenvalues f turkey ppulatins based n the mdified matrix by changing P I accrding t different levels f harvest Harvest level Al A 2, ± i ± i ± i ± i This result perhaps lks mre plausible than the previus ne. The dminant eigenvalue starts t be belw than 1 when the harvest level is 15, 15 percent turkeys f the whle ppulatin have been hunted. Future wrk This turkey mdel is very simple. The turkey ppulatin is a cmplicated system. There are t few parameters in the mdel t prperly describe the cmplexity f a real turkey ppulatin. The matrix is extensible and mre parameters can be added t the mdel. Interactins between thse parameters culd be cnsidered r the whle ppulatin may be divided int mre stages. Fr example, if the time unit is six mnths, each stage can be divided int tw grups. Then the new mdel will have twice as many as stages the ld mdel has and the parameter estimates will be mre reliable because f less envirnmental variatin. The Leslie matrix A wuld be extended t be 6*6. 36
39 In a real turkey ppulatin, the number f each grup will affect the survival and fecundity rates f the ppulatin. The influence f density may nt be ignred. If a mdel takes accunt f the effect f density dependence n the survival and fecundity rates, then the mdel will be a nnlinear mdel. Anther way t make the mdel pwerful is t allw stchasticity. The Missuri Mdel in the bk (Dicksn 1992) used a stchastic mdel. In that mdel, if the data f thse parameters which can be estimated frm a radi-marked sample f wild turkeys fr specific years is inputed, then the parameters are randmly generated frm a unifrm distributin within the 95 percent cnfidence interval arund the mean f each parameter value fr the data years. This makes the mdel stchastic. Anther type f future wrk wuld be t cnsider survival as different functins f harvest rate. The functins culd be cnstructed t reflect additivity f hunting and natural mrtality as sme degree f cmpensatin (Le. as hunting mrtality increases perhaps natural mrtality decreases t sme extent). This is a very cmplex tpic beynd the scpe f this prject. 37
40 LITERATURE Caughley, G Analysis f vertebrate ppulatins. Jhn Wiley and Sns Ltd, New Yrk, N.Y. 234pp. Cruse, DT., L. B. Crwder, and H. Caswell A stage-based ppulatin mdel fr lggerhead sea turtles and implicatins fr cnservatin. Eclgy 68: Cullen, M.R Linear mdels in bilgy. Ellis Hrwd Ltd, England. 213pp. Dicksn, J.G The wild turkey: bilgy and management. Stackple Bks. Harrisburg, PA.463pp. Ducet, P, and PB. Siep Mathematical mdeling in the life sciences. Ellis Hrwd Ltd, England. 490pp. Grenendael, J.van, H. de Krn, and H. Caswell Prjectin matrices in ppulatin bilgy. Trends in Eclgy and Evlutin 3, Lefkvitch, L. P The study f ppulatin grwth in rganisms gruped by stages. Bimetrics 21, Leslie, P.H On the use f matrices in certain ppulatin mathematics. Bimetrika 33: Luenberger, D.G Intrductin t dynamic systems: thery, mdels and applicatins. Jhn Wiley & Sns, New Yrk, N.Y. 446pp. Pielu, E. C Mathematical eclgy. Jhn Wiley & Sns, Inc., New Yrk, N Y. 385pp. Ple, R. W An intrductin t quantitative eclgy. McGraw-Hill, Inc., New Yrk, N.Y. 532pp. 38
41 Suchy, W.J., W.R. Clark, and T.W. Little Influence f simulated harvest n Iwa wild turkey ppulatins. Prc. Iwa Acad. Sci. 90: Usher, M.B A matrix mdel fr frest management. Bimetrics 25, Yu, Y A Leslie mdel, threshld functin and uncertainty fr chemical cntrl f crn earwrm. Ph.D. Thesis. Nrth Carlina State University. Raleigh, N.C. 148pp. 39
Ecology 302 Lecture III. Exponential Growth (Gotelli, Chapter 1; Ricklefs, Chapter 11, pp )
Eclgy 302 Lecture III. Expnential Grwth (Gtelli, Chapter 1; Ricklefs, Chapter 11, pp. 222-227) Apcalypse nw. The Santa Ana Watershed Prject Authrity pulls n punches in prtraying its missin in apcalyptic
More information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationBootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >
Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationThermodynamics Partial Outline of Topics
Thermdynamics Partial Outline f Tpics I. The secnd law f thermdynamics addresses the issue f spntaneity and invlves a functin called entrpy (S): If a prcess is spntaneus, then Suniverse > 0 (2 nd Law!)
More informationCHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India
CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 11: Mdeling with systems f ODEs In Petre Department f IT, Ab Akademi http://www.users.ab.fi/ipetre/cmpmd/ Mdeling with differential equatins Mdeling strategy Fcus
More informationSIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST. Mark C. Otto Statistics Research Division, Bureau of the Census Washington, D.C , U.S.A.
SIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST Mark C. Ott Statistics Research Divisin, Bureau f the Census Washingtn, D.C. 20233, U.S.A. and Kenneth H. Pllck Department f Statistics, Nrth Carlina State
More informationLecture 13: Markov Chain Monte Carlo. Gibbs sampling
Lecture 13: Markv hain Mnte arl Gibbs sampling Gibbs sampling Markv chains 1 Recall: Apprximate inference using samples Main idea: we generate samples frm ur Bayes net, then cmpute prbabilities using (weighted)
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 2: Mdeling change. In Petre Department f IT, Åb Akademi http://users.ab.fi/ipetre/cmpmd/ Cntent f the lecture Basic paradigm f mdeling change Examples Linear dynamical
More informationHypothesis Tests for One Population Mean
Hypthesis Tests fr One Ppulatin Mean Chapter 9 Ala Abdelbaki Objective Objective: T estimate the value f ne ppulatin mean Inferential statistics using statistics in rder t estimate parameters We will be
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More informationCAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank
CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal
More informationModeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function
www.ccsenet.rg/mer Mechanical Engineering Research Vl. 1, N. 1; December 011 Mdeling the Nnlinear Rhelgical Behavir f Materials with a Hyper-Expnential Type Functin Marc Delphin Mnsia Département de Physique,
More information[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )
(Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well
More informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationKinetic Model Completeness
5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins
More informationCS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007
CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is
More informationAdmissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs
Admissibility Cnditins and Asympttic Behavir f Strngly Regular Graphs VASCO MOÇO MANO Department f Mathematics University f Prt Oprt PORTUGAL vascmcman@gmailcm LUÍS ANTÓNIO DE ALMEIDA VIEIRA Department
More informationLecture 17: Free Energy of Multi-phase Solutions at Equilibrium
Lecture 17: 11.07.05 Free Energy f Multi-phase Slutins at Equilibrium Tday: LAST TIME...2 FREE ENERGY DIAGRAMS OF MULTI-PHASE SOLUTIONS 1...3 The cmmn tangent cnstructin and the lever rule...3 Practical
More informationHomology groups of disks with holes
Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationx x
Mdeling the Dynamics f Life: Calculus and Prbability fr Life Scientists Frederick R. Adler cfrederick R. Adler, Department f Mathematics and Department f Bilgy, University f Utah, Salt Lake City, Utah
More informationPhys. 344 Ch 7 Lecture 8 Fri., April. 10 th,
Phys. 344 Ch 7 Lecture 8 Fri., April. 0 th, 009 Fri. 4/0 8. Ising Mdel f Ferrmagnets HW30 66, 74 Mn. 4/3 Review Sat. 4/8 3pm Exam 3 HW Mnday: Review fr est 3. See n-line practice test lecture-prep is t
More informationNUMBERS, MATHEMATICS AND EQUATIONS
AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t
More informationChE 471: LECTURE 4 Fall 2003
ChE 47: LECTURE 4 Fall 003 IDEL RECTORS One f the key gals f chemical reactin engineering is t quantify the relatinship between prductin rate, reactr size, reactin kinetics and selected perating cnditins.
More informationSUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis
SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical mdel fr micrarray data analysis David Rssell Department f Bistatistics M.D. Andersn Cancer Center, Hustn, TX 77030, USA rsselldavid@gmail.cm
More information37 Maxwell s Equations
37 Maxwell s quatins In this chapter, the plan is t summarize much f what we knw abut electricity and magnetism in a manner similar t the way in which James Clerk Maxwell summarized what was knwn abut
More informationSequential Allocation with Minimal Switching
In Cmputing Science and Statistics 28 (1996), pp. 567 572 Sequential Allcatin with Minimal Switching Quentin F. Stut 1 Janis Hardwick 1 EECS Dept., University f Michigan Statistics Dept., Purdue University
More informationLecture 23: Lattice Models of Materials; Modeling Polymer Solutions
Lecture 23: 12.05.05 Lattice Mdels f Materials; Mdeling Plymer Slutins Tday: LAST TIME...2 The Bltzmann Factr and Partitin Functin: systems at cnstant temperature...2 A better mdel: The Debye slid...3
More informationLecture 5: Equilibrium and Oscillations
Lecture 5: Equilibrium and Oscillatins Energy and Mtin Last time, we fund that fr a system with energy cnserved, v = ± E U m ( ) ( ) One result we see immediately is that there is n slutin fr velcity if
More information(2) Even if such a value of k was possible, the neutrons multiply
CHANGE OF REACTOR Nuclear Thery - Curse 227 POWER WTH REACTVTY CHANGE n this lessn, we will cnsider hw neutrn density, neutrn flux and reactr pwer change when the multiplicatin factr, k, r the reactivity,
More informationMATHEMATICS SYLLABUS SECONDARY 5th YEAR
Eurpean Schls Office f the Secretary-General Pedaggical Develpment Unit Ref. : 011-01-D-8-en- Orig. : EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 6 perid/week curse APPROVED BY THE JOINT TEACHING COMMITTEE
More informationA Matrix Representation of Panel Data
web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins
More informationLab #3: Pendulum Period and Proportionalities
Physics 144 Chwdary Hw Things Wrk Spring 2006 Name: Partners Name(s): Intrductin Lab #3: Pendulum Perid and Prprtinalities Smetimes, it is useful t knw the dependence f ne quantity n anther, like hw the
More informationThermodynamics and Equilibrium
Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,
More informationModule 4: General Formulation of Electric Circuit Theory
Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated
More informationCHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS
CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS 1 Influential bservatins are bservatins whse presence in the data can have a distrting effect n the parameter estimates and pssibly the entire analysis,
More informationDead-beat controller design
J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable
More informationLead/Lag Compensator Frequency Domain Properties and Design Methods
Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin
More information3.4 Shrinkage Methods Prostate Cancer Data Example (Continued) Ridge Regression
3.3.4 Prstate Cancer Data Example (Cntinued) 3.4 Shrinkage Methds 61 Table 3.3 shws the cefficients frm a number f different selectin and shrinkage methds. They are best-subset selectin using an all-subsets
More informationLyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
More informationBLAST / HIDDEN MARKOV MODELS
CS262 (Winter 2015) Lecture 5 (January 20) Scribe: Kat Gregry BLAST / HIDDEN MARKOV MODELS BLAST CONTINUED HEURISTIC LOCAL ALIGNMENT Use Cmmnly used t search vast bilgical databases (n the rder f terabases/tetrabases)
More informationWe can see from the graph above that the intersection is, i.e., [ ).
MTH 111 Cllege Algebra Lecture Ntes July 2, 2014 Functin Arithmetic: With nt t much difficulty, we ntice that inputs f functins are numbers, and utputs f functins are numbers. S whatever we can d with
More informationFebruary 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA
February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal
More informationDepartment of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets
Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0
More informationBASD HIGH SCHOOL FORMAL LAB REPORT
BASD HIGH SCHOOL FORMAL LAB REPORT *WARNING: After an explanatin f what t include in each sectin, there is an example f hw the sectin might lk using a sample experiment Keep in mind, the sample lab used
More informationENSC Discrete Time Systems. Project Outline. Semester
ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding
More informationTHE LIFE OF AN OBJECT IT SYSTEMS
THE LIFE OF AN OBJECT IT SYSTEMS Persns, bjects, r cncepts frm the real wrld, which we mdel as bjects in the IT system, have "lives". Actually, they have tw lives; the riginal in the real wrld has a life,
More informationChemistry 20 Lesson 11 Electronegativity, Polarity and Shapes
Chemistry 20 Lessn 11 Electrnegativity, Plarity and Shapes In ur previus wrk we learned why atms frm cvalent bnds and hw t draw the resulting rganizatin f atms. In this lessn we will learn (a) hw the cmbinatin
More informationMath Foundations 20 Work Plan
Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant
More informationCESAR Science Case The differential rotation of the Sun and its Chromosphere. Introduction. Material that is necessary during the laboratory
Teacher s guide CESAR Science Case The differential rtatin f the Sun and its Chrmsphere Material that is necessary during the labratry CESAR Astrnmical wrd list CESAR Bklet CESAR Frmula sheet CESAR Student
More informationLeast Squares Optimal Filtering with Multirate Observations
Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical
More information4th Indian Institute of Astrophysics - PennState Astrostatistics School July, 2013 Vainu Bappu Observatory, Kavalur. Correlation and Regression
4th Indian Institute f Astrphysics - PennState Astrstatistics Schl July, 2013 Vainu Bappu Observatry, Kavalur Crrelatin and Regressin Rahul Ry Indian Statistical Institute, Delhi. Crrelatin Cnsider a tw
More informationOn Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION
Malaysian Jurnal f Mathematical Sciences 4(): 7-4 () On Huntsberger Type Shrinkage Estimatr fr the Mean f Nrmal Distributin Department f Mathematical and Physical Sciences, University f Nizwa, Sultanate
More informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More informationSAMPLING DYNAMICAL SYSTEMS
SAMPLING DYNAMICAL SYSTEMS Melvin J. Hinich Applied Research Labratries The University f Texas at Austin Austin, TX 78713-8029, USA (512) 835-3278 (Vice) 835-3259 (Fax) hinich@mail.la.utexas.edu ABSTRACT
More informationTHERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES
PREFERRED RELIABILITY PAGE 1 OF 5 PRACTICES PRACTICE NO. PT-TE-1409 THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC Practice: Perfrm all thermal envirnmental tests n electrnic spaceflight hardware in a flight-like
More informationModelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA
Mdelling f Clck Behaviur Dn Percival Applied Physics Labratry University f Washingtn Seattle, Washingtn, USA verheads and paper fr talk available at http://faculty.washingtn.edu/dbp/talks.html 1 Overview
More informationDistributions, spatial statistics and a Bayesian perspective
Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics
More informationA New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation
III-l III. A New Evaluatin Measure J. Jiner and L. Werner Abstract The prblems f evaluatin and the needed criteria f evaluatin measures in the SMART system f infrmatin retrieval are reviewed and discussed.
More informationChapter 3 Kinematics in Two Dimensions; Vectors
Chapter 3 Kinematics in Tw Dimensins; Vectrs Vectrs and Scalars Additin f Vectrs Graphical Methds (One and Tw- Dimensin) Multiplicatin f a Vectr b a Scalar Subtractin f Vectrs Graphical Methds Adding Vectrs
More informationRevision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax
.7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical
More informationmaking triangle (ie same reference angle) ). This is a standard form that will allow us all to have the X= y=
Intrductin t Vectrs I 21 Intrductin t Vectrs I 22 I. Determine the hrizntal and vertical cmpnents f the resultant vectr by cunting n the grid. X= y= J. Draw a mangle with hrizntal and vertical cmpnents
More informationSections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.
Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage
More informationYou need to be able to define the following terms and answer basic questions about them:
CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f
More informationEric Klein and Ning Sa
Week 12. Statistical Appraches t Netwrks: p1 and p* Wasserman and Faust Chapter 15: Statistical Analysis f Single Relatinal Netwrks There are fur tasks in psitinal analysis: 1) Define Equivalence 2) Measure
More informationThis section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.
Sectin 3.2: Many f yu WILL need t watch the crrespnding vides fr this sectin n MyOpenMath! This sectin is primarily fcused n tls t aid us in finding rts/zers/ -intercepts f plynmials. Essentially, ur fcus
More informationCONSTRUCTING STATECHART DIAGRAMS
CONSTRUCTING STATECHART DIAGRAMS The fllwing checklist shws the necessary steps fr cnstructing the statechart diagrams f a class. Subsequently, we will explain the individual steps further. Checklist 4.6
More informationWriting Guidelines. (Updated: November 25, 2009) Forwards
Writing Guidelines (Updated: Nvember 25, 2009) Frwards I have fund in my review f the manuscripts frm ur students and research assciates, as well as thse submitted t varius jurnals by thers that the majr
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion
.54 Neutrn Interactins and Applicatins (Spring 004) Chapter (3//04) Neutrn Diffusin References -- J. R. Lamarsh, Intrductin t Nuclear Reactr Thery (Addisn-Wesley, Reading, 966) T study neutrn diffusin
More informationResampling Methods. Chapter 5. Chapter 5 1 / 52
Resampling Methds Chapter 5 Chapter 5 1 / 52 1 51 Validatin set apprach 2 52 Crss validatin 3 53 Btstrap Chapter 5 2 / 52 Abut Resampling An imprtant statistical tl Pretending the data as ppulatin and
More informationSynchronous Motor V-Curves
Synchrnus Mtr V-Curves 1 Synchrnus Mtr V-Curves Intrductin Synchrnus mtrs are used in applicatins such as textile mills where cnstant speed peratin is critical. Mst small synchrnus mtrs cntain squirrel
More informationNGSS High School Physics Domain Model
NGSS High Schl Physics Dmain Mdel Mtin and Stability: Frces and Interactins HS-PS2-1: Students will be able t analyze data t supprt the claim that Newtn s secnd law f mtin describes the mathematical relatinship
More informationAP Statistics Practice Test Unit Three Exploring Relationships Between Variables. Name Period Date
AP Statistics Practice Test Unit Three Explring Relatinships Between Variables Name Perid Date True r False: 1. Crrelatin and regressin require explanatry and respnse variables. 1. 2. Every least squares
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationLesson Plan. Recode: They will do a graphic organizer to sequence the steps of scientific method.
Lessn Plan Reach: Ask the students if they ever ppped a bag f micrwave ppcrn and nticed hw many kernels were unppped at the bttm f the bag which made yu wnder if ther brands pp better than the ne yu are
More informationNOTE ON A CASE-STUDY IN BOX-JENKINS SEASONAL FORECASTING OF TIME SERIES BY STEFFEN L. LAURITZEN TECHNICAL REPORT NO. 16 APRIL 1974
NTE N A CASE-STUDY IN B-JENKINS SEASNAL FRECASTING F TIME SERIES BY STEFFEN L. LAURITZEN TECHNICAL REPRT N. 16 APRIL 1974 PREPARED UNDER CNTRACT N00014-67-A-0112-0030 (NR-042-034) FR THE FFICE F NAVAL
More informationSection 6-2: Simplex Method: Maximization with Problem Constraints of the Form ~
Sectin 6-2: Simplex Methd: Maximizatin with Prblem Cnstraints f the Frm ~ Nte: This methd was develped by Gerge B. Dantzig in 1947 while n assignment t the U.S. Department f the Air Frce. Definitin: Standard
More informationGeneral Chemistry II, Unit II: Study Guide (part 1)
General Chemistry II, Unit II: Study Guide (part 1) CDS Chapter 21: Reactin Equilibrium in the Gas Phase General Chemistry II Unit II Part 1 1 Intrductin Sme chemical reactins have a significant amunt
More informationDepartment of Electrical Engineering, University of Waterloo. Introduction
Sectin 4: Sequential Circuits Majr Tpics Types f sequential circuits Flip-flps Analysis f clcked sequential circuits Mre and Mealy machines Design f clcked sequential circuits State transitin design methd
More informationPattern Recognition 2014 Support Vector Machines
Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft
More information3D KMC simulation on the precipitation in the annealed ternary alloy system
3D KM simulatin n the precipitatin in the annealed ternary ally system Xuan Zhang, Mengqi Huang Abstract Kinetic Mnte arl methd is used t study the precipitatin phenmenn in binary and ternary ally system,
More informationMODULE FOUR. This module addresses functions. SC Academic Elementary Algebra Standards:
MODULE FOUR This mdule addresses functins SC Academic Standards: EA-3.1 Classify a relatinship as being either a functin r nt a functin when given data as a table, set f rdered pairs, r graph. EA-3.2 Use
More informationInternal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.
Sectin 7 Mdel Assessment This sectin is based n Stck and Watsn s Chapter 9. Internal vs. external validity Internal validity refers t whether the analysis is valid fr the ppulatin and sample being studied.
More informationLecture 7: Damped and Driven Oscillations
Lecture 7: Damped and Driven Oscillatins Last time, we fund fr underdamped scillatrs: βt x t = e A1 + A csω1t + i A1 A sinω1t A 1 and A are cmplex numbers, but ur answer must be real Implies that A 1 and
More informationChapter 3: Cluster Analysis
Chapter 3: Cluster Analysis } 3.1 Basic Cncepts f Clustering 3.1.1 Cluster Analysis 3.1. Clustering Categries } 3. Partitining Methds 3..1 The principle 3.. K-Means Methd 3..3 K-Medids Methd 3..4 CLARA
More informationAircraft Performance - Drag
Aircraft Perfrmance - Drag Classificatin f Drag Ntes: Drag Frce and Drag Cefficient Drag is the enemy f flight and its cst. One f the primary functins f aerdynamicists and aircraft designers is t reduce
More informationB. Definition of an exponential
Expnents and Lgarithms Chapter IV - Expnents and Lgarithms A. Intrductin Starting with additin and defining the ntatins fr subtractin, multiplicatin and divisin, we discvered negative numbers and fractins.
More informationDavid HORN and Irit OPHER. School of Physics and Astronomy. Raymond and Beverly Sackler Faculty of Exact Sciences
Cmplex Dynamics f Neurnal Threshlds David HORN and Irit OPHER Schl f Physics and Astrnmy Raymnd and Beverly Sackler Faculty f Exact Sciences Tel Aviv University, Tel Aviv 69978, Israel hrn@neurn.tau.ac.il
More informationUNIV1"'RSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. CUMULATIVE SUM CONTROL CHARTS FOR THE FOLDED NORMAL DISTRIBUTION
UNIV1"'RSITY OF NORTH CAROLINA Department f Statistics Chapel Hill, N. C. CUMULATIVE SUM CONTROL CHARTS FOR THE FOLDED NORMAL DISTRIBUTION by N. L. Jlmsn December 1962 Grant N. AFOSR -62..148 Methds f
More informationAP Physics Kinematic Wrap Up
AP Physics Kinematic Wrap Up S what d yu need t knw abut this mtin in tw-dimensin stuff t get a gd scre n the ld AP Physics Test? First ff, here are the equatins that yu ll have t wrk with: v v at x x
More informationI. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is
Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,
More informationActivity Guide Loops and Random Numbers
Unit 3 Lessn 7 Name(s) Perid Date Activity Guide Lps and Randm Numbers CS Cntent Lps are a relatively straightfrward idea in prgramming - yu want a certain chunk f cde t run repeatedly - but it takes a
More informationPerfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart
Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Sandy D. Balkin Dennis K. J. Lin y Pennsylvania State University, University Park, PA 16802 Sandy Balkin is a graduate student
More informationNUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION
NUROP Chinese Pinyin T Chinese Character Cnversin NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION CHIA LI SHI 1 AND LUA KIM TENG 2 Schl f Cmputing, Natinal University f Singapre 3 Science
More informationComprehensive Exam Guidelines Department of Chemical and Biomolecular Engineering, Ohio University
Cmprehensive Exam Guidelines Department f Chemical and Bimlecular Engineering, Ohi University Purpse In the Cmprehensive Exam, the student prepares an ral and a written research prpsal. The Cmprehensive
More information