Displaying 1 result from an estimated 1 matches for "k_funk_ca".
2010 Apr 08
0
khat and included polygons
...al of the maptools library.
Then I split up the SpatialPolygonsDataFrame into the included polygons by
coor <- NULL
for (i in 1:[number of polygons]) coor<- rbind(coor,map@polygons[[1]]@Polygons[[i]]@coords)
I am using now this newly generated polygon "coor" in the khat command.
k_funk_ca <- khat(as.points(cases), coor,s=0:1000)
I do not know if there is a problem using multiple polygons that way or if khat produces then only results for the "first" polygon or even completely wrong results.
Does anyone know the answer?
Thanks a lot,
Sven Schmiedel
PhD student at...