mcmcpvalue {biology} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
mcmcpvalue(samp, fac)
samp |
~~Describe samp here~~ |
fac |
~~Describe fac here~~ |
~~ If necessary, more details than the description above ~~
~Describe the value returned If it is a LIST, use
comp1 |
Description of 'comp1' |
comp2 |
Description of 'comp2' |
~~further notes~~
~Make other sections like Warning with Warning .... ~
~~who you are~~
~put references to the literature/web site here ~
~~objects to See Also as help
, ~~~
##---- Should be DIRECTLY executable !! ---- ##-- ==> Define data, use random, ##-- or do help(data=index) for the standard data sets. ## The function is currently defined as function (samp, fac) { s <- grep(fac, colnames(samp)) samp <- as.matrix(samp)[, s] std <- backsolve(chol(var(samp)), cbind(0, t(samp)) - colMeans(samp), transpose = TRUE) sqdist <- colSums(std * std) sum(sqdist[-1] > sqdist[1])/nrow(samp) }