Description Usage Arguments Details Value Note See Also Examples
This function provides a simple way to create bootstrap replications of a dataset. The replication is either non-parametrical or parametrical (for exponential or logarithmic normal data).
rpdataboot(x, b, method = c("nonp", "exp", "lnorm"))
A matrix containing numeric values
The number of bootstrap replications
The replication method applied to the data
The input matrix
x is supposed to contain (independent) observations in each column. The
bootstrap replication take this into account is done column-wise.
Depending on how the boostrap replications are further processed, the boostrap resampling should be
done either non-parametrically (
method = "nonp") or parametrically.
In the non-parametrical case, the bootstrap replications are samples drawn from the empirical distribution of the original observation, this is equivalent to drawing with replacement.
For the parametrical bootstrap replications there are currently two options: With
"exp" each bootstrap replication is a vector simulated from an exponential distribution function
whose parameter is estimated by the original observation. For
method = "lnorm" the resampling
is done by simulating from a logarithmic normal distribution whose log-mean and log-variance are
estimated from the original observation.
An array of dimension
c(dim(x), b) containing column-wise bootstrap replications of
NA's are propagated consistently. More precisely, only the non-
NA values undergo the
resampling and thus, missing values remain unchanged in the bootstrap replications.
rpdataconv for creating a suitable data matrix from a list of observation vectors, and
rpdatasim for creating such a matrix by simulation.
1 2 3 4 5 6 7 8
# Generate a data matrix of 5 samples with 10 observations each. x <- matrix(rexp(50), nrow = 10, ncol = 5) # Create (parametric) bootstrap replications x.boot.par <- bootruin:::rpdataboot(x, b = 50, method = "exp") # Create (non-parametric) bootstrap replications x.boot.nonp <- bootruin:::rpdataboot(x, b = 50, method = "nonp")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.