Dosis-Wirkungs-Probit-Analyse<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Hallo Gruppe,
Ich versuche, eine einfache Dosis-Wirkungs-Analyse zu reproduzieren, die ich in SPSS durchgeführt habe. In JMP verwende ich also das Modul „Modell anpassen“ und wähle für verallgemeinerte lineare Modelle, Binomialverteilung und Probit-Link aus.
Ich habe drei Variablen: Anzahl der durchgeführten Tests (10 Tests), Anzahl der Tests mit positivem Ergebnis und Konzentration (von 0 bis 36).
0 10 2
1 10 4
2 10 9
6 10 16
8 10 25
9 10 36
0 10 0
Ich wähle die Anzahl der durchgeführten Tests als „Y“-Variable, die Anzahl der Tests mit positivem Ergebnis als „Häufigkeit“ und die Konzentration in Protokollen als Variable aus, um die Auswirkungen zu modellieren. Aber wenn ich auf „Ausführen“ drücke, wird in den Ergebnissen diese Meldung angezeigt:
„Achtung: Hessisch nicht eindeutig. Die Singularität ist wahrscheinlich auf eine schlechte Anpassung des Modells oder auf lineare Abhängigkeiten zwischen den Modellkovariaten zurückzuführen. Passform und Ergebnis sind von fragwürdigem Wert“
Und die Analyse wird nicht durchgeführt.
Kann jemand einige Kommentare abgeben?
Danke,
Dave
Betreff: Dosis-Wirkungs-Probit-AnalyseDas ist in der Tat die richtige Reihenfolge, um die Variablen einzugeben. Danke für den Hinweis, Michael.
Dave
Betreff: Dosis-Wirkungs-Probit-AnalyseDie 0 10 0 letzte Zeile der Originaldaten verhindert aus irgendeinem Grund, dass die Analyse die richtige Reihenfolge benötigt. Wenn Sie es weglassen, ist die Reihenfolge wichtig.
Michael
Betreff: Dosis-Wirkungs-Probit-AnalyseVerfolge diesen Thread mit Interesse.
Aus den JMP-Dokumenten:
Wenn Ihre Daten als Häufigkeiten von Ereignissen und Versuchen zusammengefasst sind, geben Sie zwei fortlaufende Werte an
Säulen in dieser Reihenfolge : eine Zählung der Anzahl der Erfolge und eine Zählung der Anzahl von
Versuche. Alternativ können Sie statt der Erfolge auch die Anzahl der Fehlschläge angeben.
(meine Betonung).
Wenn Sie also die Reihenfolge ändern, in der die beiden y-Variablen eingegeben werden, werden Sie meiner Meinung nach mehr von dem sehen, was Sie erwartet haben, insbesondere wenn Sie die x-Achse so ändern, dass sie logarithmisch ist.
Allerdings kann ich nicht erklären, warum die Reihenfolge bei Ihren Originaldaten keine Rolle zu spielen scheint.
Michael
Betreff: Dosis-Wirkungs-Probit-AnalyseLieber DanO,
Es tut mir leid, dass ich Sie noch einmal störe, aber wenn ich die Zahlen auf die Demo umstelle, erhalte ich scheinbar falsche Ergebnisse. Bitte schauen Sie sich diesen Datensatz an, mit denselben Variablen, aber unterschiedlichen Zahlen. Für die drei höchsten Konzentrationen erhalte ich 100 % positive Ergebnisse
Anzahl positiver Ergebnisse Anzahl durchgeführter Tests Konzentration
0 1140 0,01
241 1140 0,1
590 1140 0,5
963 1140 1
1134 1140 5
1140 1140 10
1140 1140 50
1140 1140 100
Ich führe den Probit so aus, wie Sie es vorgeschlagen haben

Nach dem Ausführen erhalte ich ein seltsames Ergebnis, das für mich keinen Sinn ergibt. Anstatt prozentuale Korrekturen am Y vorzunehmen, scheint es etwas anderes zu tun. Was mache ich falsch?

Vielen Dank für Ihre Hilfe
Dave
Betreff: Dosis-Wirkungs-Probit-AnalyseVielen Dank, DanO, das ist eine sehr klare Antwort. Es ist interessant, dass JMP weiß, dass sie miteinander verbunden sind.
Was die Protokolltransformation betrifft, ja, ich werde das berücksichtigen (sie wurde tatsächlich in meinem ursprünglichen Beitrag erwähnt, obwohl ich zugeben muss, dass sie ziemlich vage ist!)
Am besten,
Dave
Betreff: Dosis-Wirkungs-Probit-AnalyseDave,
Ihre „Anzahl der Tests mit positivem Ergebnis“ muss ebenfalls als „Y“-Variable eingegeben werden. Die Rolle „Häufigkeit“ gibt an, wie oft eine Zeile in der Datentabelle tatsächlich vorkam. In diesem Fall erwartet das GLM-Modell die Anzahl der durchgeführten Tests und die Anzahl der positiven Ergebnisse als „separate“ Y-Werte. Die Analyse wird sie miteinander verknüpfen.

Eine Sache, die Sie nicht erwähnt haben und die bei dieser Art von Arbeit häufig durchgeführt wird, ist die Protokollierung der Konzentration. Ich habe dies im oben gezeigten Modelldialog getan. Markieren Sie den Konzentrationsmodelleffekt, klicken Sie auf die Schaltfläche „Transformieren“ und wählen Sie „Protokoll“. Die Ergebnisse werden hier veröffentlicht.

Dan