search for: chick

Dear R-users, I did some more research and I'm still not sure how to set up an ANCOVA with nestedness. Specifically I'm not sure how to express chicks nested within boxes. I will be getting Pinheiro & Bates (Mixed Effects Models in S and S-Plus) but it will not arrive for another two weeks from our interlibrary loan. The goal is to determine if there are urbanization (purban) effects on chick health (rtot) and if there are differences betw...

Hi, In the examples from the ReShape package there is a simple example of using melt followed by cast that produces a smallish amount of output about the chicks database. Here's the code: library(reshape) names(ChickWeight) <- tolower(names(ChickWeight)) chick_m <- melt(ChickWeight, id=2:4, na.rm=TRUE) DietResults <- cast(chick_m, diet + chick ~ time) DietResults My challenge is to extract an plot only a portion of this data. I woul...

Hi R experts, I am interested on the effects of two dietry compunds on the growth of chicks. Rather than extracting linear growth functions for each chick and using these in an analysis I thought using ReML might provide a neater and better way of doing this. (I have read the pdf vignette("MlmSoftRev") and "Fitting linear mixed models in R" by Douglas Bates but I am n...

according to the documentation of the cast function in the reshape function, I would expect this bit of code from the examples to calculate marginal means over only the 'diet' variable. #Chick weight example names(ChickWeight) <- tolower(names(ChickWeight)) chick_m <- melt(ChickWeight, id=2:4, na.rm=TRUE) cast(chick_m, diet + chick ~ time, mean, margins="diet") But, summary(cast(subset(chick_m, time == 0), diet + chick ~ time, mean, margins="diet")) returns:...

Dear R-experts, Quite new to R on this end, but learning fast (I hope). I am running version 2.7.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: sibli...

...down a subset of the variables to be described to describe(), but not a subset of the weights variable. describe() then searches the whatever data set is attached in order to find the weights variables, but this is in its original (i.e. not subsetted) form. Here is an example using the ChickWeight dataset that comes in the "datasets" package. data(ChickWeight) attach(ChickWeight) library(Hmisc) #this gives descriptive data on the variables "Time" and "Chick" by levels of "Diet") by(cbind(Time, Chick), Diet, describe) #trying to add weights, ho...

Dear R-users, I set up an experiment where I put up bluebird boxes across an urbanization gradient. I monitored these boxes and at some point I pulled a feather from a chick and a friend used spectral properties (rtot, a continuous var) to index chick health. There is an effect of sex that I would like to include but how would I set up a regression and look at the effect of urbanization (purban, a continuous var)) on feather properties of chicks within boxes. So th...

...y="gaussian", na.action=na.omit) where rtot = total reflectance, box = nest box (i.e., birdhouse), julian = day of the year and purbank = the proportion of urban cover in a 1 km buffer around the nest box. I'm not interested in the box effect and I've seperated males and female chicks. I've asked about nestedness before and I was given code that included "|" to indicate nestedness but this indicates a grouping does it not? I suspect that there is something wrong. In the summary I get Coefficients: Estimate Std. Error t value Pr(>|t|)...

Dear R-experts, I am running R version 2.7.1 on Windows Vista. I have a small dataset which consists of ?chick ID?, ?year (0, 1)?, ?hatching order [HO, defined as first, second and third-hatched chick]?, and the binary outcome of interest ?death (0, 1)?. So a subset of my dataset looks like this on a txt file: y ID Yr HO 1 1 1 First 0 2 1 First 0 3 1 Second 0 4 1 First 1 5 1 First 0 6 1 Third 0 7 1 First...

Estimado Oliver Nuñez Envío un ejemplo reproducible. Javier Marcuzzi # de donde tomo datos, y tiene el modelo (en el pdf) library(MCMCglmm) # librería con las funciónes que voy a usar library(nlme) datos0<-ChickWeight # creo algunos datos que agrego a los origonales Factor<-as.numeric(datos0$Chick) Factor[Factor > 0 & Factor <= 10] <- 'A' Factor[Factor > 10 & Factor <= 20] <- 'B' Factor[Factor > 20 & Factor <= 30] <- 'A' Factor[Factor >...

...the data being calculated). Currently, my approach is to create an environment, then assign things to the environement, then as.list the environment to get a list. Running the code gives, for example: > e2 <- new.env(FALSE, NULL) > assign("dude", 123, env=e2) > assign("chick", 456, env=e2) > as.list(e2) Error in as.vector(x, "list") : cannot coerce to vector Is this the best way to make a list like this? Thanks, Ben

