R help - Jan 2006 - envelopes of simulations

Hello!

I am writing you because I could not plot the confidence envelopes for
functions Jest, Jcross, Jdot, Jmulti, and L, using the Spatstat package.

I have already understood how to do that for Kest or Jest, that is:
JEnv <- plot(envelope(PPPData, Jest))
	Where PPPData is my ppp object.


However, for Jcross I must specify the two marks I want to analyse.
That is, usually I would get the Jcross doing:
Jc <- Jcross(PPPData,"Aun","Qsu")
	For marks "Aun" and "Qsu".


For L function, I can make:
K <- Kest (PPPData, correction="isotropic")
plot (K, r-sqrt(iso/pi)~r)


But I could not understand how to put all that together. 

Resuming, I do not understand how can I use the envelope function:
-	to plot Jcross (where and how do I specify types/marks ??)
-	to plot L (where do I write that I do not want K, but
r-sqrt(iso/pi)~r, instead??)
-	etc. for similar functions not so simple as Jest, Kest or Gest. 

Hope someone can please help me on this.
I have already looked at the package's manual, but I could not find any
example similar to what I am trying to do.


Best Regards,
Sara Mouro.






	[[alternative HTML version deleted]]

Sara Mouro writes:
 > Hello!
 > 
 > I am writing you because I could not plot the confidence envelopes for
 > functions Jest, Jcross, Jdot, Jmulti, and L, using the Spatstat package.

Enquiries about a package should be sent to the package maintainer
rather than R-help. 

 > I have already understood how to do that for Kest or Jest, that is:
 > JEnv <- plot(envelope(PPPData, Jest))
 > 	Where PPPData is my ppp object.
 > 
 > However, for Jcross I must specify the two marks I want to analyse.
 > That is, usually I would get the Jcross doing:
 > Jc <- Jcross(PPPData,"Aun","Qsu")
 > 	For marks "Aun" and "Qsu".

To do this, type
      envelope(PPPData, Jcross, i="Aun", j="Qsu")

Explanation:
	Looking at help(Jcross) we see that its formal syntax is 
		Jcross(X, i, j). 
	This means that when you call 
		Jcross(PPPData,"Aun","Qsu")
	the arguments are matched as 
		Jcross(X=PPPData,i="Aun",j="Qsu").

	So our problem is, how to pass the arguments 	i="Aun" and
j="Qsu"
	to the function Jcross each time it is called by `envelope'.

	The help entry for 'envelope' mentions that `envelope' will accept 
	extra arguments "..." which are then passed to the function
'fun'.
	
	So, to pass the two argument values i="Aun" and j="Qsu",
	just type
	     envelope(PPPData, Jcross, i="Aun", j="Qsu")

 > For L function, I can make:
 > K <- Kest (PPPData, correction="isotropic")
 > plot (K, r-sqrt(iso/pi)~r)

	You can compute the envelope for the K-function
	and transform it afterwards. 

	To plot envelopes of the L function discrepancy
	   LD(r) = sqrt(K(r)/pi) - r,
	just type 
		E <- envelope(PPPData, Kest, correction="isotropic")
		plot(E, sqrt(./pi) -r ~ r)

Hope this helps.
regards
Adrian Baddeley