##' Quantifies linear trends in pfam or fxn with independant variable
##'
##' @export
find_trend <- function(d, indep.var="depth", param="frac.reads", disc="pfam") {
fits <- ddply(d, disc, lm_stats, indep.var, param)
fm <- merge(d, fits, by="pfam")
fm
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.