Displaying 1 result from an estimated 1 matches for "meansqdiff".
Did you mean:
meandiff
2009 May 20
1
Comparing spatial distributions - permutation test implementation
..."Wrong number of columns in argument x")
}
} else {
dat = data.frame(x, y, var1, var2)
}
# Normalize abundances
dat$var1 = dat$var1/sum(dat$var1)
dat$var2 = dat$var2/sum(dat$var2)
# For each point (each line of dat)
# compute the squared difference in gammas from each origin
meanSqDiff = apply(dat, 1, function(d, coord, variab) {
north = (coord$x>=d[1])
east = (coord$y>=d[2])
south = (coord$x<=d[1])
west = (coord$y<=d[2])
return( mean( c(
(diff(sapply(variab[(north & east),], sum)))^2,
(diff(sapply(variab[(south & east),], sum)))^2,
(diff(sa...