GGHF {biology}R Documentation

~~function to do ... ~~

Description

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

Usage

GGHF(object)

Arguments

object ~~Describe object 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 (object) 
{
    nterms <- length(object$terms)
    GG <- function(SSPE, P) {
        p <- nrow(SSPE)
        if (p < 2) 
            return(NA)
        lambda <- eigen(SSPE %*% solve(t(P) %*% P))$values
        lambda <- lambda[lambda > 0]
        ((sum(lambda)/p)^2)/(sum(lambda^2)/p)
    }
    HF <- function(gg, error.df, p) {
        ((error.df + 1) * p * gg - 2)/(p * (error.df - p * gg))
    }
    gg.hf <- NULL
    for (term in 2:nterms) {
        SSP <- object$SSP[[term]]
        SSPE <- object$SSPE[[term]]
        P <- object$P[[term]]
        p <- ncol(P)
        gg <- GG(SSPE, P)
        error.df <- object$error.df
        hf <- HF(gg, error.df, p)
        gg.hf <- rbind(gg.hf, cbind(gg, hf))
    }
    gg.hf
  }

[Package biology version 1.0 Index]