search for: benton

Hello, I'm doing some scripting and I've noticed that R doesn't seem to have an if (cond){ do }ifelse (cond) { do } else { do } type block. Is this correct or am I missing something. THX Paul -- Research Technician Mass Spectrometry o The / o Scripps \ o Research / o Institute

Hello, Sorry to all for the lack of communication, and thanks to those with the suggestions. So I have two matrixes which are different sizes, same column number but different row number. What I would like to do is to compare A to B. In the 6 and 9th column there are ref numbers (molecular masses to be exact). I would like to check A and B by these columns. If I find that the number in A[9,x]

Thanks Tom. heatmap.2(mat, Rowv=NA, Colv=NA, col=redgreen, breaks=seq(-1,1,by=0.01), symm=TRUE, trace="none", dendrogram="none") that did the trick :) Begin forwarded message: From: "Leja, Thomas" Date: 29 March 2011 03:47:59 GMT+01:00 To: "Benton, Paul" Subject: Re: [R] comparing heatmaps Yep, you have to specify the range for colours. I usually plot the correlations from 0 to 1. Just add "breaks" option and colour without the gradient: col=redgreen, breaks=seq(0,1,by=0.1), Hope this helps :) Best, Thomas From: "Be...

...linear gradient, and I dare to ask why am I getting a linear gradient. I apologies if I haven't supplied enough information or I have made some mistake in my coding. The coding below works on some of the data but not all. Thank you for your time, Paul Benton brseq <- seq(-4.0, 4.0, by=0.1) AB<- A[,1]/A[,2] lgAB<-log(AB) freq_AB <-hist(lgAB, type="o", breaks=brseq, plot=F freq.tab <- as.data.frame(cbind(brseq, freq_AB$counts) class(fo <- (x ~ (A/(sig*sqrt(2*pi)))* exp(-1*((bin-mu)^2/(2* sig^2))))) nls.AB <- nls(f...

Umm sorry to bother everyone again but I'm having trouble with my if statement. I come from a perl background so that's probably my problem! :) So here is my code: if (any(lgAB>4) | any(lgAB<-4)){ freq_AB<-hist(lgAB, type="o", plot=F) else freq_AB<-hist(lgAB, breaks=br,type ="o", plot=F) } And I get > source("E:/R/GMDA-1.1.R") Error in

Dear all, I've been trying to find how to compare tow different heatmaps but I'm having trouble getting the colors bar to be the same. I'm doing something like the following: library(gplots) dat<-cor(matrix(rnorm(100, m=10), nrow=10)) mat<-cor(matrix(rnorm(100), nrow=10)) dev.new() heatmap.2(mat, Rowv=NA, Colv=NA, col=redgreen(75), symm=TRUE, trace="none",

Posting for a friend Begin forwarded message: From: "Geier, Florian" <florian.geier08@imperial.ac.uk<mailto:florian.geier08@imperial.ac.uk>> Subject: Fwd: readBin fails to read large files Date: September 1, 2011 4:10:53 PM GMT+01:00 To: Begin forwarded message: Date: 1 September 2011 16:01:45 GMT+01:00 Subject: readBin fails to read large files Dear all, I am trying

Hello all, I know this is a pretty easy question but I can't find it in S poetry or R help. How can I make a negative number positive. Such as -5 to be +5 I tried +(-5), but that didn't work. So no, I don't mean taking a -5^2 just to get a positive number. This is in a function so it's not just -5 it's x. :) Thanks, Paul -- Research Technician Mass Spectrometry o The

I'm having lots of trouble getting connected to samba on my little 2 computer network at home. One is a RedHat machine and the other is Windows 2000. Trying to access my just the IP of the RedHat machine from my Windows machines gives me the "Cannot find the file or item '\\192.168.20.2'" error. I've got 2 NICs in my RedHat machine so I did use the "interface

