similar to: eval(), expression()

Dear R-List, I have a dataframe X containing factor f and numeric variable x1, x2, ... I want to create a new dataframe (or possibly a matrix) that gives statistics (e.g. sum) for the variables x1, x2, ... in each group defined by factor f. What is the simplest way to do this? I tried: attach(X) Z <- data.frame(f=levels(f),x1=as.vector(lapply(split(x1,f),sum))) and stumbled on: Error

Hi! The postscript device in 0.90 does not produce symbols (pch=) as it should. E.g.: postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) plot(1:10,1:10,pch=1:5) legend(1,9,letters[1:5],pch=1:5) dev.off() The result looks quite different compared to the X-windows version. The symbol 1 (pch=1) is much to small, while the other symbols seem to big. Is it a bug? How can I

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

Dear list, I want to perform a biplot, using customized titels for the x and y axis. Setting xlab="" and ylab="" resulted in an error, e.g.: > data(USArrests) > biplot(princomp(USArrests),xlab="",ylab="") Error in biplot.default(t(t(scores[, choices])/lam), t(t(x$loadings[, : length of dimnames[1] not equal to array extent > How do I

Hello! When I know the vector of the variance of the disturbances (i.e. the structure of heteroskedasticity), say Var(u_{i})=v_{i}, what is the weights I should use as argument to lm(): M <- lm(y~x,weigths=1/v) or M <- lm(y~x,weights=1/(v^0.5)) ??? In the help pages I did not find a clear answer to this question, so please could someone help me! Thanks, Wolfgang Koller

Hello, I know this has been covered on here before, but as a complete novice, I need a little more guidance. I would like to produce boxplots with the whiskers extending to the 10 and 90th percentiles. I found this code: myboxplot.stats <- function (x, coef = NULL, do.conf = TRUE, do.out = TRUE) { nna <- !is.na(x) n <- sum(nna) stats <- quantile(x, c(.1,.25,.5,.75,.9), na.rm

Dear R-developpers, currently the function window shows the following behaviour, if start or end is out of range: > x <- ts(1:14,start=c(1950,2),frequency=4) > y <- window(x,start=c(1950,1),end=c(1953,4)) Warning messages: 1: start value not changed in: window.default(x, ...) 2: end value not changed in: window.default(x, ...) Instead of this I would like it to fill with NAs the

From: "Wolfgang Koller" <koller2@fgr.wu-wien.ac.at> Organization: University of Economics, Vienna Date: Tue, 23 Nov 1999 16:29:22 CET (+0100) MIME-Version: 1.0 Content-type: text/plain; charset=US-ASCII Content-transfer-encoding: 7BIT CC: r-bugs@biostat.ku.dk, r-help@stat.math.ethz.ch References: <153DF9230E8@fgr.wu-wien.ac.at> X-mailer: Pegasus Mail

Full_Name: Wolfgang Koller Version: 1.0.0 OS: linux Submission from: (NULL) (137.208.7.48) > data(UKLungDeaths) > ts.union(mdeaths, fdeaths,dframe=TRUE) Error in names<-.default(*tmp*, value = nmsers) : names attribute must be the same length as the vector > cbind(mdeaths, fdeaths,dframe=TRUE) Error in names<-.default(*tmp*, value = nmsers) : names attribute

Dear R-Users and Developpers, Currently R does not include functions to compute the skewness and kurtosis. I programmed it myself in the following way, but probably *real* programmers/statisticians can do that better: mykurtosis <- function(x) { m4 <- mean((x-mean(x))^4) kurt <- m4/(sd(x)^4)-3 kurt } myskewness <- function(x) { m3 <- mean((x-mean(x))^3) skew <-

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Hi! Since installation of R 0.90 there are problems with the size of text and numbers also in the x11-device (I reported about similar problems with the postscript device recently). For numbers and letters a different size is used: plot(1:10,1:10,type="n") text(2,8,"F5F5F5F5F5") text(2,7,"F5F5F5F5F5",cex=2) The same problem is also found in the postscript

Hi all, I'm just getting started with R and I would appreciate some help. I'm having trouble creating a boxplot with whiskers at the 95th and 5th percentiles instead of at 1.5 * IQR. I have read the relevant documentation, and checked existing mails on this topic. I found a small modification that should work : https://stat.ethz.ch/pipermail/r-help/2001-November/016817.html and tried to

Hi, I'm new to the whole R-thing as a replacement for Matlab, not disappointed sofar ;) I found out how to make nice looking boxplots, but i also would like the make a boxplot with 5% and 95% instead of the standard 25 and 75% quantiles. My csv input looks something like: LOCATION FILTER NR DATE VALUE MONTH Peelhorst01 1 14-Jan-94 23.07 1 Peelhorst01 1 28-Jan-94 23.68 1 Peelhorst01 1

Hi everyone I'm trying to generate a vector of the form a <- c(1,2,3,4,1,2,3,1,2,1) where n = 5 in a general way, where n can be any positive integer. I've run out of ideas. Does anyone have any suggestions? Thanks Karen Karen Kotschy Centre for Water in the Environment University of the Witwatersrand Johannesburg Tel: 011 716-2218

Hi everyone I'm trying to generate a vector of the form a <- c(1,2,3,4,1,2,3,1,2,1) where n = 5 in a general way, where n can be any positive integer. I've run out of ideas. Does anyone have any suggestions? Thanks Karen Karen Kotschy Centre for Water in the Environment University of the Witwatersrand Johannesburg Tel: 011 716-2218

Thanks for writing this great piece of code. x = rnorm(100) boxplot(x) # you shouldn't see any outliers here although sometimes yow will # lets add some outliers intentionally x = c(21, 20, 25, x) # now 10, 15 and 20 are outliers myboxplot <- boxplot(x) # now you should see your three outliers myboxplot$out # it will print the values of the outliers How does one amend

>>>>> "WoKo" == Wolfgang Koller <koller2@fgr.wu-wien.ac.at> writes: WoKo> Hi! The postscript device in 0.90 does not produce symbols WoKo> (pch=) as it should. E.g.: WoKo> postscript(file="test.ps",width=8,height=6.4,horizontal=FALSE) WoKo> plot(1:10,1:10,pch=1:5) WoKo> legend(1,9,letters[1:5],pch=1:5) WoKo>

People, It appears that there is no way of getting Boxplots to plot using Mean, SD, Max & Min - is there something else that would do what I want? I couldn't find it . . Thanks, Phil. -- Philip Rhoades GPO Box 3411 Sydney NSW 2001 Australia E-mail: phil at pricom.com.au

Dear all, first, I would like to thank R-core team for making so a great software available. Thank you all, really! Using R (1.0.0) on a Win98 platform, I need an x-server to run xgobi (through R xgobi interface). I've been using a MI/X server but trying to upgrade I noticed that MI/X is no longer available for free. Similarly, eXceed is available at evaluation version (for a limited time