search for: olau

Hi all, I'm not sure if this is a feature or a bug (and I did read the FAQ and the posting guide, but am still not sure). Some of my students have been complaining and I thought I just might ask: Let K be a vector of length k. If one types dim(K), you get NULL rather than [1] k. Is this logical? Here's the way I explain it (and maybe someone can provide a more accurate

Hi, Is there a way to make demo() print comments along with the code and output? For example, if I use something like readline("..."), demo prints both readline("...") and ...; as far as I can tell, this is also true of print() and cat(). I just want students to stop and think about every command in my demo scripts. Thanks, Olivia.

install.packages works for Linux because the directory structure is not version dependent, as it is for the Windows download. That should be fixed as soon as my sys admin consents to install R 1.9.0 so we can do the cross compile. Of course, you are right, Martin and Peter, that Zelig should be on CRAN proper. We just need to write a few .Rd files to be compliant with R check, and change our

Hi, I am working on a package which requires separate documentation (tutorial documentation, really, with a lot of beautifully latexed equations), and does not use R-help. I am trying to make it so that my help function will automatically search R-help if the help topic isn't found in my documentation. The function looks something like this: help.pkg <- function(topic) { #

Hi, I am working on a package which requires separate documentation (tutorial documentation, really, with a lot of beautifully latexed equations), and does not use R-help. I am trying to make it so that my help function will automatically search R-help if the help topic isn't found in my documentation. The function looks something like this: help.pkg <- function(topic) { #

Zelig: Everyone's Statistical Software Kosuke Imai, Gary King and Olivia Lau Version 1.0 (Available at http://gking.harvard.edu/zelig) A growing proportion of statisticians and methodologists from many disciplines are converging on R, a powerful statistics package and programming language. As an open source project, R is freely

Hi, I'm trying to find a way to take evaluate the Normal CDF over an interval and return the result on the log scale. This works, but I think it isn't numerically stable: log(pnorm(a, mean = x, sd = y) - pnorm(b, mean = x, sd = y)) Does anyone know of a single function that does the above? Or knows of a way to make it more stable? I'd really appreciate any suggestions!

Hi, I'm working on a set of demo files for a package, and I'm having a problem because I have two demo files (one called "logit", the other called "blogit") and when I type demo(logit), it gives me demo(blogit) instead. (And calls it the the demo for logit.) I've figured out that this is because the demo() function searches through all the demo directories (in

I am trying to plot two density lines on the same graph. Using the functions on the base package, I would go: plot(density(x), col = 1) lines(density(y), col = 2) And I get two distinct (one-bump) density lines. When I try to do it using lattice, I get two two-humped lines. (In other words, I think the smoothing function is taking the next set of data points and smoothing them in the same

Hi, I am trying to use boot() to refit an ordinal logit (polr in MASS) model. (A very basic bootstrap which samples from the data frame without replacement and updates the model.) I need to extract two statistics per run (the coefficients and zeta) and I tried concatenating them into a single vector after fitting, but I get the following error: Error in "[<-"(*tmp*, r, ,