Defines functions summary.rmas.risk

Documented in summary.rmas.risk

summary.rmas.risk <-
function(object, q=c(0.025,0.975),...){
  cosa <- object
  cosa.boot <- NULL 

    # extract bootstraped probabilities
  for(i in 1:length(cosa$cf.boot)){
    cosa.boot <- cbind(cosa.boot,cosa$cf.boot[[i]][,2])
    # paste together Threshold, Probanbility (cosa$cf.obs) and 95% C.I.
  tabla <- cbind(cosa$cf.obs,t(apply(cosa.boot,1,quantile, q)))
  class(tabla) <- c("summary.rmas.risk", class(tabla))
  plot(tabla, main= cosa$main)

Try the Rramas package in your browser

Any scripts or data that you put into this service are public.

Rramas documentation built on May 2, 2019, 6:52 a.m.