similar to: Cross K Ripley's function and "spatio-temporal interaction power"

Dear sirs, I'm using the Splancs package to compute standard errors of the estimates of a spatio-temporal k function. When I use as spatial and temporal distances too long vectors (respectively 60 and 80 entries) for a dataset of 1000 observations, R gives me the message Error: cannot allocate vector of size 180000 Kb Reached total allocation of 512 Mb. I ran the function memory.size()

Dear R-listers, I would like to carry out a very basic descriptive analysis of homicides rates in Italy, taking into account both the spatial dimension (103 provinces) and the temporal dimension (10 years), but no covariates. In practice, what I would like to do is to describe spatio-temporal variation of homicide rates, identifying those combinations of province-year where the homicide rate

Dear R-helpers, I have 4-dimensional atmospheric data (x,y,z,t), which I want to analyse on spatio-temporal diversities. As far as I understand there only exists the possibility to construct time series as two-dimensional matrices (mts). For the moment, I hold it in different objects: 1. a four-dimensional array for the spatial related analyses 2. a two-dimensional mts timeserie, which was

I have a dataset like this: sites years Var1 Var2 1 1960 505 3.013833 1 1961 533 4.118784 1 1962 609 14.96386 1 1963 465 -3.74409 1 1964 837 41.70164 1 1965 727 29.53478 2 1960 493 3.269235 2 1961 535 5.386015 2 1962 608 16.26244 2 1963 469 -2.09736 2 1964 830 42.01942 2 1965 715 29.92867 3 1960 489

Dear R-Users, I have a model with a latent variable for a spatio-temporal process. I would like to use EM algorithm to estimate the parameters. Does anybody know how to implement the algorithm in R? Thank you very much in advance, Michela

Dear R-Users, I'm using SSPIR package for a spatio-temporal application. Is it possible to modify the structure of the involved matrixes (Fmat, Gmat, Vmat,Wmat)? I want to create a model like this #y(t)=k*theta(t)+epsilon(t) #theta(t)=h*theta(t-1)+eta(t) #epsilon(t) N(0,V) V=sigma2*I #eta(t) N(0,W) W=sigma2_eta where the state variable theta has dimension 1(p=1) and at

Hi all, I'm using RScalapack library for parallelizing some heavy matrix operations required by MCMC methods for spatio-temporal models. The package reference manuals (dated 2005) states that the library needs LamMPI to work but we have a Linux Cluster with OpenMPI. We have found (http://cvs.fedoraproject.org/viewvc/devel/R-RScaLAPACK/) a patch for OpenMPI but we are wondering if in the

> library(SpatioTemporal) > library(plotrix) > library(maps) > palay.cov<-read.csv("C:/Users/BEDANA-PC/Desktop/STThesisWD/Thesis Data/palay.covar.csv") > palay.o<-read.table("C:/Users/BEDANA-PC/Desktop/STThesisWD/Thesis Data/palay.obs.txt") > palay.obs<-as.matrix(palay.o) > palay.stc<-read.table("C:/Users/BEDANA-PC/Desktop/STThesisWD/Thesis

Greetings I'm using R 2.8 with recent (last month) versions of the packages I need to use at present. I'm interested in examining hierarchical spatio-temporal patterns in a data set. The data consist of 94 points (X, Y, UTM coordinates) at which catch rates for a fish were recorded and there are also estimates of prey available for these fish at the same locations.

The NOAA National Centers for Coastal Ocean Science is hiring a Computational Ecologist, a statistical/computational ecologist with experience fitting advanced spatial models to marine wildlife survey data (e.g., seabirds and marine mammal transects, fisheries trawl surveys) in R and other statistical languages. This is a full-time, long-term stable contract position. We are looking for an

Hi, I've got a new dataset which I don't know how to analyze with R. My knowledge about R is limited for this kind of problem. I've tried to find a solution with some spatio temporel packages and lme/lmer functions, but didn't find any similar example. I've got 10 locations on a coast on which the number of walrus was noted. I've got a distance measure between each

I have a data assimilation problem that might be amenable to the use of GAMS, but I am not sure how feasible it is to implement. I was told the R mailing list was a great resource. My observations are spatiotemporal salinity in the San Francisco Bay at a number of instruments over a few days. The thing that I want to fit is the initial condition for a salt transport model at the beginning of this

Hi, I started using R about 3 months ago so please excuse my ignorance. I have two datasets. The first consists of 247 fixed position mosquito traps that were serviced weekly for 23 weeks. The second is dengue incidence data indexed to home address for the same time period. I would like to see if there is a space time relationship between mosquito catch rate and dengue incidence. So far I have

Dear R users, I'm currently trying to make level plots of a longitudinal study of the spatio-temporal spread of a plant disease in a field, but the results of the color key assignment isn't what I expect. Here's more info. I recorded the level of a disease on an ordinal scale from 0 (no disease) to 9 (dead plant) on a grid that has 9 rows and 14 columns. The disease level on each

Thanks to Duncan's recent work on SJava, I've got a (so far) stable platform to put together the R-Orca interface (dynamic graphics, including brushing, subseting, etc). While the first version will be looking at installation issues and study basic user interfaces for providing a few working examples, it will be important for this particular package to look at how one integrates a

Hello everybody! This is my first mail so I'll write a couple of lines of self-introduction. My name is Massimiliano, I'm from Italy and I'm studying Mathematical Engineering. I started using R in my Statistics course and have to use it to make a project which I'll discuss at the end of the course. The problem I'd like to ask you about follows. Suppose I have imported a

Dear all, It's possible to obtain the residuals with AR time series function. Thanks in advance, Massimiliano Talarico Pianificazione Commerciale - Area Crm Unicredit Xelion Banca S.p.A. Via Pirelli 32 - 20124 Milano Tel.: 02 67360 525 Fax: 02 67738 525 www.xelion.it [[alternative HTML version deleted]]

Hi all, I have problems with years of dates using "chron" package. I don't understand why R by this istruction: > dates("01/02/29",out.format="d/m/year") [1] 02/Jan/2029 > dates("01/02/30",out.format="d/m/year") [1] 02/Jan/1930 reads "29" as 2029 and "30" as 1930. How could I change to read "00" to

Dear all, I have the following R dataframe: set.seed(11) (df <- data.frame(ID=rep(1:10,1:10),a=factor(sample(1:4,55,rep=T)) )) where ID is an identification code. I need to create a new variable "b" in which I would paste the full group of "a" variable according to ID variable. For example for ID=2 "b" should be 13, for ID=3 "b" should

Hi, I have a quite "simple" setup for a particular customer that loves redundancy and failover. PDC + BDC with LDAP Passwords on two 389-ds in multimaster node + several samba member servers Actually pointing singularly on both the systems everything works great. As soon as I modify my passdb backend line from the single form to the form containing both backends that is from passdb