Displaying 1 result from an estimated 1 matches for "idngrid".
Did you mean:
idgrid
2015 Mar 03
2
Asssistance
...for (k in seq(2.5,250,by=2.5)){
i <- 1000/k
j <- 500/k
bcidata$Xgrid <- cut(bcidata$PX, breaks = i, include.lowest = T)
bcidata$Ygrid <- cut(bcidata$PY, breaks = j, include.lowest = T)
bcidata$IDgrid <- with(bcidata, interaction(Xgrid,Ygrid))
bcidata$IDNgrid <- factor(bcidata$IDgrid)
levels(bcidata$IDgrid) <- seq_along(levels(bcidata$IDgrid))
bcidata$count <- ave(bcidata$PX, bcidata$IDgrid, FUN = length)
aggregate <- aggregate(bcidata$PX,bcidata[,c("Xgrid","Ygrid","IDNgrid")],
FUN = length)...