Vincent J. Carey, Jr.
2002-Jan-02 21:05 UTC
[Rd] comparative rendering of modeling outputs
This note is to r-devel rather than r-announce because it notes an experimental package that addresses issues that intersect with broader developmental issues in R. I have posted the package cremo = Comparative REndering of Modeling Outputs for retrival at http://www.biostat.harvard.edu/~carey/cremo.html This package addresses the problem of assembling and rendering results of multiple linear model fits in formats suitable for presentation in applications contexts. An extract from the output of 'example(cremo)': cremo> data(stackloss) cremo> attach(stackloss) cremo> glm2 <- glm(stack.loss ~ Air.Flow + Water.Temp + Acid.Conc.) cremo> glm3 <- glm(stack.loss ~ Air.Flow + Water.Temp) cremo> glm4 <- glm(stack.loss ~ Air.Flow * Water.Temp) cremo> try(print(cremo(list(glm2, glm3, glm4), adeqel = c("deviance", "aic")), quote = F)) var Model 1 Model 2 Model 3 [1,] (Intercept) -39.92 -50.359 22.29 [2,] (11.896) (5.138) (35.096) [3,] Air.Flow 0.716 0.671 -0.516 [4,] (0.135) (0.127) (0.58) [5,] Water.Temp 1.295 1.295 -1.93 [6,] (0.368) (0.367) (1.58) [7,] Acid.Conc. -0.152 -- -- [8,] (0.156) -- -- [9,] Air.Flow:Water.Temp -- -- 0.052 [10,] -- -- (0.025) [11,] ----- ----- ----- ----- [12,] deviance 178.83 188.795 150.232 [13,] aic 114.576 113.714 110.916 A preliminary form of the package, which passes "R CMD check" for the 12-14 version of 1.4.0, can be obtained at http://www.biostat.harvard.edu/~carey/cremo.html Much remains to be done for the problem confronted by this package, and some avenues to be explored are mentioned at the URL cited above. Anyone interested in the problem of specifying modeling report templates is invited to get involved; the code in cremo is extremely primitive and some more sensible approaches are sketched on the cremo.html page. Distributed under GPL2. -- --- Vince Carey, PhD Ass't Prof Med (Biostatistics) Harvard Medical School Channing Laboratory - ph 6175252265 fa 6177311541 181 Longwood Ave Boston MA 02115 USA stvjc@channing.harvard.edu -----BEGIN PGP PUBLIC KEY BLOCK----- Version: PGP 6.5.8 mQCNAzqIeGUAAAEEAMJXU941vIornTS52rl6z7eo+A7wwB0km/idLnkxzIhc1uLi Qtn19OyOfG6IDSucLrtmpvwagemAnQ9jL6TVDrmlrKnqsh+FFtvUuZ37eV85L70E BsS8RZCmMYHJKfrpCwegbTVZrEkd1ByquLIN/yUwxU4IcVuHxbNQk69riQ8tAAUR tBVWaW5jZW50IEouIENhcmV5LCBKci6JAJUDBRA6iHhls1CTr2uJDy0BAdsLA/wM cCzEDsP9MqodKZfDI1s/gXW6BcCuQ6n7MdEplLgmWvyqfbvRYx4upYZ3pNp8L0zU MrlR6eCTs/eDtMO/ZbGvkqqiQO6wS2fZb1T5L/DhhtT4mEAHt0E8dNBVCj+lKr3W vYS5GqO9gY4CiT3JXFH9N19pSbUQFiNDqpmG6EbWng==DQNF -----END PGP PUBLIC KEY BLOCK----- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._