Publications by Shige
Simulating second difference using Zelig
I am trying to simulate second difference using Zelig, here is my code:——————————————————–data(turnout)# estimationz.out model = “logit”, data = turnout)summary(z.out)# first differencex.low x.high s.out s.low s.high dif # second differencex.low.low x.low.high x.high.low...
827 sym
Ask R not to create a local directory tree
I don’t like R to create a local directory tree in my home directory because new packages will automatically be installed into that directory. The way to do this is to modify the “/usr/local/lib64/R/etc/Renviron” and mark the line “R_LIBS_USER=${R_LIBS_USER-‘~/R/x86_64-unknown-linux-gnu-library/2.12’}” off. Simple as tha...
742 sym
Extension to mtable function
Here are some useful extension to the “mtable” function in the memisc package. Related To leave a comment for the author, please follow the link and comment on their blog: Shige's Research Blog. R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to po...
484 sym
Legends in ggplot2
A simple plot takes a few lines of coding:g1 g2 geom_line(aes(y=alive1, linetype=”Male”)) + scale_linetype_discrete(name = “”)g3 geom_point(aes(y=alive1, shape=”Male”)) + scale_x_continuous(“”) + scale_y_continuous(“Proportion Alive in 1982″) + opts(legend.position=c(.24, .95), legend.justification = c(1, 1)) + s...
763 sym 2 img
Ra vs. compiler package
R seems to have two byte code compilers: the Ra add-on module (and the accompanying “jit” package) and the “compiler” package came with the default installation. I wonder how they differentiate from each other and what the strengths and weaknesses of each are. Related To leave a comment for the author, please follow the li...
678 sym
Updated examples that combine Rcpp and CppBugs
This example is very informative in illustrating how to use CppBugs module with Rcpp to conduct fast MCMC simulation in R. Related To leave a comment for the author, please follow the link and comment on their blog: Shige's Research Blog. R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other to...
528 sym
GEE using Stata vs. R
I am running GEE logistic regression model for my fetal loss paper. As usual, I compare results between Stata and R and make sure they are consistent. To my surprise, the models assuming independent correlation structure give similar results but the models assuming exchangeable correlation structure give drastically different results....
1094 sym
Performance difference between Stata and R
With respect to multinomial logit model, the performance difference between the two packages are quite large, based on this post. Related To leave a comment for the author, please follow the link and comment on their blog: Shige's Research Blog. R-bloggers.com offers daily e-mail updates about R news and tutorials about learning ...
547 sym
SabreR
SabreR just released an update. It is another software package that can estimate multivariate multilevel model (other options are aML, MCMCglmm, etc.). They seem to also have a book dedicated to the software, which be worth checking out.It will be great if the author can incorporate some plotting function into the package. Related To leave a...
730 sym
Visual debugger and the debug mode of the autorun R console
The StatET team kept their promise and delivered the autorun R console with debug mode on. This, combined with the visual debugger, makes the StatET a very appealing cross-platform environment for working with R. Related To leave a comment for the author, please follow the link and comment on their blog: Shige's Research Blog. R-...
630 sym