Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Hallo alle,
Ich wäre sehr dankbar für Ratschläge zur Analyse einer Reihe von Anteilsdaten, in die ich sowohl zufällige als auch feste Faktoren einbeziehen möchte.
Die Daten stammen aus einem Experiment mit Fruchtfliegen. Jede der 3 "sozialen" Behandlungen wurde auf 3 unabhängige Fliegenpopulationen angewendet, so dass insgesamt 9 Populationen vorhanden sind. In jeder der 9 Populationen wurde dann eine Einzel- oder Hochdosisbehandlung für einzelne Fliegen durchgeführt. Jede dieser Fliegen wurde in eine separate Phiole gestellt und wetteiferte um Nachkommen. Nachkommen, die von diesen männlichen Versuchstieren gezeugt wurden, hatten rote Augen, während Nachkommen, die nicht von ihnen gezeugt wurden, weiße Augen hatten, wobei die Gesamtzahl der Nachkommen von roten und weißen Augen unter den Phiolen variierte.
Ich möchte wissen, ob der Anteil der von diesen Fliegen gezeugten Nachkommen je nach sozialer Behandlung, Diätbehandlung oder deren Interaktion unterschiedlich ist. Ich habe ein Modell mit den festen Faktoren der Ernährung (mit zwei Ebenen), der sozialen Behandlung (drei Ebenen) und der Ernährung durch Interaktion mit der Behandlung erstellt. Das Modell hat auch den zufälligen Effekt einer Bevölkerung, die innerhalb der sozialen Behandlung verschachtelt ist. Ich habe als Antwortvariable den Anteil der Nachkommen mit roten Augen in einer Durchstechflasche, gewichtet mit der Gesamtzahl der Nachkommen in der Durchstechflasche.
Beim Ausführen des Modells waren die Residuen nicht normal verteilt. Ich habe eine Arcussin-Quadrat-Wurzel-Transformation angewendet, was geholfen hat, aber die Residuen blieben nicht normal verteilt. Eine Option ist das Anwenden der Logit-Transformation, wodurch die Residuen normal verteilt werden.
Ich frage mich jedoch, ob es ein besserer Ansatz wäre, ein verallgemeinertes lineares Modell mit einer bestimmten, den Daten entsprechenden Fehlerverteilung zu verwenden. Ich kann mir vorstellen, dass dies eine Logit-Link-Funktion sein könnte. Ich würde mich sehr über Ratschläge freuen. Insbesondere darüber, ob dies besser ist als die Verwendung einer Logit-Transformation innerhalb eines linearen Modells, wie eine Fehlerverteilung für ein GLM ausgewählt wird und wie das Modell mit JMP eingerichtet wird.
Danke vielmals!
Jen
Betreff: Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Ja, JMP kann kein Modell für eine binomiale Antwort mit einem zufälligen Effekt anpassen.
Betreff: Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Danke dafür. Ich habe die Population als Zufallseffekt betrachtet, weil ich nicht an Schätzungen über die einzelnen Populationen interessiert bin. Wenn ich das Experiment erneut durchführe, würde ich neue Populationen von Fliegen aufbauen. Würden Sie etwas anderes vorschlagen?
Sieht auch so aus, als habe ich Ihren Vorschlag zur Einrichtung des Modells falsch verstanden. Ich habe eine Spalte für rote und weiße Werte und eine weitere Spalte mit Zählungen eingerichtet, so dass jeder Mann zwei Reihen hat, eine mit roten und eine andere mit weißen. Mit diesem Setup muss ich die zwei Zeilen verknüpfen, die jedem Mann zugeordnet sind und somit im Modell enthalten sind. Gibt es einen anderen Weg, um es einzurichten, ohne den Einzelnen einzuschließen?
Ist es generell so, dass ein zufälliger Faktor mit dieser Art von Modell nicht funktioniert?
Betreff: Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Ich sehe keine zufälligen Auswirkungen in Ihrer Studie. Sex ist ein fester Effekt. Bevölkerung ist ein fester Effekt. Die einzelnen Fächer sind die Residuen und müssen nicht als Zufallseffekt explizit eingegeben werden.
Betreff: Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Vielen Dank für diesen Vorschlag.
Ich habe dies so eingestellt, dass es zwei Reihen für die Anzahl der Nachkommen mit roten und weißen Augen gibt.
Beim Versuch, das Modell auszuführen, stoße ich auf die Schwierigkeit, dass der GLM keinen zufälligen Effekt zulässt.
Ich habe die zwei zufälligen Auswirkungen von einzelnen Männern und der Bevölkerung. Wenn ich sie in fixe Effekte umwandle, verliert das Modell die Freiheitsgrade, um die fixen Effekte zu testen, die mich interessieren.
Irgendwelche Ideen, wie man damit umgehen soll?
Danke noch einmal.
Betreff: Wie richte ich ein gemischtes Modell für Anteilsdaten mit entsprechender Fehlerverteilung ein?Sie können die Binomialverteilung für das Modell mit der Funktion logit link verwenden. Anstelle von Proportionen können Sie jedoch eine Antwortspalte mit verwenden rot und Weiß Werte und eine weitere Spalte mit einer Zählung für alle Bedingungen. Ordnen Sie die Spalte mit den Zählungen der Freq-Analyserolle zu.
