similar to: Is it possible to do Ripley's K and L functions analysis with R?

Hello, I am trying to conduct a bi-variate ripley's k using the k12hat function in the splancs package (found here http://rss.acs.unt.edu/Rdoc/library/splancs/html/k12hat.html). Although, I receive an error when getting to the bboxx part of the code: poly <- list(x=c(fire4$X, nar4$X), y=c(fire4$Y, nar4$Y)) plot(seq(5,80,5), sqrt(k12hat(fire4), as.points(nar4),

Hello, I have used the univariate Ripley K function in R, but does anyone know if there is a bivariate function built in? I have two species that I am dealing with. Also, how might I add error bars into the graphs (univariate and/or bivariate)? Thank you, Karin Leiderman k_leiderman at hotmail.com Graduate Student/Research Assistant Department of Mathematics Univesity of New Mexico

Hi, I have a dataset with x and y coordinates and in each point I have an identity of point, in some cases I can have more then one identity by point. My dataset is something like this: > x <- rep(c(1:4),4) > y <- rep(c(1:4),c(4,4,4,4)) > area1 <- sample(factor(rep(c("a","b","c","d"),4))) > area2 <-

In trying to streamline various optimization functions, I would like to have a scratch pad of working data that is shared across a number of functions. These can be called from different levels within some wrapper functions for maximum likelihood and other such computations. I'm sure there are other applications that could benefit from this. Below are two approaches. One uses the <<-

Dear All, I have a collections of spatial data. I have to analyze pairs of these point patterns to test their spatial interaction. I was moving towards the cross K Ripley's function. The problem, however, are the following: 1) What is the best way to get a single real value that represents the interaction "power"? 2) How to obtain a value that even allows me to rank the pairwise

I have plotted Ripley's K function for a spatial point pattern for 12 plots, with 39 Monte Carlo simulations for complete spatial randomness (CSR). I would like to analyse these data as follows: I would like to know for which plots the Ripley's K function deviates from CSR at a number of confidence intervals, namely 90%, 95% and 99%. I realise that simulation envelope I have produced

Hi Is there anyone who can help a reporter (me) with the following: I want to use the bivirate k-function to determin if burgerbars are clustered around schools. To that purpose I’ve been told to use the bivirate k-function, but my program Arcview doesn’t contain that. Im stuck on a dead end – how to I for a start get the data into R? Michael Holm DICAR – Center for Anaytical Reporting

Hi--I'm trying to run a Ripley's K analysis on a point pattern within a window of 1 square km. The maximum distance I want to use is the diagonal of the window, around 1400m, run in 50m increments. The function defaults to a 250m maximum--I thought this was simply the way I had plotted it, but when I exported the data there was nothing over 250m. I've looked into seq and ginterval,

I am trying to port the splancs library for S++ version 4.0 to R. I have modified the original FORTRAN to use dllwrap/dlltool (e.g. removed common blocks). I have been able to generate dll's (splancs.dll) using dllwrap/dlltool and Ming32 g77 with a -mdll switch. As of yet I have not been able to generate a dll that can be loaded by R. Does anyone out there have any suggestions. An example

Is it possible to avoid the loop in the following function (or make the function otherwise more efficient) and can someone point me to a possible solution? (It would be great if hours could be reduced to seconds :-). # --------------------------------------------- RanEigen=function(items=x,cases=y,sample=z) { X=matrix(rnorm(cases*items),nrow=cases,byrow=F) S=crossprod(X-rep(1,cases) %*%

I am intersted in time-space clustering of local industry. I made the point file, polygons and time table then run 'stkhat' function in Splancs, but it generates only time k function. I noticed that it has only one ks value. ... Can anybody help me?. thanks in advance. countiespt <- read.table('d:/dissertation/dailynew/jun10/countiespt_ok.txt', sep=",")

Dear All, I have a codes which calculates the result of Ripley's K function of my data. I want to repeat this process 999 times. However, i am getting an error when i use the "for i in" function. Is there any way to repeat this analysis 999 times. Here are the codes i used ; data4 <- matrix(c(sample(id),data1),203,3) a <- data4[,1] random.case=data4[a==0,]

Hi! I want to use the package "SPLANCE" to do Ripley's K analysis I have a basic question about the package. I try to find how to do it, but I cannot. I hope you can help me. I donot know how to change my data form to "point data set". my data form is: > qumo[1:5,] code species dbh x y tag status branch 223 10312 QUMO 54.5 7.83 44.80 10306 alive

The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION' Other packages can be installed without problems, except of dse_2003.4-1 with a similar error message. Why? Operating System:

Hi folks. I have some old servers running Windows and Linux and I want to create a HPC cluster. However I wanna know if Xen can distribute a single virtual machine along multiple nodes to increase it''s performance. -- Allysson Steve Mota Lacerda Laboratório de Inteligência Computacional - LABICOMP Programa de Pós-Graduação em Engenharia Elétrica – PPGEE Universidade Federal de Minas

Hello I have a XP on my Pc ,and I would like to install the splancs package on my machine what is the step for implementing them on my library? --------------------------------- [[alternative HTML version deleted]]

Hi all, I'm trying to install the splancs contributed package on a Redhat 7.3 system with R 1.5.1 installed. I've installed the readline41-4.1-10 RPM which Redhat requires for R and it works great, thanks very much. However, when I attempt to R CMD INSTALL splancs_2.01-9.tar.gz as root, I get the following error message during the install: [snip] g77 -fPIC -O2 -m486

Hi, I would need some help with the splans package in R. I am using a Shapefile (downloadable at) http://rapidshare.com/files/215206891/Redlands_Crime.zip and the following execution code setwd("C:\\Documents and Settings\\Dejan\\Desktop\\GIS\\assignment6\\DataSet_Redlands_Crime\\Redlands_Crime") library(foreign) library(splancs) auto_xy<-read.dbf("Auto_theft_98.dbf")

Hi, I have a problem that concerns combination of the package Maptools and Splancs I have 2 shapefiles that i want to manipulate (one of type point and one polygon).I import them in R using Maptools but then i can't estimate a quartic Kernel using Splancs. The package doesn't recognize the shapes (invalid points and poly argument).I don't know if this is an easy task but i have

Dear all,  Who will bail me out. Iam using R with S-Splancs. Anytime i typed in the syntex, an error will appear eg setwd("C:\\TEMP) >dat <- read.table(cheshire_fmd.cvs",header=TRUE, sep=",") > dat.<-read.table(''chesire_fmd.cvs'',header=TRUE,sep='',) Error: unexpected symbol in "dat.<-read.table(''chesire_fmd.cvs" >