sebaran spasial dua tipe titik - ipb university...ripley’s k estimator for mulitype point function...

30
Sebaran Spasial Dua Tipe Titik Praktikum Statistika Spasial | Pertemuan 3 [email protected]

Upload: others

Post on 26-Dec-2019

7 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Sebaran Spasial Dua TipeTitik

Praktikum Statistika Spasial | Pertemuan 3

[email protected]

Page 2: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

MultiType Point Patterns

• A marked point pattern in which the marks are a categorical variable is usually called a multitype point pattern. The ‘types’ are the different values or levels of the mark variable

Page 3: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Inst

alle

d D

ata

Set

Page 4: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Ants Data Set

Page 5: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Species Semut yg Diamati

Messor wasmanni Cataglyphis bicolor

Page 6: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Load Ants Data Set

Page 7: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 8: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 9: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Kuadran

Page 10: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

• CAB=Banyaknya sel yang berisi kedua himpunan PA dan PB=14

• CA0=Banyaknya sel yang hanya berisi himpunan PA =0

• C0B=Banyaknya sel yang hanya berisi himpunan PB =4

• C00=Banyaknya sel kosong =6

• CA= CAB + CA0 =14

• CB= CAB + C0B =18

• C = CAB + CA0+ C0B + C00 =24

CAB=14 C0B=4 CB=18

CA0=0 C00=6 C-CB=6

CA=14 C-CA=10 C=24

Page 11: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

CAB=14 C0B=4 CB=18

CA0=0 C00=6 C-CB=6

CA=14 C-CA=10 C=24

𝑌𝐴𝐵 =14(18)

24= 10.5 𝑌0𝐵 =

10(18)

24= 7.5

𝑌𝐴0 =14(6)

24= 3.5 𝑌00 =

10(6)

24= 2.5

𝜒ℎ𝑖𝑡𝑢𝑛𝑔2 =

14 − 10.5 2

10.5+

4 − 7.5 2

7.5+

0 − 3.5 2

3.5+

6 − 2.5 2

2.5= 11.2

𝜒𝑡𝑎𝑏𝑒𝑙2 = 𝜒0.05(1)

2 = 3.84

Hipotesis:

H0: Titik-titik pada Himpunan PA dan PB

adalah menyebar secara independen

H1: Titik-titik pada Himpunan PA dan PB

adalah menyebar secara tidak independen(Mereka berkorelasi secara spasial antarasesamanya)

𝜒ℎ𝑖𝑡𝑢𝑛𝑔2 > 𝜒𝑡𝑎𝑏𝑒𝑙

2 cukup bukti untuk menolak H0. Artinya, dapat disimpulkan bahwa pada 𝛼 = 5%, sebaran spasial

dari kedua species semut tersebut tidak saling bebas.

Page 12: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Tetangga Terdekat

The estimator Kest ignores marks. Its counterparts for multitype point patterns are Kcross, Kdot, and for general marked point patterns see Kmulti.

Page 13: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Ripley’s K estimator for Mulitype Point

Function Description

Kdot For a multitype point pattern, estimate the multitype K function which counts the expected number of other points of the process within a given distance of a point of type i.

Kcross For a multitype point pattern, estimate the multitype K function which counts the expected number of points of type j within a given distance of a point of type i.

Kmulti For a marked point pattern, estimate the multitype K function which counts the expected number of points of subset J within a given distance from a typical point in subset I.

Page 14: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Tetangga Terdekat

# diagnostic plot for independence

between Messor and other species

Page 15: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 16: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 17: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

plot(envelope(ants, Kcross,

i="Messor",

j="Cataglyphis", nsim=99,

fix.n=T, fix.marks=T))

• If the process of type i points were independentof the process of type j points, then Kij(r) wouldequal pi * r^2.

• Deviations between the empirical Kij curve andthe theoretical curve pi * r^2 may suggestdependence between the points of types i and j.

Page 18: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

mad.test(ants, Kcross, i="Messor", j="Cataglyphis", nsim=99,

fix.n=T, fix.marks=T)

Page 19: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Lansing Wood Data

• Here is the famous Lansing Woods dataset recording the positions of 2251 trees of 6 different species (hickories, maples, red oaks, white oaks, black oaks and miscellaneous trees).

Page 20: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Load the data

Page 21: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 22: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 23: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

plot(split(lansing))

Page 24: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 25: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 26: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function
Page 27: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Kuadranplot(quadratcount(split(lansing)))

Page 28: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Kuadran

Page 29: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Metode Tetangga Terdekat

Page 30: Sebaran Spasial Dua Tipe Titik - IPB University...Ripley’s K estimator for Mulitype Point Function Description Kdot For a multitype point pattern, estimate the multitype K function

Referensi

1. Baddeley, A. (2008). Analysing spatial point patterns in R.

2. Pustaka lain yang relevan