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