search for: bioinformatics_guy

Displaying 8 results from an estimated 8 matches for "bioinformatics_guy".

2008 Sep 11
4
plotCI -- multiple plots on same graph
I have a bunch of lines I want to plot using plotCI() What Id like to know is, how can I connect the points with a line and how can I print multiple lines on the same graph? -- View this message in context: http://www.nabble.com/plotCI----multiple-plots-on-same-graph-tp19435198p19435198.html Sent from the R help mailing list archive at Nabble.com.
2009 Sep 17
1
heatmap.2() problems with re-ordering of rows and columns
I have a file of the following form -11 -10 -9 -8 -10 -9 -8 NA -9 -7 NA NA -8 NA NA NA So basically a NxN matrix of log scores. I want to get a heatmap of these log scores but I'm having a problem. I'm using the following code library(gplots) data=read.table("filein.txt",header=FALSE) mat=as.matrix(data) heatmap.2(mat,dendrogram=c("none")) But on the
2008 Sep 30
3
Pattern match in R
I want to make sure this piece of code I wrote is doing what I want it to do. ll<-function(string) { grep(string,dir(),value=T) } subdir = ll("Coverage_[1-9][0-9]$") I basically wrote a little function that would grab all the files of form Coverage_[0-99] The way I wrote it, will it grab Coverage_5 or does it have to have 2 numbers (10-99)? -- View this message in context:
2009 Feb 23
2
if statement of a vector
right now I have a vector of about 1000 points. I'd like to iterate through each of these points and and test if it is greater than a certain value and if not, throw it out. x=vector y=empty vector j=0 for i (0..length[vector]) if x[i] > 10 y[j] = x[i] j++ Thats basically some pseudo code of what I want to do. Any suggestions? -- View this message
2009 Sep 21
1
Specific criteria for color palette using heatmap.2
I'm trying to display the probability space of a function and wanted to see specifically where the maximum posterior probability is. data=read.table("PosteriorData",header=F) mat=as.matrix(data) heatmap.2(mat/max(mat,na.rm=T),dendrogram=c("none"),trace=c("none"), Rowv=F,Colv=F,labRow=3*c(10:-10),labCol=3*c(-10:10),symm=T,col=rainbow(100)) As the
2008 Sep 15
1
Scripting in R -- pattern matching, logic, system calls, the works!
Im very new to R so this might be a very simple question. First I'll lay out the hierarchy of my directories, goals. I have say 5 directories of form "Coverage_(some number)" and each one of these I have text files of form "Length_(some number)" which are comprised of say 30 numbers. Each one of these Length files (which are basically incremented by 5 from 0 to 100,
2008 Oct 01
0
Linear Regression in R
I've attached below the output of a program I am running. What I want to do is, get the equation for each of the lines and determine what the overall equation is that is applicable to every scenario (I know that each line is of multiple of the proceeding line; as in the bottom's coefficient is 1x, next up is 2x, 3x, 4x and so on such that the only difference from line to line is the slope
2009 Jun 22
0
QQplots of probability vector data
I'm trying to determine if a set of data is normal from a qq plot but seem to be having a bit of difficulty. I have a file of the following form 9 36 3 37 6 38 7 39 ..... where the left column is the frequency of the number in the right column. I've found the probabilities of each number and put it in a file of the form 36 .0009 37 .0003 38