rtrf <- function(x, theta = seq(-4, 4, length = 100)){ out <- list(theta = theta) out$p <- apply(rbind(rirf(x, theta)$p), 1, sum) class(out) <- "rtrf" return(out) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.