search for: hatch

Dear all, I was wondering if anyone could help solve a problem of a missing interaction effect!! I carried out a 2 x 2 factorial experiment to see if eggs from 2 different locations (Origin = 1 or 2) had different hatching success under 2 different incubation schedules (Treat = 1 or 2). Six eggs were taken from 10 females (random = Female) at each location and split between the treatments, giving 30 eggs from each location in each treatment. Overall proportions hatching were as follows: Treat...

Dear all, I'm struggling with a plot and would value any help! I'm attempting to highlight a histogram and density plot to show a proportion of cases above a threshold value. I wanted to cross-hatch the area below the density curve. The breaks and bandwidth are deliberate integer values because of the type of data I'm looking at. I've managed to do this, but I don't think it is very good! It would be difficult, for example, to do a cross-hatch using this technique. allele.plot &l...

Hello, I'm having an error when trying to fit the next GLM: >>model<-glm(response ~ CLONE_M + CLONE_F + HATCHING +(CLONE_M*CLONE_F) + (CLONE_M*HATCHING) + (CLONE_F*HATCHING) + (CLONE_M*CLONE_F*HATCHING), family=quasipoisson) >> anova(model, test="Chi") >Error in if (dispersion == 1) Inf else object$df.residual : missing value where TRUE/FALSE needed If I fit the same model by using...

Hello, I would like to build rectangles in a plot and use color and different type of hatching for filling rectangles. I don't find the way to draw hatchings. I'm thinking to build segment by segment inside each rectangle but I'm sure that exists a better way to do that. I didn't find any documentation about that. > symbols(1,1,rectangles=cbind(1,1),bg="red",...

...inpengS1: records of temperature ? wchillpengS1: records of temperature ? minpengS2: records of temperature ? wchillpengS2: records of temperature ? nnd: nearest neighbour distance ? npd: nearest penguin distance ? eggs: numbers of eggs ? lay: laying date (julian calendar) ? hatch: hatching date (julian calendar) I have some NAs in the data. I want to test the model with all the variable then i want to remove some, but the ideal model: GLM.1 <-lmer(fledgesucc ~ +lay +hatch +elev +seadist +nnd +npd +meanterranova +minpengS1 +minpengS2 +wchillpengS1 +wchillpengS2 +(1|te...

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 the focal individual should have year == year.1st.reprod. The problem for me, I think, is that in my data set with several generations, an individual can be both offspring and pare...

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) # the year when the bird reproduced the first time (for simplicity in the dummy data, all birds reproduce the year after it was born) year.1st.reprod <- year.hatch + 1 # some status (0/1) of an i...

....)) doesn't accept lwd...) Here is the updates code I wrote: polygon2 <- function (x, y = NULL, density = NULL, angle = 45, border = NULL, col = NA, lty = par("lty"), lwd =par("lwd") ,..., fillOddEven = FALSE) { ..debug.hatch <- FALSE xy <- xy.coords(x, y) if (is.numeric(density) && all(is.na(density) | density < 0)) density <- NULL if (!is.null(angle) && !is.null(density)) { polygon.onehatch <- function(x, y, x0, y0, xd, yd, ..debug.hatch = FALSE,...

....)) doesn't accept lwd...) Here is the updates code I wrote: polygon2 <- function (x, y = NULL, density = NULL, angle = 45, border = NULL, col = NA, lty = par("lty"), lwd =par("lwd") ,..., fillOddEven = FALSE) { ..debug.hatch <- FALSE xy <- xy.coords(x, y) if (is.numeric(density) && all(is.na(density) | density < 0)) density <- NULL if (!is.null(angle) && !is.null(density)) { polygon.onehatch <- function(x, y, x0, y0, xd, yd, ..debug.hatch = FALSE,...

#I am putting a test together for an introductory biology class and I would like to put different cross hatching inside of each bar for the bar plot below color <- c("Brightly Colored", "Dull", "Neither") lizards <- c(277, 70, 3) liz.col <- data.frame(color, lizards) qplot(color, lizards, data=liz.col, geom="bar", ylab="Observed Matings", main=&qu...

I have used polygon() to mark the confidence limits of a survival curve. In another project, I have used the col parameter to fill my boxplots. The poly() description refers to filling but actually produces hatching (i.e. lines ). boxplot() does truly fill the boxes with a colour or shades of grey (e.g. col="red"). My novices perception of R graphics is: If you can hatch it, you cannot fill it. If you can fill it you cannot hatch it. Is this correct or is there a way of hatching my boxplots an...

