n Problem
and you could use the R package 'TSP' for example..
Adrian Baddeley
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
clusters="LGCP")
See the help for function 'kppm'.
It is also extensively documented in Chapter 12 of the spatstat book
(book.spatstat.org)
Adrian Baddeley
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@
rof Adrian Baddeley DSc FAA
John Curtin Distinguished Professor
Department of Mathematics and Statistics
Curtin University, Perth, Western Australia
I work Wednesdays, Thursdays and Fridays
From: Rolf Turner
Sent: Wednesday, 15 April 2020 11:55 AM
To: Mauro Pedone
://github.com/spatstat/spatstat>
Note that in the envelope() command above, you don't need the argument 'j="B"
and it is ignored.
regards
Prof Adrian Baddeley DSc FAA
John Curtin Distinguished Professor
Department of Mathematics and Statistics
Curtin University, Pert
as.linnet.SpatialLinesDataFrame() from the package 'maptools'.
Prof Adrian Baddeley DSc FAA
John Curtin Distinguished Professor
Department of Mathematics and Statistics
Curtin University, Perth, Western Australia
From: Rolf Turner
Sent: Friday, 21 June 2019 10:08 AM
To: Rolando Valdez
Cc
in(as.solist(P))
V <- lapply(P, onesplit, U=U)
Z <- Reduce(intersect.tess, V)
return(Z)
}
onesplit <- function(X, U) tess(tiles=list(A=X, NotA=setminus.owin(U,
X)), W=U)
Prof Adrian Baddeley DSc FAA
John Curtin Distinguished Pro
n you'll get NA's.
I suggest you do something like
X <- unique.cases.ppp
W <- intersect.owin(Window(X), Window(pov.f))
Y <- X[W]
m1 <- ppm(Y ~ pov.f)
Adrian Baddeley
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Then you could do
lik1 <- sum(log(lam1[X])) - integral(lam1)
lik2 <- sum(log(lam2[X])) - integral(lam2)
and compare the likelihoods.
Adrian Baddeley
From: Cenk ���Z <ci...@anadolu.edu.tr>
Sent: Friday, 15 December 2017 8:23 PM
To: Rolf
tstat package.
If you decide you need a pixel image instead, then just use 'as.im' to
convert the function
to an image.
Hope this helps.
Adrian Baddeley
Spatstat author
Prof Adrian Baddeley DSc FAA
John Curtin Distinguished Professor
Department of Mathematics and Statis
You can use the spatstat function 'edges' to extract the edges of a polygonal
window.
Example:
W <- letterR#polygonal window
E <- edges(W)
X <- runifpointOnLines(20, E)
plot(E)
plot(X, add=TRUE)
Prof Adrian Baddele
the plots
of summary functions use the "recommended range" of distance values rather than
the full "available range". To override this, when you plot the summary
function, use 'xlim' to control the range of 'r' values. (xlim should be an
argument in the 'plot' call).
Adrian
NOT convex (it's a kidney shape).
So the statement above applies to this window.
> I know this isn't quite an "r-implementation" question, but: anyone know if
> isotropic correction is unbiased for an odd shape like a country-outline
> (in this case Sweden)? Or should i use the `bo
mum distance between any two
points. Another example is the equilateral triangle.
Adrian
From: Tina Cormier <tco...@gmail.com<mailto:tco...@gmail.com>>
Date: Tuesday, 12 July 2016 8:02 pm
To: Adrian Baddeley
<adrian.badde...@curtin.edu.au<mailto:adrian.badde...@curtin.edu.au
minimal circle) is the point at the centre
of the triangle.
The circumradius is s/sqrt(3), not s/2.
Adrian
Prof Adrian Baddeley DSc FAA
Department of Mathematics and Statistics
Curtin University, Perth, Western Australia
From: Tina Cormier <tinaacorm...@gmail.com
(centre=z, radius=r)
return(w)
}
circumcircle.ppp <- function(x, ...) {
circumcircle(convexhull(x))
}
Adrian Baddeley
[[alternative HTML version deleted]]
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
not cover
geostatistical methods, but you can compute shortest-path distances in the
network and construct your own Moran¹s I.
Adrian Baddeley
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
should be equivalent as regards their fitted intensity, etc.
Adrian Baddeley
Prof Adrian Baddeley DSc FAA
Department of Mathematics and Statistics
Curtin University, Perth, Western Australia
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Please read chapter 9 of the spatstat book for a detailed explanation of what
the model formulas mean (and what they don't mean). You can download this
chapter for free from the companion website
Adrian
Prof Adrian Baddeley DSc FAA
Department of Mathematics and Statistics
Curtin University
We will continue this conversation offline
and report the results to R-Sig-Geo when done.
Adrian
Prof Adrian Baddeley FAA
Curtin University
On 25/06/15 01:39, Harold-Jeffrey Ship wrote:
I have Ubuntu 12.04 Precise. I just installed R and want to install
spatstat. This is the R information
each line by a fine grid of points, each point having a
'weight'
equal to the length of line that it replaced. Then apply density.ppp to these
points with weights.
The result is the density.psp output.
Adrian Baddeley
spatstat author
Prof Adrian Baddeley FAA
Curtin University
Adrian Baddeley FAA
Curtin University
From: Ivan Palmegiani [pan.sapiens...@gmail.com]
Sent: Wednesday, 18 February 2015 12:32 PM
To: Adrian Baddeley
Cc: r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Calculate shortest distance between points belonging
)
Prof Adrian Baddeley FAA
Curtin University
From: Ivan Palmegiani [pan.sapiens...@gmail.com]
Sent: Wednesday, 18 February 2015 10:11 PM
To: Adrian Baddeley
Cc: r-sig-geo@r-project.org; r.tur...@auckland.ac.nz; Ege Rubak
Subject: Re: [R-sig-Geo] Calculate
of
'marks' to 'factor')
Adrian Baddeley
spatstat author
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
and you want to
convert them to kilometres,
rescale(W, 1000)
or better
rescale(W, 1000, unitname=km)
The result of 'rescale' is physically equivalent to the original, but is
expressed in different units.
See help(scalardilate) and help(rescale).
Prof Adrian Baddeley FAA
University
.
In the next version of spatstat (1.37-1) we will include a partial
implementation
of the fast C code for hybrids, which will work for your example.
Prof Adrian Baddeley FAA
University of Western Australia
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
Adrian
Prof Adrian Baddeley FAA
University of Western Australia
From: Rolf Turner [r.tur...@auckland.ac.nz]
Sent: Saturday, 17 May 2014 8:29 AM
To: Hossain, Md
Cc: r-sig-geo@r-project.org; Adrian Baddeley
Subject: Re: [R-sig-Geo] Warning message: points were
if the window is a rectangle, there are only two possible marks, etc)
but no other major problems.
Thank you for communicating this question. It has spurred me to get on with
implementing a faster version of pcfcross in 'spatstat'.
Adrian Baddeley
On 28/01/14 10:39, Noam Ross wrote:
I have
, an idiom like
ttt$tiles - xxx
is dangerous, and in this case, will lead to problems, because it is no longer
true that
the union of tt$tiles is equal to tt$window. The correct way to create a
changed tessellation
is something like
tess(tlles=xxx, window=www).
Prof Adrian Baddeley FAA
. Lawrence and R. Milne and G.M. Nair},
title ={On tests of spatial pattern based on simulation
envelopes},
journal = {Ecological Monographs},
year = 2014,
note = {In press}
}
Prof Adrian Baddeley FAA
University of Western Australia
)
K1 - K[[1]]
K2 - K[[2]]
K3 - K[[3]]
Kpool - eval.fv( (K1 * nn1[1] + K2 * nn1[2] + K3 * nn1[3])/sum(nn1))
Finally
mod - lgcp.estK(Kpool)
In future, mppm() will be able to do this.
Prof Adrian Baddeley FAA
University of Western Australia
scale
so that competitors could not re-analyse the data.
Adrian Baddeley
author of 'spatstat'
Prof Adrian Baddeley FAA
University of Western Australia
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
'union.owin'
and 'dilation.owin'.
It has been fixed in spatstat 1.34-1 which has just been submitted to CRAN.
Adrian Baddeley
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
In spatstat 1.34-0 there appears to be a bug in union.owin which also affects
dilation.owin.
A fix will be distributed shortly.
Until this is fixed, the fastest way to do what you want is to convert the
window
to a pixel mask (using as.mask) before applying the dilation.
Adrian Baddeley
is performed
using the package 'polyclip' which is free open source. This means that all
uses of spatstat are now
covered by free open source licences.
Adrian
Prof Adrian Baddeley FAA
University of Western Australia /and/ CSIRO Computational Informatics /and/
University of Aarhus
... Are there any NA's in your polygon coordinates?
(Spatstat does not currently check for NA's - I will fix that!)
Can you plot the coordinates using polygon()?
If everything seems ok, please send me the data and I'll try to figure out what
is happening.
regards
Adrian Baddeley - spatstat coauthor
a larger window. This is one of the arguments
made by Diggle and Rowlingson (1994).
However, automatic methods for selecting the smoothing bandwidth (such as
'bw.relrisk' or 'bw.diggle')
will give a different answer in a larger window.
A
Prof Adrian Baddeley FAA
University of Western Australia
by fitting to the data point pattern 'X') but I'll
leave that until later.
Adrian Baddeley
Prof Adrian Baddeley FAA
University of Western Australia /and/ CSIRO Mathematics, Informatics
Statistics
Mail: cet.uwa.edu.au Skype: adrian.baddeley
want to increase the space between the panels in the plot,
use the argument 'mar.panel' described in help(plot.fasp).
A
Prof Adrian Baddeley FAA
University of Western Australia /and/ CSIRO Mathematics, Informatics
Statistics
Mail: cet.uwa.edu.au Skype: adrian.baddeley
Prof Adrian Baddeley FAA
University of Western Australia /and/ CSIRO Mathematics, Informatics
Statistics
Mail: cet.uwa.edu.au Skype: adrian.baddeley
From: Rolf Turner [r.tur...@auckland.ac.nz]
Sent: Friday, 28 June 2013 7:34 AM
To: Carsten
. my laptop will accept npoly =
25,000).
To avoid this error, you could try setting
spatstat.options(checkpolygons=FALSE)
which will simply turn off the checking. Do this if you are confident that the
polygons do not overlap.
Good luck!
Adrian Baddeley
40 matches
Mail list logo