spatialparsEst(gk, sigma.range, phi.range, spatial.covmodel,
covpars=c(), guess=FALSE)
To get a good choice of parameters, it is likely that the routine will have to be called several times in order to refine the choice of sigma.range and phi.range.