Re: [R] spatstat installation

If you don't receive a satisfactory answer here in a day or so, post on the R-sig-geo list where the experts in this sort of thing hang around. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley

Re: [R] [spatstat] Convert shapefile to pixel image

Try fasterize, converting to spatstat from raster is straightforward, happy to help. https://github.com/ecohealthalliance/fasterize Cheers, Mije On Fri, 2 Jun 2017, 07:45 Rolf Turner, wrote: > > On 02/06/17 01:17, lluis.hurt...@uv.es wrote: > > > Dear all, > > > > I

Re: [R] [spatstat] Convert shapefile to pixel image

On 02/06/17 01:17, lluis.hurt...@uv.es wrote: Dear all, I am currently working with the spatstat package, using windows and pixel images. First: My aim is to transform a shapefile (see attached) into a pixel image. My idea is to start transforming the shapefile into a Spatial Polygon

Re: [R] [spatstat] Convert shapefile to pixel image

This looks like to would be better to ask on R-sig-geo, instead of R-help. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 6/1/17, 6:17 AM, "R-help on behalf of lluis.hurt...@uv.es"

Re: [R] spatstat rmh problem

OK. Two things are going wrong. (1) There is an error in your code. You are passing the new.coef argument to density() and not to rmh(). The function density() has no such argument, but has a ... argument, so new.coef simply gets ignored. You should use:

Re: [R] spatstat rmh problem

There was indeed a bug in rmh() w.r.t. the new.coeff argument. The bug has been fixed and will not be present in the next release of spatstat. cheers, Rolf Turner On 16/09/14 16:30, Sebastian Schutte wrote: Thanks so much for your comments. Sorry for not having sent a running example from

Re: [R] spatstat rmh problem

