similar to: Quadratcount() plotting in R spatstat

Hi, I know there are older threads discussing the quadratcount function in spatstat. Unfortunately, I could not find a solution to my problem there. I'm analyzing a point pattern in an irregular polygonal window. Both the window (an entire country) and the points are projected using WGS84. When I do quadratcount with only one quadrat for the entire country it holds all my 154 points:

Hi all, I am using quadratcount in spatstat to divide a window containing a point pattern into a grid of quadrats containing the intensity of points in each quadrat. My data is in UTM co-ordinates. My window is defined as follows: >p15<-ppp(x,y,window=owin(c(341710,342100),c(3126465,3126780)),marks=NUL L, checks=TRUE) Giving me a distance of 390m in the 'x' direction and 315m in

Hello everyone, I am employing the waba method from the multilevel package for obtaining a within-group correlation (Description: http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/multilevel/man/waba.html). Does anybody know a way or a calculation for obtaining a significance value for that correlation? And another question: Does anybody know whether it is possible to save individual

Hello! I just noticed new link on R wiki on R galleries and wanted to share this info with YOU! - R graphical manuals (this is awesome page as there are all help pages of all packages on CRAN and probably even more and all graphics examples are displayed! - more than 8000 images!) http://bg9.imslab.co.jp/Rhelp/ This is a very nice addition to already existing R graph and movies galleries - R

I have photographs of plots that look like so: http://r.789695.n4.nabble.com/file/n4631960/Untitled.jpg I need to divide it up so each circle has an equal area surrounding it. So into 20 equal segments, each of which contains a circle. Quadratcount is not sufficient because if I divide it up into 36 equal quadrats, some quadrats do not contain one of the circles. I'm not even sure how to

Hi, Using the quadrats and quadratcount functions in spatstat, when I go to plot either of these, I get the quadrats coloured by their identity, i.e., using a color ramp applied to the sequence of quadrats. This only happens when the quadrats are applied to an owin which is polygonal, i.e., when I have an irregularly shaped study area. There doesn't seem to be any obvious way to

Hi, With the following code i can divides window into quadrats and counts the numbers of points in each quadrat. library(spatstat) X <- runifpoint(50) quadratcount(X) quadratcount(X, 4, 5) quadratcount(X, xbreaks=c(0, 0.3, 1), ybreaks=c(0, 0.4, 0.8, 1)) qX <-? quadratcount(X, 4, 5) plot(X) plot(qX, add=TRUE) But I want to mark each? quadrats? and select/ extract only those points by

Hello, I am using the following code to create ppp files from csv data and map shape files, but I am getting some errors which I have been unable to fix by searching them online: library(spatstat) library(maps) library(maptools) NYC2<-readShapePoly("nybb.shp") # this is a map of the NYC boroughs without waterways and no census tract divisions (but it does include lines separating

Ok guys, as requested, I will add more info so that you understand why a simple vector operation is not possible. It's not easy to explain in few words but let's see. I have a huge amount of points over a 2D space. I divide my space in a grid with a given resolution,say, 100m. The main loop that I am not sure if it's mandatory or not (any alternative is welcomed) is to go through EACH

Hi, I'm new to this mailing list so apologies if this is too basic. I have confocal images 512x512 from which I have extracted x,y positions of the coordiates of labelled cells exported from ImageJ as a.csv file. I also have images that define an underlying pattern in the tissue defined as areas of different pixel values 0 or 255 (also 512x512) I've exported these images as .txt

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Hello R Users, I am trying to do a quadrat count defined by covariate properties in spatstat. I have read my elevation raster into R (from ascii) and converted to class "im" for use in spatstat. Now I have point data of class "ppp" which window is the same extent as the elevation image. I can conveniently plot the image and the points on the image as follows:

Is there any way to add graphics (e.g. thumbnails in examples) to Rd files? A picture is worth a thousand words! Thanks John [[alternative HTML version deleted]]

I thought I came across a site that had graphics from all CRAN packages, but I can't find it. Anyone know of such a site? (I don't mean this: http://addictedtor.free.fr/graphiques/) - Tom Tom Short EPRI Solutions, Inc. -- View this message in context: http://www.nabble.com/website-with-graphics-from-all-packages--tf2666416.html#a7436096 Sent from the R help mailing list archive at

Hi everyone,   I'm trying to do an Logistic Politomic Regression in R. Because I have my resposes variables and the aswer is 0 and 1 in 3 bacterial genes. Somebody know how to do this in R in a easy way?   Thank so much,   José Bustos Facultad de Ciencias Universidad Catolica Ssma. Concepcion Chile --- El lun, 31/8/09, r-help-request@r-project.org <r-help-request@r-project..org>

I have been looking into both of these approaches to conducting a GLMM, and want to make sure I understand model specification in each. In particular - after looking at Bates' Rnews article and searching through the help archives, I am unclear on the specification of nested factors in lmer. Do the following statements specify the same mode within each approach? m1 = glmmPQL(RICH ~ ZONE,

I'm trying out a basic plot, but something about the way I subset my data leads to problems with the plot. Here is the first bit of my data set year,date,location,quadrat_juvenile,photo_location,photo_exists,genus,count,divers 2005,2005-04-30, 1 Fringing Reef,1, 1 Fringing Reef Coral Transect Pole 1-2 Quadrat 1,t,Acanthastrea,0,HP+MEM 2005,2005-04-30, 1 Fringing Reef,1, 1 Fringing Reef

I have data on the proportion of clutches experiencing different fates (e.g., 4 different sources of mortality) for 5 months . I need to test 1) if the overall proportion of these different fates is different over the entire study and 2) to see if there are monthly differences within (and among) fate types. Thus, I am pretty sure this is an RM analysis -( I measure the same quadrats each month).

I have a function (see below). This function has one object, ID. If I run the loops by themselves using a character value (ie,"VFFF1-7") instead of the function object, then the loops work fine. However, when I try to insert the character value via the function call, it doesn't work. I don't get an error, but the TotalCover.df dataframe does not update according to the loop

I have a function (see below). This function has one object, ID. If I run the loops by itself using a character value (ie,"VFFF1-7"), then the loops work fine. However, when I try to insert the character value via the function call, it doesn't work. I don't get an error, but the TotalCover.df dataframe does not update according to the loop criteria. Any obvious problems that