similar to: Random Number Testing

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

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, 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, 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)) :

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

Hi, Is there a function in R that does Cleveland's Cut-and-Stack plot (Page 190 -- 191, The Elements of Graphing Data, William S. Cleveland)? Or do I need to do it the hard way, i.e. set par(mfrow = c(m, n)) then do it one-by-one? (I have a time series data set that is almost identical to the description in Cleveland's book, hence I'm interested in trying the Cut-and-Stack plot)

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:

Hi, Suppose I have two lists. The first list is called FOO while the second is called FRED. Say FOO looks (I've simplifed it) like: [[1]] [,1] [1,] 1 [2,] 2 [[2]] [,1] [1,] 3 [2,] 4 while FRED looks like: [[1]] [,1] [1,] 5 [2,] 6 [[2]] [,1] [1,] 7 [2,] 8 Can I turn this list into a dataframe which

Hi, Is it possible to put the legend outside the plotting region? Say, I have some 24 lines within my plot, which pretty much filled the plotting region and left me no space (well not enough space) to put in a legend. So I'd like to put the legend outside the plotting region (i.e. the region enclosed by the two axis -- I hope I got my terminology right). Actually, I did the 24 lines with

Hi, This is an easy question, however I cannot remember the exact command of doing it :-( Say I have a matrix: > ff [,1] [,2] [1,] 0.1000000 116.76987 [2,] 0.6444444 108.26811 [3,] 1.1888889 95.00240 [4,] 1.7333333 112.13647 [5,] 2.2777778 107.39893 [6,] 2.8222222 115.34615 [7,] 3.3666667 97.81029 [8,] 3.9111111 105.35119 [9,] 4.4555556

Hi, Suppose I got a ts() object that looks like: Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 1983 587 498 617 502 566 495 500 522 592 574 629 400 1984 620 545 550 500 604 533 547 596 530 616 672 412 1985 625 497 624 540 645 493 604 610 505 650 633 403 and I'd like to extract out 1984 data. How can I do it? I tried to subset it as

Hi, I've been trying out postscript() when onefile is set to FALSE. However I'm having trouble understanding the documentation for postscript(). Suppose I want to do: plot(1:10) plot(1:10, type = "l") and store them in foo.ps and goo.ps. Furthermore, I'd like to do then with one postscript() command -- because this is just a simplified version of what I really need