similar to: How to do covariate adjustment in R

Hi, I am doing an analysis to see if these is tissue specific effects on the gene expression data . Our data were collected from 6 different labs (batch effects). lab 1 has tissue type 1 and tissue type 2, lab 2 has tissue 3, 4,5,6. The other labs has one tissue type each. The 'sample' data is as below:

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

Dear fellow R-help members, I hope to seek your advice on how to parse/manage a dataset with hundreds of columns. Two examples of these columns, 'cancer.problems', and 'neuro.problems' are depicted below. Essentially, I need to parse this into a useful dataset, and unfortunately, I am not familiar with perl or any such language. data <- data.frame(id=c(1:10))

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

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

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

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:

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

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 +

Dear R users, I have a data frame in the form below, on which I would like to make normality tests on the values in the ExpressionLevel column. > head(df) ID Plant Tissue Gene ExpressionLevel 1 1 p1 t1 g1 366.53 2 2 p1 t1 g2 0.57 3 3 p1 t1 g3 11.81 4 4 p1 t2 g1 498.43 5 5 p1 t2 g2 2.14 6 6 p1 t2 g3 7.85 I

My small brain is having trouble getting to grips with lme() I wonder if anyone can help me correctly set the random = argument to lme() for this kind of setup with (I think) 9 variance/covariance components ... Study.1 Study.2 ... Study.10 Treatment.A: subject: 1 2 3 4 5 6 etc. 28 29 30 Treatment.B: subject: 31

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

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

Dear r helpers, I would like to look at the interaction between two two-level factors, one between and one within participants, after accounting for any variance due to practice (31 trials in each of two blocks) in the task. It seems to require treating practice as a covariate. All the examples I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's "Practical

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 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

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 <-

Dear R users, I have a huge database and I need to adjust it somewhat. Here is a very little cut out from database: CODE NAME DATE DATA1 4813 ADVANCED TELECOM 1987 0.013 3845 ADVANCED THERAPEUTIC SYS LTD 1987 10.1 3845 ADVANCED THERAPEUTIC SYS LTD 1989 2.463 3845 ADVANCED THERAPEUTIC SYS LTD 1988 1.563 2836 ADVANCED TISSUE