R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Plot a SiZer map

plot.SiZer

R Documentation

Plot a SiZer map

Description

Plot a SiZer object that was created using SiZer().

Usage

## S3 method for class 'SiZer'.
## S3 method for class 'SiZer'
plot(x, ylab = expression(log[10](h)), 
	colorlist = c("red", "purple", "blue", "grey"), ...)

Arguments

`x`	An object created using `SiZer()`
`ylab`	What the y-axis should be labled.
`colorlist`	What colors should be used. This is a vector that corresponds to 'decreasing', 'possibley zero', 'increasing', 'insufficient data'
`...`	Any other parameters to be passed to the function `image`.

Details

The white lines in the SiZer map give a graphical representation of the bandwidth. The horizontal distance between the lines is 2h.

Value

None

Author(s)

Derek Sonderegger

References

Chaudhuri, P., and J. S. Marron. 1999. SiZer for exploration of structures in curves. Journal of the American Statistical Association 94:807-823.

Hannig, J., and J. S. Marron. 2006. Advanced distribution theory for SiZer. Journal of the American Statistical Association 101:484-499.

Sonderegger, D.L., Wang, H., Clements, W.H., and Noon, B.R. 2009. Using SiZer to detect thresholds in ecological data. Frontiers in Ecology and the Environment 7:190-195.

Examples

# data('Arkansas')
# x <- Arkansas$year
# y <- Arkansas$sqrt.mayflies

# Calculate the SiZer map for the first derivative
# SiZer.1 <- SiZer(x, y, h=c(.5,10), degree=1, derv=1)
# plot(SiZer.1)

# Calculate the SiZer map for the second derivative
# SiZer.2 <- SiZer(x, y, h=c(.5,10), degree=2, derv=2);
# plot(SiZer.2)

# By setting the grid.length larger, we get a more detailed SiZer
# map but it takes longer to compute. 
#
# SiZer.3 <- SiZer(x, y, h=c(.5,10), grid.length=100, degree=1, derv=1)
# plot(SiZer.3)

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(SiZer)
Loading required package: splines
Loading required package: boot
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/SiZer/plot.SiZer.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.SiZer
> ### Title: Plot a SiZer map
> ### Aliases: plot.SiZer
> ### Keywords: hplot smooth
> 
> ### ** Examples
> 
> # data('Arkansas')
> # x <- Arkansas$year
> # y <- Arkansas$sqrt.mayflies
> 
> # Calculate the SiZer map for the first derivative
> # SiZer.1 <- SiZer(x, y, h=c(.5,10), degree=1, derv=1)
> # plot(SiZer.1)
> 
> # Calculate the SiZer map for the second derivative
> # SiZer.2 <- SiZer(x, y, h=c(.5,10), degree=2, derv=2);
> # plot(SiZer.2)
> 
> # By setting the grid.length larger, we get a more detailed SiZer
> # map but it takes longer to compute. 
> #
> # SiZer.3 <- SiZer(x, y, h=c(.5,10), grid.length=100, degree=1, derv=1)
> # plot(SiZer.3)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>