space-time k function problem

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?. 

I think you need to coerce manufacpt to be a two-column matrix, now it is 
a data.frame, which is a different object structure.

as.matrix(manufacpt) should help, look at str(manufacpt) to see what is 
inside.
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]]

______________________________________________
R-help at stat.math.ethz.ch mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Breiviksveien 40, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 93 93
e-mail: Roger.Bivand at nhh.no

space-time k function problem

Thread (2 messages)