Exploring Pollen Data

April 2, 2012
By

(This article was first published on Noam Ross - R, and kindly contributed to R-bloggers)

I wrote a few functions to grab data and from the Global Pollen Database:

source("~/code/_Pollen/pollendatafuntions.R")
## Loading required package: stratigraph
## Loading required package: grid
## Loading required package: graphics
## Loading required package: stats
billys <- getpctAP("billys", plot = TRUE)
## Number of taxa:  105 
## Number of levels:  77 
Arboreal Pollen over time at Billy’s Lake

Arboreal Pollen over time at Billy’s Lake

Clearly, there’s are trends at different time scales. The trouble with using S-Maps on this is that there are uneven sampling intervals. One option is to use linear interpolation to approximate even intervals:

billyslin <- linpctAP(billys, ts = FALSE)
billys <- getpctAP("billys", plot = TRUE)
## Number of taxa:  105 
## Number of levels:  77 
lines(billyslin[, 1], billyslin[, 2], col = "red")
points(billyslin[, 1], billyslin[, 2], col = "red", 
    pch = 16)
With linearly interpolated evenly spaced samples

With linearly interpolated evenly spaced samples

This means introducing a new source of error into whatever model is ultimately used. What interval should be used to generate this? Here I created an equal number of points as are in the original data set. However, the sampling density is greater at more recent dates. S-Maps can handle missing data, but there need to be a good chunk of consecutive points.

The other option is start thinking about how to generalize S-Maps for data sampled at unequal intervals.

To leave a comment for the author, please follow the link and comment on his blog: Noam Ross - R.

R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...



If you got this far, why not subscribe for updates from the site? Choose your flavor: e-mail, twitter, RSS, or facebook...

Comments are closed.