/usr/lib/R/site-library/Biobase/doc/esApply.R is in r-bioc-biobase 2.14.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 | ### R code from vignette source 'esApply.Rnw'
###################################################
### code chunk number 1: R.hide
###################################################
library(Biobase)
data(sample.ExpressionSet)
###################################################
### code chunk number 2: R
###################################################
print(sample.ExpressionSet)
print(exprs(sample.ExpressionSet)[1,])
print(pData(sample.ExpressionSet)[1:2,1:3])
###################################################
### code chunk number 3: R
###################################################
print(rbind(exprs(sample.ExpressionSet[1,]),
sex <- t(pData(sample.ExpressionSet))[1,]))
###################################################
### code chunk number 4: R
###################################################
medContr <- function( y, x ) {
ys <- split(y,x)
median(ys[[1]]) - median(ys[[2]])
}
###################################################
### code chunk number 5: R
###################################################
print(apply(exprs(sample.ExpressionSet[1,,drop=F]), 1,
medContr, pData(sample.ExpressionSet)[["sex"]]))
###################################################
### code chunk number 6: R
###################################################
medContr1 <- function(y) {
ys <- split(y,sex)
median(ys[[1]]) - median(ys[[2]])
}
print(esApply( sample.ExpressionSet, 1, medContr1)[1])
###################################################
### code chunk number 7: esApply.Rnw:126-127
###################################################
sessionInfo()
|