similar to: Export kde object as shapefile

Dear all, Just wondering if someone could help me out converting my code from a for() loop into a foreach() loop or using one of the apply() function. I have a very large dataset and so I'm hoping to make use of a parallel backend to speed up the processing time. I'm having trouble getting selecting three variables in the dataset to use in the foreach() loops. My for() loop code is:

Hi R-Help, I am using R to do functional outlier detection (using PCA to reduce to 2 dimensions - the functional boxplot methodology used in the Rainbow package), and using Hscv.diag function to calculate the bandwidth matrix where this line of code is run: result <- optim(diag(Hstart), scv.mat.temp, method = "Nelder-Mead", control = list(trace = as.numeric(verbose))) Within the

Hi! I have a problem with Kernel density estimate plot for 3-dimensional data in ks-package. Here the example: # load ks, spatstat # three-dimensional kernel density of B B <- pp3(runif(300), runif(300), runif(300), box3(c(0,1))) x <- unclass(B$data)$df H <- Hpi(x) fhat <- kde(x, H=H) plot(fhat) plot(fhat, axes=FALSE, box=FALSE, drawpoints=TRUE);

It is known that function bkde2D (package "KernSmooth") returns a matrix of density estimates over the mesh induced by x1 and x2. In Details it is written that "... heights of the kernel, scaled by the bandwidths, at each datapoint are summed. This sum, after a normalization, is the corresponding fhat value in the output". There are several questions: 1) How to calculate

I'm trying to make a set of contour plots of bivariate kernel density estimates, showing three such plots overlaid, similar to this plot http://euclid.psych.yorku.ca/SCS/Private/Test/ridge-boot2.pdf except that I would like to have the contours *filled* (using transparent colors). To make this reproducible, I've saved the results of KernSmooth::bkde2D() in the following file:

Como usas la función image puedes consultar la ayuda ?image o help(image) y encontrarás el siguiente ejemplo donde se usa un diferente color Palette (mencionada por pepeceb en su respuesta). x <- 10*(1:nrow(volcano)) y <- 10*(1:ncol(volcano)) image(x, y, volcano, col = terrain.colors(100), axes = FALSE) # O puedes usar directamente el número para indicar el color image(x, y, volcano, col =

Hello there, I am presently using the locfit package in "R" and would appreciate some help here. Could anyone tell me tell me how I can obtain the various components (x and y values)of the density estimation after using the "R" command "locfit"? For example, with the command "fhat<-density(somename)" I can obtain the x and y values by simply calling

Hi, I've made several attempts to boot syslinux from the harddisk. I have a win 98 startdiskette with syslinux and copybs. I strictly follow the procedure in the syslinux documentation. Yet, the result of the boot is not that linux is starting. Instead, the Microsoft windows screen flashes for a couple of seconds and is then replaced by a Dos prompt. Syslinux is V2.04. The PC is an IBM Aptiva

Hello I wonder if anybody can help, I am using the package adehabitatHR to estimate the potential distribution of a species using the command "domain" In the example given in the AdehabitatHS manual a map containing elevation information is loaded (class= spatial pixels data frame) as well as the GPS points of the animals being tracked, these are then plotted on each other and

Hi, ? I am trying to write a shapefile with projection. I have my data in a data.frame called try and consists in xy coordinates and a numerical attribute value z1. ? Libraries loaded are: sp, rgdal, raster, maptools ? head(try) ???????? x?????? ??? y??? ? ? ? ? ? ? ??? z1 1 610237.1 ???????? 3375751 ???????????? 8.221 2 610236.1 ???????? 3375750 ???????????? 8.153 3 610236.1 ????????

I have several large shapefiles that I need to use to draw a map. I want to simplify the files and then save the output so that I can use them again. The methods that I have tried do not allow me to save the simplified output as a shapefile. I have tried: thinnedSpatialPoly in the package maptools, but I understand that you cannot save the output as a shapefile because the topographical features

[This email is either empty or too large to be displayed at this time]

I'm new to R and I'm working on point shapefiles. Is there a way that you could interpolate a shapefile via kriging in R using an attribute? All examples on the internet are using txt files and CSVs. Thanks a lot. [[alternative HTML version deleted]]

Hello I'm working with NetCDF files in an animation and am trying to superimpose a shapefile on the image as it is generates the html pages. If I take out the lines for the shapefile, it works correctly, however I'm having difficultly including the shapefile. If anyone has any ideas I would greatly appreciate the assistance. R Version 2.15.0 (2012-03-30) Platform i386-pc-mingw32/i386

My sincerest apologies, as this is a very elementary problem, but I have searched through archives and FAQs and on the web, and I am at the end of my own resources. I need to do analysis of a spatial point process occuring with four counties, using spatstat. I've had no trouble importing the shapefiles, but I can't seem to figure out what I need to do in order to use the boundary as an

Ok Thanks Richard, I'll try it ! Alain VONDRA Chargé d'exploitation des Systèmes d'Information Direction Administrative et Financière +33 1 44 39 77 76 UNICEF France 3 rue Duguay Trouin 75006 PARIS www.unicef.fr -----Message d'origine----- De : Richard W.M. Jones [mailto:rjones@redhat.com] Envoyé : mercredi 8 octobre 2014 17:23 À : VONDRA Alain Cc : libguestfs@redhat.com

Hi, I have done the interpolation for my data and I was able to create the contours in multipanel with the help of Pascal. Now, I want to clip the contour with the shapefile. I want only the portion of contour to be displayed which falls inside the boundary of the shapefile. The data mydata.csv can be found on https://www.dropbox.com/s/khi7nv0160hi68p/mydata.csv The data for shapefile can be

Hi all, who knows how to get an ESRI Shapefile for the NUTS-2 Regions of the enlarged European Union? Particularly I want to draw maps of Germany, Poland, Czech Republik, Hungary and Austria. I've found Shapefiles for the US, Russia and other countries elsewhere in the web, but for Europe it seems really difficult. With many thanks for any hint Albrecht

I would like to create a voronoi map that can show me point values like the ones in colum street_cle. I was able to do it in tileplot(latticeExtra) like the ones attached, and with the following commands: xy=list(x=x2$X, y=x2$Y) xy$street_cle <- x2$street_cle tileplot(street_cle~ x*y, xy) But this is a trellis. Is it possible to be converted to a shapefile or image so I can plot it using

Dear all, I am currently working with the spatstat package, using windows and pixel images. First: My aim is to transform a shapefile (see attached) into a pixel image. My idea is to start transforming the shapefile into a Spatial Polygon file: x <- readShapeSpatial("200001441.shp") y <- as(x, "SpatialPolygons") z <- as.owin(y) Given z, I want to identify each