tag:blogger.com,1999:blog-38987498074040053002015-01-10T11:31:03.521-08:00Fishy MathNews from the Math Bio stats group at NWFSC<br>
<a href="http://faculty.washington.edu/eeholmes">faculty.washington.edu/eeholmes</a>eehnoreply@blogger.comBlogger74125tag:blogger.com,1999:blog-3898749807404005300.post-23508438790265765832015-01-10T11:25:00.003-08:002015-01-10T11:31:03.540-08:00Winter 2015 Online Course: Applied Time Series Analysis in Fisheries and Environmental Sciences<div dir="ltr" style="text-align: left;" trbidi="on">Fish 507: Applied Time Series Analysis in Fisheries and Environmental Sciences<br />Winter 2015<br />Fisheries Dept, University of Washington<br /><br />Instructors: Eric Ward, Eli Holmes, Mark Scheuerell<br />email: eli.holmes@noaa.gov, mark.scheuerell@noaa.gov, eric.ward@noaa.gov<br /><br />Reviews current applications of univariate and multivariate time series models for biological and environmental data, emphasizing the estimation, inference, and forecasting aspects of time-series models. Explores effects of covariates and anthropogenic drivers for species that are exploited and/or of conservation concern. We taught a similar course 2 years ago. This time we are emphasizing how to fit these models in a Bayesian context with JAGS along with a MLE context with MARSS. We are recording the lectures and you can follow along with the course at: <a href="https://catalyst.uw.edu/workspace/fish203/35553/243766">https://catalyst.uw.edu/workspace/fish203/35553/243766</a></div> eehnoreply@blogger.comtag:blogger.com,1999:blog-3898749807404005300.post-73671301713154575052014-01-24T16:07:00.002-08:002014-01-24T16:07:19.662-08:00New paper on MAR modeling of community dynamics<div dir="ltr" style="text-align: left;" trbidi="on"><div><div><div class="arttitle articleTitle" style="text-align: left;"><b>Quantifying effects of abiotic and biotic drivers on community dynamics with multivariate autoregressive (MAR) models</b></div><div class="arttitle articleTitle" style="text-align: left;"><br /></div><div class="artAuthors"><strong> Stephanie E.<span class="NLM_x"> </span> Hampton<span class="NLM_x">, </span> Elizabeth E.<span class="NLM_x"> </span> Holmes<span class="NLM_x">, </span> Lindsay P.<span class="NLM_x"> </span> Scheef<span class="NLM_x">, </span> Mark D.<span class="NLM_x"> </span> Scheuerell<span class="NLM_x">, </span> Stephen L.<span class="NLM_x"> </span> Katz<span class="NLM_x">, </span> Daniel E.<span class="NLM_x"> </span> Pendleton<span class="NLM_x">, and </span> Eric J.<span class="NLM_x"> </span> Ward</strong><div class="first last"><br /></div><div class="first last">Long-term ecological data sets present opportunities for identifying drivers of community dynamics and quantifying their effects through time series analysis. Multivariate autoregressive (MAR) models are well known in many other disciplines, such as econometrics, but widespread adoption of MAR methods in ecology and natural resource management has been much slower despite some widely cited ecological examples. Here we review previous ecological applications of MAR models and highlight their ability to identify abiotic and biotic drivers of population dynamics, as well as community-level stability metrics, from long-term empirical observations. Thus far, MAR models have been used mainly with data from freshwater plankton communities; we examine the obstacles that may be hindering adoption in other systems and suggest practical modifications that will improve MAR models for broader application. Many of these modifications are already well known in other fields in which MAR models are common, although they are frequently described under different names. In an effort to make MAR models more accessible to ecologists, we include a worked example using recently developed R packages (MAR1 and MARSS), freely available and open-access software.</div></div><br /><span>Read More: <a href="http://www.esajournals.org/doi/abs/10.1890/13-0996.1">http://www.esajournals.org/doi/abs/10.1890/13-0996.1</a></span></div></div><br /><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-71185447530305278372014-01-16T10:21:00.003-08:002014-01-24T16:03:17.764-08:00Time Series (MARSS) Course offered in March in Stockholm<div dir="ltr" style="text-align: left;" trbidi="on">Mark and Eli are teaching a week-long multivariate time-series analysis course in Stockholm in March.<a href="http://timeseriescourseemb.wordpress.com/"> Course Announcement</a></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-73328103546352310182013-11-29T11:00:00.003-08:002013-11-29T11:00:51.650-08:00MARSS 3.6 up on CRAN. Significant speed increases for large models<div dir="ltr" style="text-align: left;" trbidi="on">MARSS 3.6 has been uploaded to CRAN. I fixed some inefficiencies that were causing DFA models with many time-series (n>100) and R="diagonal and unequal" to be very, very slow. My tests show 10x faster fits for n=100 and R="diagonal and unequal" for DFA models.<br /><br /><a href="http://cran.r-project.org/web/packages/MARSS/index.html" target="_blank">http://cran.r-project.org/web/<wbr></wbr>packages/MARSS/index.html</a></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-52341954040435242052013-10-08T16:39:00.000-07:002013-10-08T16:39:31.766-07:00New paper out by Jim and Eric on using Delta-GLMMs to analyze fisheries survey data<div dir="ltr" style="text-align: left;" trbidi="on">Thorson, J.T. and E.J. Ward. 2013. Accounting for space-time interactions in index<br />standardization models. Fisheries Research,147:426:433<br /><br />Scientific survey data are used to estimate abundance trends for fish populations worldwide, and are frequently analyzed using delta-generalized linear mixed models (delta-GLMMs). Delta-GLMMs incorporate information about both the probability of catch being non-zero (catch probability) and the expected value for non-zero catches (catch rates). Delta-GLMMs generally incorporate year as a main effect, and frequently account for spatial strata and/or covariates. Many existing delta-GLMMs do not account for random or systematic differences in catch probability or rates in particular combinations of spatial strata and year (i.e., space–time interactions), and do not recognize potential correlation in random space–time interactions between catch probability and catch rates. We therefore develop a Bayesian delta-GLMM that estimates correlations between catch probability and rates, and compare it with either (a) ignoring year–strata interactions, (b) modeling year–strata interactions as fixed effects, or (c) estimating year–strata interactions in catch probability or rates as independent random effects. These four models are fitted to bottom trawl survey data for 28 species off the U.S. West Coast. The posterior median of the correlation is positive for the majority (18) of species, including all five for which the posterior distribution has little overlap with zero. However, estimating this correlation has little impact on resulting abundance indices or credible intervals. We therefore conclude that the correlated random model will have a little impact on index standardization of the West Coast bottom trawl dataset. However, we propose that the correlated model can quickly identify correlations between occupancy probability and density, and provide our code to allow researchers to quickly identify whether such a correlation is likely to be significantly different from zero for their chosen data set.</div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-25994190320619109252013-07-05T15:51:00.001-07:002013-07-05T15:51:26.630-07:00Building R packages with RStudio and embedding R in your documents and reports<div dir="ltr" style="text-align: left;" trbidi="on">Building R packages with RStudio plus Embedding R in documents short-course on-line:<br /><br /><a href="http://www.iugo-cafe.org/chinook/view_node.php?id=2962" target="_blank">http://www.iugo-cafe.org/<wbr></wbr>chinook/view_node.php?id=2962</a><br /><br />Topics:<br /><ul><li>how (and why) to make an R package using RStudio </li><li>installing packages from github, git or a url to your tar.gz file</li><li>Using Sweave and RStudio to do 'reproducible research/programming'.<br /> </li><li>Using OpenOffice + R to do the same, if you don't like LaTeX</li><li>Creating web-apps that run your R code (a few links to demos)</li></ul><br /><div class=""><div class="" data-tooltip="Show trimmed content" id=":y3" role="button" tabindex="0"><img class="" src="https://mail.google.com/mail/u/1/images/cleardot.gif" /></div></div></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-35166642423418452322013-03-12T15:54:00.000-07:002013-03-12T15:54:06.060-07:00Week 10: Applied Time-Series Analysis for Fisheries and Environmental Data <div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br><br><b>Week 10: Dynamic linear models</b><br>This week, we give a brief introduction to an important class of MARSS models: dynamic linear models. These are multivariate linear regression models where the regression parameters (slope and intercepts) are treated as a AR process and thus are allowed to time evolve. We also review some of the diagnostics for MARSS models fits. <br><br><b>Lab 10</b> In the lab, you'll go through an simple example of a univariate dynamic linear model with time-varying slope and intercept. <br><b>Lecture 10</b> You can find the <a href="https://catalyst.uw.edu/workspace/file/download/3f148c16596400a18fcef983f939782221823a3e714eec63a0b70321a36efa3b">pdf of lecture</a> on the class webpage along with the link to watch a recording of the lecture. <hr></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-49288340866617691632013-03-05T14:14:00.001-08:002013-03-12T13:39:42.300-07:00Week 9: Applied Time-Series Analysis for Fisheries and Environmental Data <div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br><br><b>Week 9: Bayesian hierarchical multivariate state-space models</b><br>This week, we discuss fitting non-linear MARSS models and MARSS models with non-Gaussian errors using Bayesian methods. This is a very brief introduction and many shows some examples of how one sets up a MARSS model in JAGS and shows you what the posteriors of some models look like. <br>* Posteriors for MARSS models <br>* Intro to JAGS (the Gibbs sampler we will be using) <br>* many examples of fitting non-linear and non-Gaussian MARSS models <br><br><b>Lab 9</b> The main lab is to go through the JAGS examples shown in the lecture. <br><b>Lecture 9</b> You can find the <a href="https://catalyst.uw.edu/workspace/file/download/3f148c16596400a18fcef983f939782221823a3e714eec63a0b70321a36efa3b">pdf of lecture</a> on the class webpage along with the link to watch <a href="https://tegr.it/y/1168e">a recording of the lecture</a>. <hr></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-71949994194583205302013-02-26T18:41:00.001-08:002013-03-07T11:14:13.763-08:00Week 8: Applied Time-Series Analysis for Fisheries and Environmental Data <div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br /><b>Week 8: Estimating interactions (the B matrix)</b><br />This week, we discuss issues related to the estimation of the B matrix in the context of using it to represent species interactions in a community dynamics models. <br />* univariate discrete time Gompertz model <br />* multivariate discrete time Gompertz model <br />* including covariates <br />* spurious density dependence resulting from ignoring observation error <br />* uncertainty in B elements resulting from estimating observation variance <br />* different methods for estimating confidence intervals: bootstrapping, hessian approximation, profile likelihood <br />* diagnostics <br /><br /><b>Lab 8</b><br />The main lab is to go through case study 7 in the MARSS User Guide and the corresponding code.<br /><b>Lecture 8</b><br />You can find the pdf of lecture on the class webpage.<br /><script src="https://tegr.it/y/10mst" type="text/javascript"></script> </div><br /><hr>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-39164831709292775532013-02-26T12:31:00.000-08:002013-03-07T11:10:53.565-08:00Week 7: Applied Time-Series Analysis for Fisheries and Environmental Data <div dir="ltr" style="text-align: left;" trbidi="on"><br />Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br /><b>Week 7: Dynamic factor analysis</b><br /><br />This week, we use MARSS to do dynamic factor analysis (DFA), which allows us to look for a set of common underlying trends among a relatively large set of time series (Harvey, 1989, sec. 8.5). This is conceptually different than what we have been doing in the previous weeks. Here we are trying to explain temporal variation in a set of n observed time series using linear combinations of a set of m hidden random walks, where m << n. You can think of this as PCA for time-series data. Zuur et al. (2003) show a number of examples of DFA applied to catch data and densities of zoobenthos. <br><br /><b>Lab 7</b><br />The main lab is to go through the dynamic factor analysis chapter in the MARSS User Guide and the corresponding code.<br /><br /><b>Lecture 7</b><br />You can find the ppt of lecture on the class webpage. Technical difficulties prevented recording of the lecture.<br /><hr></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-74868252892087640572013-02-13T12:03:00.005-08:002013-03-07T11:16:23.722-08:00Week 6: Applied Time-Series Analysis for Fisheries and Environmental Data <br />Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br /><b>Week 6: Introduction to including covariates in multivariate time-series model</b><br />This week we introduce the inclusion of covariates using the framework of a multivariate autoregressive model written in state-space form. You will understand the lecture better if you read the chapter on covariates in the MARSS User Guide first. Much of the lecture is about how to include covariates in different mathematically equivalent ways. You'll want to translate the R code in the lecture into the mathematical formulas (matrix form) to see how covariates are entering the mathematical model.<br /><br /><br />Lab topic:<br />The main lab is to go through the covariate chapter and examples in the MARSS User Guide. Then we have some salmon data to play with to try different ways of including cycles (in this case driven by cohort strength) into an analysis.<br /><br /><b>Lecture 6</b><br />Click the big arrow to start. You can also find the ppt of lecture t on the class webpage.<br /><script src="https://tegr.it/y/zjha" type="text/javascript"></script><br /><br /><hr>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-41569961876689476382013-02-06T15:17:00.003-08:002013-02-06T15:55:28.649-08:00Week 5: Applied Time-Series Analysis for Fisheries and Environmental Data<div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br />Week 5: Introduction to multivariate autoregressive state-space models<br />Lecture topics:<br /><ul style="text-align: left;"><li>Review of dealing with obs error with ARIMA (from last week)</li><li>Multivariate state space models</li><li>How these are expressed mathematically</li><li>Analysis of multi-site data using this framework</li><li>Parameter estimation: Kalman filter, Newton methods and EM algorithm</li></ul>Lab topic:<br />The main lab is to go through case study 2 in the MARSS User Guide. I have a Fish 507 specific version of the case study code on the course website with questions to answer as you go through. Case study 3 and 8 are optional but going through them will help solidify your understanding of multivariate state-space models. Do go through the ARMA code as it discusses some important points about the effects of data transformation (in this case differences) on the time-series model that is appropriate for the data.<br /><br /><b>Lecture 5</b><br />Click the big arrow to start the show. You can also find just a pdf of lecture 5 on the class webpage.<br /><script type="text/javascript" src="https://tegr.it/y/yxu5"></script><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-4702373380546432852013-01-31T10:33:00.002-08:002013-01-31T10:33:35.365-08:00Week 4: Applied Time-Series Analysis for Fisheries and Environmental Data<div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br />Week 4: Introduction to univariate autoregressive state-space models<br />Topics:<br /><ul style="text-align: left;"><li>State-space models</li><li>Process versus observation error</li><li>Model Selection</li></ul><br><br /><b>Lecture 4</b><br />Click the big arrow to start the show. You can also find just the ppt of lecture 3 on the class webpage.<br /><script type="text/javascript" src="https://tegr.it/y/ybn5"></script><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-41220603213712438252013-01-22T17:59:00.003-08:002013-01-22T18:00:42.133-08:00Week 3: Applied Time-Series Analysis for Fisheries and Environmental Data<div dir="ltr" style="text-align: left;" trbidi="on"><div dir="ltr" style="text-align: left;" trbidi="on"><div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br />Week 3: Estimation, model selection, and forecasting for time series models<br />Topics:<br /><ul style="text-align: left;"><li>Summarizing ARIMA models</li><li>Estimation</li><li>Model Selection</li><li>Prediction & forecasting</li><li>Evaluating forecasts</li><li>Functions: arima(), lm(), Arima()</li></ul></div><b>Lecture 3</b><br />This is our second attempt at recording a lecture. Still much to be learned but we are getting better. Click the big arrow to start the show. You can also find just the ppt of lecture 3 on the class webpage.</div><script src="https://tegr.it/y/xfd7" type="text/javascript"></script><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-7957586085493492072013-01-15T15:27:00.002-08:002013-01-15T15:41:29.482-08:00Week 2: Applied Time-Series Analysis for Fisheries and Environmental Data<div dir="ltr" style="text-align: left;" trbidi="on">Class material: <a href="https://catalyst.uw.edu/workspace/fish203/35553/">webpage</a><br /><br />Week 2: Correlation, stationarity & stationary time-series models<br />The lecture introduces the ACF, PACF, and basic properties of AR, MA and ARMA models. The computer code section shows <a href="https://catalyst.uw.edu/workspace/file/download/983a9f707834943030ea49bc7366ea8178173f5f93ecb506bafd05474b5742df">R code</a> to analyze simulated time-series data so that participants get a feel for ACF and PACF and get a feel for AR and MA processes. The participants then move to analyzing some real time-series data using the <a href="https://catalyst.uw.edu/workspace/file/download/983a9f707834943030ea49bc7366ea81ebc5cbb4804b6cf9f657346ddcb35cdd">30+ year time-series of Lake Washington plankton</a>.<br /><br /><b>Lecture 2</b><br />This is our first attempt at recording a lecture. Ahem, there is clearly much to be learned to improve the process...Click the big arrow to start the show. You can also find just the ppt of <a href="https://catalyst.uw.edu/workspace/file/download/983a9f707834943030ea49bc7366ea81852fc591336a4267f4fa5b7b5b156182">lecture 2</a> on the class webpage.<br /><script type="text/javascript" src="https://tegr.it/y/wotn"></script><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-35488935990423956042012-11-30T13:21:00.001-08:002013-01-22T17:52:54.186-08:00Winter stats reading group starting up: Hierarchical Modeling and Analysis for Spatial Data<div dir="ltr" style="text-align: left;" trbidi="on">The <a href="http://faculty.washington.edu/eeholmes/stats_reading_group.shtml">NWFSC/SAFS stats reading group</a> is reading "Hierarchical Modeling and Analysis for Spatial Data" by Banerjee et al. this quarter. Fridays 3pm at SAFS 229 during Winter Qtr 2013. Open to interested statistical ecologists. Contact Eli.</div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-6039309235244936412012-11-12T16:49:00.003-08:002012-11-30T13:23:16.869-08:00New paper on spatial-temporal time series modeling<div dir="ltr" style="text-align: left;" trbidi="on">New paper just out by Eric Ward using Bayesian state-space time-series models.<br /><br /><i>"Applying time series models with spatial correlation to identify the scale of variation in habitat metrics related to threatened coho salmon (Oncorhynchus kisutch) in the Pacific Northwest"</i><br />Eric J. Ward, George R. Pess, Kara Anlauf-Dunn, and Chris E. Jordan<br />Canadian Journal of Fisheries and Aquatic Science (<a href="http://www.nrcresearchpress.com/doi/pdf/10.1139/f2012-096">link to paper</a>)<br /><br />Abstract: Trend analyses are common in the analysis of fisheries data, yet the majority of them ignore either observation error or spatial correlation. In this analysis, we applied a novel hierarchical Bayesian state-space time series model with spatial correlation to a 12-year data set of habitat variables related to coho salmon (Oncorhynchus kisutch) in coastal Oregon, USA. This model allowed us to estimate the degree of spatial correlation separately for each habitat variable and the importance of observation error relative to environmental stochasticity. This framework allows us to identify variables that would benefit from additional sampling and variables where sampling could be reduced. Of the eight variables included in our analysis, we found three metrics related to habitat quality correlated at large spatial scales (gradient, fine sediment, shade cover). Variables with higher observation error (pools, active channel width, fine sediment) could be made more precise with more repeat visits. Our spatio-temporal model is flexible and extendable to virtually any spatially explicit monitoring data set, even with large amounts of missing data and no repeated observations. Potential extensions include fisheries catch data, abiotic indicators, invasive species, or species of conservation concern.</div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-74223272357853329272012-08-01T15:16:00.003-07:002012-11-30T13:23:43.323-08:00Time-series analysis course winter 2012<div dir="ltr" style="text-align: left;" trbidi="on">Fish 50X: Applied Time Series Analysis in Fisheries and Environmental Sciences<br />Winter 2012<br />Fisheries Dept, University of Washington<br /><br />Instructors: Eric Ward, Eli Holmes, Mark Scheuerell<br />email: eli.holmes@noaa.gov, mark.scheuerell@noaa.gov, eric.ward@noaa.gov<br /><br />Reviews current applications of univariate and multivariate time series models for biological and environmental data, emphasizing the estimation, inference, and forecasting aspects of time-series models. Explores effects of covariates and anthropogenic drivers for species that are exploited and/or of conservation concern. Recommended: FISH 552 or prior experience with R (e.g. FISH 560), QSCI 482 or basic statistics, and at least 1 course in population dynamics (FISH 454 or 458). </div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-21151232201374628912012-08-01T15:12:00.004-07:002012-08-01T15:12:49.401-07:00R Journal article on the MARSS package<div dir="ltr" style="text-align: left;" trbidi="on">Holmes, E. E., Ward, E. J. and K. Wills. 2012. MARSS: Multivariate autoregressive state-space models for analyzing time-series data. R Journal 4: 11-19. <a href="http://journal.r-project.org/archive/2012-1/RJournal_2012-1_Holmes%7Eet%7Eal.pdf">http://journal.r-project.org/archive/2012-1/RJournal_2012-1_Holmes~et~al.pdf</a><br /></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-31712477325373675992012-08-01T15:11:00.002-07:002012-08-01T15:17:47.687-07:00MARSS 3.1 released on CRAN<div dir="ltr" style="text-align: left;" trbidi="on">MARSS 3.1 is now up on CRAN. This allows for time-varying constraints and covariates. See the updated User Guide (on CRAN). The major changes are internal and allow for us easily write customized functions for different MARSS forms (like AR-p processes and DFA models). 3.1 is considerably slower than 2.x, however this should be fixed in 3.2 or 3.3 when the Kalman filter in the KFAS package can be hooked back up to MARSS (temporarily disabled).<a href="http://cran.r-project.org/web/packages/MARSS"> http://cran.r-project.org/web/packages/MARSS</a></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-22350226607913885402012-06-23T17:36:00.000-07:002012-06-23T17:36:08.853-07:00MARSS 3.0 posted for testing<div dir="ltr" style="text-align: left;" trbidi="on">Hi MARSS users,<br /><br />A version of MARSS 3.0 is now up for testing. It should be backwards compatible with any MARSS 2.x code you have unless you use control$diffuse or control$kf.x0. diffuse now goes into your model list and kf.x0 is called tinitx and goes in the model list too. control$kf.x0="x00" is now model$tinitx=0 and control$kf.x0="x10" is model$tinitx=1.<br /><br />Big changes are<br /><ul style="text-align: left;"><li>Time-varying parameters are allowed. See the Quick_Start.pdf to get a brief intro to that feature but it should be pretty self-explanatory.</li><li>Covariates can be added in the standard way. Again see the Quick_Start.pdf for quick intro. See chapter in User Guide on estimating species interactions for an example.</li><li>There is a "form" argument in the MARSS() call that allows one to specify special types of models. Default is "marxss" which covers MARSS + covariates. The only other form now is "dfa" for Dynamic Factor Analysis. Check out the DFA chapter in the User Guide for an intro to the form="dfa" which allows you to do a standard DFA by just passing in m (number of states), data and covariates (if wanted). The dfa form is basic now. Later we will specialize its output to give loadings etc.</li><li>The AR-p models work now with method="kem" which is much, much faster than method="BFGS". See chapter in User Guide on AR-p models.</li><li>You can enter things like B=diag(list("1+2*c+3*b",0,0,"<wbr></wbr>2+3c"),2,2) in your list matrices and MARSS will know what to do, i.e. it will estimate B.c and B.b and it knows that B(1,1)=1+2c+3b . I haven't seen people want to do this..., but you can.</li><li>The print call takes a argument called "what". Use ?print.marssMLE to see how to use it. It'll make it easier to print things from your marssMLE objects (what you get back from a MARSS call). </li></ul><div id=":1bv" style="text-align: left;"> Here is are the tar.gz and .zip files for the 3.0 version. You'll find links to the User Guide and Quick_Start guide here too.<br /><br /><a href="http://fishbox.iugo-cafe.org/user/e2holmes/MARSS%20Dev%20Site" target="_blank">http://fishbox.iugo-cafe.org/<wbr></wbr>user/e2holmes/MARSS%20Dev%<wbr></wbr>20Site</a><br /> </div><div id=":1bv" style="text-align: left;">Plan is to upload to CRAN about July 1 assuming no big issues arise. Right now all the prior examples in the User Guide 2.8 and the man files work as before. I've included new examples with covariates using the new covariate code in the chapter on estimating spp interactions and I've added covariates to the DFA chapter. I've included a little more code in the AR-1 chapter on estimating those models.<br /><br />Feel free to try it out. The more real-world testing it gets before being uploaded to CRAN the better.<br /><br />Cheers,<br /><br />Eli</div></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-56268482925244507572012-03-22T17:58:00.000-07:002012-03-22T17:58:10.081-07:00Article from the group on analysis of marine plankton community structure<div dir="ltr" style="text-align: left;" trbidi="on"><a href="http://www.aslo.org/lomethods/free/2012/0054.html">Scheef, L.P., D.E. Pendleton, S.E. Hampton, S.L. Katz, E.E. Holmes, M.D. Scheuerell, and D.G. Johns. 2012. Assessing marine plankton community structure from long-term monitoring data with multivariate autoregressive (MAR) models: a comparison of fixed station vs. spatiallydistributed sampling data. Limnology & Oceanography: Methods 10: 54-64.</a><br /><br />ABSTRACT: We examined how marine plankton interaction networks, as inferred by multivariate autoregressive (MAR) analysis of time-series, differ based on data collected at a fixed sampling location (L4 station in the Western English Channel) and four similar time-series prepared by averaging Continuous Plankton Recorder (CPR) datapoints in the region surrounding the fixed station. None of the plankton community structures suggested by the MAR models generated from the CPR datasets were well correlated with the MAR model for L4, but of the four CPR models, the one most closely resembling the L4 model was that for the CPR region nearest to L4. We infer that observation error and spatial variation in plankton community dynamics influenced the model performance for the CPR datasets. A modified MAR framework in which observation error and spatial variation are explicitly incorporated could allow the analysis to better handle the diverse time-series data collected in marine environments. </div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-62002486069905837992012-03-22T17:43:00.001-07:002012-03-22T18:02:03.746-07:00Time-series analysis workshop Sat Aug 5, Portland, OR<div dir="ltr" style="text-align: left;" trbidi="on">Eli, Eric and Mark will offer their 1-day workshop on multivariate time-series analysis using the MARSS package again at the annual ESA meeting (<a href="http://www.esa.org/portland">http://www.esa.org/portland</a>). Workshop is scheduled for the Sat before the meeting. We are working on new case studies involving incorporation of covariates into analyses. MARSS 3.0 will be done by then. This is a major update that allows all parameters to incorporate time-varying covariates.</div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-1613182312261216192011-11-09T12:45:00.000-08:002011-11-09T12:45:28.095-08:00New article by Eric in Conservation Letters<div dir="ltr" style="text-align: left;" trbidi="on"><b>Integrating diet and movement data to identify hot spots of predation risk and areas of conservation concern for endangered species</b>, Eric J. Ward, Phillip S. Levin, Monique M. Lance, Steven J. Jeffries, Alejandro Acevedo-Gutiérrez<br /><i>Effective management of threatened and endangered species requires an understanding of how species of conservation concern are distributed spatially, as well as the spatial distribution of risks to the population, such as predation or human impacts (fishing, pollution, loss of habitat). Identifying high risk areas is particularly important when designing reserves or protected areas. Our novel approach incorporates data on distribution, movement, and diet of a generalist marine predator (harbor seals) to identify and map ‘hot-spots’ of predation risk for an endangered prey species (rockfish). Areas with high concentrations of seals (including some current marine reserves) are also estimated hot spots for rockfish predation. While marine reserve planning currently targets areas with good habitat and low human disturbance, our modeling suggests that future terrestrial and marine reserve design may be made more effective by incorporating other components of the food web that either directly or indirectly interact with target species.</i><br /> <br /><a href="http://onlinelibrary.wiley.com/doi/10.1111/j.1755-263X.2011.00210.x/abstract;jsessionid=7EB218B3908AAED4EC8A2E308F6C9046.d01t02">http://onlinelibrary.wiley.com/doi/10.1111/j.1755-263X.2011.00210.x/abstract;jsessionid=7EB218B3908AAED4EC8A2E308F6C9046.d01t02</a></div>eehnoreply@blogger.com0tag:blogger.com,1999:blog-3898749807404005300.post-71972768782945393132011-08-11T10:42:00.000-07:002011-08-11T10:42:04.647-07:00Post-doc opening in our research group<div dir="ltr" style="text-align: left;" trbidi="on">
<i>Time-series modeling of large-scale population and community processes</i><br />
Northwest Fisheries Science Center, NOAA Fisheries, Seattle, WA<br />
<br />
NOAA’s Northwest Fisheries Science Center (NWFSC) has a large research group using time-series modeling to study ecological dynamics. Our research interests are diverse, including estimating species-interaction strengths, inferring environmental and anthropogenic drivers of population and community dynamics, estimating stability metrics, detecting change points and regime shifts. We have numerous long-term and large-scale time-series data on fish, marine mammals, and plankton, which we use to study a wide array of basic and applied ecological and fishery questions.<br />
<br />
We are seeking a post-doctoral scientist (recent or fairly recent grad) to join our research group. Strong quantitative skills along with a background in aquatic or fisheries ecology are necessary for this position including a record of publication. The ideal candidate would also have experience and interest in one or more of the following areas:<br />
· statistical modeling, esp. hierarchical modeling<br />
· time-series analysis<br />
· population and/or community dynamics<br />
· fisheries management<br />
· salmon biology<br />
Post-doctoral positions are initially supported for 1 year with extensions up to 3 years contingent on satisfactory progress and submitted publications.<br />
<br />
<b>Why come post-doc at the NWFSC?</b> You will join a supportive, collaborative and productive team of quantitative ecologists at NWFSC who are using time-series modeling to study ecological dynamics. Our center of 300+ research scientists has a large number of post-doctoral fellows and provides a stimulating and productive environment for research. Post-docs trained in our group obtain a strong grounding in modern ecological statistics and have high success obtaining positions at both federal research and academic institutions. Close proximity to the University of Washington (a 10min walk) facilitates on-going collaborations with faculty and post-docs in the UW School for Fishery and Aquatic Sciences and other departments across campus. <br />
<br />
<b>Interested?</b> Contact one of the PIs below to discuss the position in more detail. Please attach a CV, recent publications, , and a brief statement describing your background, including any programming and modeling expertise. If you will be at the AFS meeting, we'd love to talk with you. Review of candidates will start in earnest after AFS in mid-September.<br />
<br />
PIs on this project are:<br />
<br />
Eli Holmes eli.holmes@noaa.gov <a href="http://faculty.washington.edu/eeholmes/">http://faculty.washington.edu/eeholmes/</a><br />
Mark Scheuerell mark.scheuerell@noaa.gov <a href="http://faculty.washington.edu/scheuerl/">http://faculty.washington.edu/scheuerl/</a><br />
Eric Ward eric.ward@noaa.gov <a href="http://sites.google.com/site/ericward2/">http://sites.google.com/site/ericward2/</a></div>
eehnoreply@blogger.com0