R/summ.R

`summ` <-
function (x, y) 
{
   y1 <- as.numeric(y)[order(x)]
    m.y <- tapply(y1, as.factor(sort(x)), mean)
    n <- tapply(y1, as.factor(sort(x)), length)
    sd <- tapply(y1, as.factor(sort(x)), sd)
    std.err <- sd/sqrt(length(y1))
    n.p <- table(x)
    unx <- sort(unique(x))
    y.is.u <- pava(y = m.y, w = n.p)
    y.is.d <- pava(y = m.y, w = n.p, decreasing = TRUE)


sumary <- data.frame(unx ,n, m.y ,  std.err ,y.is.u , y.is.d)
names(sumary ) <- c(" Dose", " N", "Mean  ", "S.E    ", "Isotonic Mean Up","Isotonic Mean Down")
rownames(sumary) <- NULL
return(sumary)
}

Try the IsoGeneGUI package in your browser

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

IsoGeneGUI documentation built on Nov. 1, 2018, 3:57 a.m.