similar to: print buglet with small options(width = ) (PR#759)

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

Just thought I would mention that there is a new O'Reilly book out, "Developing Bioinformatics Computer Skills", which is ok, pretty superficial about some things but at least it lets you know what is happening and where. R gets a few pages in there, 394-396, mostly nice press and accurate except that Bill Venables (hi Bill) may be surprised to find out that he is a member of the

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

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

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,

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!

Hello, I don't see how to add a method to a class. I can add "slots" (and retreive them by : myObject at slot1) but not methods, for example accessor for these slots and others functions. I've read the "methods" library documentation but I can't figure out how to do this. Thanks for any help/tips Vincent

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?

Hello, Could someone please tell me what I am thinking about incorrectly: f <- function(y) { g <- function(x) x + y g } In the following, I get what I expect based on my understanding of lexical scoping: (f(1))(3) # 4 (f(2))(3) # 5 But now, fs <- lapply(c(1, 2), f) fs[[1]](3) # 5 (Why not 4 ?) fs[[2]](3) # 5 Checking the environments of these functions, I see that

"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

Dear all, I am trying to write a package using formal methods and classes from the methods package. I have not been able to get the package to pass rcmd check without warnings, because rcmd check does not recognize my generic functions as code objects and therefore queries why they have documentation entries. I have isolated the problem in a very small trivial example which I give below. I

I have the following problem. I have an automatically generated named list with "stringified" names: a <- list("A"=..., "B"=..., "C"=..., ) then I want to refer to the elements of the list, stored as an vector of names: nn <- c("A", "B", "C"), so that I could get list elements like a$nn[1], a$nn[2], etc. Obviously it

I upgraded R-base 1.2.0 today, and then I intended to give a finishing touch to some R functions I should send out today. To my surprise, they failed in R 1.2 although they still worked fine in R 1.0.1 (I checked in another machine). The error message I got was like this: > optci(trar.g, water) Error in unique(c("AsIs", class(x))) : Object "x" not found The function

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

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

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

Hi everyone, There is no mention in ?new.env (R-1.6.0) of what the effect of setting the hash argument of new.env() actually does. What does it mean in performance terms to say that "the environment will be hashed"? Thanks, Jonathan. -- Jonathan Rougier Science Laboratories Department of Mathematical Sciences South Road University of Durham

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