fpr <- function(t) {
  f <- t[1,2]/(t[1,2]+t[2,2])
return (f)
}