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

  • alex maybury

    alex maybury

    20 Aralık 2007
  • Tek Syndicate

    Tek Syndicat

    23 Temmuz 2008
  • Truc Minh

    Truc Minh

    23 Ocak 2011