Hi, When you use the default knot choice in spm() (library(SemiPar)) a figure showing the knot locations is sent to the screen and you have to accept the knots to move on. I am trying to run simulations using this function. Is there a way to get spm() to use the default "REML" knots without needing to approve each set of knots? Here is an example: library(SemiPar) data(scallop) attach(scallop) log.catch <- log(tot.catch+1) fit <- spm(log.catch~f(longitude,latitude)) [[alternative HTML version deleted]]