similar to: appeal --- add sd to summary for univariates

Sorry, one more: on OSX, I deleted my old 2.9.2 R.app, and installed the 64 bit version of 2.9.0. I then did an "install.packages("car")" under my new 2.9.0. It seems to have worked, but alas, I still get an error that package 'car' was built under R version 2.9.2 . Where exactly does R under OSX install its packages? (is it a bug that another car is loaded?) PS:

the following code apprantely, for some grand old reason, induces a syntax error: if (seq[i] = "A") m <- trans[1,] where seq is a vector and trans is a matrix. I cannot for the life of me see why this is wrong. Syntax error is: Error: syntax error in "if (seq[i] =" Sincerely, Norman Goodacre ---------------------------------

dear R experts---now I have a case where I want to estimate very large regression models with many fixed effects---not just the mean type, but cross-fixed effects---years, months, locations, firms. Many millions of observations, a few thousand variables (most of these variables are interaction fixed effects). could someone please point me to packages, if any, that would help me estimate such

Dear R wizards: is there a clean way to assign to elements in a list? what I would like to do, in pseudo R+perl notation is f <- function(a,b) list(a+b,a-b) (c,d) <- f(1,2) and have c be assigned 1+2 and d be assigned 1-2. right now, I use the clunky x <- f(1,2) c <- x[[1]] d <- x[[2]] rm(x) which seems awful. is there a nicer syntax? regards, /iaw ---- Ivo Welch

dear R experts---I am using the coef() function to pick off the coefficients from an lm() object. alas, I also need the standard errors and I need them fast. I know I can do a "summary()" on the object and pick them off this way, but this computes other stuff I do not need. Or, I can compute (X' X)^(-1) s^2 myself. Has someone written a fast se() function? incidentally, I think

this is not an important question, but I wonder why lm returns an error, and whether this can be shut off. it would seem to me that returning NA's would make more sense in some cases---after all, the problem is clearly that coefficients cannot be computed. I know that I can trap the lm.fit() error---although I have always found this to be quite inconvenient---and this is easy if I have only

As an end-user, I wonder about Revolution R. Is the relationship between Revolution R and the R community at-large a positive one? Do the former contribute to the development efforts of the latter? Is there a competitive aspect? is their forum competitive with r-help? any other thoughts? (most of all, I simply hope that they help some of the many helpful experts on this forum, who have

Dear R wizards: I decided to take the advice in the R data import/export manual and want to learn how to work with SQL for large data sets. I am trying SQLite with the DBI and RSQLite database interfaces. Speed is nice. Alas, I am struggling to find a tutorial that is geared for the kind of standard operations that I would want in R. Simple things: * how to determine the number of rows in a

Dear community, I'm trying to install R-2.2.1 on an IRIX 6.2 (Unix System V Release 4) system without root access. Unfortunately readline is not installed in default, so I installed it locally in my home directory, more precisely in: $HOME/vol/readline-5.1, where $HOME is "/home3/fa/faga001". Afterwards I appended the path to the library with several $PATH variable, which now looks

Dear community, I'm trying to install R-2.2.1 on an IRIX 6.2 (Unix System V Release 4) system without root access. Unfortunately readline is not installed in default, so I installed it locally in my home directory, more precisely in: $HOME/vol/readline-5.1, where $HOME is "/home3/fa/faga001". Afterwards I appended the path to the library with several $PATH variable, which now looks

I'm installing Xvid 1.0.2 in wine and this is the error I receive. This aborts the installation. What could be causing this? I don't see a log file in /var for wine or else Iw ould post it, If there is any file that I should post to help diagnose this please let me know! Thanks!

Dear R wizards---I read Brian Ripley's post from 2004 which said that it was not possible to print to STDERR. Alas, I have more modest needs. I was wondering if it was possible to just send a string to STDERR with cat() while in CMD BATCH mode. Is it not possible to open STDERR in R? (Or does R use STDERR for itself and redirect it into the output stream?) This would be on a standard Unix

Dear R wizards: This must have an obvious solution, but I am stumped. I can run a linear regression giving a matrix as the independent set of variables, but if I give a data frame (which I would like to give, because it should tell the linear model the names of the variables), R does not like it. An example is: N=20; y= rnorm(N); x.m <- (matrix( nrow=N, ncol=2 )); x.m[,1]=rnorm(N);

dear R group: I have the probably fairly common problem that I would like to have one code.R file do different things if it is invoked from a symbolic link, which should be easy to uncover. $ ln -s code.R code-0.R $ ln -s code.R code-1.R $ R CMD BATCH code-1.R what needs to be in code-1.R to put code-1.r into a character vector? help appreciated. regards, /ivo welch PS : I read

dear R wizards:? I am wrestling with reshape.? I have a long data set that I want to convert into a wide data set, in which rows are firms and columns are years. > summary(rin) firm fyear sim1 Min. :1004.00 Min. :1964.0 Min. : -1.00000 1st Qu.:1010.00 1st Qu.:1979.0 1st Qu.: -0.14334 Median :1016.00 Median :1986.0 Median : 0.00116 Mean

dear R wizards: here is the strange question for the day. It seems to me that nrow() is very slow. Let me explain what I mean: ds= data.frame( NA, x=rnorm(10000) ) ## a sample data set > system.time( { for (i in 1:10000) NA } ) ## doing nothing takes virtually no time user system elapsed 0.000 0.000 0.001 ## this is something that should take time; we need to add 10,000

dear R experts: does R have "static" variables that are local to functions? I know that they are usually better avoided (although they are better than globals). However, I would like to have a function print how often it was invoked when it is invoked, or at least print its name only once to STDOUT when it is invoked many times. possible without <<- ? sincerely, /iaw

dear R experts--- t <- 1:30 f <- function(t) { cat("f for", t, "\n"); return(2*t) } g <- function(t) { cat("g for", t, "\n"); return(3*t) } s <- ifelse( t%%2==0, g(t), f(t)) shows that the ifelse function actually evaluates both f() and g() for all values first, and presumably then just picks left or right results based on t%%2.

dear r-experts: I need to speed up my monte-carlo simulations. my code is written in R (and it was also the cause of my many questions here over the last few days). my code is almost all matrix/vector algebra on panel data sets---long-difference, fixed-effects, blundell-bond, etc.. the data set is about 10MB, so 1GB per CPU core should be plenty for my operations, and with $10/GB of

I wonder if this is an intentional feature or an oversight. in some column summaries or in ifelse operations, apparently I am losing the date property of my vector. > a <- c(198012, 198101, 198102) > b <- a*100+31 > c <- as.Date( as.character(b), "%Y%m%d" ) > summary(c) Min. 1st Qu. Median Mean 3rd Qu. Max.