similar to: How to install affy package in R?

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

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

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

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

Hi! I would like to select probes (affy expression set) of genes that are "cancer-related". Conventionally the decision whether a gene is cancer-related or not is made by looking up the literature. Since this is not possible for all genes on the array I wonder if there is a way of doing this automatically? Best wishes Kristian -- _____ Dr Kristian Unger Imperial College London

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

dear all, I try to update my Rversion fro 1.7.1 to 1.8.1 on linux. I have a problem when I try to install pacakges from Bioconductors with : source("http://www.bioconductor.org/getBioC.R") and getBioC(relLevel="release") somme errors append and when I make a library(affy) for example I obtain : > library(affy) Error in setIs("character",

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

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

Dear all, The package "affy" has the following statement in file "AffyBatch.R": if (debug.affy123) cat("-->initAffyBatch\n") This is great and I would also like to use it. However, when I run my package I get the following error: Error in .local(object, ...) : object "debug.mypkg" not found Since I am not able to find the position where

Hi, i was looking into the documentation for the rma() function in affy() package, and was trying to figure out how exactly the background normalization is done. I read all three papers cited in the rma() documentation, but the most detailed explanation i could find was in Irizary et al., 2003, where they state that they compute B(PM_{ijn}) = E[s_{ijn} | PM_{ijn}] where s_{ijn} is assumed to

When I try to load cel files (hgu133a) using the ReadAffy() in R 1.8.1 command I get an error message: > x<-ReadAffy() Error: cannot allocate vector of size 102973 Kb Does anybody know what does this error mean and how to overcome it? I have tried to load the same data with R 1.7.1 and it works. There is also no error when I use R 1.8.1 to load moe430 cel files. Thanks very much for any

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 R Users, In the expresso function, which combination of these methods for data pre-processing (when using affymetrix oligo arrays) is the best: bgcorrect.metod = rma rma2 mas normalize.method = qspline quantiles loess pmcorrect.method = pmonly subtractmm mas summary.method = liwong avgdiff medianpolish mas There are many options within each method. I would appreciate a hint on the best

Estimados miembros de la lista, Estoy haciendo una análisis de supervivencia con R. Adjunto mis datos. Quiero analizar la supervivencia de 5 grupos diferentes y compararla. Para ello estoy utilizando el paquete survival. > s = Surv(c$tiempo, c$estado) > f = survfit(s ~ tratamiento, data = c) > d = survdiff(s ~ tratamiento, data = c) > d Call: survdiff(formula = s ~ tratamiento,

Hola, Manuel, No entiendo tu pregunta (la repito aqui para que sea mas explicito): hay alguna forma de comparar la probabilidad de supervivencias (en este caso anual) entre grupos sin utilizar un chi-cuadrado y un valor de P. Entiendo que lo que hace survdiff es comparar las curvas de supervivencia, pero yo quiero comparar la probabilidad de supervivencia entre grupos al final del estudio. Con

Hi!I'm doing a little data importing from .cel files, > setwd("/home/mandova/celfiles") > mydata<-ReadAffy() Error in sub("^/?([^/]*/)*", "", filenames, extended = TRUE) : unused argument(s) (extended = TRUE) Then I tried > filenames<-paste("GSM",c(seq(138597,138617,1)),".cel",sep="") >

I would like to have a function that inputs a file name, then reads that data file and does various stuff to it. Something like: doit<-function(filename) { x<-scan("~/rstuff/data/"filename) #for example, if filename is a0.dat, it will scan("~/rstuff/data/a0.dat") ... } Can someone please tell me how to do this? Thanks for any help! Bill Simpson