Simulate a collection of soil profiles based on the horizonation of a single soil profile. Now deprecated: use perturb() for perturbations of horizon thicknesses or boundaries.

sim(x, n = 1, iterations = 25, hz.sd = 2, min.thick = 2)

Arguments

x

a SoilProfileCollection object containing a single profile from which to draw simulated data

n

the number of requested simulations

iterations

sampling iterations used to determine each horizon thickness

hz.sd

standard deviation used to simulate horizon thickness, can be a vector but must divide evenly into the number of horizons found in x

min.thick

minimum horizon thickness allowed in simulation results

Value

A SoilProfileCollection object with n simulated profiles, each containing the same number of horizons and same data as x

Details

This function generates a collection of simulated soil profiles based on the horizon thickness data associated with a single "template" profile. Simulation is based on sampling from a family of Gaussian distribution with means defined by the "template" profile and standard deviation defined by the user.

Author

D. E. Beaudette

Examples


# please see documentation for perturb() for examples
#  the sim() function calls perturb() internally