similar to: weights in lmer

Hi I have data from 12 subjects. The measurement is log(expression) of a particular gene and can be assumed to be normally distributed. The 12 subjects are divided into the following groups: Infected, Vaccinated, Lesions - 3 measurements Infected, Vaccintaed, No Lesions - 2 measurements Infected, Not Vaccinated, Lesions - 4 measurements Uninfected, Not Vaccinated, No Lesions - 3 measurements

Hi everyone, I know there have been several requests regarding subsetting before, but none of them really helps with my problem: I'm trying to subset only infected individuals from the REC2 data.frame: > str(REC2) 'data.frame': 362 obs. of 7 variables: $ RINGNO : Factor w/ 370 levels "BL17546","BL17577",..: 78 81 67 41 58 66 17 $ year : Factor w/ 8

Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the same output in all 3 cases below (forgetting about the convergence problem). Would there be a reason to prefer labeling the levels one way or

Ronaldo, Further to my previous posting on your Glycogen nested aov model. Having read Douglas Bates' response and Reflected on his lmer analysis output of your aov nested model example as given.The Glycogen treatment has to be a Fixed Effect.If a 'treatment' isn't a Fixed Effect what is ? If Douglas Bates' lmer model is modified to treat Glycogen Treatment as a purely

Sad to say one of my file servers was exploited and used to run a Phishing scam. Have identified subject virus amongst other things. It appears twice in a virus scan; /sbin/z (which I assume can just be deleted) and /sys/bus/serio/drivers/atkbd/description. The latter file is also present in identical uninfected machines. I have been unable to open the file, even with root privileges, although

Dear R-help-list, I have a problem in which the explanatory variables are categorical, the response variable is a proportion, and experiment contains technical replicates (pseudoreplicates) as well as biological replicated. I am new to both generalized linear models and mixed- effects models and would greatly appreciate the advice of experienced analysts in this matter. I analyzed the

Hi, I am new to this list. I have a question regarding including both spatial and temporal random factors in lmer. These two are not nested, and an example of model I try to fit is model1<-lmer(Richness~Y+Canopy+Veg_cm+Treatment+(1|Site/Block/Plot)+(1|Year), family=poisson, REML=FALSE), where richness = integer Y & Treatment = factor Canopy & Veg_cm = numerical, continous

Hi all, i just installed MailScanner-4.60.8-1.rpm.tar.gz on centos 5 with sendmail 8.13.8. Whenever I send emails, I get bellow error. Jun 21 15:33:01 mail MailScanner[3657]: New Batch: Scanning 1 messages, 767 bytes Jun 21 15:33:02 mail MailScanner[3657]: Virus and Content Scanning: Starting Jun 21 15:33:02 mail MailScanner[3657]: ERROR: Unable to create temporary directory Jun 21 15:33:02

Hi So carrying on my use of analysis of variance to check for the effects of two factors. It's made simpler by the fact that both my factors have only two levels each, creating four unique groups. I have a highly significant interaction term. In the context of the experiment, this makes sense. I can visualise the data graphically, and sure enough I can see that both factors have different

Hello; I'm running a FreeBSD 4.10-release-p2 box and both chkrootkit 0.44 & 0.45 report that my /sbin/init file is infected. It appears as though the egrep for "UPX" in the output of "strings" triggers the infected notice. When I copy the init file from an uninfected box to this one chkrootkit continues to report it as infected. Is chkrootkit reading a copy of the

Hello list, I'm trying to figure out how exactly the specification of nested random effects works in the lmer function of lme4. To give a concrete example, consider the rat-liver dataset from the R book (rats.txt from: http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/ ). Crawley suggests to analyze this data in the following way: library(lme4) attach(rats) Treatment <-

I'm trying to install the XML package in R 1.6.1 (>install.packages("XML")). The download is OK, but during configuration the file parser.h cannot be found (it looks in libxml/ and gnome-xml/ then gives up). What is missing from my system? -- Thanks, -susan Susan J. Miller Biotechnology Computing Facility Arizona Research Laboratories Bio West 228 University of Arizona

Hello, I am a new R user. I am trying to calculate vector correlations for all pairwise comparisons in my data frame without repeats. I am familiar with the expand.grid function, but this includes repeats. Is there a way to use expand.grid and eliminate repeats? Or is there another function that can be used to do this? Thank you. Rebecca -- Rebecca Young Graduate Student Ecology &

platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 7.1 year 2003 month 06 day 16 language R I have a Dell Inspiron 4000 laptop with Windows ME. I recently installed the 1.7 version of R. The following commands are from my

Greetings, I'm in the painful process of migrating from SAS to R. In the process I've discovered that there are some basic things that I am getting hung up on. The most basic is the simple recoding of variables. Suppose I create a vector x <- rnorm(10000) and I want to recode all values of x > 1.5 to NA. How would I do that in R? TIA Cheers, Patrick

Hi folks, I am working on my student internship to USA, california and would like to attend a R Advanced Programming course while I am in USA. Where can I find such a course. I am new to this list. Trevor O'brien >From Ireland. [[alternative HTML version deleted]]

Is there is an open source ACD component for Asterisk? Likewise, is there an open source CTI component that will work with Asterisk? Regards, Jim O'Brien -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20040504/e6ee98da/attachment.htm

Hello, Can you get eigenvalues in addition to eigevectors using prcomp? If so how? I am unable to use princomp due to small sample sizes. Thank you in advance for your help! Rebecca Young -- Rebecca Young Graduate Student Ecology & Evolutionary Biology, Badyaev Lab University of Arizona 1041 E Lowell Tucson, AZ 85721-0088 Office: 425BSW rlyoung at email.arizona.edu (520) 621-4005

Hi I'm encountering some problems with coxme My data: I'm looking at the survival of animals in an experiment with 3 treatments, which came from 4 different populations, two of which were infected with a parasite and two of which were not. I'm interested if infected animals differe from uninfected ones across treatments. Factor 1: treatment (3 levels) Factor 2: infection state

Hi I?m having difficulty configuring Samba running on FreeBSD. Samba version is 2.2.8a. The scenario is that we have an NT4 domain and I am looking to migrate the printers onto Samba. All of the printers are connected via print servers, mainly HP JetDirect. As far as security is concerned I am happy for guest/anonymous logon as the server will only be used for printing but I have set security to