similar to: Data Checking

Hi, This is a very trivial question, and I do not know why I cannot remember the answer at all! How does one get the entire system information (which OS, R version....etc) using an R command? Thanks, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Postgraduate PGDipSci Student Department of Statistics University of Auckland New Zealand

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, 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, I was trying to read in the following data set by using: gas <- read.table("gas.tab", header = T) when I got the error message: Error in "row.names<-.data.frame"(*tmp*, value = row.names) : duplicate row.names are not allowed Any help are appreciated. Thanks, Ko-Kang Wang

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 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, I was just trying out an example on Page 247 of the MASS (Modern Applied Statistics with S-plus) book and saw the function "negexp.ival". It says it is supplied in the MASS library, however when I load the library in R and typed: negexp.ival it says: Error: Object "negexp.ival" not found Does it mean it only appears in the MASS library for S-plus?

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, 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 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, I'm not sure if this can be done but.. I know that with ifelse() I can do something like: ifelse(x <= 3, 1, 2) to go through each element in my vector x, and if x_i <= 3 substitute the number with 1 else with 2. Essentially I'll get a vector with 2 levels. Can I tweak it so I can get 3-levels? For example: if(x <= 3) then 1 elseif(3 < x <= 4) then 2

Hi, In the order() help file, there is an example like: a <- c(4, 3, 2, NA, 1) b <- c(4, NA, 2, 7, 1) z <- cbind(a, b) (o <- order(a, b)); z[o, ] How can I do something like "order a in ascending order, b in descending order"? And say I have a third vector c, and I'd like to add this to the previous condition "a ascending, b descending, c

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, 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, Is it possible to create a ts() object, whose data is daily based BUT measured only on working days? In other words, suppose I have a data set with 255 observations, measured from 29 June 1959 to 30 June 1960. How would I create such a data? I tried something like: ts(c(...), start(1959, 180)) but I'm not sure what to use for frequency. In other words I don't know how to

I just installed the R-patched 1.4.0 and for some reason I now have trouble compiling R packages - I could compile it a few weeks ago with R 1.4.0 released in December. This only happens if I put the --binary option in. For example: Rcmd build --binary cw returns the error: make: Entering directory `/cygdrive/c/PROGRA~1/R/RW1040/src/gnuwin32' make DLLNM= \ -C C:/MyRWork/cw PKG=cw

R 1.6.0 has gone into final feature freeze as of today. As a new feature, we'll make interim beta versions available via ftp://cran.us.r-project.org/pub/R/src/base alias http://cran.us.r-project.org/src/base (filename R-1.6.0beta_*.tar.gz, where * is the creation date). If you want to help ensure that the final 1.6.0 works satisfactorily on *your* platform, it might be a good idea to

On Thu, 19 Sep 2002 10:43:01 -0400 (EDT), you wrote in message <200209191443.KAA17404@falmouth.bwh.harvard.edu>: >Hi Duncan -- I am curious as to whether you are building/distributing >R 1.6.0 for windows? There's nothing urgent about it, I have one. >The question is how to identify a distribution URL if one exists. I will be building and distributing the final release.

Hello everybody, Could someone please send me a very simple example using Fortran from R? Say pass a value to an executable and get the result in R. Actually it seems it may be possible to call an *.f file ?? or I am wrong again? The manual is very terse on the subject. Thank you very much Stephen Elijah -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help