mcmcpvalue {biology}R Documentation

~~function to do ... ~~

Description

~~ A concise (1-5 lines) description of what the function does. ~~

Usage

mcmcpvalue(samp, fac)

Arguments

samp ~~Describe samp here~~
fac ~~Describe fac here~~

Details

~~ If necessary, more details than the description above ~~

Value

~Describe the value returned If it is a LIST, use

comp1 Description of 'comp1'
comp2 Description of 'comp2'

...

Warning

....

Note

~~further notes~~

~Make other sections like Warning with section{Warning }{....} ~

Author(s)

~~who you are~~

References

~put references to the literature/web site here ~

See Also

~~objects to See Also as help, ~~~

Examples

##---- 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)
  }

[Package biology version 1.0 Index]