icc.lme<-function(x) {
  require(nlme) # NOTE: REQUIRES version 3.1-27 or higher
  vcx<-getVarCov(x)
  d<-diag(vcx)
  d/(d+x$sigma^2)
}