Constructs a quantile-quantile plot against a given reference
distribution.
Usage
QQplot(x, a = 0.5, reference = c("normal", "exp", "student"), ...)
Arguments
x
vector, data for QQ-plot.
a
numeric, the offset fraction to be used in
ppoints(); typically in (0, 1).
reference
character, name of reference distribution.
...
ellipsis argument, passed down to quantile function of
reference distribution.
Details
Special forms like ParetoQQ plots can also be created via this function. E.g.,
to create a ParetoQQ plot, merely pass log(data) in place of data as the first
parameter and use reference = "exp" as the reference
distribution. The ParetoQQ plot should provide a linear graph when a
log transform of the data is plotted against the exponential distribution.
Value
Produces QQ-plot and returns invisibly a list of (x, y) pairs.
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(QRM)
Loading required package: gsl
Loading required package: Matrix
Loading required package: mvtnorm
Loading required package: numDeriv
Loading required package: timeSeries
Loading required package: timeDate
Attaching package: 'QRM'
The following object is masked from 'package:base':
lbeta
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/QRM/QQPlot.Rd_%03d_medium.png", width=480, height=480)
> ### Name: QQplot
> ### Title: Generic Quantile-Quantile Plot
> ### Aliases: QQplot
> ### Keywords: hplot
>
> ### ** Examples
>
> QQplot(rnorm(1000), reference = "normal")
> QQplot(rexp(1000), reference = "exp", rate = 0.3)
>
>
>
>
>
> dev.off()
null device
1
>