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]]

