geom_qq_band  R Documentation 
Draws quantilequantile confidence bands, with an additional detrend option.
geom_qq_band( mapping = NULL, data = NULL, stat = "qq_band", position = "identity", na.rm = TRUE, show.legend = NA, inherit.aes = TRUE, distribution = "norm", dparams = list(), detrend = FALSE, identity = FALSE, qtype = 7, qprobs = c(0.25, 0.75), bandType = "pointwise", B = 1000, conf = 0.95, mu = NULL, sigma = NULL, ... ) stat_qq_band( mapping = NULL, data = NULL, geom = "qq_band", position = "identity", na.rm = TRUE, show.legend = NA, inherit.aes = TRUE, distribution = "norm", dparams = list(), detrend = FALSE, identity = FALSE, qtype = 7, qprobs = c(0.25, 0.75), bandType = "pointwise", B = 1000, conf = 0.95, mu = NULL, sigma = NULL, ... )
mapping 
Set of aesthetic mappings created by 
data 
The data to be displayed in this layer. There are three options: If A A 
stat 
statistic to use to calculate confidence bands. Should be 'qq_band'. 
position 
Position adjustment, either as a string naming the adjustment
(e.g. 
na.rm 
If 
show.legend 
logical. Should this layer be included in the legends?

inherit.aes 
If 
distribution 
Character. Theoretical probability distribution function
to use. Do not provide the full distribution function name (e.g.,

dparams 
List of additional parameters passed on to the previously
chosen 
detrend 
Logical. Should the plot objects be detrended? If 
identity 
Logical. Should an identity line be used as the reference
line used to construct the confidence bands? If 
qtype 
Integer between 1 and 9. Type of the quantile algorithm to be
used by the 
qprobs 
Numeric vector of length two. Represents the quantiles used by
the 
bandType 
Character. Either 
B 
Integer. If 
conf 
Numerical. Confidence level of the bands. 
mu 
Numerical. Only used if 
sigma 
Numerical. Only used if 
... 
Other arguments passed on to 
geom 
The geometric object to use to display the data, either as a

Tailsensitive confidence bands are only implemented for Normal QQ plots. As a future update, we intend to generalize to other distributions.
Bootstrap bands are constructed based on a MLE parametric bootstrap. Hence, it is not possible to construct such bands if the sample and theoretical distributions present mismatching supports.
Thode, H. (2002), Testing for Normality. CRC Press, 1st Ed.
# generate random Normal data set.seed(0) smp < data.frame(norm = rnorm(100)) # Normal QQ plot of Normal data gg < ggplot(data = smp, mapping = aes(sample = norm)) + stat_qq_band() + stat_qq_line() + stat_qq_point() gg + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") # Normal QQ plot of Normal data with equal local levels (ell) bands bt < "ell" gg < ggplot(data = smp, mapping = aes(sample = norm)) + stat_qq_band(bandType = bt) + stat_qq_line() + stat_qq_point() + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") gg # Exponential QQ plot of mean ozone levels (airquality dataset) di < "exp" dp < list(rate = 1) gg < ggplot(data = airquality, mapping = aes(sample = Ozone)) + stat_qq_band(distribution = di, dparams = dp) + stat_qq_line(distribution = di, dparams = dp) + stat_qq_point(distribution = di, dparams = dp) + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") gg # Detrended Exponential QQ plot of mean ozone levels di < "exp" dp < list(rate = 1) de < TRUE gg < ggplot(data = airquality, mapping = aes(sample = Ozone)) + stat_qq_band(distribution = di, detrend = de) + stat_qq_line(distribution = di, detrend = de) + stat_qq_point(distribution = di, detrend = de) + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") gg ## Not run: # Normal QQ plot of Normal data with bootstrap confidence bands bt < "boot" gg < ggplot(data = smp, mapping = aes(sample = norm)) + stat_qq_band(bandType = bt) + stat_qq_line() + stat_qq_point() + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") gg # Normal QQ plot of Normal data with tailsensitive confidence bands bt < "ts" gg < ggplot(data = smp, mapping = aes(sample = norm)) + stat_qq_band(bandType = bt) + stat_qq_line() + stat_qq_point() + labs(x = "Theoretical Quantiles", y = "Sample Quantiles") gg ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.