Displaying 1 result from an estimated 1 matches for "vgrmr".
Did you mean:
varmr
2004 Nov 18
0
implementing a "loop" using by(x,x$factor,FUN)
...uot;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), l=10), seq(0,max(vgrmL$u), l=10)))
mlR<- likfit(as.geodata(LRGdata$"FALSE"),
ini=expand.grid(seq(0, max(vgrmR$v)...