similar to: Use of nesting in lmer- error in numerical expression

Dear R users, I am trying to fit a GLMM to the following dataset; tab a b c 1 1 0.6 199320100313 2 1 0.8 199427100412 3 1 0.8 199427202112 4 1 0.2 199428100611 5 1 1.0 199428101011 6 1 0.8 199428101111 7 0 0.8 199527103011 8 1 0.6 199527200711 9 0 0.8 199527202411 10 0 0.6 199529100412 11 1 0.2 199626201111 12 2 0.8 199627200612 13 1 0.4 199628100111 14 1 0.8

Dear all, I am trying to understand the output from a binomial lmer object and why the scaling of a fixed effect changes the variance components. In the model p2rec is cbind(number recruits2,number recruits 1), Pop is populations (five level factor) and ja is year (covariate running from 1955-2004). I.e. biologically I am interested to see how the proportion of recruits from the second

Hi, As some of you know, one of the R Projects ideas for GSOC 2009 is to expand the crantastic.org portal. Some ideas for what should be done is mentioned in Hadley Wickham's project description [1], others have been briefly mentioned on this mailing list before. Undoubtedly, however, there are many other good ideas out there that could be implemented. Please suggest any ideas you might have

Overview ======== This module implements the Center for Internet Security (CIS) Security Configuration Benchmark for Red Hat Enterprise Linux 6 v.1.1.0 (avilable at http://benchmarks.cisecurity.org). Each scored control has been implemented as a class or a custom fact. Installation ============ Please either: - Clone git repo from https://github.com/arildjensen/cis-puppet - Run "puppet

Hello all, If someone could take a little time to help me then I would be very grateful. I studied piping plovers last summer. I watched each chick within a brood for 5 minutes and recorded behaviour, habitat use and foraging rate. There were two Sites, the first with 4 broods and the second with 3 broods. http://r.789695.n4.nabble.com/file/n4511178/Table_PP_Maslo_et_al.png As the data within a

Hello, I hope that you can bare with me. I am new to models, but I think I have a pretty godd understanding of how to run them now, including how to use AICc and Anova. The issue is that I have many factors that I wish to compare so doing each one at a time would take forever. I came across the MuMIn package and I was so excited, however I am getting an error message and i don't know why.

Trying to find a quick/slick/easily interpretable way to collapse a data set. Suppose I have a data set that looks like this: h <- structure(list(INDEX = structure(1:6, .Label = c("1", "2", "3", "4", "5", "6"), class = "factor"), TICKS = c(0, 0, 0, 0, 0, 3 ), BROOD = structure(c(1L, 1L, 2L, 3L, 3L, 3L), .Label =

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 not

Hello, everyone. I'm hoping to prevent myself from doing 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

I'm having a little trouble with trying to install newer versions of Starcraft (updates or Brood War) on a win-free Wine installation. I can run the installer off the original Starcraft CD, and it puts all of the files in the right place, adds a couple registry entries, and then chokes on creating the entries in the Start Menu. Still, I can kill the process at that point, and I have a usable

Dear all, This question is basic but I am stumped. After running the below, I receive the message: "non-integer #successes in a binomial glm!" model1 <- glmmML(y~Brood.Size*Density+Date.Placed+Species+Placed.Emerging+Year+rate.of.parperplot, data = data, cluster= data$Patch, family=binomial(link="logit")) My response variable is sex ratio, and I have learned quickly not

Hi: I have the following dataset in R( thanks Gabor for your help) but now the problem is that all the dashes are converted to questions marks ("?") I?am trying to get the dashes back using: mydf <- data.frame(lapply(mydf,function(x) replace(mydf,"?"(x),"-"))) but isn't working. I also tried: mydf <- replace(mydf,mydf=="?","-)")

Hi: I?need?help with a dataframe(see?pic attached). is a mix of dates and text. I want to create a table either?using latex function from hmisc or xtable. I already know how to do this but the problem is getting the dataframe into R. I don't have a reproducible example?but I am hoping that the pic attachment will make it to you.?If someone is interested in helping with this task I could send

Webgen uses file paths to generate multilevel menues, bread crumbs etc. And limits the character range in paths to a-zA-Z. As the directory names are used verbatim as menu labels this is fairly limitting for non-english pages. Is there a way to map directory names to configurable UTF8 strings for menu labels? oao

Dear All I would need to perform a MANOVA with both fixed (group, sex, group*sex) and random (brood) effects. I wonder if this is at all possible and if R does that. At the moment, I only know that I can run a classic MANOVA with R. Thank you David ______________________________________________ David Costantini, PhD http://www.davidcostantini.it NERC Postdoctoral research associate

Hi all, I am trying to create images in R suitable for journal publication. I'd like to make an image that is 3.6" wide and 5.08" length, and 1200 dpi. When I create a tiff file that is 800 dpi and 4x4", it works (although doesn't look good). But when I create a tiff that is the size and dpi I want, it doesn't work. I get the following error: Error in tiff(filename =

Dear all: Does anybody know about label conflicts between xtable and Hmisc? I found a couple of e-mails similar to this problem but is not clear to me how to get around the label problem. The first table(longtable below) is generated with the latex function from Hmisc but for some reason when I try to hyperlink to it,it takes me to the top of the document. The second table created with the xtable

Hello, I have been attempting to set up a lme and have looked at numerous posts including 'R's lmer cheat-sheet' as well as reading a number of papers and other resources including R help, but I am still a little confused on how to write my model (I thought I had it). I have asked a number of questions on different forums; most of which have been resolved. My main concern right now

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

Hello! So, my problem is following. I have bird offspring growth data and I'd like to model individual growth curves (aim is to study asymptotes and inflection points) with nlme according to Pinheiro & Bates 2000: first using nlsList to generate individual curves and then nlme to study the parameters and fixed effects. The data is structured to two levels. I have broods and