Skip to content
Prev 77828 / 398502 Next

nls and na/Nan/Inf error

Use a grid search to get the starting values in which case you
will likely be close enough that you won't run into problems
even without derivatives:

attach(fldgd)
grid <- expand.grid(Vr = seq(0,.3,.1), Vm = seq(.45, 1, .05),
	alpha = seq(1,2,.25), lamda = seq(1,2,.25))
ss <- function(p) sum((Moisture - vanGen(Suction, p[1], p[2], p[3], p[4]))^2)
idx  <- which.min(apply(grid, 1, ss))
startval <- grid[idx,]
nls(Moisture ~ vanGen(Suction, Vr, Vm, alpha, lamda), start = startval)
On 9/26/05, Tony Meissner <tony.meissner at bigpond.com> wrote: