accuracy <- 
function(data=NULL, y=model$y, model, t=0.500, rate=TRUE) {
  pr.value <- model$fitted.value
  predixns <- pr.value>t
  reality  <- model$y
  acc <- sum( predixns==reality ) 
  if(rate) { acc <- acc / length(reality) }
return(acc)
}