Warning: 16 points were rejected as lying outside the specified window
Warning: data contain duplicated points
summary(sismospp)
Warning in diff(xrange) * diff(yrange): NAs produced by integer overflow
Planar point pattern: 371 points
Average intensity 3.285698e-10 points per square Metro
*Pattern contains duplicated points*
Coordinates are given to 3 decimal places
i.e. rounded to the nearest multiple of 0.001 Metros
Window: polygonal boundary
single connected closed polygon with 216 vertices
enclosing rectangle: [452302, 1801101] x [68566, 1870203] Metros
(1349000 x 1802000 Metros)
Window area = 1.12914e+12 square Metros
Unit of length: 1 Metro
Fraction of frame area: NA
*** 16 illegal points stored in attr(,"rejects") ***
Gráfico de las ubicaciones que forma el patrón observado (muestreado)
plot(sismospp, main="")
Warning in plot.ppp(sismospp, main = ""): 16 illegal points also plotted
title("Patrón de sismos en Colombia (julio-diciembre 2008)", cex.main =1, font.main=2, col.main="blue")
prueba_chi <-quadrat.test(sismospp)
Warning: Some expected counts are small; chi^2 approximation may be inaccurate
plot(quadrat.test(sismospp),cex=0.4)
Warning: Some expected counts are small; chi^2 approximation may be inaccurate
prueba_chi <-quadrat.test(sismospp)
Warning: Some expected counts are small; chi^2 approximation may be inaccurate
prueba_chi
Chi-squared test of CSR using quadrat counts
data: sismospp
X2 = 2359.8, df = 20, p-value < 2.2e-16
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Neyman modified X2 statistic NM2
data: sismospp
NM2 = 97.494, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
prueba_chi_1
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Pearson X2 statistic
data: sismospp
X2 = 2359.8, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
prueba_chi_05
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Freeman-Tukey statistic T2
data: sismospp
T2 = 1229.7, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
prueba_chi_1_3
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Cressie-Read statistic (lambda = 0.333333333333333)
data: sismospp
CR = 1315.4, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
prueba_chi_2_3
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Cressie-Read statistic (lambda = 2/3)
data: sismospp
CR = 1694, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
prueba_chi_1
Conditional Monte Carlo test of CSR using quadrat counts
Test statistic: Pearson X2 statistic
data: sismospp
X2 = 2359.8, p-value = 0.001
alternative hypothesis: two.sided
Quadrats: 21 tiles (irregular windows)
Hopkins-Skellam test of CSR
using F distribution
data: sismospp
A = 0.0048174, p-value < 2.2e-16
alternative hypothesis: two-sided
HE_sismos_clust
Hopkins-Skellam test of CSR
using F distribution
data: sismospp
A = 0.0056189, p-value < 2.2e-16
alternative hypothesis: clustered (A < 1)
HE_sismos_reg
Hopkins-Skellam test of CSR
using F distribution
data: sismospp
A = 0.0048705, p-value = 1
alternative hypothesis: regular (A > 1)
plot(density(sismospp))contour(density(sismospp), add=TRUE)title("Estimación de la intensidad de sismos en Colombia (julio-diciembre 2008)",cex.main =1,font.main=2,col.main="black")
Warning: Internal error: fvlabels truncated the function name
Warning: Internal error: fvlabels truncated the function name
Warning: Internal error: fvlabels truncated the function name
plot(pcf(K_sismos_NH))
Fest(sismospp)
Function value object (class 'fv')
for the function r -> F(r)
.....................................................................
Math.label Description
r r distance argument r
theo F[pois](r) theoretical Poisson F(r)
cs hat(F)[cs](r) Chiu-Stoyan estimate of F(r)
rs hat(F)[bord](r) border corrected estimate of F(r)
km hat(F)[km](r) Kaplan-Meier estimate of F(r)
hazard hat(h)[km](r) Kaplan-Meier estimate of hazard function h(r)
theohaz h[pois](r) theoretical Poisson hazard h(r)
.....................................................................
Default plot formula: .~r
where "." stands for 'km', 'rs', 'cs', 'theo'
Recommended range of argument r: [0, 108010]
Available range of argument r: [0, 108010]
Unit of length: 1 Metro
plot(Fest(sismospp))
Jest(sismospp)
Function value object (class 'fv')
for the function r -> J(r)
......................................................................
Math.label Description
r r distance argument r
theo J[pois](r) theoretical Poisson J(r)
rs hat(J)[rs](r) border corrected estimate of J(r)
han hat(J)[han](r) Hanisch-style estimate of J(r)
km hat(J)[km](r) Kaplan-Meier estimate of J(r)
hazard hazard(r) Kaplan-Meier estimate of derivative of log(J(r))
......................................................................
Default plot formula: .~r
where "." stands for 'km', 'han', 'rs', 'theo'
Recommended range of argument r: [0, 105610]
Available range of argument r: [0, 105610]
Unit of length: 1 Metro
plot(Jest(sismospp))
Gest(sismospp)
Function value object (class 'fv')
for the function r -> G(r)
.....................................................................
Math.label Description
r r distance argument r
theo G[pois](r) theoretical Poisson G(r)
han hat(G)[han](r) Hanisch estimate of G(r)
rs hat(G)[bord](r) border corrected estimate of G(r)
km hat(G)[km](r) Kaplan-Meier estimate of G(r)
hazard hat(h)[km](r) Kaplan-Meier estimate of hazard function h(r)
theohaz h[pois](r) theoretical Poisson hazard function h(r)
.....................................................................
Default plot formula: .~r
where "." stands for 'km', 'rs', 'han', 'theo'
Recommended range of argument r: [0, 27846]
Available range of argument r: [0, 105610]
Unit of length: 1 Metro
plot(Gest(sismospp))
Hest(sismospp)
Function value object (class 'fv')
for the function r -> H(r)
.........................................................................
Math.label Description
r r distance argument r
km hat(H)[km](r) Kaplan-Meier estimate of H(r)
hazard hat(lambda)(r) Kaplan-Meier estimate of hazard function lambda(r)
han hat(H)[han](r) Hanisch estimate of H(r)
rs hat(H)[bord](r) border corrected estimate of H(r)
.........................................................................
Default plot formula: .~r
where "." stands for 'km', 'han', 'rs'
Recommended range of argument r: [0, 272730]
Available range of argument r: [0, 1126100]
Unit of length: 1 Metro
