我很难解释relrisk
分的结果.我的数据是一个包含两个标记(两种啮齿类动物AA和RE)的多点过程,我想知道它们是否在空间上被隔离.
> summary(REkm)
Marked planar point pattern: 46 points
Average intensity 0.08101444 points per square unit
*Pattern contains duplicated points*
Coordinates are given to 3 decimal places
i.e. rounded to the nearest multiple of 0.001 units
Multitype:
frequency proportion intensity
AA 15 0.326087 0.02641775
RE 31 0.673913 0.05459669
Window: rectangle = [4, 38] x [0.3, 17] units
x 16.7 units)
Window area = 567.8 square units
relkm <- relrisk(REkm)
plot(relkm, main="Relrisk default")
The bandwidth of this relrisk
estimation is automatically selection by default(bw.relrisk
), but when I tried other numeric number using sigma
= 0.5 or 1, the results are somehow kind of weird.
How did this happened? Was it because the large proportion of blank space of my ppp?
According to chapter.14 of Spatial Point Patterns books and the previous discussion, I assume the default of relrisk
will show the ratio of intensities (case divided by control, in my case: RE divided by AA), but if I set casecontrol=FALSE
, I can get the spatially-varying probability of each type.
Then why the image of type RE in the Casecontrol=False
looks exactly same as the relrisk
estimation by default? Or they both estimate p(RE)=λRE/ λRE+λAA for each sites?
任何帮助都将不胜感激!谢谢!