diff --git a/GLiMMPS_functions.R b/GLiMMPS_functions.R
index 97f4192..8bd227f 100644
--- a/GLiMMPS_functions.R
+++ b/GLiMMPS_functions.R
@@ -93,7 +93,7 @@ glmmWald.sQTL <- function(data) {
   if ( !( inherits(testglm,"try-error")  )) { # only when it converges
     betas <-fixef(testglm)
     if (!is.na(fixef(testglm)[2])) {
-      snp.pval <-  summary(testglm)@coefs[2,4]
+      snp.pval <-  summary(testglm)$coefficients[2,4]
     }
   }
   #psi.fitted <- results$fitted.value