[R] how can I plot bar plots with all the bars (negative and positive) in the same direction????

2010-09-03 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Dear r-help mailing list, this seems stupid, but I actually don't find the solution: if I have a vector of numbers x of length n, ranging, say, from -3 to 4, if I do barplot (x) all the values below 0 go downwards, and all the positive values go upward. How can I make them all begin from the m

[R] how to create automatically names for vectors in a loop?

2010-05-27 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, I want to generate a number of vectors and store them with different names, like this: x=1 while (x<100) { vector#x# = rnorm(100) x=x+1 } where each vector has, at its hand, instead of #x# a number which goes from 1 to 99. How can I do this? Thanks Gabriele Zoppoli, MD Ph.D. F

Re: [R] "order" issue

2010-05-23 Thread Zoppoli, Gabriele (NIH/NCI) [G]
for my naive definitions...) R that now all of your values are numeric in a matrix? If you do as.numeric, everything becomes a long column of n Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy Guest Researcher, LMP, NCI, NIH

Re: [R] "order" issue

2010-05-23 Thread Zoppoli, Gabriele (NIH/NCI) [G]
1-451-8575 Mobile: 301-204-5642 Email: zoppo...@mail.nih.gov From: Ted Harding [ted.hard...@manchester.ac.uk] Sent: Sunday, May 23, 2010 6:31 PM To: Zoppoli, Gabriele (NIH/NCI) [G] Cc: R help Subject: RE: [R] "order" issue On 23-May-10 21:39:06, Zo

Re: [R] "order" issue

2010-05-23 Thread Zoppoli, Gabriele (NIH/NCI) [G]
negative number as smaller than positive ones... Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy Guest Researcher, LMP, NCI, NIH, Bethesda MD Work: 301-451-8575 Mobile: 301-204-5642 Email: zoppo...@mail.nih.gov ______________

Re: [R] "order" issue

2010-05-23 Thread Zoppoli, Gabriele (NIH/NCI) [G]
326" "Breast" 38 "LE:MOLT_4" "-3.20055" "-1.72841" "Leukemia" 36 "LE:HL_60" "-3.16745" "-3.16745" "Leukemia" 54 "OV:OVCAR_4" "-3.13137" "-0.47497" "Ovaria

[R] "order" issue

2010-05-23 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi everybody, this is a real dummy thing. I sorted a matrix based on a given column, and what I get is right, until it comes to columns of negative and positive values; than, "order" orders everything from max to min in the negative values, and then AGAIN from max to min in the positive values!

[R] memory limit

2010-05-08 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, when I try to import a microarray CEL batch, I get this error message: > myAB <- ReadAffy () Error in .Call("read_abatch", filenames, rm.mask, rm.outliers, rm.extra, : cannot allocate vector of length 1287151200 which, assuming the value is in bites, is below my RAM values (3 Gb recogni

[R] cell-based high-throughput chemical compound screens package?

2010-04-22 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi I'm looking for a package to perform quality control, normalization and analysis of high throughput cell-based chemical screens. I know that the cellHTS2 package provides this for siRNA screens. Does anybody know if something like what I'm looking for exists? Thank you! Gabriele Zoppoli,

[R] simple question

2010-04-22 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi how can I find, in a vector of characters, which is the most frequent one? Thanks Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy Guest Researcher, LMP, NCI, NIH, Bethesda MD Work: 301-451-8575 Mobile: 301-204-5642 Em

[R] cell-based high-throughput chemical compound screens package?

2010-04-22 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi I'm looking for a package to perform quality control, normalization and analysis of high throughput cell-base chemical screens. I know that the cellHTS2 package provides this for siRNA screens. Does anybody know if something like what I'm looking for exists? Thank you! Gabriele Zoppoli,

[R] exporting rotating plot into Powepoint

2010-04-06 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, I need to export a rotating plot like the one generated from the following function into a PPT slide or something that can be viewed by reviewers of a journal: library(np) library(datasets) data(faithful) f <- npudens(~eruptions+waiting,data=faithful) plot(f,view="rotate",neval=100

[R] help with Gantt chart

2010-02-27 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, I don't know to solve this error that is returned, even though I understand it: library(plotrix) Ymd.format<-"%Y/%m/%d" gantt.info<-list(labels= c("First task","Second task (1st part)","Third task (1st part)","Second task (2nd part)","Third task (2nd part)", "Fourt task","Fifth t

[R] how to do calculations in data matrices?

2010-02-13 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Please give me just a reference where I can find something useful. In summary, I need to : - find the median of each row of a matrix - create a new matrix with each value in the first matrix divided by the median of its row - if a value "a" in the second matrix is < 1, I need to substitute it w

Re: [R] logical operations with lists

2010-02-12 Thread Zoppoli, Gabriele (NIH/NCI) [G]
, LMP, NCI, NIH, Bethesda MD Work: 301-451-8575 Mobile: 301-204-5642 Email: zoppo...@mail.nih.gov From: jbreic...@gmail.com [jbreic...@gmail.com] On Behalf Of Jonathan [jonsle...@gmail.com] Sent: Friday, February 12, 2010 5:21 PM To: Zoppoli, Gabriele (NI

[R] logical operations with lists

2010-02-12 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Sorry, maybe it's easy but I haven't found anything useful: how can I obtain a list C that contains all the members in the list B that are not in list A? This are lists of nanes, not numbers! Thank you Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, Univer

[R] permutation analysis with randomly chosen subsets of a matrix

2010-01-31 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hello, here is the problem: I want to demonstrate that, on average, the Pearson's correlations of a specified subset of genes from a huge list (>18,000 columns) are higher than any randomly chosen subset of that list. I would therefore like to do a number of tests between that specified subset

[R] help on Venn's diagram

2010-01-30 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hello! I have this problem: I want to create a Venn's diagram with three lists of genes'names. The first is all the genome, the second a subset of it comprising all mitochondrial genes, and the third including all genes that correlate with a given gene. This is what I do: > library(gplots) >

[R] question about transpose

2010-01-29 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi all, if I transpose a matrix with t(data), the newly created colums do not appear to have the first row as header. How can I do to have all the newly created columns have their first row as a header? Thanks Gabriele Zoppoli, MD Ph.D. Fellow, Experimental and Clinical Oncology and Hematolog

[R] help on drawing right colors within a grouped xyplot (Lattice)

2010-01-29 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, I've lost my mind on it... I have to scatterplot two vectors, grouped by a third variable, with two different dimensions according to whether each cell line in the plot is sensitive or resistant to a given drug, and with a different color for each of 9 tissues of origin. Here's what I've

[R] how to create a histogram from a matrix?

2010-01-28 Thread Zoppoli, Gabriele (NIH/NCI) [G]
Hi, I'm trying to create a frequency histogram for all the values in a table, but when I try to do so, an error is returned, saying that I cannot create a histogram with an obkect of that class. Here's what I do: >library(lattice) >table<-read.table("C:/.../table",header=TRUE,sep="\t") >histog