[R] DISPLAYING A HISTOGRAM + finding NAs

2006-04-04 Thread mark salsburg
I have a vector of values (about 8000 obs between 0 and 1). I want to display a histogram as follows: value of vector ( bw 0 and 1)| | | | |

[R] Selecting out certain values from a MATRIX

2006-04-04 Thread mark salsburg
I have two objects, one matrix and one vector. I want to use my vector to subset certain values out of my matrix. For example: I want to tell R, to select out all rows in myMatrix into a new myMatrix2 IF that corresponding row is less than a 0.5 in myVector. So: myMatrix = a matrix of 8000 by

Re: [R] Selecting out certain values from a MATRIX

2006-04-04 Thread mark salsburg
like: myMatrix2 - myMatrix[ myMatrix[,1] myVector, ] You may find the subscripting section of chapter 1 of S Poetry to be useful. Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) mark salsburg wrote: I

[R] VARIANCE across each ROW

2006-04-01 Thread mark salsburg
I have a very large matrix. I would like to display the variance across each row. In other words, I want to output a vector containing the values of variance across row. When I use the function var(), it seems to give me the variability of each column. Any ideas?? [[alternative HTML

[R] ROWNAMES error message

2006-03-21 Thread mark salsburg
I am getting an error message, which I do not know the source to. I have a matrix SAMPLES that has preexisting rownames that I would like to change. GENE_NAMES contains these rownames. rownames(SAMPLES) = GENE_NAMES Error in dimnames-.data.frame(`*tmp*`, value = list(list(V1 = c(3843, :

[R] removing ROWS with missing values

2006-03-16 Thread mark salsburg
I am trying to find out if R can recognize specific criteria for removing rows (i.e. a prexisting function) I have a matrix myMatrix that is 12000 by 20 I would like to remove rows from myMatrix that have: -999 across all columns -999 across all columns but one -999 across all columns but two

[R] generating RANDOM ROWS from matrix

2006-03-15 Thread mark salsburg
Dear group, I would like to generate a 1000 random rows from a MATRIX with dimensions 12,000 by 20 (i.e. to generate a 1000 by 20 matrix of random rows) Does the function sample() work for this??? thank you in advance [[alternative HTML version deleted]]

[R] MEMORY limit in R

2006-03-13 Thread mark salsburg
I sent this to Bioconductor, but I think it is more of an R listserv question, because of memory allocation trouble I am having trouble getting the function heatmap() to work on the following gene expression dim(SAMPLES_log) [1] 1262620 sample1

[R] Error with write.table

2006-03-09 Thread mark salsburg
I'm having this error when I try to output out a gene expression data frame Error in file(file, ifelse(append, a, w)) : unable to open connection In addition: Warning message: cannot open file 'SAMPLES.txt' when issuing the command: write.table(SAMPLES, file = SAMPLES.txt)

[R] reading in only one column from text file

2006-03-07 Thread mark salsburg
How do I manipulate the read.table function to read in only the 2nd column??? [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] Removing Rows

2005-11-21 Thread mark salsburg
I have a data frame with the following dimensions 217 x 5 I want to create two data frames from the original. 1) One containing every tenth row of the original data frame 2) Other containing the rest of the rows. How do I do this? I've tried subset() and calling the index. thank you in

[R] Classifying values in vector

2005-08-17 Thread mark salsburg
I have a vector of size 217 called A. the values of A are not sorted and range from 0 to 1 (normalized) I am having difficulty writing a program to create a new vector B where if A's value is 0 A =0.333 then B is 0 if A's value is 0.333 A =0.666 then B is 1 if A's value is 0.666 A =1 then B is

[R] Concurrence Matrix

2005-08-17 Thread mark salsburg
Thank you for all the help. Can someone refer me to a function that can help with the creation of a concurrence matrix? I have two asymmetric vectors (C is length 217 and D is length 16063) that contain values of 0,1,2 I want to create a matrix E (where rows are D and columns are C) that

[R] tapply t.test

2005-07-26 Thread mark salsburg
I cannot find in the literature a way to conduct the following t.test on 2 objects, A and B A B col1 col2 col3 col1 col2 col3 Where col(i)'s name is identical in both A and B (they are names of tissues). How do I test (t.test) if each

[R] T test across tissues

2005-07-26 Thread mark salsburg
I cannot find in the literature a way to conduct the following t.test on 2 objects, A and B A B col1 col2 col3 col1 col2 col3 Where col(i)'s name is identical in both A and B (they are names of tissues). How do I test (t.test) if each

[R] Help on T test

2005-07-26 Thread mark salsburg
ok I created a matrix C with A B CA1 B1C1 .. . the columns contain the gene expression values.. I ran the following t.test: apply(C, 1, function(x) t.test( x[1:3], x[4,6] )$p.value ) which outputs out 16063 pvalues

Re: [R] predict.lm - standard error of predicted means?

2005-07-20 Thread mark salsburg
Can someone please refer me to a function or method that resolves this structuring issue: I have two matrices with identical colnames (89), but varying number of observations: matrix Amatrix B 217 x 89 16063 x 89 I want to creat one

[R] Combining two matrices

2005-07-20 Thread mark salsburg
Can someone please refer me to a function or method that resolves this structuring issue: I have two matrices with identical colnames (89), but varying number of observations: matrix Amatrix B 217 x 89 16063 x 89 I want to creat one

Re: [R] Combining two matrices

2005-07-20 Thread mark salsburg
The latter, I want to replicate A enough times to have the same rows as B. I realize that 16063 is not a mutliple 217, that doesn't really matter for this problem Assuming I change the colnames of one of them, how can I create a matrix C that is the composite of the replicated A with B. thank you

[R] .gct file

2005-07-19 Thread mark salsburg
I have two files to compare, one is a regular txt file that I can read in no prob. The other is a .gct file (How do I read in this one?) I tried a simple read.table(data.gct, header = T) How do you suggest reading in this file?? thank you. __

Re: [R] .gct file

2005-07-19 Thread mark salsburg
a gct. file thank you On 7/19/05, Duncan Murdoch [EMAIL PROTECTED] wrote: On 7/19/2005 12:10 PM, mark salsburg wrote: I have two files to compare, one is a regular txt file that I can read in no prob. The other is a .gct file (How do I read in this one?) I tried a simple

Re: [R] .gct file

2005-07-19 Thread mark salsburg
to read the data as tab delemited, but read in the 2nd coloumn as one value based on the quotation marks.. thanks.. On 7/19/05, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: On Tue, 2005-07-19 at 13:16 -0400, mark salsburg wrote: ok so the gct file looks like this: #1.2 (version number) 7283

[R] interactive plot showing label

2005-07-12 Thread mark salsburg
I have the following data: gene_name microarray expression A 2323 B 1983 . I have about 10,000 observations. I would like to know if there is a way to