search for: hatched

...hing was a success and not-hatching was a failure, and analysed as a binomial GLMM. I'm particularly interested in the interaction between the two factors. An expression reproducing the raw data is attached at the end of the post in case it is helpful. hatch.frame$success<-cbind(hatch.frame$Hatched,hatch.frame$Nothatched) model<-lmer(success~Origin*Treat+(1|Female),family=binomial,method="ML",data=hatch.frame) model2<-update(model,~.-Origin:Treat) anova(model,model2) Data: Models: model2: success ~ Origin + Treat + (1 | Female) model: success ~ Origin * Treat + (1 | Female)...

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

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

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", ...

Dear list, i'm a new R user, so I apologize if the topic is already being addressed by some other user. I'm trying to determine if the reproductive success of a species of bird is related to a list of covariates. These are the covariates: ? elev: elevation of nest (meters) ? seadist: distance from the sea (meters) ? meanterranova: records of temperature ? minpengS1: records

...;> >># 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...

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 ind...

Hello dear members of R-help and R-core mailing list, I am not sure if this request is a "ticket" that should be filled somewhere outside the mailing list. If so, I apologize for not doing and would like to know where I should have filled it. And to the subject matter: I would like to use a command like this: plot(c(1,8), 1:2, type="n") polygon(1:7, c(2,1,2,NA,2,1,2),

Hello dear members of R-help and R-core mailing list, I am not sure if this request is a "ticket" that should be filled somewhere outside the mailing list. If so, I apologize for not doing and would like to know where I should have filled it. And to the subject matter: I would like to use a command like this: plot(c(1,8), 1:2, type="n") polygon(1:7, c(2,1,2,NA,2,1,2),

#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",

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

Hey, I am trying to create barplot of abundances over time (in days). The period is over 171 days, so I don't want to have all labels there but only the first day of the month. I couldn't find anything like this on the forum yet. Mostrly it's 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

..."bulge" out on the left and right of the overlapped area. What I'm wanting to do is shade/hatch/color etc these areas which are unique to each density, that is, which are outside the overlapped area. My code is successful at hatching the bulging areas, but leaves me with a double hatched area for the overlap, which is distracting. If I could turn this double hatched area white, that would achieve my goal, though ideally I would like to be able to specify something like "shade only areas under the m density curve which are not also under the f density curve (and shade only u...

...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 w...

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

...tID: 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 Year ClutchSize HO Hatching SibComp Death2 Death10 1 1 1...

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:

...le) 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= "Age of Male ", ylab = "Proportion of Eggs Hatched") plot(proportion~femweight, pch=shapes[maleage], col=colors[maleage],bg="gray", xlab= "Female Weight (grams) ", ylab = " Proportion of Eggs Hatched") abline(lm(proportion~femweight)) plot(proportion~copdurationsecs, pch=shapes[maleage], col=colors[maleage],bg=&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

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"),