Displaying 1 result from an estimated 1 matches for "lrgdata".
Did you mean:
logdata
2004 Nov 18
0
implementing a "loop" using by(x,x$factor,FUN)
...quot; to
ensure convergence.
Although not the most elegant solution, the code works ... I just want to
re-code it to avoid using the "for" loop.
Any help would be greatly appreciated.
This is what I've got ....
RootDev <- 600
splits <- NULL
for (cutpoint in cutpoints) {
LRGdata <- split(gdata, gdata[,4] <= cutpoint)
vgrmL <- variog(as.geodata(LRGdata$"TRUE"), covar.col=covcol))
vgrmR <- variog(as.geodata(LRGdata$"FALSE"), covar.col=covcol))
mlL<- likfit(as.geodata(LRGdata$"TRUE"),
ini=expand.grid(seq(0, max(vgrmL$v)...