similar to: Cleveland's Cut-and-Stack Plot

Hi, I think I found the problem. It lies in my Fortran program. Is there a way, after a DO loop, to make sure it does NOT return anything? Cheers, Kevin On Mon, 15 Apr 2002, Ko-Kang Kevin Wang wrote: > Date: Mon, 15 Apr 2002 17:27:20 +1200 (NZST) > From: Ko-Kang Kevin Wang <kwan022 at stat1.stat.auckland.ac.nz> > To: R Help <r-help at stat.math.ethz.ch> > Subject:

Hi, Say I want to do something like fitting 10 different sized trees with rpart() function. The only modification I need to do is to set 10 different cp's, which I have in a vector called foo. Can I do something like: for(i in 1:10) { rpart(y ~ ., cp = foo[i], data = mydata) } My problem is, I wish to save the 10 rpart objects into 10 different names, my.rpart1 ~ my.rpart10, for

Hi, Now suppose I have just one list called FOO, which has 25 objects, e.g.: [[1]] 1 2 3 4 5 [[2]] 6 7 8 9 10 . . . And I want to do something like: FRED <- data.frame(cbind(unlist(FOO[[1]]), unlist(FOO[[2]]), # ... for all 25 subsets )) Is it possible to do this, without doing unlist(FOO[[i]]) 25

Hi, Is it possible to draw barplot with x-axis being shown? I looked up the help file and I couldn't seem to find it. For example, I tried to do x <- 1:9 p <- log10(1 + 1/x) barplot(p, xlab = "Digit d", ylab = "Probability", ylim = c(0, 0.35), axes = F, main = "Benford's Law Probability") axis(1, 1:9) axis(2, seq(0, 0.35, by = 0.05),

Hi, Is it possible to add a subtitle that appears directly below the main title? I tried the "sub" parameter, but it adds sub-title to the bottom of the plot. Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Postgraduate PGDipSci Student Department of Statistics University of Auckland New Zealand Homepage:

I've just completed(?) an R Guide for Windows users. It was extended from my "Compile R for Windows" that I put up a few months ago. To get it go to http://www.stat.auckland.ac.nz/~kwan022/rinfo.php Here is the table of contents: 1 Introduction 2 Installation 2.1 Installing R Base 2.2 Installing packages 3 Running R 3.1 Rgui 3.2 Rcmd 3.3 Rterm 4 Comile

Hi, Suppose I have a function: myfunc <- function(x, y) { ... } And within the function I want to print out the name of the x, y vectors. For example, if I do: > myfunc(foo, goo) [1] "foo" "goo" It shall return "foo", "goo" (with or without quotes is fine), where foo and goo are two vectors with numbers. I know this sounds strange, but I'd

Hi, I'm learning Fortran and trying to load a Fortran subroutine into R. I've done: R SHLIB Fibonacci.f and it compiled fine. Then I went into R and done: > dyn.load("Fibonacci.so") > Fib <- function(n) { + .Fortran("Fibonacci", + as.integer(n))[[1]] + } > Fib(5) Error in .Fortran("Fibonacci", as.integer(n)) :

Hi, Suppose I have created something like this in R: foo <- "myfoo" and I want to find out the number of character in foo (in other words, R should return 5 since "myfoo" has 5 charactors. How can I do it? I tried: length(foo) but it returned 1. Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang

Hi, I've been playing with: > joinDate <- format(strptime(as.vector(forum[,2]), "%d-%b-%y"), + "%d-%b-%Y") > today <- format(strptime(as.vector("14-Jun-04"), "%d-%b-%Y"), + "%d-%b-%Y") > joinDate [1] "04-Feb-2004" "13-Feb-2004" "26-Feb-2004"

Hi, Suppose I have a for() loop that draws 6 boxplots as follows: par(mfrow = c(2, 3)) for(i in 2:length(spam.sample)) { boxplot(split(spam.sample[,i], yesno)) } Where spam.sample is a data frame with 7 columns, and I'm interested in plotting column 2 ~ 7 against column 1 (yesno). The boxplots appeared fine, however I'm trying to add a meaningful title, x and y labels to them.

Hi, Say I did something like: par(mfrow = c(1, 2)) plot(1:10) plot(1:10) and I'd like to have an overall title, i.e. a title that would appear, centered, on the top of both plots, rather than the top of the last plot. I tried: title("FOO", outer = T) but it doesn't appear right. Half of the FOO is outside the figure region and can't be seen... Is there another

Hi, Does anyone know how to manually configure the number of digits printed out from xtable()? For example, I'm exporting a data frame through xtable() into a LaTeX table, I only have two columns in the data frame so by default I only get two decimal places. But I'd like at least 5 decimal places. I had a look at ?xtable() but can't seem to find an example. Cheers, Kevin

Hi, Suppose I have: 459 1789 23590 and I'd like to round them to: 400 1700 24000 On the other hand, say if I have: 232 1234 23120 that need to be rounded to: 300 1300 24000 I tried the round(), floor() or ceiling() and can't get what I want. Is there any tricks I can use to achieve this goal? Cheers, Kevin

I have only just discovered the joys of the dotchart (since I am reading William Cleveland's -- Sean Carmody The Stubborn Mule http://www.stubbornmule.net http://twitter.com/seancarmody [[alternative HTML version deleted]]

I would rather use cleveland dot plots than bar charts to display my study results. I have not been able to find (or figure out) an R package that is capable of producing the publication quality dot charts Im looking for. I have either not been able to get error bars (lattice), cannot order the data display properly (latticeExtra), or cannot make adjustments to axes. Does anyone have a quick

Hi, Suppose I have the following lines at the end of a function: answer <- c(2, 1, 0, 4, 5) # In fact, answer will be generate in my # function print(answer) # Print the answer # Now, find the best fitted n degree polynomial print(paste("The best fit is with", which.min(answer) - 1, "-degree polynomial")) this will return:

I am trying to understand the nuances of STL (seasonal trend decomposition with loess) based on William Cleveland's (and others?) original development. I do not understand the specification or use of "frequency components" or equivalent "low-pass filter" components in the stl() function. I have run the stl() function on a standard example data (co2) in both S-Plus and

Hi, This is actually a theory question. I'm a bit confused by the s.window parameter in the stl() function (which is in the ts package). For example, in the stl documentation it uses the nottem data, and then: plot(stl(nottem, s.win = 4, t.win = 50, t.jump = 1)) What does it mean by s.win = 4? Is it because a year has 4 seasons (namely Spring, Summer, Autumn and Winter)? If so will it

Hi, Is there a package that contains testing algorithms (e.g. runs test, serial test, autocorrelation test...) for random number generators ? Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Postgraduate PGDipSci Student Department of Statistics University of Auckland New Zealand Homepage: http://www.stat.auckland.ac.nz/~kwan022