Als Titel muss ich so etwas zeichnen:
Kann ggplot oder andere Pakete, wenn ggplot nicht in der Lage ist, verwendet werden, um so etwas zu zeichnen?
r
data-visualization
ggplot2
funnel-plot
lokheart
quelle
quelle
stat_quantile()
, bedingte Quantile auf ein Streudiagramm zu setzen. Mit dem Formelparameter können Sie dann die Funktionsform der Quantilregression steuern. Ich würde Dinge wie Formel = vorschlageny~ns(x,4)
, um eine glatte Passform zu erhalten.Antworten:
Obwohl es Raum für Verbesserungen gibt, ist hier ein kleiner Versuch mit simulierten (heteroskedastischen) Daten:
quelle
Wenn Sie nach dieser Art von Trichterdiagramm (Metaanalyse) suchen , ist möglicherweise Folgendes der Ausgangspunkt:
quelle
linetype=2
Arguments in denaes()
Klammern - das Zeichnen der 99% -Zeilen - führt zu einem Fehler "Kontinuierliche Variable kann nicht dem Linientyp zugeordnet werden" mit dem aktuellen ggplot2 (0.9.3.1). Änderunggeom_line(aes(x = number.seq, y = number.ll999, linetype = 2), data = dfCI)
zugeom_line(aes(x = number.seq, y = number.ll999), linetype = 2, data = dfCI)
Werken für mich. Fühlen Sie sich frei, die ursprüngliche Antwort zu ändern und diese zu verlieren.Siehe auch das Cran-Paket berryFunctions, das einen Trichter-Plot für Proportionen ohne Verwendung von ggplot2 enthält, falls jemand ihn in Basisgrafiken benötigt. http://cran.r-project.org/web/packages/berryFunctions/index.html
Es gibt auch das Paket extfunnel, das ich nicht angeschaut habe.
quelle
Der Code von Bernd Weiss ist sehr hilfreich. Ich habe unten einige Änderungen vorgenommen, um einige Funktionen zu ändern / hinzuzufügen:
geom_segment
anstellegeom_line
der Linie verwendet, die das metaanalytische Mittel abgrenzt, sodass es dieselbe Höhe wie die Linien hat, die die 95% - und 99% -Konfidenzbereiche abgrenzenMein Code verwendet als Beispiel einen metaanalytischen Mittelwert von 0,0892 (se = 0,0035), aber Sie können Ihre eigenen Werte ersetzen.
quelle