similar to: Possible memory leak in loop.

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 wondering if there is a way to get read in a single keystroke at a time in R as a string, akin to ncurses-style interfaces. I looked into readLines, readChar, etc. using stdin, but these all require the use of an end of line. Has anyone ever had need to do this or have any ideas on how to do this? Thanks, Jon PS I apologize if this double-sends, but I am having mail client issues.

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, I am a stata user and I am moving my first steps in R. I am dealing with a silly issue concerning looping over variables. I read previous posts on similar topics in the R help archive, but I did not find a solution. Here is my case: I run a simple bivariate linear regression saving the output in objects called: pd.memb.0, pd.memb.1, pd.memb.2, pd.memb.3 pd.memb.0 <- lm(E.fs.memb ~

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

What fun. I don't often stoke up the boiler and actually compile open-source code, but here I am at the beginning of just such a joyous and wondrous adventure. I'm running the latest level of AIX (V5.3...V4.3 has been obsolete for at least, what, 5 years) and using the IBM compiler (xlc is similarly obsolete, from back when the compiler came with the O/S...oh, for the halcyon days of

Dear all, I am trying to optimize a logistic function using optim, inside the following functions: #Estimating a and b from thetas and outcomes by ML IRT.estimate.abFromThetaX <- function(t, X, inits, lw=c(-Inf,-Inf), up=rep(Inf,2)){ optRes <- optim(inits, method="L-BFGS-B", fn=IRT.llZetaLambdaCorrNan, gr=IRT.gradZL, lower=lw, upper=up, t=t, X=X)

> Can you explain why you were unable to add a small unit test to the testsuite? This is the case I have. ****** r3 = memb(r1+#1) r4 = memb(r1+ #0) .. .. { p0 = r3 /* Should almost never emit this. */ p1 = r4 /* Should almost never emit this. */ .. .. p0 = or(p1, p0) ***** The CHECK tags I have are ; CHECK: [[IntReg0:r[0-9]+]]

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