similar to: pairing data using combn with criteria

All, I have: sales <- c(2300,900,155,102,42,10) names(sales) <- c("Christmas","Valentine's Day", "Mother's Day","Father's Day", "Thanksgiving","New Year's Day") barplot(sales,ylim=c(0,2500)) But it doesn't place all of the name labels on the plot. So I tried: sales <-

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to

[cc'ed back to r-help] I've started to take a look, and there's nothing immediately obvious about the problem with the fit (the warnings and errors are about a "non-positive-definite Hessian", which usually means an overfitted/poorly identified model) -- still working on whether there's a way to get more useful information. As it turns out, glmmADMB's default

Hi everyone! I have some doubts about mixed effect models and I hope someone could help me. I´m trying to analyze a dataset coming from samples of dung beetles in the same forest fragments along 3 consecutive years (1994, 1995 and 1996) and 14 years after (2010). I sampled dung beetles in 18 different fragments with different sizes and different degrees of isolation. My aim is to determine

Dear Colleagues, I hope to know how ordinal logistic regression with a mixed model is made in R. We (My colleague and I) are studying the behavior of a beetle. The attraction of beetles to a stimulus are recorded: the response is Slow, Mid, or Fast. They are based on the time after the presentation of the stimulus to the beetles. Because we do not observe the behavior continuously but do

Hi Everyone, I would like to add errorbars to a lattice plot. I already have a function that adds error bars to a plot (originally written by George Gilchrist) which takes as arguments the coordinates of the point and then the size of the error bar. I would like to integrate it in my lattice code but I seem to be stumped. I can put the error bar function inside my code as if it were a panel

Bradley Taylor shot me this *goldmine* of ugly ugly ugly dogs: http://www.sonoma-marinfair.org/uglydogvote.shtml I *have* to use one of these for the Mongrel BUGS Mascot. He''ll go on our bug list page and replace the little beetles on the left. Pick the dog you think best represents a lovable but defective pooch and reply to this with your +1. Let the voting begin! -- Zed A. Shaw

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 geom_smooth() in ggplot2 package. ggplot2 seems to use only lm, glm, gam, loess and rlm as smoothing methods. Any advice on the functions I should use to accomplish this will be very helpful. Thank you very much.

Hello, I am using function lmodel2 to calculate RMA regression between life-history traits in ladybirds beetles. It works well but I am not able to plot an RMA regression line on the scatterplot of my data. I am of course unable to plot the confidence intervals. For ordinary least square regression abline(lm(x~y))works well but for RMA regression, abline(lmodel2(x~y))does not do it. If somebody

Dear all I was wondering whether such a long post could be fortune-ed. What do you think? Regards Liviu On Sat, Aug 21, 2010 at 9:33 PM, Sharpie <chuck at sharpsteen.net> wrote: > Well, I can think of three ways it can go down: > > > 1. ?You want a shiny new pony. > > You ask about it on the mailing list and it seems that everyone else in the > world responds

Hi, I have a couple of samples that produce interesting artifacts when encoded with the CVS snapshot of 2001-02-13. Both are about a meg. ftp://slumber.dhs.org/tmp/4.wav.bz2 When encoded with oggenc -b 128, there is a sort of stereo separation in the drums, while in the original they are 'solidly' positioned in the stereo field. This also occurs with -b 160 and is barely audible with

I have some data I would like to model which involves choice of food by dung beetles. There are a number of experiments where in each case, there are five choices. Overall there are more than 5 different foods being compared (including a placebo) and different experiments use different comparisons. The problem is a generalization of Bradley-Terry but it differs from some generalizations in

Dear All, I'm new with mapping and using shapefiles in R. I use the version: R 3.2.1 GUI 1.66 with Snow Leopard Mac OS (6956) I want to build a map of the countries of Africa, in order to add points of collecting species of beetles, for later publishing. library(maptools) library(mapdata) When I put the extreme coordinates of Africa, the map began only at longitude 0? and thus lacks the

Hello! I am building a package, which includes also one Fortran subroutine, which works fine if I compile it as a shared library and load it into R via dyn.load(). However, when I launch R CMD check it doesn't stop with checking examples. It's just doing and doing ... I pasted the whole output from R CMD check. Does anyone have any suggestions? I'm still using R 2.0.1.

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(x){ storage <- matrix(nrow=length(x),ncol=1) storage[1,] <- x[1] / max(x) # as there is no previous value for( i in 2:length(x)){ p <- i - 1

Hello everyone Not strictly an R question but close... hopefully someone will be able to help. I wish to compare the switchpoints in two switchpoint regressions. The switchpoints were estimated using the segmented library running in R, and I have standard errors for the estimates. I initially thought I could just bootstrap confidence intervals for the difference between the switchpoints,

Hello all I'm trying to do a simulation that involves identifying the minimum point between two peaks of a (usually) bimodal distribution. I can do this easily if there are only two peaks: CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys to Xs for(p in 1:512) ifelse(CnBdens$y[p]>CnBdens$y[p-1],peak1<-p,break) #identifies first peak in probability

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 values for each variable are not the same - for example, some parasitoids were not released on days that collections from the trap took place, whilst sometimes releases did occur on the same

Hi, i am working on a problem where i need to compute the products of all possible combinations of size m of the elements of a vector. I know that this can be achieved using the function combn(), e.g.: > vector <- 1:6 > combn(x = vector, m = 3, FUN = function(y) prod(y)) In my case the vector has 2000 elements and i need to compute the values specified above for m = 32. Using combn() i

-------- Original Message -------- Subject: [AGDG-LIST:428] Summer Course in Guelph Date: Fri, 22 Dec 2006 09:12:24 -0500 From: Larry Schaeffer <lrs at uoguelph.ca> Reply-To: lrs at uoguelph.ca To: Animal Geneticist's Discussion <agdg-list at colostate.edu> The Centre for Genetic Improvement of Livestock at the University of Guelph is pleased to announce a one week summer course