Dear R-users,
I'd like to use a corStruct matrix in fitting a binomial model, but I
need to define it using a pre-calculated distance matrix (instead of
feeding it with spatial coordinates). Do you have any idea on how to do
this in a simple manner?
So far the only way I could imagine was to defi
On Fri, 26 Aug 2011, Roger Bivand wrote:
On Thu, 25 Aug 2011, Karl Ove Hufthammer wrote:
Dear list members,
Is it possible to define which FID values should be used when saving a
SpatialPointsDataFrame to a GML files using writeOGR (rgdal package)?
The FID values are stored as (and F1, F2
Mathieu,
as stplot for STIDF's is still pretty "bare", or simple, you might as
well work directly with xyplot from lattice. Having your data, it turns
out you only have a limited number of dates, the following seems to work:
library(lattice)
xyplot(y_sref~x_sref|time, as.data.frame(ciotat.t),
On Thu, 25 Aug 2011, Karl Ove Hufthammer wrote:
Dear list members,
Is it possible to define which FID values should be used when
saving a SpatialPointsDataFrame to a GML files using writeOGR
(rgdal package)?
The FID values are stored as (and F1,
F2, …). If I insert a new row in the SPDF,
No worries. good luck!
On 8/26/2011 10:45 AM, Michał Kwieciński wrote:
Learning about existence of "band1" was what I actually needed. Apart
from one bigger spatial econometrics project, I haven't done much work
with R yet, so I still lack some knowledge on basic operations - I
will try to reme
Hi Reza,
I wrote a simple function that creates choropleths in R. You can check out
the code on my github repos and some documentation on this post.
Github repos:
https://github.com/darribas/darribas-R-tools/blob/master/choropleth.r
Blog post w/ docs:
http://dreamessenceonthemove.wordpress.com/20
Learning about existence of "band1" was what I actually needed. Apart
from one bigger spatial econometrics project, I haven't done much work
with R yet, so I still lack some knowledge on basic operations - I
will try to remember about the names() function next time. I thought
the problem would lie
If you want to see more info of the dataset, you can also view it with
str(map@data)
or
head(map@data)
Matt
On 8/26/2011 9:57 AM, Forrest Stevens wrote:
As you said, readGDAL will give you a SpatialGridDataFrame object, which
behaves much like a normal data.framet. Your best bet is to run
Dear Users,
I am fitting a spatial model a multilevel data (mouth-tooth-surface) of very
small prevalence. Therefore in my case every person will act as a map. I am
looking at a binary outcome of a tooth surface which is nested in the tooth and
nested in the mouth. So I have a surface-tooth-mou
As you said, readGDAL will give you a SpatialGridDataFrame object, which
behaves much like a normal data.framet. Your best bet is to run it through
names to see what data are available to work with:
names(map)
Like a regular data.frame this will give you the bands or data read in from
your ASCII
Hello,
I use the raster package for asci files
it works great and has a very hands on access to get statistics or do
any kind of raster operations
hth
a
MSc. Alexander Arpaci
Assistant researcher
Department of Forest - and Soil Sciences
Institute of Silviculture
Tel. +43-1-47654-4081
e-mail: al
Hi,
I am having difficulty in reading data for further use that is in
ASCII grid format (it's the UN population data from
http://sedac.ciesin.columbia.edu/gpw). I used readGDAL function which
leaves me with SpatialGridDataFrame class object, but apart from
creating simple image with
map = readGDA
Thank you,
this is what I need!
Van: adrian.badde...@csiro.au [adrian.badde...@csiro.au]
Verzonden: vrijdag 26 augustus 2011 10:42
Aan: r-sig-geo@r-project.org
CC: Quets Jan; adrian.badde...@csiro.au; r.tur...@auckland.ac.nz
Onderwerp: [R-sig-Geo] using 'fo
Hi,
I have a question about stplot.
Let's say I'd like to plot via multi-panel a series of plots separated by a
given time interval. How can I accomplish that?
Also, let's say this time interval is specified between thow given time
instances, for example: every week between 14/07/11 and 28/08/1
Dear Edzer,
my.STSDF <- as(my.STFDF, "STSDF") does exactly what it should. I have a
STSDF object with several columns of data that have NAs for one variable but
not for other variables. I wanted to keep only those observations in space
and time that had data (no NAs) for all variables (see complet
Jan Quets wrote:
> env is a envelope from spatstat where the argument 'savefuns' is set to 'TRUE'
> I would like to put the code at the end of the mail in a for loop like:
> for (i in 1:999){
> sims[i][[1]] = attr(env,"simfuns")$???
> }
According to `help(envelope)' in the package 'spatstat',
16 matches
Mail list logo