similar to: missing data and completed missing data

Hi, How can I plot two scatterplots on the same, one panel? I have two times series (price data for two goods) for the same period. Many thanks, -R [[alternative HTML version deleted]]

Hi everyone, I have tried to solve a simple problem for days but I can't figure out how to run it properly. If someone could give me a hint, this would be really great. Basically, I want to run a standard economist's fixed, and random effects regression (corresponds to xtreg in STATA) but with _variable_ weights (they correspond to changing industry shares in the market). Here is

I have a question on the output generated by randomForest in classification mode, specifically, the confusion matrix. The confusion matrix lists the various classes and how the forest classified each one, plus the classification error. Are these numbers essentially averages over all the trees in the forest? If so, is there a way I can get the standard deviation values out of the randomForest,

Dear R users, I was wondering if anyone knows how to obtain(subset) the first and/or the last record of a subject in a longitudinal setup. Normally in SAS one uses first.variable1 and last.variable1. So my question is that is there an R way of doing this. Regards, -- Luwis Diya, Phd student (Biostatistics), Biostatistical Center, School Of Public Health, Catholic University of Leuven, U.Z. St

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

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

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

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

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

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

I have the feeling that this is a dumb question but oh well ... When I have an object x with several attributes, y <- x will make a copy of x with all the attributes included. What if I just want y to get the value assigned but nothing else (no class, names, dim or other attributes)? Johann -- Phone/FAX/Adr/PGPkey ... at http://www.ai.univie.ac.at/~johann

Consider the following: > A <- 1:10 > A [1] 1 2 3 4 5 6 7 8 9 10 > dim(A) NULL > dim(A) <- c(2,5) > A [,1] [,2] [,3] [,4] [,5] [1,] 1 3 5 7 9 [2,] 2 4 6 8 10 > dim(A) [1] 2 5 > dim(A) <- 10 > A [1] 1 2 3 4 5 6 7 8 9 10 > dim(A) [1] 10 Would it not make sense to have dim(A) = length(A) for all vectors?

To whom it may concern: Given that my sample size is n, my mean is 100, and my sd is 10, I need to use a random number generator (which I believe is the function rnorm(5,100,10)), but I need to repeat it a large number of times, and then plot the sampling distributions of the sample means, sd's, and variances of those generated sets. I'm having a real hard time trying to figure out how

Hi, What is the best text editor for programming in R? I am using JEdit as the text editor, however, it does not have anything specific for R. It will be nice to have a developing environment where the keywords are highlighted, plus some other debugging functions. Yi-Xiong [[alternative HTML version deleted]]

[Using R 2.2.0 on Windows XP; OK, OK, I will update soon!] I have noticed some undesirable behaviour when applying ifelse to a data frame. Here is my code: A <- scan() 1.000000 0.000000 0.000000 0 0.00000 0.027702 0.972045 0.000253 0 0.00000 A <- matrix(A,nrow=2,ncol=5,byrow=T) A == 0 ifelse(A==0,0,-A*log(A)) A <- as.data.frame(A) ifelse(A==0,0,-A*log(A)) and this is the output:

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