SORU
7 NİSAN 2011, PERŞEMBE


p-değeri çekin ve r-kare bir doğrusal regresyon

Nasıl p-değeri (tek açıklayıcı değişken sıfır olmayan katsayısının önemi) çekin ve R-Kare regresyon modelinden değeri musunuz? Örnek...

x = cumsum(c(0, runif(100, -1,  1)))
y = cumsum(c(0, runif(100, -1,  1)))
fit = lm(y ~ x)
summary(fit)

summary(fit) biliyorumgörüntüler-p ve R-Kare değeri, ancak diğer değişkenler bu sokmak mümkün olmak istiyorum.

CEVAP
7 NİSAN 2011, PERŞEMBE


R-Kare doğrudan özeti nesneden değer summary(fit)$r.squared iade edebilirsiniz. Doğrudan özü size tüm öğeleri listesi için names(summary(fit)) bkz.

This blog post p-değeri: dönmek için bir işlev özetliyor

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(fit)
[1] 1.622665e-05

Alternatif olarak, Özet nesnesine benzer bir şekilde anova(fit) nesneden p-değeri yukarıda tut.

Bunu Paylaş:
  • Google+
  • E-Posta
Etiketler:
R

YORUMLAR

SPONSOR VİDEO

Rastgele Yazarlar

  • Feel The Electricity!

    Feel The Ele

    20 ŞUBAT 2010
  • National Geographic

    National Geo

    7 Mayıs 2006
  • Vintendo Power

    Vintendo Pow

    2 Ocak 2007