similar to: spatstat 1.11-0

A new version 1.9-5 of the package `spatstat' has been uploaded to CRAN. What it is: ---------- spatstat is a package for analysing spatial data, mainly Spatial Point Patterns. What's in it: ------------ Functions for exploratory data analysis, model-fitting, simulation, spatial sampling, model diagnostics, and formal inference. Data types include point patterns, line segment patterns,

A new version 1.9-5 of the package `spatstat' has been uploaded to CRAN. What it is: ---------- spatstat is a package for analysing spatial data, mainly Spatial Point Patterns. What's in it: ------------ Functions for exploratory data analysis, model-fitting, simulation, spatial sampling, model diagnostics, and formal inference. Data types include point patterns, line segment patterns,

This is a cry for help. My apologies for taking up bandwidth with an issue that is not really on topic. But I really do want to acquire the requested information. In the course of preparing my PBRF (performance based research funds) portfolio for the upcoming PBRF exercise here in New Zealand, I have been advised that I should provide some estimate of the number of "active users" of

Hello R users, I am fitting spatial point process models in spatstat using ppm and my models actually contain spatial covariates. The results of the models are returned with some levels of the covariates absent. eg when I fit a model with "aspect" using: ppm(points, ~aspect, covariates=list(aspect=aspect)) the results give coefficients for "North"," South",

Dear all, I'm working with two point patterns (ppp) in spatstat. I turned one of them into a spatial covariate (im) object. After that, I used this im object to fit a Poisson model for the second point pattern, using the covariate layer from the first one. In R, the whole thing looks somewhat like this: my_first.im <- as.im(my_first.ppp) test.ppm <- ppm (my_second.ppp, ~my_first,

Dear r-help, I would like to use the 'ppm' function of the 'spatstat' package to fit a Strauss inhibition model. I understand that I can specify a parametric model for the "background" trend, but how would I specify a trend which is estimated using a Kernel density smoother? In particular, I would like to use the 'kde' function of the 'ks' package to

Hello, I have fitted a Poisson Process model in spatstat using >fit1<-ppm(points, ~elevation, covariates=list(elevation=elevation.im)) This far, everything went well, but I try to simulate the fitted model using the function: >sim1<-rmh(fit1) But I get the error message: "Extracting model information...Evaluating trend...done. Checking arguments..determining simulation

Hello all, I'm attempting to conduct spatial analysis of trees within a plot. I want to see if the trees are spatially correlated to soil characteristics, say pH, or moisture content. I think one way to do it is with mpl, however, my soil characteristics were not taken at exactly the same locations as my trees and further, the vectors aren't the same length. I'm getting the

Recently Brian Ripley had occasion to mock my inability to see a comment in the code for arima(), in the stats package. After considerable dredging around in the r-news archives I found reference to keep.source() and keep.source.pkgs(), which I conjectured just ***might*** possibly be the ``obvious'' resource to which Prof. Ripley thought I ought to have resorted. However, after

Hi Everyone, I am trying to figure out the spatstat package for the first time and am having some trouble. Unfortunately, I can't post my data set but I'll hopefully post enough details for some help. I want to model the intensity of a spatial point process using 2 covariates from my data. After reading through the documentation, I have successfully created 2 "ppp" objects. The

Hello, I am trying to plot a quadratcount object over a ppp object in the spatstat package. I ultimately want to get something like this http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/spatstat/man/images/big_quadratcount_001.png http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/spatstat/man/images/big_quadratcount_001.png See

I upgraded to R version 2.6.0 yesterday; now I find that I can't load spatstat. I re-installed spatstat (that part seemed to work without complaint) but when I execute > library(spatstat) I get the error message Error in dyn.load(file, ...) : unable to load shared library '/Users/rturner/Rlib/spatstat/libs/ i386/spatstat.so':

Rainer M Krug writes: > Thorsten Wiegand used in his paper Wiegand T., and K. A. Moloney 2004. > Rings, circles and null-models for point pattern analysis in ecology. > Oikos 104: 209-229 a statistic he called O-Ring statistic which is > similar to Ripley's K, only that it uses rings instead of circles. > > http://www.oesa.ufz.de/towi/towi_programita.html#ring

I'm testing some different formulations of pairwise interaction point processes in Spatstat (version 1.11-6) using R 2.5.0 on a Windows platform and I wish to simulate them using the Metropolis-Hastings algorithm implemented with Spatstat. Spatstat utilizes Fortran77 code with the preprocessor RatFor to do the Metropolis-Hastings MCMC, but the Makefile is more complicated than any I have

A full set of course notes on 'Analysing spatial point patterns in R' is now available at <http://www.csiro.au/resources/pf16h.html> The course is based on the package 'spatstat'. It includes a brief introduction to R, a detailed introduction to the 'spatstat' package, and a discussion of statistical methodology. Adrian Baddeley

Hi I'm working with the function Kest in the package SpatStat (under LINUX with R 2.1.0). In order to evaluate the statistical significance of my point pattern I'm doing 999 Montecarlo replications. The script that use the Kest function runs OK for most of the different point patterns that I have but for a particular point pattern, which have only 17 points, it runs until the 34th

Dear R developers A question about the class 'listof', defined in package 'stats'. Other than its definition and use in the code for 'anova', we can't see that the class 'listof' is used for anything else (in recommended packages, or elsewhere). In the spatstat package we have been using a 'listof' to represent a list of spatial objects of the same

Greetings - Is it possible, inside .First.lib, to find out the version number of the package that is being loaded? If only one version of the package has been installed, we could scan the DESCRIPTION file, something like .First.lib <- function(lib, pkg) { library.dynam("spatstat", pkg, lib) dfile <- system.file("DESCRIPTION", package="spatstat")

Ranjan Maitra writes: > Does anyone have a suggestion (or better still) code for sampling > from the uniform distribution over the convex hull of a set of > points? This is implemented in library 'spatstat'. If x and y are vectors of coordinates of your initial set of points, library(spatstat) W <- convexhull.xy(x, y) P <- runifpoint(42, W) will compute

Hi all, I am using quadratcount in spatstat to divide a window containing a point pattern into a grid of quadrats containing the intensity of points in each quadrat. My data is in UTM co-ordinates. My window is defined as follows: >p15<-ppp(x,y,window=owin(c(341710,342100),c(3126465,3126780)),marks=NUL L, checks=TRUE) Giving me a distance of 390m in the 'x' direction and 315m in