similar to: janka data

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

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

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

I'd like to find the dataset files cloud.data, janka.data, house.data ... used in Appendix C of "Notes on R" by Bill Venables and Dave Smith. Can someone give me a pointer to a download site, please? With thanks, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk> Date: 17-Jun-98

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

Hello: I sent this question the other day with the wrong subject heading and couple typos, with no response. So, here I go again, having made those corrections. I would like to estimate, for lack of a better description, a partially additive non-parametric model with the following structure: z~ f(x,y):w1 + g(x,y):w2 + e In other words, I'd like to estimate the marginals with respect to

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

I was reading Chapter 2 of the R manual that came with 64.1, but couldn't get the behavior out of expressions as described in the manual. Is there any example on the subject? For example, "x sub i" is supposed to be produced by "x[i]", but where to put it? I tried mtext(text=x[i]) and R simply interprets that as a regular expression. Thanks in advance for any clue that

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.

I am embarrassed to have to ask this but can anyone tell me of a Tukey multiple comparisons procedure available for R? I have looked through the search page, through the FAQ, and in the index of V&R (1999), and I still can't find such a thing. I see there is a ptukey function and a qtukey function but that is as far as I got. Do I need to roll my own? -- Douglas Bates

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

I am just about working through Faraways excellent tutorial "practical regression and ANOVA using R" on page 24 he makes the x matrix: x <- cbind(1,gala[,-c(1,2)]) how can I understand this gala[,-c(1,2)])... I couldn't find an explanation of such "c-like" abbreviations anywhere. thanks for a hint. another problem: I couldn't load the faraway library, using the