similar to: RE: Comparison of correlation coefficients - Details

Dear all I apologize for cross-posting, but first it is accepted custom to thank the repliers and give a summary, and second I have still the feeling that this problem might be a general statistical problem and not necessarily related to microarrays only, but I might be wrong. First, I want to thank Robert Gentleman, Mark Kimpel and Mark Reiners for their kind replies. Robert Gentleman kindly

Dear Spencer Thank you for this extensive explanation of the problem. I was just curious. Best regards Christian ============================================== Christian Stratowa, PhD Boehringer Ingelheim Austria Dept NCE Lead Discovery - Bioinformatics Dr. Boehringergasse 5-11 A-1121 Vienna, Austria Tel.: ++43-1-80105-2470 Fax: ++43-1-80105-2683 email: christian.stratowa at

Dear R-experts I would like to read all files from a directory, the files have names "myname0001.txt" etc. I paste the directory plus file names and use "read.delim()". My problem is that some file names are missing, so I get an error and my program stops. Is there a way to check for a null pointer analogous to C, so that I can simply skip non-existent filenames? Please do

Dear all, First I apologize for cross-posting, but I think that this could be of interest to BioC users, too. For DNA copy-number analysis I have downloaded PLASQ500K from: http://genome.dfci.harvard.edu/~tlaframb/PLASQ/ After creating sub-directories SND and STD containing 3 Sty Mapping arrays each, I tried to compute parent-specific copy number: > library(PLASQ500K) > psCN <-

Dear all Just for fun, I have just downloaded the paper mentioned below and checked it with R-1.6.1. Everything is ok with exception of Table 2b, where I get always 1 instead of 0.5: > pbinom(1e15,2e15,0.5) [1] 1 Which value should be correct? Best regards Christian Stratowa ============================================== Christian Stratowa, PhD Boehringer Ingelheim Austria Dept NCE Lead

Hello, I'm fitting linear mixed models to gene-expression data from microarrays, in a data set where 4608 genes are studied. For a sample of 5 subjects and for each gene we observe the expression level (Intensity) in four different tissues: N, Tp, Tx and M. I want to test whether the expression level is different accross tissues. Between-subject variability is modeled with a random

Hey. Is ist possible to guarantee CPU Performance to virtual machines? E.g.: I have a 4 CPU host machine. 3 guests are running. Each has 2 CPU. When VM1 needs more resscources than the other 2, can i guarantee that VM1 gets the most Power of the system at this time? Thomas Diederich ************************************************** * Boehringer Ingelheim Pharma GmbH & Co.KG * A

Does anyone know if Xen is currently supporting Itanium on 64 Bit!? Thomas Diederich ************************************************** * Boehringer Ingelheim Pharma GmbH & Co.KG * A Informationsverarbeitung / Diplomant Systemtechnik * * Mail: diederit@ing.boehringer-Ingelheim.com <mailto:diederit@ing.boehringer-Ingelheim.com> * Phone: +49 (0)6132/77-98151

Hey, just pulled the new updates from the repository, and did a male dist, make install and an ./install.sh Then i made an initrd an bootet the system. But when I do an lsmod, it shows me: NOTHING! What''s wrong there?! Thomas Diederich ************************************************** * Boehringer Ingelheim Pharma GmbH & Co.KG * A Informationsverarbeitung / Diplomant

Hi, I have a problem executing "library(ROracle)" in R: OS/Software: Redhat 8.0, all available patches applied Oracle 9i v9.2.0.1.0 R v1.7.0 ROracle v0.5-0 DBI v0.1-5 The compilation and installation of ROracle went fine. However when I try to load ROracle I get the following: >> library(ROracle) >Error in dyn.load(x, as.logical(local), as.logical(now)) : > unable

Cd /xen/xen-unstable/linux/2.6.12-xen0/ Then: make ARCH=xen xconfig Thomas Diederich ************************************************** * Boehringer Ingelheim Pharma GmbH & Co.KG * A Informationsverarbeitung / Diplomant Systemtechnik * * Mail: diederit@ing.boehringer-Ingelheim.com * Phone: +49 (0)6132/77-98151 -----Ursprüngliche Nachricht----- Von:

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

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

> Dear expeRts, > > Currently, my colors are as follows: > mycol <- > c("blue1","blue2","blue3","blue4","black","yellow4","yellow3","yellow2","y > ellow1") > heatmap(snp, Rowv=NA, Colv=NA, col=mycol) > > However, I would like to have the following colors: > bright blue ->

dear members, i would like to write a variable in a plot title (main="") but i don't know the right syntax:(...i tried a lot of different ways without success. here my example: y=30 z=33 for (i in 10:length(tissue)) { png(filename = tissues[i], width = 1024, height = 768, pointsize = 12, bg = "white") gene.graph("ENSG00000115252", rma.affy, gps=list(1:3,

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

Hello, I will start out with the caveat that I'm not a statistician by training, but have a fairly decent understanding of probability and likelihood. Nevertheless, I'm trying to fit a nonlinear model to a dataset which has two main factors using nlme. Within the dataset there are two Type categories and four Tissue categories, thus giving me 8 datasets in total. The dataset is in

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