I am intersted in time-space clustering of local industry.

I made the point file, polygons and time table then run 'stkhat' function in Splancs,

but it generates only time k function. I noticed that it has only one ks value. ...



Can anybody help me?. 

thanks in advance.





countiespt <- read.table('d:/dissertation/dailynew/jun10/countiespt_ok.txt', sep=",")

 polymap(countiespt)

 

 manufacpt <- read.table('d:/dissertation/dailynew/jun10/manufacpt.txt', sep=",")

 names(manufacpt)<-c("x", "y")

 

 

 manufac23<-read.table('d:/dissertation/dailynew/jun10/manufac23.txt', sep=",")

 names(manufac23)<-c("x", "y", "t")



manufac23<-read.table('d:/dissertation/dailynew/jun10/manufac23.txt', sep=",")
names(manufac23)<-c("x", "y", "t")
bur2 <- stkhat(manufacpt, manufac23$t, countiespt, c(34060, 37712),  seq(1,60,60), seq(91,1005,60))
oldpar <- par(mfrow=c(2,1))
plot(bur2$s, bur1$ks, type="l", xlab="distance", ylab="Estimated K",
  main="spatial K function")
plot(bur2$t, bur2$kt, type="l", xlab="time", ylab="Estimated K",
  main="temporal K function")
par(oldpar)



	[[alternative HTML version deleted]]


