Search Archives

> But why not just > c(x[1], diff(x))/max(x) So simple! Thank you ever so much Berend. Best wishes, Ben -- View this message in context: http://r.789695.n4.nabble.com/error-in-self-made-function-cannot-deal...

The below works for one trap. Surely there is a quick way to scale it to all of the traps? > temp.ACAP1<-outer(release.days[Trap=="ACAP1"],collection.days.2[Trap=="ACAP1"],'-') > temp.ACAP1<-ifelse(temp.ACAP1>=0,NA,temp...

Hi all, I?m trying to plot confidence intervals for the fitted values I get with my lme model in R. Is there any way I can plot this in the form of a shaded band, like the output of...

Hi all, My code: > temp<-outer(release.days,collection.days,'-') > temp<-ifelse(temp>=0,NA,temp) > release.diff<-apply(temp,2,max,na.rm=TRUE) works for one trap and does what I want. That is, it determines the time...

Dear All, I have a dataframe made up of individual beetles consisting of individual number, family number, mother's family number, father's family number, and sex of the beetle. I would like to pair up the individuals for breeding...

Hi Jon, Sorry about the lack of information. I'll give the data from one of my traps: (nb. all of the data for each trap is referenced to the start day of the trap. ie ACAP1 started collecting at...

I have a function to calculate the rate of increase (the difference between the value and the previous value divided by the total number of eggs in a year) of egg production over the course of a year: rate <- function...

I am attempting to emulate the VLOOKUP function from Excel in R. I want to compare one column (coll.minus.release) with another (release.days) to get the number of parasitoid released at that time (TotalParasitoids). for example: coll.minus...

Hi all, I have scoured the archives of this forum but nothing quite seems to fit the bill... I would like to plot a graph displaying two variables (y axes) that share date as the x axis. However, the date...