similar to: lapply and kernelUD (adehabitatHR package): Home Range kernel estimation for a list of individuals

Hi,  I am looking to perform a kernel density analysis in R, following the script that a friend of mine is using successfully. However, unlike for her, when I enter >kud<-kernelUD(data, h="LSCV"), R crashes and closes every time. This does not happen when I perform a different kernel analysis (>kud<-kernelUD(data)).. I have tried to un-install and re-install R, but it still

Hello, I am working through examples of generating Home Range using the adehabitatHR package. Everything is going fine until I load the code for the kernelUD function, it looks like this; > data(puechabonsp) > kud <- kernelUD(puechabonsp$relocs[,1], h="href") > kud I then get an error of; Error: could not find function "kernelUD" Any thoughts on what I am

Dear all, I'm trying to get a Brownian bridge kernel (kernelbb) for each combination of two consecutive animal locations (see commands below) and put them, with a loop, inside a list. It works well at the beginning but after 42 runs, it appears the following warning : >Error in seq.default(yli[1], yli[2], by = diff(xg[1:2])) :  >  invalid (to - from)/by in seq(.) I looked at the

Dear all, I have just uploaded 4 new packages on CRAN, which are on the long term designed to replace the "old" package adehabitat: * adehabitatMA: functions to perform spatial operations (morphology, buffer, etc.) * adehabitatHS: functions for the analysis of habitat selection by wildlife * adehabitatHR: functions for home range estimation of animals * adehabitatLT: functions for

Dear all, I have just uploaded 4 new packages on CRAN, which are on the long term designed to replace the "old" package adehabitat: * adehabitatMA: functions to perform spatial operations (morphology, buffer, etc.) * adehabitatHS: functions for the analysis of habitat selection by wildlife * adehabitatHR: functions for home range estimation of animals * adehabitatLT: functions for

Thank you for providing the example code... for the request of running it multiple times it would have helped if you could have confirmed that the example ran through without errors... there were a lot of mistakes in it. Look into using the reprex package to check your example next time. I don't do this kind of analysis... I really don't know what to expect from the functions. The

Dear List, It's the first time that I try to request your support. I intent do use the BRB (Biased Random Bridge) function with a boundary. I believe I have made sure to fulfill all the rules for boundaries (segment length, angle of junction) When calculating the UD I get the following Error: Error in adehabitatHR:::.fbboun(sor2, boundary, sigg, max(so$h)) : non convenient sign Can

I am having difficulty loading a points shapefile (1600 records) I generated from ArcGIS as a PointsSpatialDataFrame using the maptools readShapePoints function. library(maptools) surveypts <- readShapePoints(system.file("C:/temp/survey.shp", package="maptools")[1] Error in getinfo.shape(filen) : Error opening SHP file I've tried numerous re-exports and renames of the

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org Reproducible Rterm crash. I drew a complicated graph, then resized the graphics window. R reported > Insufficient memory for resize. Killing device

Hello, I'm new to R and am trying to use the adehabitat package for home range and habitat selection analysis for some animal radiotelemetry data sets I have, but I can't get the home range functions (mcp & kernelUD) to work with my data. I think that my problem is that I don't properly understand how the package uses data frames and factors, and I may not be defining the id

Dear all, I have a data similar to this: myframe<- data.frame (ID=c("Ernie", "Ernie","Ernie","Ernie"), Timestamp=c("24.09.2012 08:00", "24.09.2012 09:00", "24.09.2012 10:00", "25.09.2012 10:00"), Longitude=c("8.481","8.482","8.483","8.481"),

I am attempting to write what MATLAB does as : [S,A]= shaperead() which returns an N by 1 geographic array structure S containing geometric information, and a paralell N by 1 attribute structure array, A containing feature attribute information. How can this be done with R code I tried using shp<-readShapePoints() but this doesnt return the same structure that I am looking for. Please

Dear list members I am calculating kernel home-ranges and would like to visualize them with the function image(). The xy-values consist of longitude and latitude data. example:       x        y 36154.97 355143.0   R-code: xy<-hares[,c("x","y")] id<-hares[,"date"]   (ud1 <- kernelUD(xy, id)) ud2 <- getvolumeUD(ud1) image(ud2)   My problem is that the image

Dear all, I have uploaded to CRAN the version 1.7 of the package 'adehabitat'. Significant changes are listed below: * The Brownian bridge kernel estimation algorithm has been greatly improved. It now takes more than 80% less time than the previous version. A new function "liker" has also been added, which estimates the one of the two smoothing parameters of the bridge

Dear all, I have uploaded to CRAN the version 1.7 of the package 'adehabitat'. Significant changes are listed below: * The Brownian bridge kernel estimation algorithm has been greatly improved. It now takes more than 80% less time than the previous version. A new function "liker" has also been added, which estimates the one of the two smoothing parameters of the bridge

Hi all, I'm faced with a problem applying the sp package: The projection argument in readShapePoly(Shapefile,proj4string="CRS class argument") e.g.: CRS("+proj=aea +lat_1=46 +lat_2=73 +lat_0=60 +lon_0=84 +x_0=0 +y_0=0 +ellps=clrk66 +units=m +no_defs") doesn't have any impact on the plotted object. I also tested the simple example: xy = cbind(x = 2 * runif(100) -

Hi ALL, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2",

Hi, Doing more search i've discovered package RGDAL that can write a geotiff file with projection. I saved a geotiff file in UTM projection and if i read the file back in R and check the projection seems that everything is OK. But if i load the file in ArcGIS (ESRI product) i get the warning that the file is missing spatial reference so it cannot be projected, but it is displayed correctly.

Hi I'm trying to analyse some data and need to set the geographic coordinate system before I can do the analysis. I've been trying to use the project command in rgdal but keep getting an error message saying: Error in project(locationsMatrix, PROJECTION.OUT) : latitude or longitude exceeded limits ( PROJECTION.OUT <- "+proj=aea +lat_1=-18 +lat_2=-36 +lat_0=0 +lon_0=132

summary: I can display a lon-lat map on a lattice::layerplot, and I can display a Lambert conformal conic (LCC) map on a spam::image, but I can't display an LCC map on a lattice::layerplot. Example follows. What am I doing wrong? details: I've been using `lattice` (via `rasterVis`) successfully to display global atmospheric data, which works well enough (though I am definitely intrigued