Generates random variates from the normal distribution by inversion. Optionally graphs the population cumulative distribution function and associated random variates, the population probability density function and a histogram of the random variates, and the empirical cumulative distribution function versus the population cumulative distribution function.
inorm(u = runif(1), mean = 0, sd = 1,
minPlotQuantile = 0.01,
maxPlotQuantile = 0.99,
plot = TRUE,
showCDF = TRUE,
showPDF = FALSE,
showECDF = FALSE,
show = NULL,
plotDelay = 0,
maxPlotTime = 30,
showTitle = TRUE,
respectLayout = TRUE)
vector of uniform(0,1) random numbers, or NULL
mean of the distribution
standard deviation of the distribution
controls the minimum quantile to appear in the plots
controls the maximum quantile to appear in the plots
logical; if TRUE
(default), one or more plots will appear (see parameters below); otherwise no plots appear
logical; if TRUE
(default), cdf plot appears, otherwise cdf plot is suppressed
logical; if FALSE
(default), pdf plot is suppressed, otherwise pdf plot appears
logical; if FALSE
(default), ecdf plot is suppressed, otherwise ecdf plot appears
shortcut way of denoting plots to display; either a binary vector of length three or an integer in [0,7] (see "Details" below)
delay, in seconds, between generation of the random variates
maximum time, in seconds, to plot all of the random variates
logical; if TRUE
(default), displays a title in the first of any displayed plots
logical; if TRUE
(default), respects existing settings for device layout (see "Details" below)
A vector of normal random variates
Generates random variates from the normal distribution, and optionally, illustrates
the use of the inverse-cdf technique,
the effect of random sampling variability in relation to the pdf and cdf.
When all of the graphics are requested,
the first graph illustrates the use of the inverse-cdf technique by graphing the population cdf and the transformation of the random numbers to random variates,
the second graph illustrates the effect of random sampling variability by graphing the population pdf and the histogram associated with the random variates, and
the third graph illustrates effect of random sampling variability by graphing the population cdf and the empirical cdf associated with the random variates.
All aspects of the random variate generation algorithm are output in red. All aspects of the population distribution are output in black.
The normal distribution has density
for
The algorithm for generating random variates from the normal distribution is synchronized (one random variate for each random number) and monotone in u. This means that the variates generated here might be useful in some variance reduction techniques used in Monte Carlo and discrete-event simulation.
Values from the u vector are plotted in the cdf plot along the vertical axis as red dots. A horizontal, dashed, red line extends from the red dot to the population cdf. At the intersection, a vertical, dashed red line extends downward to the horizontal axis, where a second red dot, denoting the associated gamma random variate is plotted.
This is not a particularly fast variate generation algorithm because it uses
the base R qnorm
function to invert the values contained in u
.
All of the elements of the u
vector must be between 0 and 1.
Alternatively, u
can be NULL
in which case plot(s) of the
theoretical pdf and cdf are displayed according to plotting parameter values
(defaulting to display of both the pdf and cdf).
The show
parameter can be used as a shortcut way to denote plots to
display. The argument to show
can be either:
a binary vector of length three, where the entries from left to right
correspond to showCDF
, showPDF
, and showECDF
,
respectively. For each entry, a 1 indicates the plot should be
displayed, and a 0 indicates the plot should be suppressed.
an integer in [0,7] interpreted similar to Unix's chmod command. That is, the integer's binary representation can be transformed into a length-three vector discussed above (e.g., 6 corresponds to c(1,1,0)). See examples.
Any valid value for show
takes precedence over existing individual
values for showCDF
, showPDF
, and showECDF
.
If respectLayout
is TRUE
, the function respects existing
settings for device layout. Note, however, that if the number of plots
requested (either via show
or via showCDF
, showPMF
, and
showECDF
) exceeds the number of plots available in the current layout
(as determined by prod(par("mfrow"))
), the function will display all
requested plots but will also display a warning message indicating that the
current layout does not permit simultaneous viewing of all requested plots.
If respectLayout
is FALSE
, any existing user settings for device
layout are ignored. That is, the function uses par
to explicitly set
mfrow
sufficient to show all requested plots stacked vertically to
align their horizontal axes, and then resets row, column, and margin settings
to R default values on exit.
The minPlotQuantile
and maxPlotQuantile
arguments are present in
order to compress the plots horizontally. The random variates generated are
not impacted by these two arguments. Vertical, dotted, black lines are
plotted at the associated quantiles on the plots.
The plotDelay
and maxPlotTime
arguments can be used to slow down
the variate generation for classroom explanation.
In the plot associated with the pdf, the maximum plotting height is associated with 125% of the maximum height of pdf. Any histogram cell that extends above this limit will have three black dots appearing above it.
# NOT RUN {
inorm(0.7, mean = 3, sd = 1)
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and PDF plots
set.seed(8675309)
inorm(runif(10), 3, 1, showPDF = TRUE)
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and ECDF plots
set.seed(8675309)
inorm(runif(10), 3, 1, showECDF = TRUE)
par(mfrow = c(3,1)) # 3 rows, 1 col for CDF, PDF, and ECDF plots
set.seed(8675309)
inorm(runif(10), 3, 1, showPDF = TRUE, showECDF = TRUE)
par(mfrow = c(1,1)) # 1 row, 1 col for PDF only plot
set.seed(8675309)
inorm(runif(10), 3, 1, showPDF = TRUE, showCDF = FALSE)
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and PDF plots
inorm(runif(120), 3, 1, showPDF = TRUE, minPlotQuantile = 0.02, maxPlotQuantile = 0.98)
# overlay visual exploration of ks.test results
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and ECDF plots
set.seed(54321)
vals <- inorm(runif(10), 3, 1, showECDF = TRUE)
D <- as.numeric(ks.test(vals, "pnorm", 3, 1)$statistic)
for (x in seq(2.5, 3.5, by = 0.05)) {
y <- pnorm(x, 3, 1)
segments(x, y, x, y + D, col = "darkgreen", lwd = 2, xpd = NA)
}
# plot the PDF and CDF without any variates
par(mfrow = c(2,1)) # 2 rows, 1 col for PDF and CDF plots
inorm(NULL, 3, 1)
# plot CDF with inversion and PDF using show
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and PDF plots
inorm(runif(10), 3, 1, show = c(1,1,0))
inorm(runif(10), 3, 1, show = 6)
# plot CDF with inversion and ECDF using show
par(mfrow = c(2,1)) # 2 rows, 1 col for CDF and ECDF plots
inorm(runif(10), 3, 1, show = c(1,0,1))
inorm(runif(10), 3, 1, show = 5)
# plot CDF with inversion, PDF, and ECDF using show
par(mfrow = c(3,1)) # 3 rows, 1 col for CDF, PDF, and ECDF plots
inorm(runif(10), 3, 1, show = c(1,1,1))
inorm(runif(10), 3, 1, show = 7)
# plot three different CDF+PDF+ECDF vertical displays
par(mfcol = c(3,3)) # 3 rows, 3 cols, filling columns before rows
set.seed(8675309)
inorm(runif(20), 3, 1, show = 7)
inorm(runif(20), 3, 1, show = 7)
inorm(runif(20), 3, 1, show = 7)
# plot three different CDF+PDF+ECDF horizontal displays, with title only
# on the first display
par(mfrow = c(3,3)) # 3 rows, 3 cols, filling rows before columns
set.seed(8675309)
inorm(runif(20), 3, 1, show = 7)
inorm(runif(20), 3, 1, show = 7, showTitle = FALSE)
inorm(runif(20), 3, 1, show = 7, showTitle = FALSE)
# exhibit use of the respectLayout = FALSE option
par(mfrow = c(3,3)) # this will be ignored below since respectLayout = FALSE
set.seed(8675309)
inorm(runif(20), 3, 1, show = 7, respectLayout = FALSE)
par("mfrow") # will have been reset c(1,1)
par("mar") # will have been reset to c(5.1, 4.1, 4.1, 2.1)
# }
Run the code above in your browser using DataLab