Re: [R-sig-Geo] variance estimation in spgwr

2008-02-16 Thread Takatsugu Kobayashi
Hi, I am curious if spgwr or some other spatial packages estimate local variances. I am particularly interested if gwr residuals are locally siginificant or not. In this case, I believe I should estimate local variances. I though an alternative to spgwr is gam in mgcv. Suppose I have xy coor

Re: [R-sig-Geo] Changing bar scale: spplot

2008-02-11 Thread Takatsugu Kobayashi
Hi, I am reading sp and trellis manuals, but how can I change a bar scale? I am looking for something like zlim=c() Thank you! tk ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] while loop help: SOLVED

2008-01-25 Thread Takatsugu Kobayashi
I figured it out. Thanks for taking a look. thanks Takatsugu Kobayashi wrote: > Hi, I must apologize both for this lengthy thread and for posting a > thread based on the same example I have been using. > > I successfully created a code that re-arrange 10*10 grid cells of the >

Re: [R-sig-Geo] while loop help

2008-01-24 Thread Takatsugu Kobayashi
Hi, I must apologize both for this lengthy thread and for posting a thread based on the same example I have been using. I successfully created a code that re-arrange 10*10 grid cells of the same size into new grid cells of varying sizes. Now I am attempting to create grid cells of the same si

Re: [R-sig-Geo] bitmap plots in loop

2008-01-21 Thread Takatsugu Kobayashi
hi, I have read a couple of threads about png and jpg not being supported in R CMD BATCH. So I decided to use bitmap() instead. I would like to repeat printing jobs 1000 times. I only got empty png files. Doesn't this loop open bitmap, plot, close bitmap everytime? rm(list=ls(all=T)) ### Incl

Re: [R-sig-Geo] Polygon perimeter calculation in PBSmapping? [SOLVED]

2008-01-21 Thread Takatsugu Kobayashi
I figured out by reading vignette('sp') and a PBSmapping manual. Thanks. Takatsugu Kobayashi wrote: > Hi, > > I created a grid of cells of varying sizes using makeGrid() in sp. I > could calculated the areas of these cells of varying sizes using > calcArea() and combin

Re: [R-sig-Geo] overlay points and polygons and show maps

2008-01-19 Thread Takatsugu Kobayashi
Hi, I am trying to show grid cells with colors. Colors depends on the number of points falling within grid cells. Here is what I did; rm(list=ls(all=T)) library(sp) library(maptools) library(PBSmapping) ### Define parameters for a grid start.point <- -5 cco <- c(start.point,start.point) csize

Re: [R-sig-Geo] Polygon perimeter calculation in PBSmapping?

2008-01-18 Thread Takatsugu Kobayashi
Hi, I created a grid of cells of varying sizes using makeGrid() in sp. I could calculated the areas of these cells of varying sizes using calcArea() and combinePolys() in PBSmapping. I wonder if there is a function with which to compute the perimeter length. Or should I convert polygons to pol

Re: [R-sig-Geo] construct neighbor list from polygon list with user defined distance thresholds, as found in GeoDa?

2008-01-12 Thread Takatsugu Kobayashi
Also, look at this thread: https://stat.ethz.ch/pipermail/r-sig-geo/2005-July/000502.html Debarchana Ghosh wrote: > Hi, > > After reading a ESRI shape file I want to create a neighbor list (class nb) > from the polygon list by defining a distance thresholds for example 2500m. > This is a option i

Re: [R-sig-Geo] construct neighbor list from polygon list with user defined distance thresholds, as found in GeoDa?

2008-01-12 Thread Takatsugu Kobayashi
?dnearneigh in spdep package is the one I can think of. Taka Debarchana Ghosh wrote: > Hi, > > After reading a ESRI shape file I want to create a neighbor list (class nb) > from the polygon list by defining a distance thresholds for example 2500m. > This is a option in GeoDa weights creation. How

Re: [R-sig-Geo] Creating another grid from n*n grid

