R-ohjelmointi.org

Tilastotieteellistä ohjelmointia R-kielellä

lokakuu, 2011

Logistisen regression selitysaste

Tein erääseen julkaisuun logistisia regressiomalleja, ja käsikirjoituksen vertaisarvioija halusi nähdä mallien max rescaled R squared [sic] -arvon. Onneksi on Google, joka tiesi kertoa, että se on SAS:n logistiselle regressiolle ilmoittama pseudo-selitysaste, mutta sen laskeminen R:ssä olikin toinen juttu. Laajennuspaketista pscl löytyi nopeasti komento pR2(), joka laskee erilaisia pseudo-selitysasteita: n <- 100 x <- c(rnorm(n), 1+rnorm(n)) […]

Posted in Aiheeton | No Comments »

lmer() ja p-arvot

R:ään on saatavilla useampiakin laajennuspaketteja, joilla voidaan sovittaa aineistoon lineaarisia sekamalleja. Eräs ensimmäisiä paketteja oli nlme, josta on sittemmin jatkokehitetty paketti lme4. Monia käyttäjiä harmittava piirre on, ettei lme4-paketin funktio lmer() anna mallin parametrien estimaateille p-arvoja. Syynä on, ettei paketin kehittäjä Douglas Bates luota p-arvojen laskemiseksi tarvittavien vapausasteiden nykyisiin määrittämismenetelmiin, eikä siten ole halunnut tuoda […]

Posted in Aiheeton | No Comments »

Category