Kysymys:
Kysymyksiä siitä, miten satunnaiset vaikutukset määritetään lmerissä
alwin hoff
2012-07-03 18:58:55 UTC
view on stackexchange narkive permalink

Mittasin äskettäin, kuinka uuden sanan merkitys saavutetaan toistuvien altistusten aikana (käytäntö: päivät 1 - päivä 10) mittaamalla ERP: t (EEG), kun sanaa tarkasteltiin eri tilanteissa. Hallitsin myös kontekstin ominaisuuksia, esimerkiksi sen hyödyllisyyttä uuden sanan merkityksen löytämisessä (korkea vs. matala). Olen erityisen kiinnostunut käytännön vaikutuksista (päivät). Koska yksittäiset ERP-tallenteet ovat meluisia, ERP-komponenttien arvot saadaan keskiarvolla tietyn tilan kokeiden aikana. Toiminnolla lmer käytin seuraavaa kaavaa:

  lmer (ERPindex ~ käytäntö * konteksti + (1 | osallistujat), data = base)  

ja

  lmer (ERPindex ~ käytäntö * konteksti + (1 + käytäntö | osallistujat), data = base) 

I Olemme myös nähneet vastaavanlaisia ​​seuraavista satunnaisista vaikutuksista kirjallisuudessa:

  lmer (ERPindex ~ käytäntö * konteksti + (harjoittelu | osallistujat) + (käytäntö | osallistujat: konteksti), data = pohja ) 

Mitä saavutetaan käyttämällä satunnaiskerrointa muodossa osallistuja: konteksti ? Onko olemassa hyvää lähdettä, joka antaisi matriisi-algebran pelkkää tietoa tuntevan henkilön ymmärtää tarkalleen, mitä satunnaiskertoimet tekevät lineaarisissa sekamalleissa ja miten ne tulisi valita?

Kolme vastused:
Macro
2012-07-04 19:26:28 UTC
view on stackexchange narkive permalink

Kuvaan, mikä malli kukin lmer () -puhelusi sopii ja miten ne eroavat, ja vastaan ​​sitten viimeiseen kysymykseesi satunnaisefektien valitsemisesta.

Kummassakin kolmessa mallissasi on kiinteät tehosteet käytännölle , context ja näiden kahden vuorovaikutukselle. Satunnaiset vaikutukset vaihtelevat mallien välillä.

  lmer (ERPindex ~ käytäntö * konteksti + (1 | osallistujat), data = pohja) 

sisältää satunnaisen sieppauksen, jonka jakavat henkilöt, joilla on sama arvo osallistujille . Toisin sanoen jokaisen osallistujan regressioviivaa siirretään ylös / alas satunnaisella summalla, jonka keskiarvo on 0 $.

  lmer (ERPindex ~ käytäntö * konteksti + (1 + harjoittelu | osallistujat), data = base) 

Tämä malli satunnaisen sieppauksen lisäksi , sisältää myös satunnaisen kaltevuuden käytännössä . Tämä tarkoittaa sitä, että yksilöiden oppimisnopeus on erilainen. Jos yksilöllä on positiivinen satunnainen vaikutus, niin se kasvaa keskimääräistä nopeammin harjoituksen aikana, kun taas negatiivinen satunnainen vaikutus osoittaa, että hän oppii keskimäärin vähemmän nopeasti harjoituksen aikana tai mahdollisesti huononee harjoittelun aikana satunnaisen varianssin mukaan vaikutus (tämä olettaa, että käytännön kiinteä vaikutus on positiivinen).

  lmer (ERPindex ~ käytäntö * konteksti + (harjoittelu | osallistujat) + (harjoittelu | osallistujat: konteksti), data = pohja)  

Tämä malli sopii satunnaiselle kaltevuudelle ja sieppaukselle käytännössä (sinun on tehtävä (käytäntö-1 | ...) estämään sieppaus), aivan kuten edellinen malli teki, mutta nyt olet myös lisännyt satunnaisen kaltevuuden ja siepannut tekijä osallistuja: konteksti , joka on uusi tekijä, jonka tasot ovat jokainen yhdistelmä osallistujissa ja context ja vastaavat satunnaisvaikutukset jaetaan havainnoilla, joilla on sama arvo sekä osallistujilla että context : lla. Tämän mallin sovittamiseksi sinulla on oltava useita havaintoja, joilla on samat arvot sekä osallistujille että context : lle, muuten mallia ei voida arvioida. Monissa tilanteissa tämän vuorovaikutusmuuttujan luomat ryhmät ovat hyvin harvinaisia ​​ja johtavat erittäin meluisiin / vaikeasti sovitettaviin satunnaisvaikutelmamalleihin, joten haluat olla varovainen, kun käytät vuorovaikutustekijää ryhmittelymuuttujana.

