if (FALSE) {
# Load the dataset (choose 10k or 300k version)
data('CrystalData10k')
# Fit a GeDS/GeDSboost model and compare how well the intensity peaks are captured
Gmod <- NGeDS(FQ ~ f(Q), data = CrystalData10k, phi = 0.999, q = 3)
Gmodboost <- NGeDSboost(FQ ~ f(Q), data = CrystalData10k, phi = 0.999, q = 3)
par(mfrow = c(1,2))
plot(Gmod, n = 2)
plot(Gmodboost, n = 2)
}
Run the code above in your browser using DataLab