Nothing
SLR <- function(data, numer = NA, denom = NA, weight = TRUE)
{
data <- as.matrix(data/apply(data, 1, sum))
if (!weight[1])
weights <- rep(1/ncol(data), ncol(data))
if (weight[1])
weights <- apply(data, 2, mean)
if (length(weight) == ncol(data)) {
if (sum(weight <= 0) > 0)
stop("Error: some weights zero or negative")
if (sum(weight) != 1)
print("Sum of weights not exactly 1, but are rescaled")
weights <- weights/sum(weights)
}
if (is.na(numer[1]) | is.na(denom[1]))
stop("Numerator and/or denominator not specified")
if (length(intersect(denom, numer)) > 0)
stop("Error: numerator and denominator intersect")
if (length(numer) == 1)
num <- data[, numer]
if (length(numer) > 1)
num <- apply(data[, numer], 1, sum)
if(sum(num==0)>0)
stop("Error: some numerators are zero")
if (length(denom) == 1)
den <- data[, denom]
if (length(denom) > 1)
den <- apply(data[, denom], 1, sum)
if(sum(den==0)>0)
stop("Error: some denominators are zero")
num.wt <- sum(weights[numer])
den.wt <- sum(weights[denom])
amalglr <- log(num/den)
amalglr.weight <- num.wt * den.wt
names(amalglr.weight) <- paste(paste(colnames(data)[numer],
collapse = "&"), paste(colnames(data)[denom], collapse = "&"),
sep = "/")
return(list(LR = amalglr, LR.wt = as.numeric(amalglr.weight)))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.