Derivative to be plotted. Default is 0 - which plots the smoothed function.
CI.method
What method should be used to calculate the confidence interval about the estimated line. The
methods are from Hannig and Marron (2006), where 1 is the point-wise estimate, and 2 is the row-wise estimate.
alpha
The CI has a 1-alpha/2 level of significance.
use.ess
ESS stands for the estimated sample size. If at any point along the x-axis, the ESS is too small, then
we will not plot unless use.ess=FALSE.
draw.points
Should the data points be included in the graph?
...
Additional arguments to be passed to the graphing functions.
Author(s)
Derek Sonderegger
References
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.
See Also
locally.weighted.polynomial
Examples
data('Arkansas')
x <- Arkansas$year
y <- Arkansas$sqrt.mayflies
model <- locally.weighted.polynomial(x,y)
plot(model)
model <- locally.weighted.polynomial(x,y,degree=2)
plot(model, derv=1)
plot(model, derv=2)
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.LocallyWeightedPolynomial.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.LocallyWeightedPolynomial
> ### Title: Plot a LocallyWeightedPolynomial object
> ### Aliases: plot.LocallyWeightedPolynomial
> ### Keywords: hplot
>
> ### ** Examples
>
> data('Arkansas')
> x <- Arkansas$year
> y <- Arkansas$sqrt.mayflies
>
> model <- locally.weighted.polynomial(x,y)
> plot(model)
>
> model <- locally.weighted.polynomial(x,y,degree=2)
> plot(model, derv=1)
> plot(model, derv=2)
>
>
>
>
>
> dev.off()
null device
1
>