korreloivat toistensa kanssa (kun taas yksilöitä, joita ei pidä korreloida) - satunnaiset vaikutukset toteuttavat tämän. Jos luulet havainnot, jotka jakavat sekä osallistujien että context tasot, ovat enemmän samanlaisia ​​kuin näiden kahden osan summa, "vuorovaikutuksen" satunnaisvaikutuksen sisällyttäminen voi olla tarkoituksenmukaista.

Muokkaa: Kuten @Henrik kommenteissa mainitsee, sopivat mallit, esim .:

  lmer (ERPindex ~ käytäntö * konteksti + ( 1 + harjoitus | osallistujat), data = base)  

tee niin, että satunnainen kaltevuus ja satunnainen leikkaus korreloivat keskenään ja että malli arvioi korrelaation. Rajoittaaksesi mallia siten, että satunnainen kaltevuus ja satunnainen leikkaus ovat korreloimattomia (ja siksi riippumattomia, koska ne ovat normaalisti jakautuneita), sovitat sen sijaan mallin: konteksti + (1 | osallistujat) + (käytäntö-1 | osallistujat), data = pohja)

Valinnan näiden kahden välillä tulisi perustua siihen, ajatteletko esimerkiksi osallistujilla s, joiden lähtötaso on keskimääräistä korkeampi (ts. positiivinen satunnainen sieppaus), todennäköisesti myös muutosnopeus on keskimääräistä suurempi (ts. positiivinen satunnainen kaltevuus). Jos näin on, sallit näiden kahden korreloinnin, jos ei, rajoitat heidän olevan riippumattomia. (Tässäkin esimerkissä oletetaan, että kiinteän vaikutuksen kaltevuus on positiivinen).

