## load example data set
data(sample, envir = environment())
## sieve sample (in phi units)
sample_sieved <- prepare_Sieving(
sample = sample,
interval = c(5, 6))
## plot results
plot(density(
x = sample$grainsize,
from = -1,
to = 11))
lines(density(
x = sample_sieved$grainsize,
from = -1,
to = 11),
col = 2)
Run the code above in your browser using DataLab