similar to: mathematics in R graphics

It was a real surprise, but a student in my class found that the function eigen is buggy. He traced to the problem from his inability of getting principal component analysis to work on his data. Chong Gu Here is a matrix I generated through X'X, where X is 2x3. > jj [,1] [,2] [,3] [1,] 0.8288469 -1.269783 -0.7533517 [2,] -1.2697829 2.162132 2.0262917 [3,]

I am using R to teach, and here are a couple of things that I thought would work didn't work. 1. I noticed the utility data(***,package=***) recently and like it very much, but unless I type in the whole word "package" I'll get an error in 1.4.0. For example, data(cats,package=MASS) works fine but data(cats,pac=MASS) doesn't. 2. drop1 doesn't seem to be as smart as

I am not sure if I should call this a bug or just a design imperfection. I have noticed that the data.frame function does not preserve the integrity of multivariate components. Here is a simple illustration: > x <- matrix(1:6,3,2) > y <- 1:3 > z <- data.frame(x=x,y=y) > z$x # NULL > z$y # [1] 1 2 3 One can however get by that by using model.frame: > zz <-

survexp in survival5 doesn't seem to work for me. see below: > library(survival5) Attaching Package "package:survival5": The following object(s) are masked from package:base : sort.list > library(chron) > data(ratetables) > survexp(~ratetable(year=julian(6,1,1991), + sex=1,age=35*365.24),times=(0:30)/6*365.24) Error in as.character(as.date(c(min(R[, 3]),

I was trying to add points to filled.contour plots, but the points(...) function apparently uses a different coordinate system than the filled.contour(...) function; it appears that points(...) paints points where they would be if the contours were stretched over the whole plotting frame including the legend bar. Am I missing something or is it a "misalignment" in the code? BTW, I have

I'd like to create a data frame with components > jk$x1 [1] 2 > jk$x2 [,1] [,2] [1,] 0 0 I used to be able to do it with > jk <- data.frame(x1=2,x2=I(matrix(0,1,2))) But now I get a error message. Can I still do what I want? Thanks for any help. Chong Gu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

I needed to generate some data from the Inverse Gaussian distribution but it is not in R. A google search took me to http://www.maths.uq.edu.au/~gks/s/invgauss.html which contains the documentation for the d, p, r, q functions with a download link of the functions. The functions work fine in R. I am curious why the facilities are not included in R. Is it for legal reasons or just the

Can someone confirm that the EISPACK routines for eigenvalues of symmetric matrix are in base R. They seem to be, but I can't seem to locate where they are in the src tree. Thanks. Chong Gu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi Once more I have doubts about references. I want to give a reference of the packages I've used (gss, MASS, geoR). For the gss package I'm using @TechReport{gss00, Author = {Chong Gu}, Title = {The gss package}, Institution = {CRAN}, year = 2000, month = {22 } # dec, } and something similar for the others, but this doesn't look

Two questions: 1. inverse.gaussian is up there as one of the glm families, but do people ever use it? There is no inverse.gaussian in the R distribution family, and when I checked McCullagh & Nelder, it only appeared twice in the book (according to subject index), once in the table on p. 30 and once on p. 38 in a passing sentence. Is there a good reference on this distribution? 2. When I

Folks, I am not sure if it's a feature or a "bug". The same is observed in Splus. Suppose I have Poisson counts, and I would like to estimate the parameter using glm. I would assume I can feed it the individual counts, or I can feed it the distinctive counts with the frequency as the weights, and I would get the same results. I do, but the deviance df are returned differently.

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

I am not sure whether it is a design decision or just an oversight. When I ask for the standard errors of the predictions with predict(budwm.lgt,se=TRUE) where budwm.lgt is a logistic fit of the budworm data in MASS, I got Error in match.arg(type) : ARG should be one of response, terms If one is to construct a CI for the fitted binomial probability, wouldn't it be more natural to do

I was hoping that the following specification postscript(file="fig.ps",height=2,width=4) par(mar=c(5,5,2,2)+.1,mfrow=c(1,2),mex=.6) would yield the same scaling effect as postscript(file="fig.ps",height=2,width=6) par(mar=c(5,5,2,2)+.1,mfrow=c(1,3),mex=.6) but the fact is that the 'mfrow=c(1,2)' yields smaller plotting frame and bigger font size.

Dear R-users, I am working to install Chong Gu's new gss library on my Linux system (Redhat 5.1/i386), and have run into the following roadblock: Dr. Gu's R code is a front end for a package of Fortran routines called RKPACK. While I have no trouble compiling the individual bits of Fortran code with fort77 in Linux, I don't know how to link the resulting *.o files into a shared

I can save to png, TIFF and jpg but get an error "Error: invalid graphics state" when trying to save as pdf and I have to restart R. This happens when I add mtext lines. There are a few other questions e.g. I want to move the mtext on side 1 to the left, but that is the main issue for now. I'm not using the ylab inside gplot.hexbin because it writes over the numbers on the axis.

Hi, due to some FTP problems in Auckland and me not being in office yesterday R 0.64.2 took some time to travel from Denmark to New Zealand and then back to Europe ... but now it is finally available at the CRAN master site http://www.ci.tuwien.ac.at/R here in Vienna and it should propagate to the various miroors within the next 24 hours or so. Best, Fritz

Hi, due to some FTP problems in Auckland and me not being in office yesterday R 0.64.2 took some time to travel from Denmark to New Zealand and then back to Europe ... but now it is finally available at the CRAN master site http://www.ci.tuwien.ac.at/R here in Vienna and it should propagate to the various miroors within the next 24 hours or so. Best, Fritz

Hi, I wanted to plot 2 lines on a single graph. Each graph has one axis that can be common. The code that I'm using is: ------------------------------------------------------- par(mfrow=c(1, 1)) x1 <- c(2, 4, 6, 8, 10, 12) x2 <- c(10, 20, 30, 40, 50, 60) y1 <- c(10,12,15,22,34,21) y2 <- c(40, 130, 150, 145, 40, 30) par(las=1, mar=c(4, 4, 2, 4)) plot.new()

I have been playing with the postscript graphics in 64.1 and don't seem to be able to get the character size right. It doesn't seem to respond to size options such as cex. Am I missing something that is specific to R? Chong Gu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send