similar to: R-help Digest, Vol 31, Issue 30

You've omitted a comma. races2000 is a data frame, which for purposes of extracting rows behaves like a 2-dimenional object. The following works fine: hills2000 <- races2000[races2000$type == 'hill', ] Additionally, you might like to ponder > type <- races2000[names(races2000)=="type"] > type[1:4] Error in "[.data.frame"(type, 1:4) :

Hello, how do I obtain standard errors of variances and covariances of the random effects with LME comparable to those of for example MlWin? I know you shouldn't use them because the distribution of the estimator isn't symmetric blablabla, but I need a measure of the variance of those estimates for pooling my multiple imputation results. Regards, Roel.

I'd use mcmcsamp() to examine the posterior distribution, under a relatively uninformative prior, of of the parameter estimates. For estimates that are based on four or five or more "degrees of freedom", I'd surmise that the prior will not matter too much. With estimates where the number of "degrees of freedom" is one or two or three, the posterior distribution may vary

expression() accepts multiple expressions as arguments, thus: plot(1:2, 1:2) legend("topleft", expression(y == a * x^b, "where "* paste(y=="wood; ", x=="dbh"))) Is there a way to do this when values are to be substituted for a and b? i.e., the first element of the legend argument to legend() becomes,

I am puzzled by the warning message in the output below. It appears whether or not I fit the seasonal term (but the precise point of doing this was to fit what is effectively a second seasonal term). Is there some deep reason why AR parameters ("Warning message: some AR parameters were fixed: ...") should somehow intrude into the fitting of a model that has only MA terms? >

I have posted, at http://wwwmaths.anu.edu.au/~john/r/newsets/ image (.rda) files, and first stabs at .Rd files for various data on deaths in London from 1629 to 1939. (There are of course gaps.) The sources (Guy 1882 & Stocks 1942) are documented in the .Rd files: (1) poxetc: measles, smallpox & total deaths: 1629-1881 [I have deliberately left several inconsistencies that were in

Under help(matrix) it is written: 'is.matrix' tests if its argument is a (strict) matrix. It is generic: you can write methods to handle specific classes of objects, see InternalMethods. Further down, under "Details", the meaning of "strict" is explained more explicitly: 'is.matrix' returns 'TRUE' if 'x' is a matrix

Typing in statements from the command line transfers to the computer a macro that was perhaps in the user's head. Putting together small sequences of carefully thought out code, in which key components have been carefully thought out and tested, has to be, for most of us, better than trying to make it all up on the run. An exception is necessary for those unusual people who who (akin to the

The following failed to make it to r-help when Murray Jorgensen sent his message a little time ago. [It seems to have been an issue of the machine from which my message originated.] Typing in statements from the command line transfers to the computer a macro that was perhaps in the user's head. Putting together small sequences of carefully thought out code, in which key components have been

I'd like to make the following change to termplot(): Add panel=points as an extra parameter. Replace if (partial.resid) points(xx, pres[, i], cex = cex, pch = pch, col = col.res, ...) by if (partial.resid) panel(xx, pres[, i], cex = cex, pch = pch, col = col.res, ...) Also, I think, add col.smooth, span and iter as named parameters; see below.

The following works fine with the x11 device, though it may well be that an initial plot is overwritten. With a pdf or postscript device, I get two plots, the first of which still has the red border from having the focus, while the second is the plot that I want. library(lattice); library(grid) plt <- xyplot(uptake ~ conc, groups=Plant, data=CO2) print(plt)

Although base graphics does not mix with lattice in the one graph, I've found that print.trellis(position=..., ) and the use of par(fig=...) to put regular and trellis graphics on the one graphics page works like a treat, at least in version 2.0.1 of R. [Base graphics functions that are themselves inconsistent with par(fig=...) are obviously disallowed.] I am wondering whether there are

The following is from version 1.6.2 of R under Windows, or 1.6.1 under Mac OSX/X11 > density(1) Error in if (!(lo <- min(hi, IQR(x)/1.34))) (lo <- hi) || (lo <- abs(x[1])) || : missing value where logical needed I am not sure how this should be handled. I encountered it in connection with densityplot(). In that connection, it might be enough to modify density() so that it

Is there any way, after use of print.trellis(), to obtain the co-ordinates of the plot region, e.g., in what are then the native co-ordinates? e.g. library(DAAG) library(lattice); library(grid) data(cuckoos) pushViewport(viewport(layout=grid.layout(2, 1))) pushViewport(viewport(layout.pos.row=1)) cuckoos.strip <- stripplot(species ~ length, data=cuckoos) print(cuckoos.strip,

This message is aimed at anyone who may be using exercises from my book with John Braun for teaching purposes. I am using this channel of communication in the absence of any other obvious effective channel. I ask the forbearance of list members. Our intention is to post solutions to selected exercises (the more challenging exercises) on the web, via a link from

In co-operation with Markus Hegland and myself, Zhongwen Ding has written a package, based partly on Markus Hegland's code, that provides a parallel processing interface to a remote multi-processsor system. Pyro (Python Remote Objects) and R must both be installed, both on the client machine and on the remote server. The system uses rsync, with an ssh protocol, to handle file transfer. Once

Peter Dalgaard writes: > Bill.Venables at cmis.csiro.au writes: > ........ > > [WNV] Now there's an insight! In fact I think that once you get used to > > the style you get to appreciate it. The commercial world substitutes > > politeness for candour. True candour can be disconcerting but it does get > > the message across much quicker. > Also, this

<soap-box> While I agree with John about abuse of methodology, I can't subscribe to the "failure to read the manual" proposition. I don't mean to lecture, but it would seem to me that people who, like me, will probably gain far more from this list than we will be able to contribute (at least in the near term) owe it to those who will likely contribute far more than they

There are answers that could and should be applied in specific situations. At least in academia and in substantial research teams, statisticians ought to have a prominent part in many of the research teams. Senior statisticians should have a prominent role in deciding the teams to which this applies. why should it be ok to do combine high levels of chemical expertise with truly appalling

termplot() does not carry subsetting over to carriers that are in the environment but not in the data frame. This generates a "subscript out of bounds" error. > data(ToothGrowth) > logdose <- log(ToothGrowth$dose) > tooth.lm <- lm(len ~ logdose, data=ToothGrowth) > termplot(tooth.lm) ## Works fine > toothVC2.lm <- lm(len ~ poly(dose,2),