similar to: a question related to clustering

Dear all. I need to draw a scatter plot of 23 chromosome copy numbers (y axes) against chromosome and physical location within each chromosome in one plot. The data matrix looks as below: chr location copy_num 1 118345 1.320118 1 3776202 1.133879 1 4798845 0.989997 1 5350951 1.100967 . more data here . . 2 118345 2.459119 2 157739 1.915919 2 1530065 1.924372 2

Dear all. I am trying to buy a hp server to run R and to complete some other tasks with limited bugets. The r-project.org site recommended that R will run on hppa-hp-hpux. However this system is out of our buget and ML system from hp is much cheaper. Is there anybody running R on ML370 or ML 570 systems? If you can provide me with some other information related that will be appreciated.

Hello all, I have a data frame containing 170 columns and would like to generate a new data frame containing all rows but only 67 columns selected according to column names. How can I do this? Thanks in advance, AG Lee [[alternative HTML version deleted]]

Hi, I have built R (current development version) and BioConductor 1.7 with portland group compiler on a AMD Opteron. When I ran qc assessment on Affymetrix latin square data set, I got the following output, Loading required package: affy Loading required package: Biobase Loading required package: tools Welcome to Bioconductor Vignettes contain introductory material. To view,

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 am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

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 am getting the following error in my script. I am very very new to R and have obtained this script from another person. #read file in (dummy data) starburst.plot<-function(affy.fold, affy.FDR)(ifelse( ((affy.fold) >=0), -1*log10(affy.FDR), 1*log10(affy.FDR))) starburst.plot<-function(meth.fold, meth.FDR)(ifelse( ((meth.fold) >=0), -1*log10(meth.FDR), 1*log10(affy.FDR))) At my next

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

Hi, I decided to upgrade to 1.8 today... :-) Anyway, we are writing our own package that is dependent on a bioconductor library - 'affy'. I've checked and when I fire up R, library(affy) behaves as expected... so it all seems to be installed and OK... In the DESCRIPTION file in my package source I have the line: Depends: affy When I run R CMD check simpleaffy I get to: ... *

Hello all, I am trying to make a genome association plot for p-values related to SNPs and was fortunate to find that R contains a package that produces Manhattan plots which is what's preferred for my current project. The function mhtplot() is found in the 'gap' package which I installed in R 2.11.1 on Windows. I thought I'd test out the function first with the examples they

I have run into a problem loading a just saved R object using R-devel. I have been saving and loading this particular type of R object for a long while and never ran into this problem. I save, then immediately reload (to test save) and get "ReadItem: unnknown type 65". This error is reproducible after logout from server and restart of emacs and R. Below is my output and

I have run into a problem loading a just saved R object using R-devel. I have been saving and loading this particular type of R object for a long while and never ran into this problem. I save, then immediately reload (to test save) and get "ReadItem: unnknown type 65". This error is reproducible after logout from server and restart of emacs and R. Below is my output and

I am trying to upload into R 143 Affymetrix chips onto using R on the NIH Nimbus server. I can load 10 chips without a problem, however, when I try to load 143 I receive a error message: cannot create a vector of 523263 KB. I have expanded the memory of R as follows: R --min-vsize=10M --max-vsize=2500M --min-nsize=10M -max-nsize=50M (as specified in help in R). After running this command the

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 R-ers, I am struggling with my x-axis in a association plot. What I would like is to place the labels of the x-axis between the tick markers and normally the labels are printed at the place where the tick marker is placed. I don???t want to move the tick marker (it gives the switch between one chromosome and the next) but I just want to put the chromosome number in between the

Hi, all. Something strange happened to my ESS. I use Mac OS X. When R is used in stand alone mode, it recognizes external libraries like "affy", "e1071" etc. However, when R is invoked through ESS in Emacs, it produces errors like: Error in library(affy) : There is no package called 'affy' Interestingly, I remember I successfully used these libraries through ESS in

Dear list, I?m using R embedded in another program (coded in tcl/tk) under Windows Vista. In this context I don?t launch Rgui.exe but rather Rtem.exe. Now I have a problem finding libraries not specifically installed as administrator (and which are not in the path ?program files? but in Contacts\Documents of the current user). To be precise, the user launching my tcl/tk program can?t find

Hi, I am looking for a way to find out the version information of installed R packages. ex: >library(affy) >SOME_COMMAND_FOR_VERSION(affy) I know I can do either getRversion() or R.Version() could give me the version of R. I assume there must be a way for me to get version information for the add-on paclages,right? Thanks! Yupu