Blog Archives

Simulating species abundance data with coenocliner

July 31, 2014
By
Simulating species abundance data with coenocliner

Coenoclines are, according to the Oxford Dictionary of Ecology (Allaby 1998), “gradients of communities (e.g. in a transect from the summit to the base of a hill), reflecting the changing importance, frequency, or other appropriate measure of different species populations”. In much ecological research, and that of related fields, data on these coenoclines are collected and analyzed in a variety...

Read more »

Simultaneous confidence intervals for derivatives of splines in GAMs

June 16, 2014
By
Simultaneous confidence intervals for derivatives of splines in GAMs

Last time out I looked at one of the complications of time series modelling with smoothers; you have a non-linear trend which may be statistically significant but it may not be increasing or decreasing everywhere. How do we identify where in the series the data are changing? In that post I explained how we can use the first...

Read more »

Identifying periods of change in time series with GAMs

May 15, 2014
By
Identifying periods of change in time series with GAMs

In previous posts (here and here) I looked at how generalized additive models (GAMs) can be used to model non-linear trends in time series data. In my previous post I extended the modelling approach to deal with seasonal data where we model both the within year (seasonal) and between year (trend) variation with separate smooth functions....

Read more »

Modelling seasonal data with GAMs

May 9, 2014
By
Modelling seasonal data with GAMs

In previous posts I have looked at how generalized additive models (GAMs) can be used to model non-linear trends in time series data. At the time a number of readers commented that they were interested in modelling data that had more than just a trend component; how do you model data collected throughout the year over many years with...

Read more »

Summarising multivariate palaeoenvironmental data

January 9, 2014
By
Summarising multivariate palaeoenvironmental data

The horseshoe effect is a well known and discussed issue with principal component analysis (PCA) (e.g. Goodall 1954; Swan 1970; Noy-Meir & Austin 1970). Similar geometric artefacts also affect correspondence analysis (CA). In part 1 of this series I looked at the implications of these “artefacts” for the recovery of temporal or single dominant gradients from multivariate palaeoecological data....

Read more »

Decluttering ordination plots part 4: orditkplot()

December 31, 2013
By
Decluttering ordination plots part 4: orditkplot()

Earlier in this series I looked at the ordilabel() and then the orditorp() functions, and most recently the ordipointlabel() function in the vegan package as means to improve labelling in ordination plots. In this, the fourth and final post in the series I take a look at orditkplot(). If you’ve created ordination diagrams before or...

Read more »

Summarising multivariate palaeoenvironmental data

December 28, 2013
By
Summarising multivariate palaeoenvironmental data

Ordination methods that yield orthogonal axes of variation are often used to summarise the multivariate data obtained from sediment cores. Usually the first or, less often, the first few ordination axes are taken as directions of change or the main patterns of variance in the multivariate data. There is an oft-overlooked issue with this approach that has the potential...

Read more »

New version of permute on CRAN

December 17, 2013
By

After some time brewing on my machines, I’m happy to have released a new version of my permute package for R. This release took quite a while to polish and get right as there was a lot of back-and-forth between vegan and permute as I tried to get the latter working nicely for both useRs and developers, whilst Jari...

Read more »

New version of analogue on CRAN

December 14, 2013
By

It has been almost a year since the last release of the analogue package. At lot has happened in the intervening period and although I’ve been busy with a new job in a new country and coding on several other R packages, activity on analogue has also progressed a pace. As the version 0.12-0 of the package hits a...

Read more »

Time series plots in R

October 23, 2013
By
Time series plots in R

I recently coauthored a couple of papers on trends in environmental data (Curtis & Simpson in press; Monteith et al. in press), which we estimated using GAMs. Both papers included plots like the one shown below wherein we show the estimated trend and associated point-wise 95% confidence interval, plus some other markings. The coloured sections show...

Read more »