similar to: (no subject)

I have a need to find the order of number to get a scaling parameter as a power of 10. I have a function that works *so far*, but it is ugly and probably buggy. In the interest of avoiding code-based outliers in my data, I thought I would ask if anyone here has a better way. > scl <- function(x){ + length(charToRaw(format(trunc(x), scientific = F)))-1} > a <- 123456789 > b <-

I was trying to use memory.size() to determine whether a code loop I am executing created a memory leak, since one replicate of the simulation takes 670.98 seconds according to proc.time(), while 5 replicates takes 170762 seconds. So I set it up as: memA <- memory.size() looping code... memB <- memory.size() memA returns as 9.3, and memB returns 11.3. I'm not familiar with

Dear R Users, In the expresso function, which combination of these methods for data pre-processing (when using affymetrix oligo arrays) is the best: bgcorrect.metod = rma rma2 mas normalize.method = qspline quantiles loess pmcorrect.method = pmonly subtractmm mas summary.method = liwong avgdiff medianpolish mas There are many options within each method. I would appreciate a hint on the best

Tinn-R (http://www.sciviews.org/Tinn-R/) is one of the topmost suggestions when googling an R-(text-)editor for Windows. However, to me it appears dissappointing that Tinn-R does not handle utf-8 (mac-roman, or any other) encoded R-scripts or, in general, text files. Besides Emacs and the R built-in editor, could you recommend a good editor for Windows, even some commmercial for a small

Dear All, How to cite the PDF user's guide for the LIMMA package? This is not about how to cite the LIMMA package. Roger Roger L. Vallejo, Ph.D. Computational Biologist & Geneticist U.S. Department of Agriculture, ARS National Center for Cool & Cold Water Aquaculture 11861 Leetown Road Kearneysville, WV 25430 Voice: (304) 724-8340 Ext. 2141 Email: roger.vallejo@ars.usda.gov

Hi Martin, In using the Cluster Package, I have results for PAM and DIANA clustering algorithms (below "part" and "hier" objects): part <- pam(trout, bestk) # PAM results hier <- diana(trout) # DIANA results GeneNames <- show(RG$genes) # Gene Names are in this object But

Dear R-Help, I am using the function "topTable" from the LIMMA package. To estimate adjusted P-values there are several options (adjust="fdr" , adjust="BH") as shown below: topTable(fit, number = 10, adjust = "BH", fit$Name) I guess any of these options (fdr, BH, etc.) is using a default of FDR=0.05 which is quite conservative (i.e., very

Dear R users, I am using the beadarray package. I am trying to upload raw bead-level data using these commands: ######################################################## library(beadarray) datadir <- ("C:/Computer_programs/R/beadarray/cecilia") targets = read.table("targets.txt", sep = "\t", header = TRUE, as.is = TRUE) BLData = readIllumina(arrayNames =NULL,

In using the NLME package (R 2.6.1 for Windows), I am having a problem in running an R script that used to run with no problems using a Linux OS in 2004. So I am wondering if during these last ~3 yrs we had major changes in the syntax of the NLME package that I am not aware. This is the R script: library(nlme) treat=as.factor(c(1,2,1,2,1,2,1,2)) mouse=as.factor(c(1,1,2,2,3,3,4,4))

Dear Masters, I'm driving crazy with the lowess() function.... my intent is smoothing confidence intervals for predicted y values in a linear model lm() setting since in the predict() function there exists an option for predicting NA values, I instead encounter problems when I fit a missing values x variable to the predicted terms.....,impossible task!!! I've been banging on my head

Is it possible to send data from an executing R script, to external R script files, on linux in this case, get that external R script to execute, and pass results back to the central script? If so what commands should I be looking at? I've googled for R and external, stuff like that but no luck. The reason for the question is that I need to generate simulation data, then perform several

Hi I am new to R and am trying to set up a two-dimensional array/matrix with the elements defined by the function similar to below. Been trying to use outer with apply but can't seem to get the indexing quite right. Is their a simple way of accomplishing this task ?? ----- / | 1 x < 0.5 & y < 0.5 |

Hi, I'm running R 2.11 Does anyone know if it possible to transform one character vector to one character string ? Many thanks Benoit -- Benoit Wastine Laboratoire des Sciences du Climat et de l?Environnement (LSCE/IPSL) CEA-CNRS-UVSQ CE Saclay Orme des merisiers B?t 703 - Pte 13A 91191 Gif sur Yvette Cedex France Tel : 33 (0)1 69 08 21 97 Fax : 33 (0)1 69 08 77 16

Hi Everyone, Can anyone help me resolve a problem that i'm having with nlminb. The problem is that it stops after just one iteration and returns the same values as "start" ones. Thank you very much for your help. Sincerely. -- Kamel Gaanoun (+33) (0)6.76.04.65.77 [[alternative HTML version deleted]]

The guide of R Code style recommends to write script with a maximal length of 80 characters? how to break a line longer than 80 characters like filename_name <- ("/blaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa/"/ blaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa/data.csv") I did not find any answer while googling. Thank you!

Dear R users Is there a way to obtain the residuals from a model fitted by quantile regression? Thank you. Thanaset -- View this message in context: http://r.789695.n4.nabble.com/Quantile-Regression-Extracting-Residuals-tp3225423p3225423.html Sent from the R help mailing list archive at Nabble.com.

Dear colleagues, I am trying to get someone to use R on MS Windows with the browser disabled. My question is how does he/she get to use R-help which goes off the browser (and correspondingly complains about the inability to start firefox, etc). In linux, which is what I use, this is not a problem: what needs to be set in Windows to replicate this behavior? Sorry if I am not clear: please feel

R folks, I am trying to create a mhtplot plot that will alter the point size via the plot point value. The size of the point should change via a set criteria such as: if point value <10^-6 size 3x else if 10^-3> point value >10^-6 size 2x else point value >10^-3 size is x I consider myself an R newbie, and I am assuming that I could do this via some sort of function to change

Which one is more efficient? x2=c() for (i in 1:length(x)) { x2=c(x2,func(x[i])) } or x2=x for (i in 1:length(x)) { x2=func(x[i]) } where func is any function? Dirk -- View this message in context: http://r.789695.n4.nabble.com/appending-to-a-vector-tp3449109p3449109.html Sent from the R help mailing list archive at Nabble.com.

HI, Dear R community, I try to create 100 dummy variables like the following: ack$id_1 <- (ack$ID==1)*1 ack$id_2 <- (ack$ID==2)*1 .. . ack$id_100 <- (ack$ID==100)*1 I used the following codes: for(i in 1:100){ ack$id_[i] <- (ack$ID==i)*1 } But only one column is created, can anyone help me? Thanks a lot! -- Sincerely, Changbin --