SHARE
TWEET

Untitled

nzcoops May 22nd, 2012 198 Never
  1. coef.lme <- function(mod){
  2.     res <- data.frame(
  3.         "Beta.CI" = paste(round(summary(mod)$coefficients$fixed, 3), " (",round(summary(mod)$coefficients$fixed-1.96*sqrt(diag(mod$varFix)), 2), ",", round(summary(mod)$coefficients$fixed+1.96*sqrt(diag(mod$varFix)), 2),")", sep=""),
  4.         "P.value" = round(2 * pt(-abs(summary(mod)$coefficients$fixed/sqrt(diag(mod$varFix))), summary(mod)$fixDF[[1]]), 3)
  5.     )
  6.     return(res)
  7. }
RAW Paste Data
Top