similar to: Suppressing (or changing) the colours when using spatstat plot.quadratcounts

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

Hi, I know there are older threads discussing the quadratcount function in spatstat. Unfortunately, I could not find a solution to my problem there. I'm analyzing a point pattern in an irregular polygonal window. Both the window (an entire country) and the points are projected using WGS84. When I do quadratcount with only one quadrat for the entire country it holds all my 154 points:

A few general comments about stepwiseAIC and a suggestion of how to select models a) Apart form the problem, that stepwise selection is not a garanty to get the best model, you need to have a lot of data to avoid overfitting if your model includes 7 parameter plus interactions (> 10 observations per parameter is what you are ideally looking for). b) Have a look at Anderson and Burnham's

Hi, With the following code i can divides window into quadrats and counts the numbers of points in each quadrat. library(spatstat) X <- runifpoint(50) quadratcount(X) quadratcount(X, 4, 5) quadratcount(X, xbreaks=c(0, 0.3, 1), ybreaks=c(0, 0.4, 0.8, 1)) qX <-? quadratcount(X, 4, 5) plot(X) plot(qX, add=TRUE) But I want to mark each? quadrats? and select/ extract only those points by

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 have photographs of plots that look like so: http://r.789695.n4.nabble.com/file/n4631960/Untitled.jpg I need to divide it up so each circle has an equal area surrounding it. So into 20 equal segments, each of which contains a circle. Quadratcount is not sufficient because if I divide it up into 36 equal quadrats, some quadrats do not contain one of the circles. I'm not even sure how to

Hi, Is there any equivalent of type="n" when constructing barplots which will still construct the axes (plot=F, as it says doesn' plot anything at all). Alternatively I tried setting col="white" and border="white" but the border command does not seem to be operational. True?? Any other ideas? What I'm actually trying to do is construct vertical abline()'s

Hi, I'm new to this mailing list so apologies if this is too basic. I have confocal images 512x512 from which I have extracted x,y positions of the coordiates of labelled cells exported from ImageJ as a.csv file. I also have images that define an underlying pattern in the tissue defined as areas of different pixel values 0 or 255 (also 512x512) I've exported these images as .txt

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

http://bugzilla.mindrot.org/show_bug.cgi?id=1054 Summary: Nmap Causing SSH Session to Prematurely End Product: Portable OpenSSH Version: 3.8.1p1 Platform: All OS/Version: Mac OS X Status: NEW Severity: normal Priority: P2 Component: ssh AssignedTo: bitbucket at mindrot.org ReportedBy:

Dear list members, I'd like to perform a glm analysis with a hierarchically nested design. In particular, I have one fixed factor ("Land Use Classes") with three levels and a random factor ("quadrat") nested within Land Use Classes with different levels per classes (class artificial = 1 quadrat; class crops = 67 quadrats; and class seminatural = 30 quadrats). I have four

Hi all, When I detect the spatial point pattern, I want to use the Cramer-von Mises statistic to assess the curve-wise significance of deviations from null hypotheses. Who can tell me which function in R package "Spatstat" can do this work? Thanks a lot Jeff [[alternative HTML version deleted]]

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':

Hi, I tried to install spatstat using: install.packages("spatstat", dependencies = TRUE, lib = "./R/i486-pc-linux-gnu-library/2.8") but got the following msg: --- Please select a CRAN mirror for use in this session --- <here i selected one mirror, and tried several mirrors) Loading Tcl/Tk interface ... done Warning message: In install.packages("spatstat",

Hello I am currently studying spatial distribution of lizards in Spain and Germany. For this I am using the spatstat package of R to make a Kolmogorov-Smirnov test of CSR of unmarked points data. I have problems transforming the data to ppp class, and I would like to know if you can help me. I have already made the test with the data of 'japanesepines' that comes with the package, and

Dear R list, I'm trying to create an object of class "owin" (observation window) in the package spatstat from GIS mapping data. Here's an example of my problem. Everything goes well until the last line of code. I get the error message shown at the bottom: library(spatstat) library(sp) library(maptools) con <- url("http://gadm.org/data/rda/CHE_adm1.RData") class(con)

Hello, Somehow I cannot control my envelope plots on Spatstat package. I would like to plot an envelope with no legend and also I would like to label the plot as such: >plot(envelope(data, nsim=39, main="My K Plot", legend=FALSE, xlab="Distance") But somehow the "main", "legend" and "xlab" calls do not have any effect on the plot. Can someone

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",

Hi I want to extract the points from an object of type ppp. How can I do this? Rainer