[R] trying to use standard notation

From: Andrew Gelman <gelman_at_stat.columbia.edu>
Date: Wed 10 May 2006 - 00:46:03 EST

Hi, all. In setting up my package for post-processing regression models, I am trying to use standard notation as much as possible: thus, I use coef() to access estimated coefficients. I wrote a function called se.coef() to grab standard errors, and se.fixef() and se.ranef() to grab se's from coefficients estimated from lmer().

I also need a function to access sigma-hat (the residual sd of the regression). I've called this function sigma.hat(), but I'm open to using a better name if there's something more standard.

I also have a function called sim() that makes simulations using lm() or glm() output (based on the curvature of the likelihood function or, to put it another way, approximate Bayesian inference assuming a flat prior distribution). sim() returns a list: with two components: (1) beta, a matrix of simulated coefficients (with dimension n.sims by k), and (2) sigma, a vector of simulated sigma's (with length n.sims).

I'll be doing more, but these are the basics for now. I need to be able to pull these objects out as necessary. And I just was wondering if anyone had suggested names for the functions that are more consistent with R naming conventions.


Andrew Gelman
Professor, Department of Statistics
Professor, Department of Political Science

Statistics department office:
  Social Work Bldg (Amsterdam Ave at 122 St), Room 1016
Political Science department office:
  International Affairs Bldg (Amsterdam Ave at 118 St), Room 731

Mailing address:
  1255 Amsterdam Ave, Room 1016
  Columbia University
  New York, NY 10027-5904
  (fax) 212-851-2164

R-help@stat.math.ethz.ch mailing list
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Wed May 10 00:51:18 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 10 May 2006 - 06:10:04 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.