similar to: a matrix in a seprate window?

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]

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

Hello, I hope this doesn't turn into a statistics question but here I go. I am using the nls function with a Gaussian distribution, see coding below. When I run the nls I get an error back saying that I have a linear gradient. I then, of course am unable to do anything else. The data that I am using are intensity values from some mass spectrometry data. Is there something I can

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

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

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 '

HI Madhu, I guess you got your solution from Rui: ?dat1<-data.frame(x=c(1,1,2,2,2,3,4,4,4),y=c(23,45,87,46,78,12,87,79,76)) s<-sample(unique(dat1[,1]),length(unique(dat1[,1]))*0.8) ?s #[1] 3 4 2 You can have a list containing both the dataframes list1<-list(dat1[dat1$x%in%s,],dat1[!dat1$x%in%s,]) list1 [[1]] #? x? y #3 2 87 #4 2 46 #5 2 78 #6 3 12 #7 4 87 #8 4 79 #9 4 76 #[[2]] ?# x? y

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

Dear all, Thank you for any help. I have a data.frame and would like to plot it in 3D. I have tried wireframe() and cloud(), I got scatterplot3d(xs) Error: could not find function "scatterplot3d" > wireframe(xs) Error in wireframe(xs) : no applicable method for "wireframe" > persp(x=x, y=y, z=xs) Error in persp.default(x = x, y = y, z = xs) : (list)

Hello all, I consider myself pretty new to the whole OO based programming so I'm sorry if I'm doing something stupid. > xml<-read.metlin(url) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "read.metlin", for signature "url" >read.metlin standardGeneric for "read.metlin" defined from package

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))*

hi, i've setup a LDAP server with account information, and compiled samba with ldap support. everything works great, except for the password changes i still have to run two seprate commands ( passwd, smbpasswd ) to change a users password. i've tried to put the pam_smbpasswd.so module into system-auth, but that does work? any pointers? thanks adriaan putter

Hello list, I posted this a few days ago on opensolaris-discuss@ list I am posting here, because there my be too much noise on other lists I have been without this zfs set for a week now. My main concern at this point,is it even possible to recover this zpool. How does the metadata work? what tool could is use to rebuild the corrupted parts or even find out what parts are corrupted. most but

Hi. I tried the following: F<-function(x,y,b1,b2,b3) b1/(1+b2*(x+b3*y)); t<-data.frame(ExM=c(T$ExM,F(T$S1,T$S2,0.2,0.2,51.28205)),S1=rep(T$S1,2),S2=rep(T$S2,2),gr=c(rep(1,length(T$ExM)),rep(2,length(T$ExM)))); library(lattice); wireframe(ExM ~ S1 * S2, data = t, groups=gr); where T is da dataframe like T<-data.frame(S1=...,S2=...,ExM=...) But the resulting response

Good morning all, I have been slowly tweaking the wiki document ref my XPS M1530 and think it's at a stage where I can remove the Draft status. However, before I do this I'd like opinions just in case I have missed something. Anyone care to add anything? On a seprate note I have been unable to get Bluetooth working under Linux so I will try and get more info as I can. At present I just

I am having trouble locating webvmail on my * server. Is this a seprate porgram or does it come with *. I am running version asterick*CLI> show version Asterisk CVS-03/26/04-17:08:20 built by root@localhost.localdomain on a i686 running Linux asterick*CLI> Thanks Kurt __________________________________ Do you Yahoo!? Yahoo! Photos: High-quality 4x6 digital prints for 25¢

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