similar to: principal component analysis in affy

Hi, The file I am reading is a text file, whose contents are a matrix that has 15 rows and 58 columns. The first row has column names, and the first column has row names, so the format is correct as far as using read.table is concerned. The other values in the table are all float values (numeric). So when I read in the file using data1 <- read.table("HAL001_HAL0015_Signals.txt"), it

I have been using R and Perl. When I read in a text file using the read.table option, and I try to mathematically manipulate the individual elements in the table, I keep getting an "object is not subsettable" error. If I try to use a different method, it works, but takes too much time(basically, I then need to read in values individually into R instead of as a 2D array, so the number of

Hi, I've just started using R and RSPerl. I have some code as follows: &R::initR("--no-save"); &R::call("t.test", (\@array1, \@array2)); where @array1 and @array2 are both 1-dimensional arrays in Perl having 54675 elements each. On execution the output is as follows: Calling R function name `t.test', # arguments: 3 1) Arg type 3 Got a reference to a

I have passwords turned off, and require keys to match. The zombie armies swarming outside are trying brute force attacks that in part involve guessing login NAMES. If they guess the wrong NAME, this is logged in syslog. If they guess a working user name, then the attack has PARTIALLY SUCCEEDED, but this information is IGNORED. That is, it is not logged. If the zombie army has tell when it

Hi, I have an advance search page and its code is as follows - I have written my code like this <div class="container"> <%= form_tag search_index_path, method: :get do %> <%= radio_button_tag ''user_type'', ''customer'' %><p>Customer</p> <%= radio_button_tag ''user_type'',

Hi all Maybe someone knows a way to solve this anomaly in sample(): I like to compute a sample (n=100) with replications from a population of 2500 units but if I draw repeated samples from it I dont get what seems to be a representative sample if I look at other partitions of the population. Enclosed is the population g99 with 4 columns: (units, partition 1 (site), partition 2 (type), weights);

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

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

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,

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

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

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

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

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

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]

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

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