similar to: Reading polygon shapefiles in splancs

Hi, I have a problem that concerns combination of the package Maptools and Splancs I have 2 shapefiles that i want to manipulate (one of type point and one polygon).I import them in R using Maptools but then i can't estimate a quartic Kernel using Splancs. The package doesn't recognize the shapes (invalid points and poly argument).I don't know if this is an easy task but i have

Dear forum, I would like to use the kernel2d or spkernel2d in the Splancs-package, but it does not recognize my polygon data. "Error in kernel2d(as.points(ptsbin), polygonprov, h0 = 2, nx = 100, : " is the error message. Invalid poly argument The data are defined as follows: polgonprov<-list(x=polyprov$X, y=polyprov$Y) with X and Y coordinates in the Lambert1972

Dear all, I had to place points at random, one in each of larger number of polygons (actually in objects of class 'SpatialPolygonsDataFrame' , see sp library), and tried first to do it using spsample (from sp). Surprisingly, every 5-15 trials, the output was a NULL value. The doc says that ' this may occur when trying to hit a small and awkwardly shaped polygon in a large

Hi, I would need some help with the splans package in R. I am using a Shapefile (downloadable at) http://rapidshare.com/files/215206891/Redlands_Crime.zip and the following execution code setwd("C:\\Documents and Settings\\Dejan\\Desktop\\GIS\\assignment6\\DataSet_Redlands_Crime\\Redlands_Crime") library(foreign) library(splancs) auto_xy<-read.dbf("Auto_theft_98.dbf")

Hello, I am trying to conduct a bi-variate ripley's k using the k12hat function in the splancs package (found here http://rss.acs.unt.edu/Rdoc/library/splancs/html/k12hat.html). Although, I receive an error when getting to the bboxx part of the code: poly <- list(x=c(fire4$X, nar4$X), y=c(fire4$Y, nar4$Y)) plot(seq(5,80,5), sqrt(k12hat(fire4), as.points(nar4),

>From version 2.01-17, splancs, a package for spatial point pattern analysis, depends on package sp. Because both packages have defined a bounding box function bbox(), the bbox() method in sp has been given priority, and the default object for dispatch is an array with two or more columns (as in splancs until now). However, the deprecated splancs function returned a non-closed rectangle of

>From version 2.01-17, splancs, a package for spatial point pattern analysis, depends on package sp. Because both packages have defined a bounding box function bbox(), the bbox() method in sp has been given priority, and the default object for dispatch is an array with two or more columns (as in splancs until now). However, the deprecated splancs function returned a non-closed rectangle of

Greetings: are there any simple ways to convert an sp class object to a splancs ppp class object, outside of reading the coordinates and computing a bounding box? cheers, -- Dylan Beaudette Soils and Biogeochemistry Graduate Group University of California at Davis 530.754.7341

Dear R Users; I have been using the space-time K-function analysis to detect the existence of clustering in both space and time. I wish to ask and know what the minimum number of cases should be for such analyses to be valid. I have data with 23 cases and don't know how careful i should be in the interpreting the results. Secondly, i wish to find out if there is any bond between the

I am trying to port the splancs library for S++ version 4.0 to R. I have modified the original FORTRAN to use dllwrap/dlltool (e.g. removed common blocks). I have been able to generate dll's (splancs.dll) using dllwrap/dlltool and Ming32 g77 with a -mdll switch. As of yet I have not been able to generate a dll that can be loaded by R. Does anyone out there have any suggestions. An example

Howdy Gurus I am try to compute envelope of Khat from simulations of complete spatial randomness using R package of "Splancs." I got the following error: > UL.khat <- Kenv.csr(length(X_coord), bnd, nsim=100, s) Doing simulation 1 Error in runif(n, min, max) : invalid arguments In addition: Warning message: NAs introduced by coercion I think that I have a problem with

Hi, I previously received help in extract data from a shapefile and now my question is about kernel density estimation. My objective is to have 3 kernel density plots; 2 for the each set of cases and the 3rd is the difference in kernel densities between the 2 sets of cases. Previously, I used the density function from the stats package, which worked but I wanted finer control of the bandwidth.

Dear all, I want my PhD thesis which I hand in tomorrow to look even nicer: Does polymap in the splancs library not allow horizontal plotting of y-labels? I have tried polymap(studyarea, xlab = "x (m)", ylab = "y (m)", las = 1) but it doesn't change the labels? Mayby some function in library(spatstat) support las? Thanks! Sincerely, Tord

I have a problem with function 'inout()' in package 'splancs' on CRAN-R Version 1.2.3 under FreeBSD4.3-STABLE. The following script produces and draws points and a polygon-surrounded area. Repeating the same script many times shows, that points on the polygon-line often, but not always, are outside of the polygon-area. library(splancs) # dataset with polygon (convex

I installed R-1.4.0 and then updated all the packages I have, using update.packages(). I got an error for splancs and sgeostat. Trying to install them 'manually' with R CMD INSTALL, I got the following errors. 1) for splancs: Tue<15:26>library[82]#R CMD INSTALL splancs_2.01-5.tar.gz Installing *source* package `splancs' ... libs f95 -PIC -xO5 -dalign -xlic_lib=sunperf -c

I am not sure a previous e-mail reached the list (no mail aknowledgement from R-boundle etc.). The question was how to write polygon or segment coordinates into a shapefile set or any other ArcGIS supported format. The library shapefiles seems to do something but the documentation is a bit beyond of my mind.... and I cannot get the meaning of the functions write**** and its application to the case

The main limitation of the shapefiles package that I put together is that it does not create shapefiles from R objects - rather it only writes shapefiles that have been read into R and manipulated within the constraints of the existing file structure. By this I mean that for example you can change the coordinates of points and write them back out. Or you can add a bunch of blank columns in the

Dear List, I've been running a numerical simulation and I found this odd error in my code where the which command could not identify which rows of a column of data.frame were corresponding to the value 0.3. There are 7 unique values in this column (0.01,0.05,0.1,0.2,0.3,0.4,0.5), and this does not work only for 0.3. So I looked at the column and manually tried to use the which() command, and

Dear List, I have a set of distributions recorded at an equal interval of time and I would like to plot them as series of horizontal histograms (with the x-axis representing time, and y-axis representing the bins) since the distribution shifts from unimodal to multimodal in several occasions. What I would like to see is something close to a violinplot, but I do not want a kernel density

Dear List, I am struggling with the following problem. Suppose I have the following data.frame: coord<-expand.grid(x=1:10,y=1:10) and I want to extract the row numbers of those matching the criteria defined by the following data.frame: choice<-data.frame(x=c(2,1,2),y=c(1,2,2)) the result should return in this case: 2,11,12 All I managed to do is the following: which(data