similar to: plot.table() ?

Dear R core team Using the following code produces an empty plot (similar to col = NA): > plot(1:9, col = factor(rep(1:3,3), labels = c("red", "blue", "black"))) My question: Shouldn't one get at least a warning (or an error) if one tries to use a factor as col argument? Thanks for an answer. Regards, Christoph Buser

I downloaded the latest version of R (1.5.0) and issued a read.table() command that used to work with (1.3.1). Unfortunately I erased 1.3.1 from my hard drive after I downloaded R 1.5.0. Here's the command and the error message I get under 1.5.0: > golub1<-read.table("C:/Microarrays/code/data_set_ALL_AML_train.txt",sep="\t",q

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

At a snail's pace I keep on translating an introduction to R into italian; I have reached the section describing the glm() function, in which some example code is presented. The very last line of code, before the beginning of the section on Poisson models is: ldp <- ld50(coef(fmp)); ldl <- ld50(coef(fmp)); c(ldp, ldl) which of course gives results 43.663 and 43.663; the correct code

At a snail's pace I keep on translating an introduction to R into italian; I have reached the section describing the glm() function, in which some example code is presented. The very last line of code, before the beginning of the section on Poisson models is: ldp <- ld50(coef(fmp)); ldl <- ld50(coef(fmp)); c(ldp, ldl) which of course gives results 43.663 and 43.663; the correct code

Hello, I thought some of you might like to know about a GNU project that is complementary to R in some ways, namely gretl http://ricardo.ecn.wfu.edu/gretl (GNU Regression, Econometrics and Time-series Library). gretl (a library with cli and gui clients, the gui using GTK) is designed to be very user-friendly, and suitable for teaching econometrics. It has a fairly wide variety of least-squares

Hi, I am having trouble with the survival library, particualrily the coxph function. the following works coxph(jtree9$cph.call,z,rep(1,dim(z)[1])) Call: coxph(formula = jtree9$cph.call, data = z, weights = rep(1, dim(z)[1])) coef exp(coef) se(coef) z p SM 0.2574 1.294 0.0786 3.274 1.1e-03 Sex -0.1283 0.880 0.1809 -0.709

I'm trying to do a Windows build of Ray Brownrigg's port of the maps library, and running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c

I'm trying to do a Windows build of Ray Brownrigg's port of the maps library, and running into some problems. There's a makefile for the shared library. I suspect that I just need to change that to get it right, and everything will build properly, but I'm running into trouble. Here's the original makefile, along with my comments about changes: >SRC= mapget.c mapproject.c

The point estimates produced by wilcox.test are perverse (not wrong, just brain damaged). The Hodges-Lehmann estimator that goes with the signed rank test is the median of the Walsh averages. The Hodges-Lehmann estimator that goes with the rank sum test is the median of the pairwise differences. wilcox.test agrees except that it uses the following very peculiar definition of "sample

I'm trying to build the 1.4.0 development version of R for Windows. I'll be posting problems as they arise, until I give up (or succeed), or someone asks me to stop. 1. I don't have zlib, the compression library. The source file gnuwin32/unzip/WHERE says to get it from ftp://ftp.info-zip.org/pub/infozip/zlib/zlib.html; that web site says it has been closed and redirects me to

Hi Duncan, I've been looking around at GUI toolkits to use with Python (and by extension R). I initially tried Tk, but was disappointed because there is no consistent set of high level widgets. I've subsequently done some looking at QT, wxWindows, and GTK. QT is out because I am a *commercial* user, and I don't want to have a license hassle. GTK seemed to have a lower level of

Hi, Kurt Hornik and I will start to publish a quarterly electronic ``R Newsletter'' beginning with January 2001 (if all goes well). The idea is to have a mixture of articles describing 1) new features in R itself 2) contributed add-on packages 3) nice applications 4) more general statistical computing issues related to R 5) or ... OK, now comes my christmas wish to you: If you are

I have compiled R version 1.3.1 on a Linux system. No errors or warnings are reported. However, "make check" fails thus: ===== make[5]: Leaving directory `/proj/informatics-build/GML/R/R-1.3.1/src/library' running code in base-Ex.R ... ../../bin/R --vanilla < base-Ex.R > base-Ex.Rout make[4]: *** [base-Ex.Rout] Error 139 ===== The base-Ex.Rout being produced at this

Hi, Can someone set me straight as to how to write formulas in R to indicate: complete independence [A][B][C] joint independence [AB][C] conditional independence [AC][BC] nway interaction [AB][AC][BC] ? For example, if I have 4 factors: hair colour, eye colour, age, sex does > mosaicplot( frequency ~ hair + eye + age + sex) mean that the model fitted is of complete independence of

This dataset below is one sample answer the questioner from our customer. > testbank <- read.table("testbank.txt", header=T) > testbank age married income gender ownhouse class 1 20-30 no high female yes 1st 2 30-40 no high female yes 1st 3 40-50 no low female yes 1st 4 50-60 no high female yes 1st 5 60-70

For a variety of displays (mosaicplots, barplots, ...) one often wants to either abbreviate or wrap long labels, particularly when these are made up of several words. In general, it would be nice to have a function, abbreviate.or.wrap <- function(x, maxlength=10, maxlines=2, split=" ") { } that would take a character vector or a list of vectors, x, and try to abbreviate or wrap

I have been tinkering with mosaicplot() and friends as a way of learning R. As part of this, I've written a pair.table() method for mosaic matrices, and would like to extend mosaicplot to work with loglin and logln (MASS) objects. I'm using R 1.4.0 on Win 98. I've been trying to figure out the formula interface, and think there's a bug, but not sure how to find it, yet alone fix

Hello. I've been using R for a couple of months and enjoying it a lot. This is my first post to R-help. I'm using the vcd package to make mosaic plots with labels on the tiles indicating the number of items in each cell. For example, I've made this plot: > allmorph<-structure(c(10, 26, 17, 100, 70, 97, 253, 430, 185, 177, > 25, 1), .Dim = as.integer(c(6, 2)), .Dimnames

Hello Pardon for the elementary question, I did try searching the archives with various terms, but to no avail. I am using R1.9.1 on a windows machine One of the great advantages of R (to me, anyway) is being able to see the code for a function , e.g. by typing sd one sees the code for getting a standard deviation. However, for many functions this only provides info. including UseMethod, eg.