Displaying 1 result from an estimated 1 matches for "olddiff".
Did you mean:
olddf
2006 Apr 18
2
Unfound objects in function
...how I did or
how I can cancel it!
Thanks in advance,
Colin
library (spatstat)
kernelEst <- function (data, max.r = 10, edge = TRUE) {
if (!is.ppp (data)) stop ("data must be a point process object")
i <- 0.1
diff <- 1
smo <- 1
while (diff > 0) {
olddiff <- diff
oldsmo <- smo
smo <- density.ppp (data, sigma = 2 * i, edge = edge)
lambda <- smo[data]
Ki <- Kinhom (data[,data$window], lambda = lambda, r = seq (0,
max.r, length = 100), correction = "trans")
diff <- sum (Ki$theo - Ki$trans)...