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... [Read more...]
I am trying to simulate second difference using Zelig, here is my code:--------------------------------------------------------data(turnout)# estimationz.out [Read more...]
Stata's new "margins" command is as powerful as "Zelig" package for R for the purpose of conducting post-estimation manipulation of the results (e.g. producing predicted probabilities, confidence intervals, etc.). While "margins" uses delta method, "Ze... [Read more...]
I am doing a simple comparison of different estimation procedures in dealing with a simple binomial model. Here is where I got started:---------------------------------------------library(INLA)library(npmlreg)library(MCMCglmm)library(DPpackage)data(See... [Read more...]
The increasing number of R-oriented Bayesian computational tools such as MCMCpack, MCMCglmm, DPpackage, R-INLA, spBayes, have made BUGS less and less crucial for day to day Bayesian computation. Honestly, I cannot figure out a single analysis that BUGS... [Read more...]
Both lattice and ggplot2 seem really interesting and worthy of learning. But I only have time to learn one of them, and the choice is not an easy one.Here is an awesome reference; this blog is generally very interesting; and here is something... [Read more...]
Two things that are crucial for a wider use of R among applied researchers. The first one is data manipulation/reshaping tool. I think the package "reshape" and "reshape2" have done good job and have largely removed the barrier. The second one is ... [Read more...]
The tikzDevice package is quite amazing. Here are two graphs I just made, with (lower) and without (upper) using the tikzDevice package. The difference in quality is huge.Since the LaTeX source file for the figure is quite large in size and may take si... [Read more...]
I requested an academic copy of the Revolution R Enterprise today. Since they only have Windows or Redhat version, I installed it on an old spare machine dual-booting Linux and Windows. I like what they did with the IDE (the debugger is nice) and ... [Read more...]
According to landroni, R can also be easily integrated with Geany:Download 0.19 (non-optional), then you need to set send_selection_unsafe=true in geany.conf and assign "Send selection to terminal" to a ctrl+r keybinding or similar. As long as you... [Read more...]
I am running a mixture logit model on a sample of 10,000 cases. The model takes a long time (3-4 hours) to run, which makes it nearly impossible to try different specifications and such. It probably makes sense to use a subsample for exploratory analys...
In my attempt to estimate a mixture regression model for the long-term impact of prenatal famine exposure, I compared results obtained from NPMLREG and Flexmix using the "rainfall" data set came with the package "forward" and was able to achieve identi...
This is another R package for latent variable analysis. I hope it can soon include multilevel and mixture modeling capacity. Another promising package is OpenMx, which is also under active development and very promising. At this moment, OpenMx seems to...
MCMCglmm does not produce more than one chain. So the way to check convergence is to run the model twice, each with different starting values, then conduct the diagnostic with:gelman.diag(m1$Sol, m2$Sol)