GGHF {biology} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
GGHF(object)
object |
~~Describe object 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 section{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 (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 }