search for: rtot

...y 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 between sexes (sex) and the effect of being in the same clutch (box). The model is rtot = sex + purban + (chick)box. I've loaded the package lme4. And the code I have so far is bb <- read.csv("C:\\eabl\\eabl_feather04.csv", header=TRUE) bb$sex <...

...d we quantified certain aspects of the color of those feathers. Since I often have more than one sample from a nest, I thought I should use a nested design. Here's the code I've been using and I'd appreciate if someone could look it over and see if it was correct. bb.glm1 <- glm(rtot ~ box/(julian +purbank), data=bbmale, family="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 eff...

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 the model should look something like rtot = se...

...ad.table(file=fn,header=FALSE,skip=2) #Substitute the header: names(DF) <- c("year","month","day","hour","FN","Punc","Pcor","Pmelt","ETP","ETR","RS","RI","RB","Rtot","Ssnow","SI","SSM","SUZ","SLZ") #Create datetime vector library(chron) DF$datetime <- with(DF,chron(paste(month,day,year,sep="/"))+hour/24) #Try to convert DF into ts - object DFts DFts <- as.ts(DF) # works, but gives back...

...wing example from the help menu ## Compute row and column sums for a matrix: x <- cbind(x1 = 3, x2 = c(4:1, 2:5)) dimnames(x)[[1]] <- letters[1:8] apply(x, 2, mean, trim = .2) col.sums <- apply(x, 2, sum) row.sums <- apply(x, 1, sum) rbind(cbind(x, Rtot = row.sums), Ctot = c(col.sums, sum(col.sums))) My questions are: 1. I don't understand what values I would need to enter for my 2 variables (data is included at end of this message) in the 'cbind' command 2. For the second row, what modification is required? 3. If I wanted to comput...