similar to: Using optim function for logistic model simulation

Hello! I have a dataset with the following two vectors: year<-c(1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009)

I'm trying to test a method of identifying individuals (birds) based on measured data (their calls). I have test data from known individual birds, and I am using the Mclust package to see if the program can correctly identify which calls come from different birds. So far, mclust has correctly ID'd the number of birds in the test data set (i.e., the correct # of clusters). However I

Dear Community, I have a data set with two columns, bird number and mass. Individual birds were captured 1-13 times and weighed each time. I would like to remove those individuals that were captured only once, so that I can assess mass variability per bird. I¹ve tried many approaches with no success. Can anyone recommend a way to remove individuals that were captured only once? Thanks, Ray

Dear all, I have received some data on birds that looks sth like this: # a unique id for each individual id <- c(1,1,1,2,2,2,3,3,3,4,4,5,6) # the year the bird was measured year <- c(1995, 1996, 1997, 1995, 1996, 1997, 1996, 1997, 1998, 1996, 1997, 1997, 1998) # the year the bird was hatched year.hatch <- c(1995, 1995, 1995, 1995, 1995, 1995, 1996, 1996, 1996, 1996, 1996, 1997, 1998)

Hello! I am doing an analysis on a questionnaire of hunters taken in 4 different districts of some mysterious foreign country. The aim of the study was to gather info on the factors that determine the hunting success of a peculiarly beautiful bird in that area. All variables are factors, i.e. they are variables such as "Use of Guns - yes / no", "Use of Dogs - yes / no" and

I'm fairly new to R and still learning how to use it. I could really use some help with the following problem. I have a huge .csv file containing thousands of measurements on 34 different birds. Measurements include longitude, latitude, altitude, speed, time, etc. All birds have a different number (ranging from 121 to 542). All measurements have a tripID (1 for the first trip of every bird, 2

Sorry that I was unclear. For an individual to qualify for my analysis I want both of the following two criteria to be fulfilled: First, I want to select measurement taken at a certain age: for the focal individual the year of measurement (year) should be the same as year.hatch Second, I want the focal individual to be born by a mother that reproduces for the first time. So the /parents /of

Dear all, I am experiencing problems using the glmmPQL function in the MASS package (Venables & Ripley 2002) to model binomial data with spatial autocorrelation. My question - is the presence of birds affected by various hydrological parameters? Presence/absence data were collected from 83 sites and coupled against hydrological data from the same site. The bird survey sampling effort

Dear R-experts, I have a table with following collumns: State, SamplePlot, Species and BodySize. I sampled bird species at 34 SamplePlots and 5 States (regions) monthly during two years. On each bird record I measured bodysize and identified the species. So I have many records of each species (about 150 species) at each SamplePlot and each Region (State). Now I would like bootstrap

Are there any plans to integrate Virtual Earth''s Bird''s Eye view to mapstraction? I tried an "integration hack" but then found out that there are too many issues involved (check whether birds eye is available for the current map area, disabling zoom controls,...). Has anyone else tried it? Franz

Hi, I was reading the book on "Mixed Effects Models and Extensions in Ecology with R" by Zuur et al. In Section 6.2, an example is discussed where a gamm-model is fitted, with a smoother for time, which differs for each value of ID (4 different bird species). In earlier versions of R, the following code was used BM2<-gamm(Birds~Rain+ID+

Dear List, I am running constrained correspondence analysis for abundance data of 7 birds. However, I would like to check which bird prefers which environment gradient by showing the species with different colors of the dots in cca plot (package vegan). Please kindly help and thank you Elaine [[alternative HTML version deleted]]

Hello, We are attempting to use nlme to fit a linear mixed model to explain bird abundance as a function of habitat: lme(abundance~habitat-1,data=data,method="ML",random=~1|sampleunit) The data consist of repeated counts of birds in sample units across multiple years, and we have two questions: 1) Is it necessary (and, if so, how) to specify the repeated measure (years)? As written,

Hi all, I'm currently working through the "Beginner's Guide to R" (Zurr et al.) and I'm wondering about the first exercise in chapter 3: I imported the data from BirdFluCases.txt and executed the 'names' and 'str' functions as follows: Bird = read.table(file='C:\\rbook\\BirdFluCases2.txt', header=TRUE) names(Bird) [1] "Year"

Hi all, I'm currently working through the "Beginner's Guide to R" (Zurr et al.) and I'm wondering about the first exercise in chapter 3: I imported the data from BirdFluCases.txt and executed the 'names' and 'str' functions as follows: Bird = read.table(file='C:\\rbook\\BirdFluCases2.txt', header=TRUE) names(Bird) [1] "Year"

Hello r-help, I'm trying to fit birds counts over years using glm. In fact I'm trying to reproduce an analysis already perform with genstat (attach document). I have done (with Estate and year as factors): Model1 <- glm(Females~Estate+Year+offset = log(area)), family = quasipoisson(link = log), na.action = "na.exclude") After I have calculated the prediction using: Pred1

Hello, I’m trying to fit the following model: Dependent variable: MAXDEPTH (the maximum depth reached by a penguin during a given dive) Fixed effects: SUCCESSMN (an index of the “individual quality” of a bird), STUDYDAY (the day of the study, from -5 to 20, with 0=Dec 20), and the interaction SUCCESSMN*STUDYDAY Random effect: BIRD (the bird id, as each bird is performing several dives)

Hi, I am trying to write a simulation of the movements of four animals between six patches. The movement between patches is based on a first-order Markov chain so that the next patch they visit depends on the patch they were in before. I have written code that allows me to simulate the movement of one animal but when I add more there seems to be a problem and all chains come back as N/A. I can

I am trying to make it so two columns with similar data use the same internal numbers for same factors, here is the example: > read.csv("test.csv",header =FALSE,sep=",") V1 V2 V3 1 sun moon stars 2 stars moon sun 3 cat dog catdog 4 dog moon sun 5 bird plane superman 6 1000 dog 2000 > data <-

I am using package ks() to build 3D representations of bird territories and calculate territory volume from spatial data (simply x, y, and z coordinates). What I want to do is determine at what sample size (# locations collected) does the territory volume stop increasing. This should give me an idea of the number of points needed for future seasons. So I have a couple of birds each with 200