“Como extrair valor p de LM em R” Respostas de código

Como extrair valor p de LM em R

lmp <- function (modelobject) {
    if (class(modelobject) != "lm") stop("Not an object of class 'lm' ")
    f <- summary(modelobject)$fstatistic
    p <- pf(f[1],f[2],f[3],lower.tail=F)
    attributes(p) <- NULL
    return(p)
}

> lmp(modelobject)
[1] 1.622665e-05