similar to: comparing poisson distributions

Hi, I need to convert numbers to characters in order to index an array when I encountered the number 100000 which turned into "1e+05". Does anyone know how to get around this problem? Lana Schaffer Biostatistics/Informatics The Scripps Research Institute DNA Array Core Facility La Jolla, CA 92037 (858) 784-2263 (858) 784-2994 schaffer at scripps.edu

Hi, I am using this subset statement and it works outside a function. LIS[[i]]<- lapply(LI, subset, select=cov[[i]]) However, wrapped inside a function this statement produces the same values for every LIS[[1]] which is only the first subset of LI. Does anyone know why is not working correctly inside a function? ff = factor(covariate) nLev <- nlevels(ff) cov <-

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

Hi, I am trying to compile the R-dev version on a unix Suse machine and got errors. Would someone be able to help me determine what to do to fix these errors: make[1]: Entering directory `/lustre/people/schaffer/R-devel/m4' make[1]: Nothing to be done for `R'. make[1]: Leaving directory `/lustre/people/schaffer/R-devel/m4' make[1]: Entering directory

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

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

I am having difficulty getting Netatalk to authenticate using NT Domain via Winbind. I have Samba successfully installed, and winbindd is working and running. (Tested via wbinfo) I have placed pam_winbind.so in the /usr/lib/security directory, and have installed netatalk with the --with-pam option, but I still cannot seem to get it working. Any help would be appreciated. John Helt Systems

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

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

Hi, I would like to view matrices I am working with in a clean, easy to read, separate window. A friend showed me how to do something like I want with edit(). I can view the matrix in the 'R Data Editor': For a sample matrix: > mat=matrix(1:15,ncol=3) > mat [,1] [,2] [,3] [1,] 1 6 11 [2,] 2 7 12 [3,] 3 8 13 [4,] 4 9 14 [5,] 5 10 15

All, I was having trouble trying to create a new class of data and pass it on to splom (in the lattice library). I mentioned this to Martin Morgan after a talk he gave. Following is not so much a question, but rather an answer from Morgan that might be useful to others. Here is the edited part of an email conversation with him: On Thursday, August 26, 2010 1:36 PM, Martin Morgan wrote:

OK. This is becoming very frustrating!! I have Samba 2.2.3a installed on a Sun running Solaris 8. Winbind is installed and seems to be working fine. My PC users can map the Samba share via their NT domain account information (PDC is an NT4 box). As soon as I copy nsswitch/libnss_winbind.so to /lib and create the symbolic links (nss_winbind.so.1 & nss_winbind.so.2) my PC users can no

Have made a stacked area plot, but now want to manipulate the x axis: make in even increments (50) in order to suppress the tick marks (forming solid bar under plot). However, the plot functions do not seem to work, and I cannot find documentation for use of xat in stackpoly. This is a time series of data covering 579 years, and has been successfully converted from zoo to matrix. Data plotting

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

Hi. I'm trying to get 'mu' to show up as a Greek symbol but, despite trying every example I could find, can't get it to work. Any insights would be welcome. This is what I'm using that works, but displays mu with the letter u. plotTimeXMastPAR <- qplot(DT,MastPAR, data=A, xlab = "", ylab = quote(PAR (uE ~m^-2 ~s^-1)), geom="line") +

Thanks again everyone, I have a problem trying to install RMySQL on Suse 10.0. I have used both the install.packages() and the 'R CMD INSTALL RMySQL..." neither worked. I get the error message telling me that it cannot find mysql. So I have set the flags using export PKG_CPPFLAGS="-I</usr/bin/>" export PKG_LIBS="-L</usr/lib/mysql>" After that I run R

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

Hi, I need to add on a plot text containing italic Greek characters using the function mtext (i.e. I cannot use Hershey vectors). The characters are nicely displayed when the file is saved as png but not when saved as eps. See code below as example: #postscript('test.eps') png('test.png') plot(1,1, type = 'n') mtext(side = 3, line = 2,