bryce benton wrote: > number_to_currency(1234567890.506) => $1,234,567,890.51 BTW agile book does mention this .. will this number_to_currency(x) be internationalized ? or i can set it to Indonesian for example (or to any other languages rather than just English) ? -- Salam hangat, <%= Arie %>...

Hello all, I'm looking for a way to be able to read a text file into R. It's a csv file but when I do "txt <-read.table("F00.csv", header=T, sep=",")" It doesn't read the file properly, and I only get 2 columns. If I open it up in OOc or Excel it open right with 7 columns. What I would really like to do is read the file as text and then split it and

Hello, I apologise if someone has already answered this but I searched and googled but didn't find anything. I have a matrix which gives me the similarity of each item to each other. I would like to turn this matrix into something like what they have in the graph package with the nodes and edges. http://cran.r-project.org/doc/packages/graph.pdf . However I cannot find a method to

...1. another diablo II prob (shemm@gmx.de) > 2. Re: VXDs under linux? (Chris) > 3. Is there a list of games/applications working with Wine? (Christian Hedin) > 4. Re: VXDs under linux? (Zac) > 5. Drive Y (Rusmir =?ISO-8859-2?Q?Du=B9ko?=) > 6. wine.conf newbie problem (Benton) > 7. Re: Drive Y (g.patel@wanadoo.fr) > 8. Re: Is there a list of games/applications working with Wine? (Benjamin FRANCOIS) > 9. Re: wine.conf newbie problem (g.patel@wanadoo.fr) > 10. Re: VXDs under linux? (Daniel Foesch) > 11. Re: Not enough system resources to create...

Dear all, Just a quick one, hopefully. I have a histogram made from the method 'hist()'. How do I get labels on the bars? Such that the bars will have the x axis on the bar, not the frequency of the point but the number of the point itself. To make a quick summary, I want the the numbers from 'losses' (below) to be on the bar's. Thanks, Paul

So thanks for the help, I have a matrix (AB) which in the first column has my bin numbers so -4 - +4 in 0.1 bin units. Then I have in the second column the frequency from some data. I have plotted them and they look roughly Gaussian. So I want to fit them/ find/optimize "mu", "sigma", and "A". So I call the nls function : nls_AB <- nls(x ~ (A/sig*sqrt(2*pi))*

After noticing a commit to the HG repository, http://hg.dovecot.org/dovecot-1.0/rev/f88b1e225be5 about, " We didn't update last_output timestamp with long fetches, causing Dovecot to idle-disconnect the client." I chekced out the latest from the DOVECOT-1.0 repo, and built again. Unfortunately, this problem still exists. Login simply disconnects. Any help, anyone?

hello all, Quick question, how do I get the p value out of the anova? Thanks, Paul > pb<-aov(as.numeric(diff[5,16:33]) ~ grF) > summary(pb) Df Sum Sq Mean Sq F value Pr(>F) grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 * Residuals 14 3.0783e+10 2.1988e+09 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 '

Hello all, I'm having a spot of trouble. I have made a very nice distance matrix from some comparison calculations. However, the matrix is 3096 x 3096, so it's relatively large. Currently I'm using the svn package as it seems to work for my data. The graph package errors and says the "adjacency matrix must be symmetric for undirected graphs". I though 3096 x 3096 was

Dear All, I have a matrix with data that is not organised. I would like to go through this and extract it. Each feature has 2 vectors which express the data. I also have an index of the places where the data should be cut. eg. >class(cc) "matrix" >cc [,1] [,2] [1,] 1 26 [2,] 2 27 [3,] 3 28 [4,] 4 29 [5,] 5 30 [6,] 6 31 [7,] 7 32

Hello all and thank you, I will try to make this as clear as I can. I have a matrix, at the top of the matrix are 4 time points, the first col describes the sample, then under each time point I have the intensity values of that sample. So : Samples 0 10 30 120 <- time A 200 300 400 1000 <- data B 100 400 200 100 C 500 400 200 1000 ^ | samples I would like to plot these however, all I