R/NormalizeByPatient.R

normalize.by.patient<-function(x,pat){
  f<-try(lme(x~1,random=~1|pat))
  out<-x;
  if (!inherits(f, "try-error")) {f<-lme(x~1,random=~1|pat); out<-x-f$fitted[,2]}
  return(out)
}
mssm-msf-2019/BiostatsALL documentation built on May 22, 2019, 12:16 p.m.