...n categories only i want to search which are flagged as 0. please help me with the R script. *Menu Name* 9\ bobbie" 9\ chz steak" 9\ tuna" provolone 20\ bobbie" bottled soda 20oz cran-slam ww american small chips medium drink 9\ meatball" capriotti's water 20'' chicken cheese steak 9\ veg turkey" medium chips 9\ capastrami" 12\ bobbie" 12'' chicken cheese steak cookie 12\ chz steak" 9\ cole turkey" kid grilled cheese white 12\ italian" 12\ meatball" 12\ capastrami" turkey sand w 20\ slaw be jo" swiss 12\ cole...

...shootdowns. The answer to that question is `no''. There appears to be very little, if any additional cost to the guest of using persistent grants. There is perhaps a small saving, from the reduced number of hypercalls performed in granting, and ending foreign access. Signed-off-by: Oliver Chick <oliver.chick@citrix.com> --- drivers/block/xen-blkback/blkback.c | 149 +++++++++++++++++++++++++--- drivers/block/xen-blkback/common.h | 16 +++ drivers/block/xen-blkback/xenbus.c | 21 +++- drivers/block/xen-blkfront.c | 186 ++++++++++++++++++++++++++++++----- include/xen/...

Dear R crew: I think I am in the right mailing list. I have a very simple dataset consisting of two variables: cestode intensity and chick size (defined as CAPI). Intensity is clearly overdispersed, with way too many zeroes. I'm interested in looking at the association between these two variables, i.e. how well does chick size predict tape intensity? I fit a zero inflated negat. binomial model using the "pscl" package....

...seen below in the code. > I have also tried to use predict, with exactly the same result: the > offset is ignored. This applies to both lms and glms. > > Am I missing something here? > Thank you > Samuel Riou > > CODE > #no offset lmA<-lm(MassChange24h~DATEN1, subset(Chicks1, Year==2007? & > AGE>10), na.action=na.exclude) > summary(lmA) > > #linear offset > lmAO<-lm(MassChange24h~DATEN1+offset(-0.37356*AGE), subset(Chicks1, > Year==2007? & AGE>10), na.action=na.exclude) > summary(lmAO) > > > print(Chicks$DATEN1[Year==2...

An excerpt from dataset ChickWeight: ??? weight Time Chick Diet 1?????? 42??? 0???? 1??? 1 2?????? 51??? 2???? 1??? 1 3?????? 59??? 4???? 1??? 1 I am interested in the residuals of the dataset.? Specifically in saving them to another format. I have been creating text files with sink. CW.lm <- lm(weight ~ Diet, ChickWeight)...

Dear R crew: I am sorry this question has been posted before, but I can't seem to solve this problem yet. I have a simple dataset consisting of two variables: cestode intensity and chick size (defined as CAPI). Intensity is a count and clearly overdispersed, with way too many zeroes. I'm interested in looking at the association between these two variables, i.e. how well does chick size predict tape intensity? Since I have a small sample size, I fit a zero inflated negat. Bi...

...SMTP; 16 Jan 1997 19:48:39 -0000 Received: from dhp.com (dhp.com [199.245.105.1]) by dhp.com (8.8.4/8.6.12) with SMTP id OAA18650; Thu, 16 Jan 1997 14:48:14 -0500 Date: Thu, 16 Jan 1997 14:48:13 -0500 (EST) From: Matt <panzer@dhp.com> Approved: R.E.Wolff@BitWizard.nl To: Christopher Hicks <chicks@chicks.net> cc: linux-security@redhat.com Subject: Re: [linux-security] Re: logwatching In-Reply-To: <Pine.LNX.3.94.970116124053.13190B-100000@yakko.chicks.net> Message-ID: <Pine.LNX.3.95.970116144453.18239A-100000@dhp.com> MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCI...

Dear R-help subscribers Would you tell me how to apply SSfpl with binary data as below? Unfortunately, there is not the EXAMPLE in help(SSfpl) for binary data but for quantitative data(Chick). V1: dose V2: log-transformed dose V3: response (rate) V1 V2 V3 1 0.775 -0.2548922 0.1666667 2 5.000 1.6094379 0.8148148 3 10.000 2.3025851 0.5000000 4 20.000 2.9957323 0.6428571 I also attach the error message in R 1.5.1 in following. > targan V1 V...

Estimados Tengo un error que me desconcierta, es un código que simplifiqué de otro trabajo donde no hay problemas, sin embargo me da un error. Una diferencia es que en mi otro trabajo uso spline y ahora polinomio, este es de segundo grado y se encuentra tanto en efectos fijos como aleatorios, el modelo es correcto, corre con MCMCglmm pero no con nlme. grid <- expand.grid(Tiempo=4:6,