search for: kesting

Displaying 17 results from an estimated 17 matches for "kesting".

Did you mean: testing
2005 Sep 30
2
Creating an array [0 to 101] of a list
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
2006 Feb 27
2
about"Riply's K function"and "envelope"
hi! I did Riply's K function and envelope in the package "SPATSTAT". When the species number is less, it can work well and it is quickly. But when the species number is more(example:2000), it told me"memory limit". So I change "memory limit". But the speed is very slow. I took 13 hours to run the programe. Perhaps there are other methods that can do it. Or just
2005 Sep 01
2
SpatStat Kest - Error Message help
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
2013 Feb 07
2
Spatstat - Kest, envelope problem
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
2005 Oct 10
3
Vectorizing loop
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() }
2003 Oct 20
0
calculating K-function (Kest) in spatstat
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
2006 Jan 12
1
envelopes of simulations
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
2008 Mar 04
1
Ripley's K function within envelope--how to change the maximum distance?
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,
2005 Oct 10
5
Show Progress in loop
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
2004 Dec 20
11
winbind problems
Hello, I am running a Samba server (3.0.7) on a Suse 9.2 box. I have connected this server successfully to a Windows 2000 Active Directory (mixed mode). I have nsswitch.conf, krb5.conf configured and winbind seems to be running properly for the most part. With wbinfo I can get all of my user and group information. Problem is, it seems that at random times, the samba server just stops
2009 Sep 29
1
How to parsing data like this in R
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
2007 Nov 19
1
Convert from ppp to spp objects
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
2010 Apr 14
1
envelope in spatstat
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
2005 Jun 01
0
cor.test on spatial data
I have spatial data on various species of trees within plots. What I'd like to know is if their patterns are positively associated or not. Could I use cor.test on the values generated by Kest or Gest functions? or would it be more powerful to do a G i to j and then G j to i analysis instead? Thanks, Suzie __________________________________________________________________ Switch to Netscape
2009 Feb 20
0
Random labeling hypothesis in spatstat
Hi, I am trying to test the random labeling hypothesis with my data. I have a point pattern of bird nests and each nest has a mark that it is either predated or not-predated. I want to test the null hypothesis that the K distribution of predation events is equal to the K distribution of all events (nests). This is the code I am using: Kdif <- function(X,...,i) { Kidot <-
2012 Jul 02
1
SPATSTAT: Minimum points for a Ripley K to be sensible?
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
2004 Dec 21
1
Winbind problem revisited
Okay, I started over from scratch with my samba server rebuild, but I am still getting some weird issues. Here are my config files of importance: --------------/etc/samba/smb.conf-------------------- # Samba Configuration File [global] workgroup = WAYNE realm = WAYNE.LOCAL server string = Samba Server security = ADS password server =