similar to: new book

Dear R-help, I think I have kinda isolated the problem I had to the following: i <- 0 while( {i <- i + 1} < 5) { if(i < 3) next print(i) } This seems to go into an infinite loop. After I break the execution, i has the value 1. At the R prompt, if I start from i <- 0 and keep typing {i <- i + 1} < 5, it eventually evaluate to TRUE. So why does the while loop not work?

Prof. Tierney, Thanks very much for the info. Why does the loop work if I move the assignment out of the condition? E.g., the following works: i <- 0 while(i < 5) { i <- i + 1 if(i < 3) next print(i) } Regards, Andy > -----Original Message----- > From: Luke Tierney [mailto:luke at nokomis.stat.umn.edu] > Sent: Tuesday, October 02, 2001 3:36 PM > To: Liaw, Andy

I have a, no doubt, simple question. I wish to write a function such that a <- 9 b <- 10 changer _ function(x,y) { if (y>x){ x <<- Y+1}} Of course there are easier ways to accomplish the task above, but I am more interested in how to have the "x <<- Y+1" part of the function to change x in place for purposes of a much larger function. I have been wrestling with

Is there an implementation of a reasonable k-nearest neighbor finder already in one of the packages? -- +---------------------------------------------------------------------------+ | Robert Gentleman phone : (617) 632-5250 | | Associate Professor fax: (617) 632-2444 | | Department of Biostatistics office: not yet

Dear all, Does anyone knows how I could extract the p-value in: > survdiff(Surv(tempo,status) ~ grupo,data=dados1,rho=1) Call: survdiff(formula = Surv(tempo, status) ~ grupo, data = dados1, rho = 1) N Observed Expected (O-E)^2/E (O-E)^2/V grupo=1 21 5.12 12.00 3.94 14.5 grupo=2 21 14.55 7.68 6.16 14.5 Chisq= 14.5 on 1 degrees of freedom,

Byron Ellis has been making some progress on a hdf5 library for microarray data (and any other thing you want to put in there). In doing so some issues have arisen that are of more general interest. 1) hdf5 supports annotation (through comments) so it would be nice if the comment function in R became generic. I think this is backward compatible and basically not really an issue of any

> From: Robert Gentleman [mailto:rgentlem@jimmy.harvard.edu] > > Greg, > you might be better off writing a wrapper, > > R->foo where foo has less than 65 args and does nothing more than > unpack them and calls bar (the original) with as many args as you > want, not much code change and not much cost... Hi Rob, [BTW, sorry about the crosspost

It was my fault, as you have seen. Y matrix must contain successes and failures, and not successes and trials. Sorry for my ignorance and thanks a lot for your help. Oscar -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20001220/341c2519/attachment.html

Prof. Gentleman (and R-help), Here's an example of what didn't work. I still don't understand why. Function: trycon <- function(file, n) { f.con <- file(file, open="rt") on.exit(close(f.con)) i <- 0 while( length(readln <- readLines(f.con, 1)) > 0 ) { x <- unlist(strsplit(readln, " ")) if(length(x) <= 6 && x[3] ==

"An Introduction to R" says that arithmetic operators coerce logical vectors to numeric vectors. This doesn't seem to be true for unary "-" however: > x <- 1:10 > p1 <- (x > 5) > p2 <- (-p1) > is.logical(p2) [1] TRUE > is.numeric(p2) [1] FALSE Since "==" always does binary comparison, this can lead to some bizarre behavior: > p1

The standard R boxplot appears to use quartiles to determine the height of the rectangles and a range parameter - RNG - (default=1.5 I think) that determines the length of the whiskers as <= RNG x Interquartile Range. Is it possible to instead define the range as extending to the 95th percentile? If so, how would this be done? nb, I'm plotting multiple boxplots on a single chart so a

Hi, The help for `assign' says In all the assignment operator expressions, `x' can be a name or an expression defining a part of an object to be replaced (e.g., `z[[1]]'). The name does not need to be quoted, though it can be. So I would expect that by a<-1:4 assign("a[1]",2) the first component of a gets the value 2, but the object "a[1]" is

I realize this isn't quite the correct forum (but then I don't think that there is one). Anyways, I've just been given enough money to hire two programmer for at least two years to work on the following project: (for experienced programmers only, please, salary in the 70k range) BRIEF DESCRIPTION FOR JOB POSTING To develop a system for the manipulation and analysis of microarray

Hello all, Since I got no answer to my first mail, I suspect I did not formulate my question adequately, so I try again. In R, one can debug or trace a function 'foo' with debug(foo) or trace(foo), respectively. This leads to a special treatment of the function 'foo' until one enter undebug(foo) or untrace(foo). I would like to know if there is a convenient way to know at any

Robert Gentleman wrote: > On Tue, Oct 16, 2001 at 04:40:34PM +0100, Rita Ribeiro wrote: > > Hi, > > > > In order to avoid deep copies by passing large arguments to functions or > > > > returning values, I'm trying to do the assignment of variables in a > > given environment. The problem is when I try to assign a structure: a > > list for example.

Environment variables set with Sys.putenv() disappear (i.e. become "") after a while, especially after heavy-duty I/O. Example: R> x <- matrix(1., 3000, 3000) R> save(x, file="myx.RData") R> Sys.putenv(HOME="/tmp") R> while (Sys.getenv("HOME") != "") {cat("ok\n"); load("myx.RData")} The loop prints

Full_Name: Felix Hernandez-Campos Version: R-1.5.0 OS: IRIX, FreeBSD, Solaris Submission from: (NULL) (208.30.174.32) [R-1.5.0] The installation of the package cluster fails in FreeBSD (4.2) and IRIX (6.5) (but it works fine for Solaris). In FreeBSD, the installation hangs right after f77 -fPIC -g -O2 -c clara.f -o clara.o (the first FORTRAN compile). In IRIX, the program also hangs at

All: I am getting an error code from the optimization function. The code is Error in optim(p,fun.LLike, lower=low, upper = up, method = "L-BFGS-B", : non-finite finite-difference value [0] If I add a trace=6 option to my control list the last message before this error is: At X0, 0 variables are exactly at the bounds Any ideas on where I should start would be

This is in R 1.1.1 and "R-devel" : > options(width = 10) ; pi [1] [1] 3.141593 Note the extra "[1]". -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject

Hello! I don't know if there is some R-package able to associate ontology to a long list of GeneBank Name (a txt-tab file or an XML file), i.e. I would as output a formatted file with 4 columns (1:GeneBank Name 2,3,4:ontology). I know that I have to perform a mapping of genes, I got a look on AnnBuilder pkg, but I 've not idea from where to start. Some suggestion? Thanks in advance!