Drei Formen eines „gemischten Modells“ interpretieren

Es gibt eine Unterscheidung, die mich mit gemischten Modellen auslöst, und ich frage mich, ob ich Klarheit darüber bekommen könnte. Nehmen wir an, Sie haben ein gemischtes Modell der Zähldaten. Es gibt eine Variable, von der Sie wissen, dass Sie sie als festen Effekt haben möchten (A) und eine andere Variable für die Zeit (T), die beispielsweise als "Site" -Variable gruppiert ist.

So wie ich es verstehe:

glmer(counts ~ A + T, data=data, family="Poisson") ist ein Modell mit festen Effekten.

glmer(counts ~ (A + T | Site), data=data, family="Poisson") ist ein Zufallseffektmodell.

Meine Frage ist, wenn Sie etwas haben wie:

glmer(counts ~ A + T + (T | Site), data=data, family="Poisson")was ist T? Ist es ein zufälliger Effekt? Ein fester Effekt? Was wird eigentlich erreicht, wenn man T an beiden Stellen setzt?

Wann sollte etwas nur im Abschnitt für zufällige Effekte der Modellformel erscheinen?

Dies wird möglicherweise deutlicher, wenn Sie die Modellformel für jedes dieser drei Modelle aufschreiben. Sei die Beobachtung für die Person in Site in jedem Modell und definiere analog, um sich auf die Variablen in Ihrem Modell zu beziehen. i j A i j , T i $Y_{ij}$$i$$j$$A_{ij}, T_{ij}$

glmer(counts ~ A + T, data=data, family="Poisson") ist das Modell

Das ist nur ein gewöhnliches Poisson-Regressionsmodell.

glmer(counts ~ (A + T|Site), data=data, family="Poisson") ist das Modell

Dabei sind zufällige Effekte, die bei jeder Beobachtung von Personen von Standort . Diese zufälligen Effekte dürfen in dem von Ihnen angegebenen Modell frei korreliert werden (dh es werden keine Einschränkungen für ). Um Unabhängigkeit zu erzwingen, muss man sie in verschiedene Klammern setzen, zB würde es tun. Dieses Modell geht davon aus, dass ist für alle Seiten , aber jede Seite hat eine zufällige Offset ( ) und hat eine zufällige lineare Beziehung mit sowohl .j ≤ log ( E ( Y i j ) ) α 0 η j 0 A i j , T i $\eta_{j} = (\eta_{j0}, \eta_{j1}, \eta_{j2}) \sim N(0, \Sigma)$$j$$\Sigma$(A-1|Site) + (T-1|Site) + (1|Site)$\log \big( E(Y_{ij}) \big)$$\alpha_0$$\eta_{j0}$$A_{ij}, T_{ij}$

glmer(counts ~ A + T + (T|Site), data=data, family="Poisson") ist das Modell

Also hat eine "durchschnittliche" Beziehung zu , gegeben durch die festen Effekte aber diese Beziehung ist für jede Site unterschiedlich und diese Unterschiede werden durch die Zufallseffekte . Das heißt, die Basislinie ist zufällig verschoben, und die Steigungen der beiden Variablen sind zufällig verschoben, und alle Benutzer derselben Site haben dieselbe zufällige Verschiebung. A i j , T i j θ 0 , θ 1 , & thgr; 2 γ j 0 , γ j 1 , γ j $\log \big( E(Y_{ij}) \big)$$A_{ij}, T_{ij}$$\theta_0, \theta_1, \theta_2$$\gamma_{j0}, \gamma_{j1}, \gamma_{j2}$

was ist T? Ist es ein zufälliger Effekt? Ein fester Effekt? Was wird eigentlich erreicht, wenn man T an beiden Stellen setzt?