...about year to year data. This is the relevant part of my table: datum month abundance 26/03/11 March 1 27/03/11 March 0 28/03/11 March 1 29/03/11 March 0 30/03/11 March 0 31/03/11 March 0 1/04/11 April 0 2/04/11 April 0 3/04/11 April 0 ... 9/09/11 September 1 I tried making a sequence: hatch$axis.ticks2=seq(from=hatch$datum[1],to=hatch$datum[171], by=30) But this sequence is too short to use the names.arg function. Next, I tried more ideas from this forum that might work. But the best I got was by just replacing the day by a month and reducing the size of the tickmarks. Or by blanking...

Dear Listers, I am drawing a plot of two density curves, for male and female incomes. I would like to shade/hatch/color (whatever) the areas under the curves which are distinctive for each gender. This is the code I have tried so far: m <- density(topmal.d$y, bw = "sj") f <- density(topfem.d$y, bw = "sj") par(mfrow = c(1,1)) plot(x = c(0,400), y = c(0,0.02), type = "n",...

...ing a lot of pointing and clicking in Excel. I have a dataframe of bird nest check observations, in which I know the date of the first check, the date of the second check (both currently in Julian date format), the status of the nest at the second check (alive or failed), and the date that the nest hatched (i.e. changed from Incubation stages to Brood-rearing stage). Many nests have more than one record, as there were several nest checks throughout the duration of the nesting attempt. What I want to do is assign a nest Stage variable, either Incubation or Brood-rearing. It's very easy to do so...

M. Camanm posted in Jul 1999 the following message: " Is there any way to fill the bars in a barplot() with solid lines for postscript output, i.e. cross hatch or parallel lines, or a halftone gray rather than (semi) continuous-tone gray produced by gray()? S allows this, or at least used to, via the angle and density arguments to barplot(). The objective of course, is to produce camera ready black and white plots on a laser printer. The gray values produ...

....1 on Windows Vista. I have small dataset which consists of: # NestID: nest indicator for each chicken. Siblings sharing the same nest have the same nest indicator. # Chick: chick indicator consisting of a unique ID for each single chick. # Year: 1, 2. # ClutchSize: 1-, 2- , 3-eggs. # HO: hatching order within each clutch (1, 2, 3 [first, second and third-hatched chick]). # SibComp: sibling competence: present/ absent (0, 1) # Death2: death at two days post-hatch (0, 1) # Death10: death at ten days post-hatch (0, 1) So a subset of my dataset looks something like this: NestID Chick Y...

I've been trying to run AOPA's Real Time Flight Planner (RTFP) from Jeppesen Sanderson Inc. It's almost works, but I get the following error when it tries to draw the route: fixme:gdi:ExtCreatePen PS_USERSTYLE not handled fixme:gdi:ExtCreatePen Hatches not implemented X Error of failed request: BadValue (integer parameter out of range for operation) Major opcode of failed request: 58 (X_SetDashes) Value in failed request: 0x0 Serial number of failed request: 63199 Current serial number in output stream: 70195 I don't actually...

...an add the legend later using imaging software such as paint but i lose resolution by doing this. Is it possible to place the legend out side of figure axes on the white area around all three panels? Any help would be greatly appreciated. Thank you Grant single<-read.table("C:\\Grant\\hatchfirstonly.txt",header=T) attach(single) names(single) par(mfrow=c(1,3)) maleage<- factor(maleage, labels = c("Old","Young")) colors <- c('gray', "navyblue") shapes <- c(17,19) proportion<-(eggs1-unhatch1)/eggs1 plot(proportion~maleage, xlab= &qu...

I'm trying to convert some S-Plus code which generates barplots and other shaded area plots to R. If I specify that I want hatching using the angle and density arguments, the messages 1: argument `density' is not used (yet) in: .NotYetUsed("density", error = FALSE) 2: argument `angle' is not used (yet) in: .NotYetUsed("angle", error = FALSE) so apparently these features have yet to be implemen...

Hi R experts, I have been using ReML as follows... model<-lmer(late.growth~mtf+year+treat+hatch.day+hatch.day:year+hatch.day:treat+ mtf:treat+ treat:year+ year:treat:mtf+(1|fybrood), data = A) then I wanted to plot the results of the three way interaction using lines() as follows... tmp<-as.vector(fixef(model)) graph1<-plot(mtf,fitted(f2), xlab=list("Brood Size"), ylab=list(&...