Displaying 1 result from an estimated 1 matches for "nggis".
Did you mean:
nggi
2012 Mar 03
2
contour for plotting confidence interval on scatter plot of bivariate normal distribution
...the existing scatterplot
par(new=TRUE)
contour(den, levels=confidence.border, col = "red", add = TRUE)
Unfortunately I doubt very much this is correct :( In fact I am sure
this is wrong, because the border for probs=0.05 is drawn outside the
values.... So please help and check.
Pascal H?nggis code seems to work, but I don't understand the magic he
does with
pp <- array()
for (i in 1:1000){
z.x <- max(which(den$x < x[i]))
z.y <- max(which(den$y < y[i]))
pp[i] <- den$z[z.x, z.y]
}
before doing the very same as I did above:
confidencebound &...