T γ j 1 T log ( E ( Y i j ) $T$ ist eine deiner Kovariaten. Es ist kein zufälliger Effekt - es Siteist ein zufälliger Effekt. Es gibt einen festen Effekt von , der je nach dem zufälligen Effekt, den - im obigen Modell verleiht, unterschiedlich ist . Durch Einbeziehen dieses Zufallseffekts wird eine Heterogenität zwischen Standorten in der Beziehung zwischen und .$T$Site$\gamma_{j1}$$T$$\log \big( E(Y_{ij}) \big)$

Wann sollte in der Modellformel nur etwas im Bereich für zufällige Effekte erscheinen?

Dies ist eine Frage dessen, was im Kontext der Anwendung sinnvoll ist.

In Bezug auf den Achsenabschnitt - Sie sollten den festen Achsenabschnitt aus vielen Gründen dort belassen (siehe z. B. hier ). Betreff: Der zufällige Abschnitt erster Linie eine Korrelation zwischen Beobachtungen, die an derselben Stelle gemacht wurden. Wenn eine solche Korrelation keinen Sinn ergibt, sollte der Zufallseffekt ausgeschlossen werden.$\gamma_{j0}$

In Bezug auf die zufälligen Steigungen spiegelt ein Modell mit nur zufälligen Steigungen und keinen festen Steigungen die Annahme wider, dass für jede Site eine Beziehung zwischen und Ihren Kovariaten für jede Site besteht Wenn Sie diese Effekte jedoch über alle Websites mitteln, besteht keine Beziehung. Wenn Sie beispielsweise eine zufällige Steigung in aber keine feste Steigung hatten, würde dies bedeuten, dass die Zeit im Durchschnitt keine Auswirkung hat (z. B. keine säkularen Trends in den Daten), aber jede in eine zufällige Richtung im Zeitverlauf weist. das könnte Sinn machen. Auch hier kommt es auf die Anwendung an.$\log \big( E(Y_{ij}) \big)$$T$Site

Beachten Sie, dass Sie das Modell mit und ohne zufällige Effekte anpassen können, um festzustellen, ob dies der Fall ist. Im festen Modell sollten keine, im nachfolgenden Modell jedoch signifikante zufällige Effekte zu sehen sein. Ich muss Sie darauf hinweisen, dass Entscheidungen wie diese häufig besser auf der Grundlage eines Verständnisses der Anwendung als durch Modellauswahl getroffen werden.

Sie sollten beachten, dass Tes sich bei keinem der Begriffe Ihres Modells um zufällige Effekte handelt, sondern um einen festen Effekt. Zufällige Effekte sind nur die Effekte, die nach dem |in einer lmerFormel erscheinen!

Eine ausführlichere Beschreibung dieser Spezifikation finden Sie in dieser früheren FAQ-Frage .

Aus diesen Fragen sollte Ihr Modell Folgendes ergeben (für Ihren festen Effekt T):

• Ein globaler Hang
• Ein Ausdruck mit zufälligen Steigungen, der die Abweichung von der Gesamtsteigung für jede Stufe von angibt Site
• Die Korrelation zwischen den zufälligen Steigungen.

Und wie von @ mark999 gesagt, ist dies tatsächlich eine übliche Spezifikation. Bei Entwürfen mit wiederholten Messungen möchten Sie im Allgemeinen zufällige Steigungen und Korrelationen für alle Faktoren mit wiederholten Messungen (innerhalb von Subjekten) haben.

Im folgenden Artikel finden Sie einige Beispiele (die ich hier immer wieder zitiere):

Judd, CM, Westfall, J. & amp; Kenny, DA (2012). Stimuli als Zufallsfaktor in der Sozialpsychologie behandeln: Eine neue und umfassende Lösung für ein allgegenwärtiges, aber weitgehend ignoriertes Problem. Zeitschrift für Persönlichkeits- und Sozialpsychologie , 103 (1), 54–69. doi: 10.1037 / a0028347

Etwas sollte nur im zufälligen Teil erscheinen, wenn Sie nicht besonders an seinem Parameter interessiert sind, sondern ihn einschließen müssen, um abhängige Daten zu vermeiden. Wenn beispielsweise Kinder in Klassen verschachtelt sind, möchten Sie Kinder normalerweise nur als zufälligen Effekt.