similar to: select affy probes of cancer-related genes

Hi there I have a problem with using Sweave in combination with the option driver = cacheSweave. Whichever code I try to run - when it comes to converting the tex file into pdf it comes up with the same errors (\csname \endcsname errors). Does anybody have an idea what it going wrong? > Sweave("pgfSweave-example.Rnw",driver = cacheSweaveDriver) Writing to file

Hi there, I recently switched to Mac and I wonder if there is any way to get the R console to autocomplete names of list objects or slots of objects. The command line version allows to type e.g. list$ and two times the tab button and then comes up with the names of the list objects. Same for slots of e.g. Eset objects like eset. Unfortunately this does not work in the Mac R console. Is there a

Hi there, I wonder if there is a way of efficiently generating a correlation matrix of two expression matrices. I want to correlate miRNA and mRNA expression and used the following code: ##dat.mi miRNA expression matrix, dat.m mRNA expression matrix nc <- nrow(dat.mi) cor.mat <- data.frame(rep(NA,nrow(dat.m))) pval.mat <- data.frame(rep(NA,nrow(dat.m))) for(i in 1:nc) { cr <- vector()

Hi, I have been using the prcomp function to perform PCA on my example microarray data, (stored in metric text files) which looks like this: 1a 1b 1c 1d 1e 1f ...................................................4r 4s 4t g1 1.2705 1.2766 ...........................................................2.0298 g2 0.1631

Hi, I have a question concerning the analysis of some affymetrix chips. I downloaded some of the data from GEO GSE11324 (see below). In doing so I'm stuck after I identified the probesets with significant changes. I have problems in assigning probeset specific gene names as well as getting the genomic coordinates. Furthermore I have no clue how to deal with the fact, that most genes have

For any given pre-specified gene or short list of genes, yes the Cox model works fine. Two important caveats: 1. Remeber the rule of thumb for a Cox model of 20 events per variable (not n=20). Many microarray studies will have very marginal sample size. 2. If you are looking at many genes then a completely different strategy is required. There is a large and growing literature; I like Newton

Hi there is it possible that pdfs generated using the pdf() function with default settings leads to loss of information? I was plotting copy number changes from Agilent 180k data in form of rectangles (rect()) while each rectangle represents one region of copy number change. When plotting into a pdf I noticed that some very small rectangles do not appear (even after extensive zooming) in the pdf

Hi all, I tried to do normalization of affymetrix data with bioconductor on a Linux server. When I read in the cel files all seemed ok. But the next step caused an error. With Win XP all works fine. Did anyone experience similar problems? Thanks, Thomas > PI <- ReadAffy() > PI AffyBatch object size of arrays=712x712 features (14 kb) cdf=ATH1-121501 (??? affyids) number of

Hello, I am trying to install affy package in R as follow: >R CMD INSTALL -l lib ~/rstuffs/affy_1.4.32.tar.gz Then I get an error at the end: Warning message: There is no package called 'Biobase' in: library(package, character.only = TRUE, logical = TRUE, warn.conflicts = warn.conflicts, [1] "ProgressBarText" [1]

Hello, I have a problem with read.affy.mixed function. I want to read in together a set of CEL files from chip types Affymettrix HGU133A_2 and HGU133_Plus_2. I have my files to be read in in one directory together with a white space delimited file describing them (covdesc). In this directory I give a command: > merge <- read.affy.mixed() Error in merged[[i]] : subscript out of

Hi I try parallelising some code using the snow package and the following lines: cl <- makeSOCKcluster(8) pfunc <- function (x) (if(x <= (-th)) 1 else 0) ###correlation coefficient clusterExport(cl,c("pfunc","th")) cor.c.f <- parApply(cl,tms,c(1,2),FUN=pfunc) The parApply results in the error message: > cor.c.f <- parApply(cl,tms,c(1,2),FUN=pfunc) Error

Dear All, There is a question I met when using Affy package in Bioconductor. I asked it in BioC and didn't get any responses. Sorry to post again: Could anyone tell me how to draw a deep-blue Affymetrix image through Image() function in Affy package? The default settings of image() draw me a black-white image and if I modify it to 256 colors, I get a somehow yellowish image. The reason

Does anyone have nice quality controlls for affy arrays,.... Can't find any tools as are being used for 2 dye arrays..... cheers, marinus This e-mail and its contents are subject to the DISCLAIMER at http://www.tno.nl/disclaimer/email.html [[alternative HTML version deleted]]

I would like to build a package for the HT HG-U133+ PM arrays from affy, but I can't find any good documentation on how to go about it. Naively using makecdfenv's make.cdf.package() causes R to seg-fault. I'm unfamiliar with the CDF format as such, but I'm guessing that it's changed somewhat because the PM arrays no longer have P/A and mismatches. I'm looking to build

Hi: I have clustered microarray gene expression data and trying to map between microarray probe, gene, pathway, gene ontology, and homology for a set of (affy) microarray probes. Is there any package in R which facilitates this? I am looking at bioconductor, but till now could not find a solution. A link to some worked example would be appreciated. Thanks and regards. John [[alternative HTML

Hello All, I need your help. I am analysing affymetrix data and have to select the probe-set that has median expression among all the probe-sets for same gene. This way I want to remove the redundancy by keeping the analysis to single gene entry level. I am fully aware that it is not a nice thing to do but I just have to do it. To do so, I came across 'findLargest' function of

Dear all, I was trying to normalize a subset of affy data those transcribe are either P or M (called pm_filter). I am able to normalize pm_filter subset by using RMA method, however MAS5 is not working. For RMA method, I used the following commend: est<-rma(affydata, subset=pm_filter). Could any help me; how do I do this by using MAS5 method? Regards, Anup Som -- Dr. Anup Som,

Hi , I was trying read in an affy array/matrix and then to calculate the mean value for all the duplicates. Is there a function in the R package that would do this? I tried the help function and searched for 'duplicate'. Although it provides a list of functions that would eliminate the duplicate probe ids, I couldn't find one that would calculate the mean for the duplicates. Any

hi all i started recently using R and i found myself stuck when i try to analyze microarray data. i use the "affy" package to obtain the intensities of the probes, i have two CTRs and two treated. HG.U133A.Experiment1.CEL HG.U133A.Experiment2.CEL HG.U133A_Control1.CEL HG.U133A_Control2.CEL 1007_s_at 2156.23115 467.75615 364.60615 362.11865

I'm analyzing data using Random Forest Regression. For some of the species I am analyzing, the percent variation explained is negative. Could you please explain to me what that means? If you need more information, please let me know. Thank you. Sincerely, Rachel Unger [[alternative HTML version deleted]]