Displaying 1 result from an estimated 1 matches for "kernsum".
Did you mean:
kernsur
2004 Nov 18
0
Memory problems...leaky?
...probability <- matrix (NA, rowsize, colsize)
for (p in 1:rowsize) {
or (q in 1:colsize){
distance <- sqrt ((p - i)^2 + (q - j)^2)
probability [p,q] <- (2 / (bflow * pie * (1 + (distance /
bflow)^2)))
}
}
kernsum <- sum ((probability * input / input), na =
TRUE)
intmean <- sum (input * probability, na = TRUE) /
kernsum
result [i,j] <- sum (((input - intmean)^2 * probability), na =
TRUE) / kernsum
}...