Ist das ein echter Effekt? Übereinstimmende Paare - Bivariate - RegressionHallo,
bitte helfen sie einem neuling hier - ist das ein echter effekt.
Ich teste eine neue Methode zum Testen (Mittelwert 10) einer neuen Methode nach einem Goldstandard (paco2). Ich untersuche auch die Auswirkungen einer anderen Variablen, um ein verbessertes Modell zu testen.
Ich habe mit einer Bland-Altman-Analyse (gemischte Paare) begonnen, und meine Methoden zeigen eine gute Korrelation von 0,69 (kleines CI). Ich teste auch die unabhängige Variable auf meinen Goldstandard und selbst hier finde ich einen "Effekt".
Also begann ich eine bivariate Analyse für beide (Mittelwert 10) - nicht gezeigt - und für SaO2. Der Mittelwert 10 zeigt eine gute Korrelation zu paCO2, obwohl das r-Quadrat nur 0,48 beträgt. Auch meine andere Variable (SaO2) zeigt eine gewisse Wirkung in der bivariaten Analyse mit p <0,001.
Keine Analyse dieser beiden Effekte (Bild rechts) - beide scheinen einen signifikanten Effekt zu haben. Und das Modell verbesserte sich mit einem Quadrat von r von 0,67.
Aber ist das ein wirklicher Effekt? Oder mache ich meine Analyse hier falsch?
Alle Kommentare und Unterstützung sehr geschätzt! Vielen Dank, Mark
Re: Ist das ein echter Effekt? Übereinstimmende Paare - Bivariate - RegressionPeter,
vielen dank für den weisen rat. Es ist auch ein überraschendes physiologisches Ergebnis für mich, aber am Ende musste ich meinen Berechnungen vertrauen. Vielen Dank, Marc
Re: Ist das ein echter Effekt? Übereinstimmende Paare - Bivariate - RegressionKennzeichen,
Vielen Dank für Ihre umfassende Beratung! Das war sehr nützlich.
Ja, ich sah mir die Residuen an und bekam kein passendes Muster, sondern normal verteilt. (siehe Abbildung 1)
Ich habe mir auch unsere Datenbank genauer angesehen, um einige Echtzeitvariablen zu erstellen. Tatsächlich ist die Anpassung von paCO2 bei der Implementierung der zweiten Variablen besser. In einer grafischen Analyse ist die Vorhersageformel für paCO2 basierend auf Mean und SpO2 (ich habe SaO2 durch SpO2 ersetzt) perfekt und gibt mir einen perfekten Korrelationskoeffizienten g (Abbildung 2).
Kann ich die Zeilen speichern, die außerhalb der 90% -Dichte-Ellipse liegen?
Nochmals vielen Dank für Ihre Hilfe! Marc
fig 1
fig 2
Re: Ist das ein echter Effekt? Übereinstimmende Paare - Bivariate - RegressionSie fragen, "... ist die Wirkung wirklich?" Mein Kollege Mark Bailey hat kluge Ratschläge gegeben. Ich werde einfach auf seinen vorletzten Satz setzen, um empirisch zu bestätigen / zu validieren und ein wenig komplementär zu betrachten.
Was sagt Ihnen Ihr Wissen über den Prozess, dass eine Wirkung möglich ist? Einer der Punkte, die ich immer an diejenigen heranziehe, die in der Welt der statistischen Inferenz neu sind, ist, dass, wenn Ihre p-Werte, t-Werte, F-Verhältnisse und jede andere Folgerungsstatistik Ihnen sagen, dass Sie jetzt bewiesen haben 'und beachten Sie die Zitate, denn statistische Inferenzen beweisen wirklich nichts (aber das ist eine andere Geschichte für später), ein Effekt, der' Wasser läuft bergauf 'ähnelt. Na dann stimmt etwas nicht. Vielleicht lauerten Variablen im Experiment. Vielleicht gab es ein Messsystem-Bias. Vielleicht war da (Leerzeichen ausfüllen).
Wenn Ihr Wissen über Physik, Chemie, Biologie oder sozioökonomisches Verhalten sagt, dass der Effekt unwahrscheinlich oder unmöglich ist ... dann übertreffen in meinem Buch die p-Werte und ihre Brüder jeden Tag.
Re: Ist das ein echter Effekt? Übereinstimmende Paare - Bivariate - RegressionIch bin nicht sicher über den ersten gezeigten Vergleich, paCO2 gegen SaO2 , aber die Handlung zeigt, dass die Übereinstimmende Paare Analyse ist nicht angemessen . Warum? Bei dieser Analyse wird davon ausgegangen, dass der einzige Unterschied im Durchschnittsniveau liegt. Das heißt, es kann eine feste Differenz sein, aber keine proportionale Differenz. Die Grafik sollte ein horizontales Datenband zeigen. Stattdessen können Sie sehen, dass es von der oberen linken Ecke zur unteren rechten Ecke geneigt ist. Die Berechnungen sind korrekt, aber die Ergebnisse sind ohne Bedeutung.
Der zweite Vergleich ist Mittelwert 10 gegen paCO2 . Dieser Fall ist angemessen für Übereinstimmende Paare weil die Datenpunkte als horizontales Band erscheinen. Das Konfidenzintervall und der zweiseitige t-Test für den Unterschied sind bei alpha = 0,05 beide signifikant. Die Punktschätzung zeigt das an Mittelwert 10 ist um 9 Einheiten niedriger.
Der dritte Vergleich ändert den Fokus. Jetzt ist es paCO2 gegen SaO2 . Ich kann nur den oberen Teil des Berichts zur Analyse der Abweichung sehen. Es ist wie hoch F-Verhältnis ,> 500, was wahrscheinlich bei alpha = 0,05 signifikant ist. Ich kann den Parameter Estimates-Bericht nicht sehen. Diese beiden Variablen scheinen linear zu sein und ignorieren andere Variablen.
Der vierte Vergleich ist eine mehrfache Regression von paCO2 gegen Mittelwert 10 und SaO2 . Diese Analyse ist auch ein anderer Fokus. Vergleichst du? Mittelwert 10 zu paCO2 oder nicht? Die Analyse zeigt, dass eine enge Beziehung besteht paCO2 und beide Prädiktoren. Das Mittelwert 10 und SaO2 sind jedoch nicht ausreichend korreliert, um die Signifikanz zu beeinflussen. Mittelwert 10 hat eine Steigung von 0,9 (vorausgesetzt, dass der gesamte Fehler in ist Mittelwert 10 ). SaO2 hat eine Steigung von -0,3, was eine negative Vorspannung an anzeigt paCO2 .
Ist dieser Effekt echt? (Ich gehe davon aus, dass du meinst SaO2 Befangenheit, wenn Sie nach dem Effekt fragen.) Ich würde zunächst die Regressionsanalyse mit Diagnosefunktionen fortsetzen, damit Sie den Ergebnissen soweit vertrauen können. Haben Sie die Residuen überhaupt untersucht? Haben Sie den Varianzinflationsfaktor für die Parameterschätzungen geprüft? Schließlich müsste ein empirisches Modell empirisch bestätigt / verifiziert / validiert werden. Können Sie Proben mit den gewünschten Werten machen Mittelwert 10 und SaO2 zur Bestätigung der Modellvorhersage für paCO2 ?