Blog Archives

Convergence diagnostics with MCMCglmm

July 5, 2010
By
Convergence diagnostics with MCMCglmm

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)

Read more »

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 »

Sponsors

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)