2008-01-12 Thread Takatsugu Kobayashi
Hi, (As an exercise)I am trying to create a SpatialPolygons object from n*n grid. By reading the material (from vignette{'sp'), I managed to convert grid to polygons. Now I would like plot this polygons, and here is what I have done so far: Step1: Create a 9*9 grid n<-9 # number of cells in ea

Re: [R-sig-Geo] GDA: set up

2008-01-09 Thread Takatsugu Kobayashi
Hi, I am trying to install rgdal package on my local directory in my school server. I uploaded gdal-1.5.0, built and compiled it on my local directory. After ./configure I got: GDAL is now configured for x86_64-unknown-linux-gnu Installation directory:/usr/local C compiler:

Re: [R-sig-Geo] adding new values to readShapePoly() object

2008-01-01 Thread Takatsugu Kobayashi
Hi Please forgive me if I am asking a simple question. I would like to add values in a column vector to readShapePoly() object and show a map based on these values using spplot() or addPoly(). Basically, my raw data are employment census data with the coordinates (X,Y). tmp1 <- read.csv("tmp1

Re: [R-sig-Geo] Create points to polygons

2007-12-31 Thread Takatsugu Kobayashi
/GIS on R. I am basically reading [R], [R-Geo], and R News... Thank you very much. Takatsugu Kobayashi PhD student Indiana University Bloomignton ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] Shortest network path distances

2007-09-15 Thread Takatsugu Kobayashi
; T. +44(0)1865 318842 > F. +44(0)1865 311383 > [EMAIL PROTECTED] > www.earthwatch.org > > > On Thu, 13 Sep 2007, Takatsugu Kobayashi wrote: > > >> Hi, >> >> I would like to know if there are packages to compute shortest network >> path distances, some

Re: [R-sig-Geo] Shortest network path distances

2007-09-13 Thread Takatsugu Kobayashi
Hi, I would like to know if there are packages to compute shortest network path distances, something that can be down on TransCAD and ArcGIS network analyst. I have cencus tract centroids and network line shapefiles for 126 US urbanized areas. Basically, I would like to move the census tract c

[R-sig-Geo] re; Any good book for R spatial stats?

2007-04-21 Thread Takatsugu Kobayashi
Hi I wonder if there are nice books for spatial analysis in R, in addition to reading R news, Journal of Computers and Geosciences, and Gnu manuals, etc. Eventually I would like to be capable of doing research using R-spatial packages and R-GRASS. My PC platform is RHEL Client 5. ArcGis is still s

Re: [R-sig-Geo] Map projectin for Europe...

2007-04-20 Thread Takatsugu Kobayashi
ct structure... Sorry for this R-geo newbie question Thank you very much while you must be busy at AAG. Best regards, Taka Roger Bivand wrote: > On Fri, 20 Apr 2007, Takatsugu Kobayashi wrote: > > >> Hi, >> >> I has a meteorology project in which I have to map

Re: [R-sig-Geo] Map projectin for Europe...

2007-04-19 Thread Takatsugu Kobayashi
Hi, I has a meteorology project in which I have to map average wind speeds over EU, as well as European countries' boundary. I have been looking online for free European maps like the one I can find on US Census for N. America. But it seems like google doesn't hit anything useful. I wonder if

Re: [R-sig-Geo] Directional statistics and weighted mean directions

2007-03-04 Thread Takatsugu Kobayashi
R users: This is more of a general question. I am trying to conduct hypothesis testing about the difference of two angles extending from the same point. I am reading directional statistics by Mardia, but it seems to talk about only frequency distribution on the circumference drawn by the unit

Re: [R-sig-Geo] How can I view a list of neighbors from the outcome knn2nb

2007-03-04 Thread Takatsugu Kobayashi
Hi I am trying to create a list of nearest neighbors using knearneigh and dinearneigh. To understand these I have been coding what these functions do: For example, for dnearneigh, suppose I have x-y coordinates and the corresponding z values # calculate distances between observations data<-cbi

Re: [R-sig-Geo] Question about dnearneigh

2007-03-01 Thread Takatsugu Kobayashi
Hi, I have been doing exploratory analysis of the link between global Moran's I coefficients and directional biases of frequencies. I created a 10*10 matrix that contains random numbers using randomNumbers(). This went well till I had to define a list of nearest neighbor of each cell. I define

Re: [R-sig-Geo] How to list neighbors IDs and corresponding frequencies

2007-02-23 Thread Takatsugu Kobayashi
Rusers: I am trying to compute a Moran's I coefficient for census tracts for Indianapolis in R. I have a shapefile for Indianapolis and uploaded it into R. I understand how to compute a Moran's I coefficient using knn2nb command. How can I compute a Moran's I based on a edge/border contiguity

Re: [R-sig-Geo] geographically weighted negative binomial regression model?

2006-10-29 Thread Takatsugu Kobayashi
Hi, I wonder if anyone know whether there is a R package for implementing a geographically weighted negative binomial regression model. Or is it possible to do this model by using a mgcv package with a family = negative.binomial? I am not sure how to determine the adaptive bandwidth in a gam pa

Re: [R-sig-Geo] create a grid data

2006-09-18 Thread Takatsugu Kobayashi
Dear Dr. Bivand: I do have both datasets. I will try interp() for an irregularly spaced data. On Fri, 25 Aug 2006 [EMAIL PROTECTED] wrote: > Dear R users- > > > I have a dataset that has X, Y and other values (Zs) attached to points > (X,Y). To show a map of this dataset, a matrix of Zs shou