similar to: plotting multiple animal tracks against Date/Time

Dear list,   Dear list, I am using read.csv.sql() from the sqldf package to read individual-based data from a csv file that is too large for R. My original file contains a column called “ID” that identifies the individual. I would like to read in data for only one individual at a time, for example "Bobby". read.csv.sql("filename",sql = 'select * from file where ID =

Dear Jim and Henrique, thank you both for your help. I have done this but run into another problem: ? In the example?below "loc1,loc2" occurs in the (now correct, thanks to your advice) "list" twice. ? ? trips=("loc1,loc2","loc2,loc3") ? DF$listoftrips=("loc1,loc2", "loc1,loc3", "loc2,loc3","loc1,loc2"). ? I?am

Dear list, ? I am using match() to match pairs of locations, e.g. trip="loc1,loc2" from a list of such pairs, e.g. list=("loc1,loc2", "loc1,loc3", "loc2,loc3","loc2,loc1"). ? In this example match() will match "trip" with the first element of "list", but not the 4th, because the order is reversed. ? How can I get a match with

Dear list,   I am looking for an R-only implementation of a Nelder-Mead process that can find local maxima of a spatially distributed variable, e.g. height, on a spatial grid, and outputs the coordinates of the new point during each evaluation. I have found two previous threads about this topic, and was wondering if something similar has been implemented since those messages were posted.   Thank

Dear members, How can I read and plot an ArcGIS raster file into R ? The file has extension .aux and contains floating point bathymetry data. The purpose is to create a spatial model in R that uses ArcGIS map data. I have managed to read and plot various shape files into my R project, but I am stuck with this now. I am new to this list and also to R, so any help would be much appreciated. Many

Dear list, ? would anybody be able to tell me why the statement ? Tripstatistics=aggregate(TripsData[2:3],by=list(Trip=Tripmatch),FUN="mean") ? seems to work well with TripsData 1 but not with TripsData 2 ? ? With TripsData 2 it yields ? Error in FUN(X[[1L]], ...) : arguments must have same length I can't see a difference in the two data sets. Could someone shed light on the error

Dear list, ? I am trying to program semi-random movement within a circle, with no particles leaving the circle. I would like them to bounce back when they come to close to the wall, but?I don't seem to be able to get this right. ?Would somebody?be able to?give me a hint ? This is my code so far, the particle starts at some point and moves towards the wall, but I don't get the

Dear list, I am trying to program semi-random movement within a circle, with no particles leaving the circle. I would like them to bounce back when they come to close to the wall, but I don't seem to be able to get this right. Would somebody kindly give me a hint ? This is my code so far, the particle starts at some point and moves towards the wall, but I don't get the "bouncing

Hi, I'm new to R and new to this forum. I'm struggling with trying to extract certain rows of data from my data.frame. The data.frame has eleven columns. Among those columns are "FISH_ID" and "DATE_TIME". FISH_ID is a factor. For each of my 21 unique FISH_IDs (levels) I have a few to a few thousand rows, each row with a unique DATE_TIME value. I would like to obtain,

Dear list, ? I am intersted in estimating movement based kernel densities for fish that were relocated at fixed receivers positioned along the coast. These data tend to display both a drift movment between receivers and a random movement component that can be estimated from the mean and the variance of the transit time between receivers. ?If I obtain an estimate of the diffusion coefficient from

Hi, I am looking to do some quantitative genetic analyses using animal models and was wondering if someone could suggest an appropriate package in R. It would help if it was similar to the ASReml genetic analyses software. Thanks, Deepa Senapathi Deepa Senapathi Centre for Agri-Environmental Research (CAER) School of Agriculture, Policy & Developement. University of Reading

Hello list; I am working with a dataset containing animal locations over time for a large number of individuals. I would like to compute a center of activity by finding the minimum harmonic mean on a grid overlaying the points of an animal (a "standard" to express center of activity in animal ecology). I have searched the archives at Jonathan Barons server but have not been able

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 R-help, I have fitted a glm logistic function to dichotomous forced choices responses varying according to time interval between two stimulus. x values are time separation in miliseconds, and the y values are proportion responses for one of the stimulus. Now I am trying to extrapolate x values for the y value (proportion) at .25, .5, and .75. I have tried several predict parameters, and they

Hi, You should look at Pinheiro and Bates (2000) Mixed-effects models in S and S-Plus. It describes how to format the correlation matrix to pass to functions lme and gls. Basically, the correlation matrix has to be one of the corStruct classes, probably corSymm for your example. So in the call to lme (or gls if you really have no random effects), use something like:

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

Hi, I am trying to simulate animal movement in a gridded landscape made up of cells. At each time step (iteration), the animal moves from one cell to another in a random fashion. This is how I am simulating movement, where a and b are the x,y co-ordinates of the animal at the previous time step: for (i in 1:no.of.steps){ direction <- sample(1:8, 1) if(direction == 1){ a <- a b <- b -

I have a data frame containing the Id, Mother, Father and Sex from about 10,000 animals in our colony. I am interested in graphing simple family trees for a given subject or small number of subjects. The basic idea is: start with data frame from entire colony and list of index animals. I need to identify all immediate relatives of these index animals and plot the pedigree for them. We're

Dear All, I have two data like this : $cat main.csv name,id,memory,storage mohan,1,100.20,1.10 ram,1,200,100 kumar,1,400,50 xxx,1,100,40 aaa,1,800,45 mount,1,200,80 > main <- read.csv(file='main.csv',sep=',' , header=TRUE) > main name id memory storage 1 mohan 1 100.2 10 2 ram 1 200.0 100 3 kumar 1 400.0 50 4 xxx 1