similar to: about"Riply's K function"and "envelope"

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,

Hello, I am writing cause I have a problem when try to create confidence envelopes in spatstat. I have a point data representing firms in my study area. When I try to create confidence envelopes the error message turns up: [> p <- envelope(point, fun=Kest, nsim=99) Error in marks.ppp(Y, dfok = FALSE) : Sorry, not implemented when the marks are a data frame] Could you write me how can I

Hello! I am writing you because I could not plot the confidence envelopes for functions Jest, Jcross, Jdot, Jmulti, and L, using the Spatstat package. I have already understood how to do that for Kest or Jest, that is: JEnv <- plot(envelope(PPPData, Jest)) Where PPPData is my ppp object. However, for Jcross I must specify the two marks I want to analyse. That is, usually I would get the

Hi R users, This query is regarding the use of the 'envelope' function in Spatstat. My data can be represented as a point process with CONTINUOUS marks: points <- ppp(x=x,y=y, marks=m, window= wind) However the marks are alignments (lines), and so have to be treated differently to normal scalar marks. Hence to create a mcf object with the appropriate test function for alignment

Hi I'm working with the function Kest in the package SpatStat (under LINUX with R 2.1.0). In order to evaluate the statistical significance of my point pattern I'm doing 999 Montecarlo replications. The script that use the Kest function runs OK for most of the different point patterns that I have but for a particular point pattern, which have only 17 points, it runs until the 34th

Hi I looked, but I didn't find it: I need an array [0 to 101] where each element is a list (the result of Kest in spatstat). I.e. I want to do: A[2] <- Kest(pp1) A[3] <- Kest(pp2) ... A[101] <- Kest(pp100) How can I create A ? Rainer -- Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Department of Conservation Ecology University of Stellenbosch Matieland

Hi I have the following loop and would like to vectorize it. Any ideas if it is possible? Thanks, Rainer Tha Loop: for (i in 2:Result$NoSims) { ppp <- runifpoint(Result$NoPlants) K <- Kest(ppp) Result$LSim[i,] <- sqrt(K$iso / pi) - K$r CM <- (Result$LSim[i,] * Result$LSim[i,]) / abs(K$r[2] - K$r[1]) Result$SigCM[i] <- sum(CM, na.rm=TRUE) print(i) flush.console() }

Hello, What are the minimum number of points in a point pattern before a clustering analysis using a Ripley K function loses any meaning? I haven't been able to find much comment in the literature about this --- instead case examples typically use 100-200 points as a minimum. Regards, Sebastian Pucilowski

Dear R-helpers, After putting my somewhat sparse (21 points) data into a 'ppp' format, I am unable to calculate the K function using Kest. I get the error "some 'x' not counted; maybe 'breaks' do not span range of 'x' ". So I tried to calculate the K function using the swedishpines dataset. Of course Kest(swedishpines) works fine but when i break it

Hi All! As a new useR, I am currently working in R 2.5.1 to produce Ripley's L (the linear version of Ripley's K) for a spatial point dataset (originally a ESRI shapefile). I have coerced the data to a ppp object (because I was originally using 'kest' in the spstat package to produced Ripley's K) but I need to convert my spatial data to a spp object to be able to use

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 have a loop which is doing time consuming calculations and I would like to be able to have some feedback on where it is in it's calculations. I tried to simply show the counter variable in the loop, but id doesn't work as all display seems to be delayed until the loop is completed. Is there any way of displaying the progress of a loop? Rainer The loop: for (i in

Hi, R-users, I met a problem: Items:[Anna 'moi =) akku loppu joskus 4ltä. Kestää kauan nää..'\tAmer, Tuusula (0:20)\t20\t12\t16\t00\t00\t11]/Anne 'Ei jakoa,uus päivä muistio et 4n niin peruin. Hups'\t (0:16)\t0\t12\t18\t00\t00\t11/Elina 'Konsertissa. En tod. vastaa teille'\tEtu-Töölö, Helsinki (2:40)\t24\t12\t18\t00\t00\t11 I want to parsing the above data into the

Hello, Somehow I cannot control my envelope plots on Spatstat package. I would like to plot an envelope with no legend and also I would like to label the plot as such: >plot(envelope(data, nsim=39, main="My K Plot", legend=FALSE, xlab="Distance") But somehow the "main", "legend" and "xlab" calls do not have any effect on the plot. Can someone

Hi, Can anyone please tell me how to get ".probe" function in "struct xenbus_driver" called? I want to do one simple test of event channel communication between peer modules of domU and dom0. I wrote one backend and one frontend driver (I tried registered as misc/blk/input device) where I put in all my xenbus routines, but cannot get .probe called after insmod my module, I

Hello. There are so many sound files in /var/lib/asterisk/en. Is there an easy way to let me play back all of them one by one while I am watching CLI to see the current file name? Thanks for help. -- Jian Gao IT Technician SJ Geophysics Ltd. <http://www.sjgeophysics.com> jian.gao at sjgeophysics.com <mailto:jian.gao at sjgeophysics.com> Tel: (604)582-1100

Our computing services people have recently reconfigured the Sun compilers structure on our network, and I can no longer install, on our Sun/Solaris machine, a certain package (spatstat) with which I am heavily involved. (The same problem may well pertain to other packages as well; I haven't yet experimented except with spatstat.) Strictly speaking I can ***install*** it, but when I try to

Hi all, I am using "spatstat" to investigate the spatial structure of some plant populations, and I want to detect these patters with IPP and a Thomas process based on pair-correlation function. I know the function "pcfinhom" is available to characterize the IPP, but I have no idea about how to use the pcf with Thomas process? Additionally, generating simulation envelopes

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 sure what to do! Richard * Installing *source* package ?deldir? ... ** libs gfortran

Hello,   I have a dataframe: dd <- data.frame(b = c("chr2", "chr1", "chr15", "chr13"),         x = c("A", "D", "A", "C"), y = c(8, 3, 9, 9),         z = c(1, 1, 1, 2))   >dd       b x y z 1  chr2 A 8 1 2  chr1 D 3 1 3 chr15 A 9 1 4 chr13 C 9 2 Now I want to sort them according column "b", but only its