[R] as.character on OS X tiger

2005-09-22 Thread Marcus Davy
Hi, When I run this example code on a G5 the last digit is dropped off on a binary installation of R (R.app), whereas on a linux machine the digit is not removed. # Code: a - seq(0,1, length=1000)[12] a as.character(a) as.character(as.numeric(as.character(a))) # Results I get: a - seq(0,1,

RE: [R] How to extract x rows to get x pvalues using t.test

2005-03-16 Thread Marcus Davy
Hi, if a striped down version of t.test is required for speed, before implementing a rewrite I would check out the multtest Bioconductor package. It takes a fraction of a second to do 56000 t-tests on a 2.4 Ghz PIV. dim(X) [1] 56000 4 pcols - 4 replicates - rep(0:1, each=2)

[R] Re: [BioC] limma, FDR, and p.adjust

2004-12-19 Thread Marcus Davy
Mark, there is a fdr website link via Yoav Benjamini's homepage which is: http://www.math.tau.ac.il/%7Eroee/index.htm On it you can download an S-Plus function (under the downloads link) which calculates the false discovery rate threshold alpha level using stepup, stepdown, dependence methods

Re: [R] Farmating printed Numbers and Text

2004-12-13 Thread Marcus Davy
Or ?formatC, e.g. tmp - c(1,10,100,1000,1) formatC(tmp, digits=4, format=d) [1] 110 100 1000 1 sapply(tmp, function(x){sprintf(%5.0f, as.double(x))}) [1] 110 100 1000 1 marcus Andrew Robinson [EMAIL PROTECTED] 14/12/2004 1:51:51 PM Or ?sprintf Andrew

Re: [R] Is k equivalent to k:k ?

2004-12-12 Thread Marcus Davy
Havent seen anything on R-devel yet on this topic, but if the elseif statement on line 43 of seq.default was else if (length.out == 0) double(0) instead of else if (length.out == 0) numeric(0) then that *might* satisfy Richard A. O'Keefe's comment for length=0:n cases

Re: [R] a question about swap space, memory and read.table()

2004-12-09 Thread Marcus Davy
On a 32-bit windows (standard install or R) machine you cannot allocate more than 2Gigs of memory, R --max-mem-size=2G or R --max-mem-size=2000M If you try to allocate more than 2G (eg 4000 Meg) I suspect it then defaults back to 1Gig. Check you memory limit with memory.limit(). There is lots

Re: [R] Is k equivalent to k:k ?

2004-12-08 Thread Marcus Davy
Bringing up an old topic on a small technicality. In the help documentation on seq. Value: The result is of 'mode' 'integer' if 'from' is (numerically equal to an) integer and 'by' is not specified. The arguement in seq specifying length also creates double which is not obvious in

Re: [R] printing tabular data nicely

2004-06-03 Thread Marcus Davy
Hi, If you choose to use LaTeX and xtable, you can also enhance the visual presentation of the output tables using the booktab LaTeX package. Have a look at: http://www.ecs.soton.ac.uk/~srg/softwaretools/document/start/booktabs.pdf You can change the horizontal line tags \hline to \toprule,

Re: [R] Accessing more than 2GB memory in Windows

2004-05-19 Thread Marcus Davy
Hi, I am also interested in the same topic, that is windows 32bit restriction to 4Gigs of virtual address space partitioning into User and System space. For your reconfigured server 3:1 User to System space what does memory.limit() say is available? For those interested a windows FAQ link on

[R] floor of n observations in number generators

2003-12-30 Thread Marcus Davy
)) [1] 99 length(runif(m*(1-pi0),min=0, max=1)) [1] 99 marcus Marcus Davy Bioinformatics __ The contents of this e-mail are privileged and/or confidenti...{{dropped}} __ [EMAIL PROTECTED] mailing