<div dir="ltr"><p class="MsoNormal"><span lang="EN-US">Hi Rolf and
Julia</span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">Many thanks
for your help, I appreciate that. I am sorry for the inconvenient, I am new in
these things.</span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">Just to
clarify, my window is complex because is created from a "shapefile"
that delimit my study area. This study
area (called "W" in the analysis) have the spatial locations of my
data ("ficus" and "tienda").</span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">With the
purpose of reproduce the analysis and try to find the problem, in the
attachment documents you can find the "shapefile" of my study area, a
"txt" document of my data and my R script (Kcrossinhom(DV)). </span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">Again, I
really appreciate your help. Thank you.</span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p>
<p class="MsoNormal"><span lang="EN-US">Cheers</span></p></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/10/21 Julia Chacon Labella <span dir="ltr"><<a href="mailto:juliachacon@gmail.com" target="_blank">juliachacon@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Yes, I tried to reproduced your example, and I didn´t have any problem. Have you check how are your lambdas?? Why are you using "at points" instead of the default "by pixel"??<div>
<br></div>
<div>It seems that could be the lambda values what is wrong. Try to visualized them, and try with the pixel image.<div><br><div>Other option is to built a new window for your pattern. How did you built it?<br><div><br></div>
</div></div></div></div><div class="gmail_extra"><div><div class="h5"><br><br><div class="gmail_quote">2013/10/22 Rolf Turner <span dir="ltr"><<a href="mailto:r.turner@auckland.ac.nz" target="_blank">r.turner@auckland.ac.nz</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Again, your problem is not reproducible. We do not have the objects<br>
"north", "west", "W" and "tipo" which are necessary to duplicate your<br>
steps.<br>
<br>
I tried to reproduce your problem with an artificial example:<br>
<br>
W <- owin(c(0,3350),c(0,6150))<br>
set.seed(42)<br>
Z <- runifpoint(201,win=W)<br>
m <- sample(rep(c("ficus","tienda")<u></u>,c(131,70)),201)<br>
marks(Z) <- factor(m)<br>
ff <- with(split(Z),bw.diggle(ficus)<u></u>)<br>
tt <- with(split(Z),bw.diggle(<u></u>tienda))<br>
lambdaf <- with(split(Z),density(ficus,<u></u>ff,at="points"))<br>
lambdat <- with(split(Z),density(tienda,<u></u>tt,at="points"))<div><br>
K <- Kcross.inhom(Z,"ficus","<u></u>tienda",lambdaf,lambdat)<br>
<br></div>
No problema.<br>
<br>
It is possible that the error has something to do with the rather intricate<br>
(1716 vertices) boundary of your observation window. However since I don't<br>
have that window I can't tell.<br>
<br>
cheers,<br>
<br>
Rolf Turner<div><div><br>
<br>
<br>
On 10/22/13 10:16, david.villalobos wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Sure Julia,<br>
<br>
First:<br>
<br>
I have the following object called "Z"<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Z <- ppp(north, west,window=W,marks=tipo)<br>
summary(Z)<br>
</blockquote>
Marked planar point pattern: 201 points<br>
Average intensity 4.46e-05 points per square unit<br>
<br>
Coordinates are integers<br>
i.e. rounded to the nearest 1 unit<br>
<br>
Multitype:<br>
frequency proportion intensity<br>
ficus 131 0.652 2.90e-05<br>
tienda 70 0.348 1.55e-05<br>
<br>
Window: polygonal boundary<br>
single connected closed polygon with 1716 vertices<br>
enclosing rectangle: [484546.8, 487897.4] x [1147521.3, 1153682.9] units<br>
Window area = 4509680 square units<br>
<br>
Second:<br>
<br>
I calculated the bandwidth sigma for the kernel estimator of point process<br>
intensity with a Cross Validated Bandwidth Selection for Kernel Density<br>
(bw.diggle)<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
attach(split(Z))<br>
f<-bw.diggle(ficus)<br>
t<-bw.diggle(tienda)<br>
</blockquote>
Third:<br>
<br>
I calculated with the Kernel Smoothed Intensity of Point Pattern<br>
(density.ppp) the values of the the estimated intensity of the sub-process<br>
of points of type i and j:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
lambdaf<-density.ppp(ficus,f,<u></u>at="points")<br>
lambdat<-density.ppp(tienda,t,<u></u>at="points")<br>
</blockquote>
and Finally:<br>
<br>
I ran the script for Kcross.inhom for the method (1): estimate intensities<br>
by nonparametric smoothing, but I had the error.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
K <- Kcross.inhom(Z, "ficus", "tienda", lambdaf, lambdat)<br>
</blockquote>
Error in Kwtsum(dclose, bI, weight, b, 1/lambdaI, breaks) :<br>
Weights in K-function were infinite or NA<br>
<br>
I tryed the method (2): leave-one-out, but I had the same problems<br>
<br>
Any help is appreciated. A lot of thanks.<br>
<br>
<br>
<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://r-sig-geo.2731867.n2.nabble.com/Problems-with-Kcross-inhom-tp7584914p7584920.html" target="_blank">http://r-sig-geo.2731867.n2.<u></u>nabble.com/Problems-with-<u></u>Kcross-inhom-<u></u>tp7584914p7584920.html</a><br>
Sent from the R-sig-geo mailing list archive at Nabble.com.<br>
<br>
______________________________<u></u>_________________<br>
R-sig-Geo mailing list<br>
<a href="mailto:R-sig-Geo@r-project.org" target="_blank">R-sig-Geo@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/<u></u>listinfo/r-sig-geo</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
R-sig-Geo mailing list<br>
<a href="mailto:R-sig-Geo@r-project.org" target="_blank">R-sig-Geo@r-project.org</a><br>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/<u></u>listinfo/r-sig-geo</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br>Julia<br>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr"><div><br></div><div>Saludos</div><div><br></div><div>David Villalobos Chaves<br></div><div>Maestría Académica en Biología</div><div>Universidad de Costa Rica</div>
</div>
</div>