"f2fpds" <-
function(f, rate=NA) {
if(is.na(rate)) return(f)
if(! check.fs(f)) return(FALSE)
if(rate <= 0) {
warning("invalid rate ? <= 0 parameter")
return(rep(NA,length(f)))
}
g <- (log(f) + rate) / rate
g[g < 0] <- NA
g[g > 1] <- NA
return(g)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.