You should use: plot(density(rmh(mod,new.coef=c(1,200 Sorry, my bad, typo in the example code. (2) However, even when the correct call is given you still wind up with identical densities!!! Hm. I think this may be a bug; I'll will check with the other authors of spatstat and

Re: [R] spatstat rmh problem

Your example is not reproducible. We don't have cshape or im.pop (and are possibly lacking other bits and pieces; I didn't check the details since the example fails to run from the get-go). Please provide a *reproducible* example. Also I am puzzled by the line mod - ppm (ppp, ~ pop ,

Re: [R] spatstat rmh problem

Thanks so much for your comments. Sorry for not having sent a running example from the start. Here it is: library (spatstat) #Load example data
data(demopat) #Generate a random point pattern within the polygon set.seed(12345) pdat - rpoint(200,win=demopat$window) #Generate a distmap, which

Re: [R] spatstat package: Simulating from fitted Matern cluster process model fails

There appears to be a small bug in the code. What is happening is that occasionally there are simulated points that lie inside your triangular window but do not lie inside any pixel of the image created from your distorigin() function. This will be fixed in a future release of spatstat.

Re: [R] spatstat output

On 07/07/13 22:12, catalin roibu wrote: Dear R users, Is there a possibility to extract only the r, CI's envelope and L function from the output of spatstat? I use this code E - alltypes(df1, Kest, nsim = 100, envelope = TRUE,savepatterns=TRUE,correction=isotropic) And second question, is there

Re: [R] spatstat error

Why do you first say min(Datos$x) and then give the numeric value of this minimum? Just delete all numerical values: danta=ppp(Datos$x, Datos$y, c(min(Datos$x), max(Datos$x)), c(min(Datos$y), max(Datos$y))) or if you really want to type the numeric constants (probably a bad idea)

Re: [R] spatstat error

Questions about spatstat should be directed to the R-Sig-Geo list, or better still, to the package maintainers. To add a bit to what Blaser Nello has already told you: (1) Your syntax in respect of the use of the min() function is indeed incomprehensibly bizarre. (2) Blaser Nello has

Re: [R] spatstat error

Questions about spatstat should be directed to the R-Sig-Geo list, or better still, to the package maintainers. To add a bit to what Blaser Nello has already told you: (1) Your syntax in respect of the use of the min() function is indeed incomprehensibly bizarre. (2) Blaser Nello has indicated

Re: [R] spatstat kstest

Hello, Try the following. maintainer(spatstat) [1] Adrian Baddeley adrian.badde...@csiro.au Hope this helps, Rui Barradas Em 07-02-2013 02:05, Hiroshi Saito escreveu: Dear sir, Which mailing list is appropriate to ask for kstest in spatstat? Regards, Hiroshi Saito

Re: [R] Spatstat - Kest, envelope problem

You can solve this problem by actually *reading* the error message. Evidently the object point is a marked point pattern whose marks are in the form of a data frame. So either: (1) Do not create point as a marked point pattern --- i.e. don't add a marks component when you create it,

Re: [R] Spatstat - Kest, envelope problem

P. S. I just noticed that you referred to the envelopes produced as confidence envelopes. They are ***NOT*** confidence envelopes!!! They are critical envelopes. This is a very different concept. Repeat after me, 50 times: critical envelopes, critical envelopes, critical envelopes,

Re: [R] spatstat kstest

It is also acceptable (and often useful) to ask on the R-Sig-Geo list. Bug reports should of course be sent to Prof. Baddeley. cheers, Rolf Turner On 02/08/2013 12:46 AM, Rui Barradas wrote: Hello, Try the following. maintainer(spatstat) [1] Adrian Baddeley

Re: [R] spatstat kstest

Thanks. On Thu, 07 Feb 2013 11:46:36 + Rui Barradas ruipbarra...@sapo.pt wrote: Hello, Try the following. maintainer(spatstat) [1] Adrian Baddeley adrian.badde...@csiro.au Hope this helps, Rui Barradas Em 07-02-2013 02:05, Hiroshi Saito escreveu: Dear sir, Which

Re: [R] Spatstat: Mark correlation function

On 24/11/12 06:36, AMFTom wrote: I normally use the following code to create a figure displaying the mark correlation function for the point pattern process A: M-markcorr(A) plot(M) I have now started to use the following code to perform 1000 Monte Carlo simulations of Complete Spatial

Re: [R] SPATSTAT: Minimum points for a Ripley K to be sensible?

Sebastian Pucilowski s.pucilow...@student.unimelb.edu.au writes: What are the minimum number of points in a point pattern before a clustering analysis using a Ripley K function loses any meaning? It depends what are your definition of `meaningful'. The K-function doesn't become meaningless

Re: [R] Spatstat users.

Yes I am an active spatstat use On Wednesday, 28 March 2012, Rolf Turner wrote: 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

Re: [R] Spatstat users.

No I am not an active spatstat user. David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Spatstat users.

Further to my cri de coeur of yesterday, asking people to indicate whether they are active spatstat users. I've had quite a few replies but I'd like a lot more. Please answer; it'll only take you a few seconds. Just reply to this email with a Yes (yes I am an active spatstat user) or No (no I

Re: [R] Spatstat - coordinates in observation window

On 16/03/12 23:02, Lucie V wrote: Dear R users, I wish to run spatial point pattern analysis (e.g. pair correlation function, mark correlation function) for which I need to create an observation window (window=owin) from which the spatial analysis is generated. The command I used to create

Re: [R] Spatstat - problem with which.marks and as.pp

Basically the problem is that you have a *VERY* old version of spatstat. Upgrade! Comments interpolated below. On 05/12/11 21:52, karajamu wrote: I forgot to change the header, so I guess no one read my mail. That's why I'm trying it again... Hello everbody, I am new to this mailing

Re: [R] spatstat = owin + image

Hej, Thanks for the idea. I will try and see if I can make progress. + I appologise for the shortcut explanations. I'll dedicate myself to greater accuracy in future. tkd -- View this message in context: http://r.789695.n4.nabble.com/spatstat-owin-image-tp3837023p3843284.html Sent from the R

Re: [R] spatstat = owin + image

On 24/09/11 02:32, tkdweber wrote: Dear Community I am at my wits end and seek advice. My wish is to plot coordinates (x,y in WGS84_UMTS for the ones interested) of sampling points. This I can do by the standard spatstat prodcedure via owin. You appear to be rather confused from the very

Re: [R] spatstat regression troubles

On 16/04/11 15:50, Gregory Ryslik wrote: 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

Re: [R] spatstat regression troubles

Hi Mr. Turner, You are correct that I am confused a bit by the RCode. Basically, I have 3772 observations of data and only about 500 of them correspond to where people exist. For the other observations, I just have the covariate values so I thought that this was appropriate. Thus, where

Re: [R] spatstat regression troubles

On 17/04/11 02:17, Gregory Ryslik wrote: Hi Mr. Turner, You are correct that I am confused a bit by the RCode. Basically, I have 3772 observations of data and only about 500 500? You said 944 previously. Doesn't really matter, but. of them correspond to where people exist. For the

Re: [R] spatstat regression troubles

On 17/04/11 02:17, Gregory Ryslik wrote: Hi Mr. Turner, You are correct that I am confused a bit by the RCode. Basically, I have 3772 observations of data and only about 500 500? You said 944 previously. Doesn't really matter, but. of them correspond to where people exist. For the

Re: [R] spatstat sub-selection

You can't specify it in rthin, which I can understand for the simple reason that this violates the independence of the sampling. If you put a fixed limit on your resampling, the chance of being selected depends on when some point is selected. You could do it by hacking the object : rthin.exact -

Re: [R] Spatstat, markcorr, max. radius limited??

On 21/11/2009, at 12:01 AM, rudi1...@gmx.de wrote: Hello, could please somebody help me. I want to apply the mark correlation function but for radii up to 75 meters (in 75 individual 1m steps). Unfortunately, There is a sensible default for the values of the argument r at which the mark

Re: [R] Spatstat: xy binary data into mask type to use in owin(mask=)

If 'X' is a data frame containing columns 'x', 'y' and 'value, try m - with(X, tapply(value, list(y,x), all)) z - im(m, xcol=sort(unique(X$x)), yrow=sort(unique(X$y))) w - as.owin(z) This will only work if the x, y values form a rectangular grid in some order. Adrian

Re: [R] Spatstat: xy binary data into mask type to use in owin(mask=)

In the first instance, questions about contributed packages should be addressed to the package maintainers rather than to the R-help list. On 20/10/2009, at 2:06 AM, Javier PB wrote: Dear users, I am trying to export polygons from Arcmap into Spatstat to run some simulations using functions

Re: [R] Spatstat

Richard Chirgwin wrote: Hi all, I am trying to install Spatstat on OpenSUSE 11.1. install.packages(spatstat, dependencies = TRUE) fails on the basis of various compiler packages (full message below). I have gcc version 4.3.2, which should include gfortran and g++ - so I'm not

Re: [R] Spatstat

Peter Dalgaard wrote: Richard Chirgwin wrote: Hi all, I am trying to install Spatstat on OpenSUSE 11.1. install.packages(spatstat, dependencies = TRUE) fails on the basis of various compiler packages (full message below). I have gcc version 4.3.2, which should include

Re: [R] Spatstat - K2 index

On 2/12/2008, at 4:30 AM, Gough Lauren wrote: Hi all, I'm using spatstat to investigate the spatial structure of an arid shrub population. The first-order intensity of my data does not appear to be homogenous, so I would like to use inhomogeneous techniques. I realise there is a

Re: [R] Spatstat - K2 index

Gough Lauren wrote: I'm using spatstat to investigate the spatial structure of an arid shrub population. The first-order intensity of my data does not appear to be homogenous, so I would like to use inhomogeneous techniques. I realise there is a inhomogeneous K-function available in spatstat,

Re: [R] Spatstat help - quadratcount query

2008 05:34 To: R-help Forum Cc: Gough Lauren; Rolf Turner Subject: Re: [R] Spatstat help - quadratcount query Gough, Lauren wrote: 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

Re: [R] Spatstat help - quadratcount query

Gough, Lauren wrote: 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. However, when I look at the data for quadrat counts it seems the function is not keeping the size of the

Re: [R] Spatstat - Several density plots using the same scale

Arthur Weiss wrote: Hi everyone, I am using the package spatstat for ploting kernel maps of my data. It is a marked point pattern, the result of mosquito surveillance in a area in a week. For each trap, the number of individuals captured is the mark of the point. plot(density(X,

Re: [R] Spatstat - Several density plots using the same scale

On 10/10/2008, at 4:48 AM, Arthur Weiss wrote: Hi everyone, I am using the package spatstat for ploting kernel maps of my data. It is a marked point pattern, the result of mosquito surveillance in a area in a week. For each trap, the number of individuals captured is the mark of the point.