similar to: lme4: How to specify nested factors, meaning of : and %in%

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

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

Hi, I make this model using lme m.lme <- lme(Glycogen~Treatment,random=~1|rTrt/Liver) How to make this using lmer? I try > m.lmer <- lmer(Glycogen~Treatment+(1|rTrt/Liver)) Erro em lmer(Glycogen ~ Treatment + (1 | rTrt/Liver)) : entry 0 in matrix[0,0] has row 2147483647 and column 2147483647 Al??m disso: Mensagem de aviso: / not meaningful for factors in: Ops.factor(rTrt, Liver)

Hi, I've been doing an experiment, measuring the dead-zone-diameters of bacteria, when they've been grown with paper diffusion disks of antimicrobial. There are two groups, or treatments - one is bacteria that have been cultured in said antimicrobial for the past year, the other group is of the same species, but lab stock and has not gone had any prior contact with the antimicrobial.

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

Hi, I''m trying to understand the lme output and procedure. I''m using the Crawley''s book. I''m try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :125.0 Min. :1 Min. :1.0 Min. :1 1st Qu.:135.8

Hi All, I would like to ask a question on fixed and random effecti in lme. I am fiddlying around Mick Crawley dataset "rats" : http://www.bio.ic.ac.uk/research/mjcraw/statcomp/data/ The advantage is that most work is already done in Crawley's book (page 361 onwards) so I can check what I am doing. I am tryg to reproduce the nested analysis on page 368:

Hi! How to calculate the correct SE of mean in a nested or spliplot anova? Nested example: --------------------- m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > m Call: aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver)) Grand Mean: 142.2222 Stratum 1: Treatment Terms: Treatment Sum of Squares 1557.556 Deg. of Freedom 2 Estimated

Dear All, I'm not sure of the interpretation of interactions with contrasts. Can anyone help? I do an ANCOVA, dryweight is covariate, block and treatment are factors, c4 the response variable. model<-aov(log(c4+1)~dryweight+treatment+block+treatment:block) summary(model); Df Sum Sq Mean Sq F value Pr(>F) dryweight 1 3.947 3.947 6.6268 0.01076 *

Perhaps you can see somethign I can't - or perhaps there is a better way for me to get information for you ? Let me know if there is as this server is not live yet... All the best, Noel NB ATTACHMENTS REMOVED FOR LIST POSTING Domain/network info: Domain = UK Win2000 DC (192.168.5.4) = BRAIN Live Samba server 2.2.3a (192.168.5.5) = BELLY New Samba server 2.2.3a

Hi all! I've got this error while running example(Glm) library("rms") > example(Glm) Glm> ## Dobson (1990) Page 93: Randomized Controlled Trial : Glm> counts <- c(18,17,15,20,10,20,25,13,12) Glm> outcome <- gl(3,1,9) Glm> treatment <- gl(3,3) Glm> f <- glm(counts ~ outcome + treatment, family=poisson()) Glm> f Call: glm(formula = counts ~

Hi everybody: I?m trying to rewrite some routines originally written for SAS?s PROC NLMIXED into LME4's glmer. These examples came from a paper by Nelson et al. (Use of the Probability Integral Transformation to Fit Nonlinear Mixed-Models with Nonnormal Random Effects - 2006). Firstly the authors fit a Poisson model with canonical link and a single normal random effect bi ~ N(0;Sigma^2).The

Hello, unfortunately I have I big problem I can't solve. I have to analyse if a gene is tissue specific. For example for the gene xyz I have following expression values: Heart Liver Brain 8.998497 10.013561 12.277407 9.743556 10.137574 11.033957 For every tissue I have two values from two different experiments. Now I want to test if Heart is significant higher

Hi group, I'm trying to do a Tukey test to compare the means of a factor ("treatment") with three levels in an lme model that also contains the factors "site" and "time": model = response ~ treatment * (site + time) When I enter this model in csimtest, it takes all but the main factor "treatment" as covariables, not as factors (see below). Is it

Hello all! I am once again analyzing patient survival data with chronic liver disease. The severity of the liver disease is given by a number which is continuously variable. I have referred to this number as "meld"--model for end stage liver disease--which is the result of a mathematical calculation on underlying laboratory values. So, for example, I can generate a Kaplan-Meier plot

Hi, I'm reading van Belle et al "Biostatistics" and trying to run a cox test using a dataset from: http://faculty.washington.edu/~heagerty/Books/Biostatistics/chapter16.html (Primary Biliary Cirrhosis data link at top of the page), I'm using the following code: --------------- start of code library(survival) liver <-

Hi, I have a experiment with block, plots, sub-plots, and sub-sub-plots with repeated measures and 3 factors (factorial design) when we have been observed diameter (mm), high (cm) and leaves number (count). However, we don't have one treatment in one factor, so, my design is unbalanced. On a previous message here, a friend tell me that "It appears to me that your design is a split-split

Hello, I am just a beginner of R 1.9.0. I try to construct a predictive score for the development of liver cancer in cirrhotic patients. So dependant variable is binanry (cancer yes or no). Independant variables are biological data. The aim is to find out a cut-off value which differentiate (theoratically) from normal to pathological state for each biological data. How can I step in procedue to

Hi, I have a question about how to do covariate adjustment. I have two sets of 'gene expression' data. They are from two different tissue types, 'liver' and 'brain', respectively. The purpose of my analysis is to compare the pattern of the whole genome 'gene expression' between the two tissue types. I have 'age' and 'sex' as covariates. Since

Dear All I am trying to do a repeated measures analysis using lmer and have a number of issues. I have non-orthogonal, unbalanced data. Count data was obtained over 10 months for three treatments, which were arranged into 6 blocks. Treatment is not nested in Block but crossed, as I originally designed an orthogonal, balanced experiment but subsequently lost a treatment from 2 blocks. My