center <- function(x,grouping.var) { mean.table <- summarize(x,by=llist(as.factor(grouping.var)),FUN=mean,na.rm=T,stat.name="Group.Mean") y <- x - mean.table$Group.Mean[as.factor(grouping.var)] return(as.vector(y))}