Displaying 1 result from an estimated 1 matches for "myowin".
Did you mean:
mowin
2004 May 29
1
Rhelp: Need help interpreting plots in spatstat
...t the empirical Ghat and the simulation envelope...but that
first plot stumps me.
You can see a pdf of the resulting plot if you download it from my
website : www.bol.ucla.edu/~lschwei/sample.pdf WARNING: this is
immediate download.
THE CODE: ######
ghat.env <- function(n, s, r, win=owin(myOwin){
hold <- matrix(0, s, length(r))
for(i in 1:s){
hold[i,] <- Gest(runifpoint(n, win=myOwin), r=r)$raw
}
mn <- apply(hold, 2, mean)
Up <- apply(hold, 2, max)
Down <- apply(hold, 2, min)
return(data.frame(mn, Up, Down))
}
tsds.ghat <- Gest(tsdspoints)
tsds.ghat.short <- tsds.ghat...