Blog Archives

A quetion

June 21, 2010
By
A quetion

I posted a question to the R-Mixed-Model mailing list but have not had any responses yet:--------------------------------------------------------------------------------------- Dear All,I wonder if anybody has tried to make glmmBUGS work with JAGS. Mya...

Read more »

R-INLA package

June 19, 2010
By
R-INLA package

Another R package for mixed effect modeling. Looks promising.

Read more »

Sweave with Emacs and ESS, problem solved!

May 9, 2010
By
Sweave with Emacs and ESS, problem solved!

With help from Seb, the nagging problem I had when trying to sweaving with Emacs and ESS (http://old.nabble.com/ESS-and-Sweave-td28339734.html) has been solved.First of all, apply the patch for ess-swv.el, provided by Seb. Second, set the default ...

Read more »

Sweave vs. pgfSweave

April 28, 2010
By
Sweave vs. pgfSweave

I think pgfSweave should replace the default of Sweave system because it brings significant improvements over the original Sweave. I certain have made it my default Sweave tool on my system.------------------------------- EDITED ON MAY 1, 2010 --------...

Read more »

Trouble with ESS and Sweave

April 23, 2010
By
Trouble with ESS and Sweave

Last time I tried to sweave a document from with Emacs+ESS, I was using an earlier version of ESS (the current version is 5.8), and things seemed to be fine. Today when I tried to sweave a simple document and produced PDF output, I got error message of...

Read more »

GLMM using DPpackage

April 9, 2010
By
GLMM using DPpackage

I was able to fit a semi-parametric Bayesian GLMM model using DPpackage. It took me many hours to sample from the posterior distribution (DPM prior):MCMC scan 1000 of 5000 (CPU time: 18950.080 s)MCMC scan 2000 of 5000 (CPU time: 22510.100 s)M...

Read more »

Measuring the length of time to run a function

March 17, 2010
By
Measuring the length of time to run a function

This post describes how to time the run time of a R function.

Read more »

GLMM revisted

March 5, 2010
By
GLMM revisted

A short while ago, I reported some discrepancies between the results produced by "lme4" and other R packages as well as Stata. Today I upgraded to the most recent version of "lme4a" and re-ran my model. The error of false convergence disappea...

Read more »

An interesting paper

February 27, 2010
By
An interesting paper

Ben Bolker has an interesting paper (outline of a paper) comparing different approaches to estimate GLMM in R environment, which is very helpful to what I am doing right now.The paper pointed out the following options to fit GLMM using R:glmerglmmMLglm...

Read more »

Generalized linear mixed effect model problem

February 16, 2010
By
Generalized linear mixed effect model problem

I am trying to compare cohort difference in infant mortality using generalized linear mixed model. I first estimated the model in Stata:xi:xtlogit inftmort i.cohort, i(code)which converged nicely:Fitting comparison model:Iteration 0:   log likelih...

Read more »