R: Mit "segment" - Befehl zu setzen, vertikale Linien zwischen den gezeichneten Punkte

Habe ich eine kleine dataframe, von dem bin ich Plotten 3 Spalten, um ein Risiko abzuschätzen, und die 95% Konfidenzintervalle. Jetzt habe ich diese 3 Sätze von Vektoren dargestellt als Punkte, aber ich möchte Sie verbinden mit "segment".

Hier ist ein Beispiel der dataframe gezeichnet:

                      Diagnosis  age.group X..change X..lower X..upper
1                    Dysrythmia All adults        16        0       35
2                 Heart failure All adults        -4      -20       14
3                        Asthma All adults        10       -5       28

Hier ist mein plot-code:

plot(dt[,4], pch="-", ylim=c(-20, 50), axes=F, ann=F, cex=1.5)
abline(h=0, col=1, lty=2)
points(dt[,3], pch=16, col="black", bg="black" )
points(dt[,5], pch="-", cex=1.5)
axis(1, at=1:10, lab=dt[,1], las=3, lwd=0, cex.axis=0.7, pos=-22)
axis(2, at=5*-20:54, las=1, cex.axis=0.7, cex.lab=0.7, col=1)
title(main="Risk estimates: All Adults", col.main="black", font.main=1)
title(ylab="Increase in risk (%)", col.lab=rgb(0,0.5,0))
box()

Punkte sind die Schätzungen und die Striche sind die Konfidenzintervalle. Ich will schließen Sie diese drei Punkte für jede Diagnose. Ich habe die R-notation, aber es hilft mir nicht herauszufinden, wie man zu sagen R, die xy "Daten" ich will zeichnen Sie die Segmente anschließen, denn ich habe Vektoren verwendet hier anstelle der Werte? Kann mir jemand helfen, schreiben Sie ein segment ein Zeile code? Danke

Was über das hinzufügen von Zeilen(dt[,3], pch=16, col="black" bg="black")? Oder muss es "segment"?
Sie sollten mit GGPLOT. Es ermöglicht Ihnen das erstellen von CI-Linien/error bars ganz einfach: R Kochbuch/ggplot2

InformationsquelleAutor mEvans | 2012-05-09

  1. 1

    OK, so dass Sie im Grunde zeichnen möchten GUS-Staaten. Die Pfeile Befehl ist wahrscheinlich besser, für das, was Sie wollen. Hier ist ein kurzes Beispiel, arbeitet mit dem, was Sie haben, für den unteren CI.

    plot(dt[,3], pch="•", ylim=c(-20, 50), axes=F, ann=F, cex=1.5, bty = 'o')
abline(h=0, col=1, lty=2)
arrows(1:3, dt[,3], 1:3, dt[,4], angle = 90, length = 0.08)

    Können Sie sich daher verlassen sich die extra-Punkte Befehle und neue Segmente Befehle. Es ist alles viel übersichtlicher.

    Aber, wenn Sie darauf bestehen, hinzufügen von Segmenten zu dem, was Sie haben, es ist nur...

    segments(1:3, dt[,3], 1:3, dt[,5])
segments(1:3, dt[,3], 1:3, dt[,4])
    Wenn John ' s segments Ansatz ok ist, warum dann nicht for (k in 1:nrows(dt)) lines(dt[k,3:5]) ?

    InformationsquelleAutor John

Schreibe einen Kommentar