Clam Proxy Data

Description

Proxy data from the Orknay Islands

Usage

data(CLAM1)

Format

The format is: List of 2 $ x: num [1:91] 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 ... $ y: num [1:91] 2.14 2.11 2.4 2.21 2.32 2.44 2.85 2.78 2.27 2.05 ...

Details

Duplicate data has been removed.

Source

One-year data interval (4.68-9.31 mm) of delta 18O record of an archaeological limpet Patella vulgata (specimen QG2-1064-1) from Orkney, Scotland (Surge and Barrett, 2012).

References

Wang, T., Surge, D., and Lees, J. M., (2015) ClamR: A Statistical Evaluation of Isotopic and Temperature Records in Sclerochronologic Studies. Palaeogeography, Palaeoclimatology, Palaeoecology, doi:10.1016/j.palaeo.2015.07.008.

Examples

data(CLAM1)
## maybe str(CLAM1) ; plot(CLAM1) ...
plot(CLAM1$x, CLAM1$y, type="b", xlab="Distance", ylab="d18O" )

## Not run: 
###### this is an example from Wang et al.:
####  it takes too long to run on CRAN, but should work

shellx=CLAM1$x[38:70]
shelly=CLAM1$y[38:70]

window_shell=windowsize(shellx,shelly,1.8,9.4,0.2)

#the window size is 5mm, and make all the plots together
gout_shell = proxyJK(shellx, shelly, 5)

par(mfrow=c(3,2))
plot(shellx,shelly,type="b", xlab="Distance from Margin (mm)",
     ylab=expression(delta*"18O(ppm VPDB)"),  
     xlim=c(4,10), ylim = c(1.5,4))
plot((window_shell$win)/2,window_shell$error,xlab="Window Size (mm)", 
     ylab="Error",  xlim=c(1.6/2,9.4/2), ylim=c(0,0.5))
abline(v=4.63/2, lty="dotdash",col="black")
abline(v=5/2, col="black")
plotproxy1(shellx, shelly, gout_shell, xlim=c(4,10), ylim = c(1.5,4),
          xlab="Distance from Margin (mm)",
          ylab=expression(delta*"18O(ppm VPDB)"), main="")
plotproxy.all(gout_shell,YAXstyle=1, xlim=c(4,10), ylim1=c(0,4), 
             ylim2=c(-15,5))
plotproxy.error(shellx, shelly, gout_shell, type = 1, xlim=c(4,10), 
               ylim = c(1.5,4), xlab="Distance from Margin (mm)",
               ylab=expression(delta*"18O(ppm VPDB)"))
plotproxy.error(shellx, shelly, gout_shell, type = 2, xlim=c(4,10),
               ylim = c(1.5,4), xlab="Distance from Margin (mm)",
               ylab=expression(delta*"18O(ppm VPDB)"))




## End(Not run)

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(ClamR)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ClamR/CLAM1.Rd_%03d_medium.png", width=480, height=480)
> ### Name: CLAM1
> ### Title: Clam Proxy Data
> ### Aliases: CLAM1
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(CLAM1)
> ## maybe str(CLAM1) ; plot(CLAM1) ...
> plot(CLAM1$x, CLAM1$y, type="b", xlab="Distance", ylab="d18O" )
> 
> ## Not run: 
> ##D ###### this is an example from Wang et al.:
> ##D ####  it takes too long to run on CRAN, but should work
> ##D 
> ##D shellx=CLAM1$x[38:70]
> ##D shelly=CLAM1$y[38:70]
> ##D 
> ##D window_shell=windowsize(shellx,shelly,1.8,9.4,0.2)
> ##D 
> ##D #the window size is 5mm, and make all the plots together
> ##D gout_shell = proxyJK(shellx, shelly, 5)
> ##D 
> ##D par(mfrow=c(3,2))
> ##D plot(shellx,shelly,type="b", xlab="Distance from Margin (mm)",
> ##D      ylab=expression(delta*"18O(ppm VPDB)"),  
> ##D      xlim=c(4,10), ylim = c(1.5,4))
> ##D plot((window_shell$win)/2,window_shell$error,xlab="Window Size (mm)", 
> ##D      ylab="Error",  xlim=c(1.6/2,9.4/2), ylim=c(0,0.5))
> ##D abline(v=4.63/2, lty="dotdash",col="black")
> ##D abline(v=5/2, col="black")
> ##D plotproxy1(shellx, shelly, gout_shell, xlim=c(4,10), ylim = c(1.5,4),
> ##D           xlab="Distance from Margin (mm)",
> ##D           ylab=expression(delta*"18O(ppm VPDB)"), main="")
> ##D plotproxy.all(gout_shell,YAXstyle=1, xlim=c(4,10), ylim1=c(0,4), 
> ##D              ylim2=c(-15,5))
> ##D plotproxy.error(shellx, shelly, gout_shell, type = 1, xlim=c(4,10), 
> ##D                ylim = c(1.5,4), xlab="Distance from Margin (mm)",
> ##D                ylab=expression(delta*"18O(ppm VPDB)"))
> ##D plotproxy.error(shellx, shelly, gout_shell, type = 2, xlim=c(4,10),
> ##D                ylim = c(1.5,4), xlab="Distance from Margin (mm)",
> ##D                ylab=expression(delta*"18O(ppm VPDB)"))
> ##D 
> ##D 
> ##D 
> ## End(Not run)
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>