similar to: Macros versus functions

How do you stop R from putting up a dialog box when you quit Rgui? (I use Windows and I never save workspaces that way) Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz Fax 7 838 4155 Phone +64 7 838 4773 wk +64 7 849 6486 home

I want to interlace two vectors. This I can do: > x <- 1:4 > z <- x+0.5 > as.vector(t(cbind(x,z))) [1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 but this seems rather inelegant. Any suggestions? Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz

Greetings all. Any help with the following would be appreciated. I want to create a data frame for each file in a directory. The following code does not work but it may show what I am trying to do: carmakes <- c('BMW','Chrysler','Citroen','Fiat','Ford','Holden','Honda',

I got a shock a few days ago when I accidentally visited www.r-project.com . I thought that the r-project site had been hacked until I realised my mistake. There is also a site www.r-project.net. Both of these sites appear to be Japanese. Does anyone know anything about them? I suppose that it is not unusual for names close to those of popular sites to be used. It is good that they use a

Well I don't know if anyone can help with this but it will be interesting to know if others have had the same problem. I can't start R at home on my laptop [ I'm using 1.8.1 under Windows XP]. When I click on the shortcut I get the usual Windows box for when an application needs to close. A couple of clicks down it displays the following: Error signature AppName: rgui.exe

I found myself wanting to average a vector [vec] within each level of a factor [Fac], returning a vector of the same length as vec. After a while I realised that lm1 <- lm(vec ~ Fac) fitted(lm1) did what I want. But there must be another way to do this, and it would be good to be able to apply other functions than mean() in this way. Cheers, Murray -- Dr Murray Jorgensen

I am exploring the result of clustering a large multivariate data set into a number of groups, represented, say, by a factor G. I wrote a function to see how categorical variables vary between groups: > ddisp <- function(dvar) { + csqt <- chisq.test(G,dvar) + print(csqt$statistic) + print(csqt$observed) + print(round(csqt$expected)) + round(csqt$residuals) + } > > x

I am introducing the scan() function to my class. Consider the following file (Scanexamp.Rnw ) \documentclass[12pt]{article} \begin{document} <<>>= height = scan() 64 62 66 65 62 69 72 72 70 part = scan(what = character(0)) "Soprano" "Soprano" "Soprano" "Alto" "Alto" "Tenor" "Tenor" "Bass"

I'm going to try my hand at converting some Fortran programs to R. Does anyone know of any good articles giving hints at such tasks? I will post a selective summary of my gleanings. Cheers, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz

I have just noticed that quite a few columns of a data frame that I am working on are numeric factors. For summary() purposes I want them to be vectors. I tried, for example > indx <- c(1:18,21:37,40,41) > i <- 0 > i <- i+1 > summary(as.vector(sflows1[indx[i]])) Length Class Mode min.pkt.sz 3000 factor numeric but this does not give the five-number

G'day list! I'm thinking about programming likelihood and score intervals for generalized linear models in R based on the paper "On the computation of likelihood ratio and score test based confidence intervals in generalized linear models" by Juha Alho (1992) (Statistics in Medicine, 11, 923-930). Being lazy, I thought that I would ask if anyone else on the list has

Does anyone else miss email-friendly character graphics such as the following example, produced using Minitab? Histogram of C6 N = 478 N* = 21 Each * represents 2 observation(s) Midpoint Count -12 16 ******** -11 53 *************************** -10 63 ******************************** -9 83

I'm running R 1.2.2 under windows 98 on a Pentium 133 laptop. I can't seem to retrieve the package MASS data sets: > library(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) > data(abbey) Warning message: Data set `abbey' not found in: data(abbey) And yet all the .rda files for the MASS datasets are in D:\Program

I wanted a sort of matrix product of an array and a matrix. As there does not seem to be any array multiplication apart from outer() I proceeded as follows: lambda <- array(0, c(n,m,d)) # stuff omitted # zed is an n by m matrix # # \lamb.star_{ik} lamb.star <- matrix(0, nrow=n, ncol=d) for (i in 1:n) { for (k in 1:d) { for (j in 1:m) { lamb.star[i,k] = lamb.star[i,k] +

Has anyone got a cute way of compiling a table of the available plotting symbols? I've been repeatedly pasting in cno <- cno + 1 cno plot(x,y,pch=cno) which is tedious, and besides I have to note down the symbols by hand. Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at

Having just update to R 2.6.2 on my old Windows laptop I notice that the number of packages is growing exponentially and my usual approach of get-em-all may not be viable much longer. Has any thought been given to dividing "contributed" binaries into a recommended set, perhaps a couple of hundred, and the remained. That way one could install the recommended ones routinely and add in

Suppose you fit a linear model > model.1 ~ lm(v1 ~ ..., data=myframe) and v2 is some other column of myframe typically not in the model. You will often want to try > plot(v2, model.1$residuals) but this will fail if there are NAs in the response v1 as model.1$residuals has length equal to the number of nonmissing values in v1. I suppose > plot(v2[!is.na(v1)],

Here is the file minimal.Snw: \documentclass[a4paper]{article} \title{R tips and tricks} \author{Murray Jorgensen} \usepackage{Sweave} \begin{document} \maketitle \section*{Entering data from a single variable} The following data are transformed tensile strength measurements on polyester fibres. They may be found on the file \texttt{TENSILE.DAT}. We may enter this data into R using the

I want to paste in the following numbers into a scan: 0.023 0.032 0.054 0.069 0.081 0.094 0.105 0.127 0.148 0.169 0.188 0.216 they are separated by tabs alone, unless my mailer has done something to the tabs. Now have a look at this: > scan() 1: 0.0230.0320.0540.0690.0810.094 1: 0.1050.1270.1480.1690.1880.216 Error in scan() : "scan" expected a real, got

I can ask my question using and example from Chapter 1 of Pinheiro & Bates. > # 1.4 An Analysis of Covariance Model > > OrthoFem <- Orthodont[ Orthodont$Sex == "Female", ] > fm1OrthF <- + lme( distance ~ age, data = OrthoFem, random = ~ 1 | Subject ) > summary( fm1OrthF ) Linear mixed-effects model fit by REML Data: OrthoFem AIC BIC