similar to: difficulty setting the random = argument to lme()

Hello, I'm using aov() to analyse changes in brain volume between males and females. For every subject (there are 331 in total) I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). The data looks like this: Subject Sex Lobe Tissue Volume subect1 1 F g 262374 subect1 1 F w 173758 subect1 1 O g 67155 subect1 1 O w 30067 subect1 1 P g 117981

I'm working with a nested model (mixed). I have four factors: Patients, Tissue, sex, and tissue_stage. Totally I have 10 patients, for each patient, there are 2 tissues (Cancer vs. Normal). I think Tissue and sex are fixed. Patient is nested in sex,Tissue is nested in patient, and tissue_stage is nested in Tissue. I tried aov and lme as the following, > aov(gene ~ tissue + gender +

Hi Dear all, I have some gene expression data samples from different tissue types ----------------------------------------------- - 120 samples from blood (B) - 20 samples from Liver (L) - 15 samples from Kidney (K) - 6 samples from heart (H) ----------------------------------------------- All the samples are from different individuals, so there are in total 161 individuals from which the DNA was

Hi all (really probably just Deepayan): In the plot below I want to add text on either side of each violin plot that indicates the number of observations that are either positive or negative. I'm trying to do this with ltext() and I've also monkeyed about with panel.text(). The code below is generally what I want but my calls to ltext() are wrong and I'm not sure how to fix them.

Hi, There are 20 subjects grouped by Gender, each subject has 2 tissues (normal vs. cancer). In fact, it is a 2-way anova (factors: Gender and tissue) with tissue nested in subject. I've tried the following: Model 1: lme(response ~ tissue*Gender, random = ~1|subject) Model 2: response ~ tissue*Gender + subject Model 3: response ~ tissue*Gender It seems like Model 1 is the correct one

Hi I'm having difficulty working out how to get what I think is the appropriate partitioning of variability in a repeated measures setup. I have G=5 treatment-groups, each containing n=6 subjects, and a response is measured on each subject on t=4 occasions. I think the anova degrees of freedom should partition as follows - Between-subjects: G*n-1=29 [ between-groups: g-1 = 4 ,

this is a clearer (I hope) version of an earlier post - My problem is formulating the random = argument to give estimates of all 9 random components for this kind of setup where there are (I think) 9 variance/covariance components ... Study.1 Study.2 ... Study.5 Treatment T1: subject: 1 2 3 4 5 6 ... 13 14 15 Treatment T2: subject: 16 17

Hi, I'm analysing a dataset in which the same 5 subjects (male.pair) were subjected to two treatments (treatment) and were measured for 12 successive days within each treatment (layingday). Overall 5*2*12=120 observations. I want to test the effect of treatment, time (layingday) and their interaction. I have done so through the ANOVA below: >

Dear friends - We have 25 rats, 14 of these subjected to partial removal of kidney tissue, 11 to sham operation, and then followed for 6 weeks. So far we have data on 26 urine metabolites measured by NMR 7 times during the observation. I have smoothed the measurements by b.splines in fda including a roughness penalty, and inspecting the mean curves for nephrectomized and sham animals indicate

Hello, I've used this command to analyse changes in brain volume: mod1<-aov(Volume~Sex*Lobe*Tissue+Error(Subject/(Lobe*Tissue)),data.vslt) I'm comparing males/females. For every subject I have 8 volume measurements (4 different brain lobes and 2 different tissues (grey/white matter)). As aov() provides only type I anovas, I would like to use lmer() with type II, however, I have

Hello I am trying to fit a REML to some soil mineral data which has been collected over the time period 1999 - 2004. I want to know if the 19 different treatments imposed, differ in terms of their soil mineral content. A tree model of the data has shown differences between the treatments can be attributed to the Magnesium, Potassium and organic matter content of the soil, with Magnesium being the

Dear R list: I modelled changes in a variable (mconc) over time (d) for individuals (replicate) given one of three treatments (treatment) using: mconc.lme <- lme(mconc~treatment*poly(d,2), random=~poly(d,2)|replicate, data=my.data) summary(mconc.lme) shows that the linear coefficient of one of the treatments is significantly different to zero, viz. Value Std.Error

Hi, I posted the message below last week, but no answers, so I'm giving it another attempt in case somebody who would be able to help might have missed it and it has now dropped off the end of the list of mails. I am fairly new to R and still trying to figure out how it all works, and I have run into a few issues. I apologize in advance if my questions are a bit basic, I'm also no

Hi, I have searched the archives and can't quite confirm the answer to this. I appreciate your time... I have 4 treatments (fixed) and I would like to know if there is a significant difference in metal volume (metal) between the treatments. The experiment has 5 blocks (random) in each treatment and no block is repeated across treatments. Within each plot there are varying numbers of

Dear all, I would like to know if it is possible to fit in R a Cox ph model with time-dependent covariates and to account for hierarchical effects at the same time. Additionally, I'd like also to know if it would be possible to perform any feature selection on this model fit. I have a data set that is composed by multiple marker measurements (and hundreds of covariates) at different time

Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a lack of understanding. We have growth curves of new dental bone that can well be modeled by a linear growth curve, for two different treatments and several subjects as random parameter. By definition, newbone is zero at t=0, so I tried to force the

Hello fellow R users, I am having a problem finding the estimates for some overall treatment effects for my mixed models using 'lme' (package nlme). I hope someone can help. Firstly then, the model: The data: Plant biomass (log transformed) Fixed Factors: Treatment(x3 Dry, Wet, Control) Year(x8 2002-2009) Random Factors: 5 plots per treatment, 5 quadrats per plot (N=594 (3*5*5*8)

Dear R-Community! For example I have a study with 4 treatment groups (10 subjects per group) and 4 visits. Additionally, the gender is taken into account. I think - and hope this is a goog idea (!) - this data can be analysed using lme as below. In a balanced design everything is fine, but in an unbalanced design there are differences depending on fitting y~visit*treat*gender or

Hi, I need help creating a code for a multiple BACI design (Before-After Control-Impact) ANOVA. I'm new to R and basically need to run a complex mixed model ANOVA that treats location as a random factor. Data are from a fire experiment, run 2001-2004 (2 years pre, 2 years post). Response is bird abundance. 4 Treatments had 3 replicates each (forest stands): 1. Control, 2. Prescribed fire

Hi I am running some linear and non-linear mixed effect models and would like to do some planned contrasts (a priori contrasts) I have looked in the help and in many forums and it seems possible to do so but don't understand how to write the function and I couldn't find an example in Pinheiro and Bates. lme {nlme} has a contrasts argument but I can't understand how to code it.