En halua olla nirso, mutta eikö toisessa mallissa ole myös korrelaatiota leikkauksen ja rinteiden välillä. Pelkästään rinteiden lisäämisen tulisi olla: `lmer (ERPindex ~ käytäntö * konteksti + (1 | osallistujat) + (0 + harjoittelu | osallistujat, data = pohja)` Vai olenko väärässä? (Ei yhteyttä: Anteeksi pieni muokkaukseni viestissäsi. Jos olet eri mieltä selvennyksestä, vaihda se takaisin)
@Henrik, kyllä, olet oikeassa, että se arvioi myös kahden satunnaisen vaikutuksen korrelaation. Kirjoittaessani tätä vastausta yritin antaa "kokonaiskuvan" käsityksen näiden mallien tapahtumista, johon ei sisältynyt mainintaa satunnaisvaikutusten välisestä korrelaatiosta, jolla ei ole yksinkertaista "kahden sentin" kuvausta. kaltevuus ja sieppaus tekevät :) Joka tapauksessa, en usko, että tämä puute tekee tulkinnasta, jonka tein vastauksessani väärin. BTW, kiitos muokkauksesta.
@Henrik, Lisäsin huomautuksen erosta satunnaisvaikutusten korreloinnin ja korreloinnin välillä, mikä mielestäni parantaa vastausta - kiitos, kun osoitit sen.
Kiitos. Yritän päästä mukaan sekoitettuun mallintamiseen ja kamppailen myös kysymyksen kanssa siitä, miten ja milloin mitä satunnaisefektirakennetta käytetään, jonka halusin vain varmistaa. Yhteenvetona, hyvä vastaus (+1).
On tärkeää huomata, että jos käytät satunnaisia ​​tehosteita, älä käytä niitä eristämään tietojoukossa esiintyviä epähomogeenisuuden taskuja, jotka ovat erottamattomia kiinnostuksistasi. Esimerkiksi ei pidä yrittää sovittaa satunnaisleikkausta aiheittain tekijälle, joka vaihtelee aiheittain.
@rpierce, satunnaisen leikkauksen varianssia, joka vastaa tekijää, joka eroaa kaikkien havaintojen välillä, ei tunnisteta jäännösvarianssista. Olet oikeassa, että sinun ei pitäisi yrittää tehdä tätä, mutta vaikka tekisitkin, malli ei olisi arvioitavissa.
@Macro, kiitos tästä upeasta vastauksesta. Haluaisin kiinnittää huomiota siihen, että lisäsit korreloimattomia satunnaisia ​​vaikutuksia. Ehdotettu malli: lmerillä (ERPindex ~ käytäntö * konteksti + (1 | osallistujat) + (käytäntö-1 | osallistujat), data = pohja) on itse asiassa enemmän parametreja kuin siinä, jolla on korreloiva leikkaus ja kaltevuus. Tämä näyttää tapahtuvan, koska korrelaation arvioinnin lisäksi sillä on nyt kaksi satunnaista sieppausta. Seuraava vaihtoehtoinen muotoilu on sama asia. lmer (ERPindex ~ käytäntö * konteksti + (1 | osallistujat) + (0 + käytäntö | osallistujat), data = pohja) Minulla ei valitettavasti ole tunnusta
@pom, kiitos kohteliaisuudesta. Re: kommenttisi, olen testannut tätä simuloiduilla tiedoilla ja mielestäni sinulla on se taaksepäin. Toisessa muokkaamassani mallissa on yksi parametri vähemmän kuin ensimmäisessä. Tämä johtuu siitä, että toinen malli rajoittaa korrelaation kahden satunnaisen vaikutuksen välillä nollaksi. Muut kuin mallit ovat samat. En ole varma, mitä kohtaat, mutta toistettava esimerkki auttaisi. Tässä on minun: `x <-normi (1000); id <- rep (1: 100, kukin = 10); y <- normaali (1000); g <- lmer (y ~ (1 + x | id)); g2 <- lmer (y ~ (1 | id) + (x-1 | id)); attr (logLik (g), "df"); attr (logLik (g2), "df"); "
gung - Reinstate Monica
2013-08-11 23:53:13 UTC
view on stackexchange narkive permalink

@Macro on antanut täällä hyvän vastauksen, haluan lisätä vain yhden pienen pisteen. Jos jotkut tilanteestasi käyttävät:

  lmer (ERPindex ~ käytäntö * konteksti + (harjoittele | osallistujat) + (harjoittele | osallistujat: konteksti), data = pohja)  

Epäilen heidän tekevän virheen. Harkitse: (käytäntö | osallistujat) tarkoittaa, että jokaiselle osallistujalle on satunnainen kaltevuus (ja leikkaus) käytännön vaikutukselle, kun taas (käytäntö | osallistujat: asiayhteys) tarkoittaa, että jokaiselle -käyttäjälle kontekstin mukaan yhdistelmä . Tämä on hieno, jos he haluavat sitä, mutta epäilen, että he haluavat (käytäntö: context | osallistujia) , mikä tarkoittaa, että vuorovaikutusvaikutukselle käytännöstä kontekstin mukaan jokaiselle osallistujalle .

Michael R. Chernick
2012-07-03 23:24:20 UTC
view on stackexchange narkive permalink

Satunnaisten vaikutusten tai sekoitettujen vaikutusten mallissa satunnaista vaikutusta käytetään, kun haluat käsitellä havaittua vaikutusta ikään kuin se johtuisi jonkinlaisesta vaikutusten todennäköisyysjakaumasta.

Yksi parhaista Esimerkkejä voin antaa, kun mallinnetaan kliinisten tutkimusten tietoja monikeskisestä kliinisestä tutkimuksesta. Sivustovaikutus mallinnetaan usein satunnaisena vaikutuksena. Tämä tapahtuu, koska tutkimuksessa tosiasiallisesti käytetyt noin 20 sivustoa otettiin paljon suuremmasta potentiaalisten sivustojen ryhmästä. Käytännössä valinta ei ehkä ole ollut sattumanvarainen, mutta silti voi olla hyödyllistä käsitellä sitä ikään kuin se olisi.

Vaikka sivuston vaikutus olisi voitu mallintaa kiinteäksi vaikutukseksi, olisi vaikea yleistää tuloksia suuremmalle väestölle, jos emme ottaisi huomioon sitä tosiseikkaa, että eri valitun 20 paikkaa olisi erilainen. Käsittelemällä sitä satunnaisena vaikutuksena voimme ottaa sen huomioon tällä tavalla.

-1, koska tämä vastaus ei käsittele todellisia kysymyksiä tässä.


Tämä Q & A käännettiin automaattisesti englanniksi.Alkuperäinen sisältö on saatavilla stackexchange-palvelussa, jota kiitämme cc by-sa 3.0-lisenssistä, jolla sitä jaetaan.
Loading...