hinzufügen x=y-Linie zu scatterplot -

Ich bin mit der scatterplot-Funktion aus dem car Paket erzeugen Sie ein Streudiagramm.
Ich möchte in der Lage zum erzeugen einer Referenz-Linie in der Handlung werden sollte, dass x=y ist.
Ich habe versucht, mit abline und es tut fügen Sie eine Zeile, aber es ist nicht der x=y-Linie. Kann mir da jemand helfen?

Mein code ist wie folgt:

scatterplot(phenos$P1~pheno$P0, data=pheno,spread=FALSE,ylab="6 month timepoint", xlab="Baseline Timepoint", jitter=list(x=1, y=1))
abline(0,1)

Dank.

Ich nehme an, Sie erstellen Sie einen boxplot anstelle eines Punktdiagramms. Richtig?
Was bekommen Sie? Welche Art von Daten P1 und P0?
abline(0,1) fügt eine Linie Steigung=1 und intercept=0, so sollte es Sie geben, was Sie brauchen. Wenn Sie ein paar Daten zu Plotten (so können wir sehen, was Sie sehen, wenn Sie führen Sie das Streudiagramm und abline Befehle), können wir mehr helfen.
zu bekommen ist x=y , lines(pheno$P0,pheno$P0)
car::scatterplot macht viel fancy stuff intern, einschließlich der Verwendung der layout() Funktion, wird es viel schwieriger zu überlagern, etwas einfaches wie eine 1:1-Linie ...

InformationsquelleAutor CadisEtRama | 2013-12-06

line plot r scatter-plot

1

Dies ist tatsächlich ziemlich schwierig/hackish, weil scatterplot() intern verwendet layout, das macht es schwer, die Kontrolle der Nebenhandlung, die momentan von der Grafik-Treiber. (Update: es ist schwieriger als ich dachte -- Einstellung par("mfg") muss gearbeitet haben, mehr oder weniger durch Zufall.)

Make-up-Daten (update: die Verwendung von Daten mit dem Mittelwert x und y nicht gleich null, und nicht einander gleich, als es illustrieren die Schwierigkeiten, die naiv mit abline() mehr klar)
```
set.seed(1)
d <- data.frame(x=rnorm(10,mean=10),y=rnorm(10,mean=12))
library(car)
```
Versuchen meine alte Strategie (die ja eigentlich nicht arbeiten, oder arbeiten nur unvorhersehbar):
```
scatterplot(y~x,data=d,reset.par=FALSE)
k <- 1              
for (i in 1:2) {
   for (j in 1:2) {
      par(mfg=c(i,j,2,2))
        abline(0,1,lwd=3,col=k)
        k <- k+1
  }
```
}

Je nachdem, wie ich das mache, entweder ich bekomme Warnungen und Fehler oder falsche Antworten. Es scheint egal, ob ich die scatterplot() Aufruf innerhalb einer Funktion ... ??

Zweiten Versuch, mehr konservativ: rekonstruieren Sie das layout von Grund auf neu.
```
 scatterplot(y~x,data=d)
 uu <- par("usr")
 ## mimic layout frolm car:::scatterplot.default.  Would be different if we were drawing only one
 ## of x-boxes or y-boxes
 layout(matrix(c(1, 0, 3, 2), 2, 2), widths = c(5, 95), 
        heights = c(95, 5))
 oldmar <- par(mar=rep(0,4))  ## zero out margins so we can plot in sub-boxes without errors
 ## now skip through the first two sub-plots
 par(new=TRUE); plot.new(); par(new=TRUE); plot.new()
 par(oldmar)  ## reset margins
 ## blank plot with user limits set and 'interior' axis calculations
 plot(0:1,0:1,xlab="",ylab="",xlim=uu[1:2],ylim=uu[3:4],xaxs="i",yaxs="i")
 ## add annotation
 abline(a=0,b=1,col=4,lwd=3)              
```
Angesichts der Höhe der Aufwand und Fragilität dieser Lösung, es könnte in der Tat sein bestes, um hack scatterplot Optional ermöglichen abline() festgelegt werden zusätzlich, oder Fragen Sie den maintainer, für die Fähigkeit ...
- Ben, ich fand, dass ich nicht brauchen, um einen der parameter-Einstellung zu bekommen y=x Linie in den plot richtig. Ich habe gerade scatterplot(y~x,data=d) und dann abline(0,1) - und der y=x Linie wurde korrekt aufgetragen. Tun Sie das gleiche Verhalten?
- Nachtrag: ich habe versucht, die gleiche Sache mit verschiedenen random-Daten-sets und festgestellt, dass manchmal abline gibt eine Zeile mit y-Achsenabschnitt = 0 und manchmal ist der y-Achsenabschnitt ist verschoben sehr leicht unter null (aber ich bekomme immer die richtige Steigung=1). Also ich stehe korrigiert.
- Ich denke es wird davon abhängen, ob der plot setup kommt vor zu entsprechen, etwas, wo (0,0) liegt in der Nähe der Mitte des Diagramms. Versuchen Sie es mit Wege von (0,0) und sehen, was passiert ...
- wenn ich versuche die 'abline(0,1)', bekomme ich eine Linie, die ein wenig unter null. Es ist fast, "es" aber nicht wirklich. Ben, dass der code nicht funktioniert für mich, ich habe immer Fehler, ich fand es schwer zu Folgen, wie Sie dies getan haben, sorry! Ich hauptsächlich verwenden wollte 'scatterplot()' für die boxplots an den Achsen. Ich kann haben, zu versuchen, wieder auf Basis-Graphiken in R erzeugen die x=y-Linie. Haben Sie Jungs wissen, wie man hinzufügen box-plots auf die Achsen mit base Grafik dann?
- geben Sie uns eine reproduzierbar Beispiel von etwas, das nicht funktioniert, bitte.
- wie vorhergesagt, die y-interecept bewegt sich leicht, je nachdem mittels der gezeichneten Punkte.
- das ist schön (ich nehme an, du meinst, wenn Sie nur mit abline() ohne die fancy stuff in meiner Antwort), aber ... ich bin nicht überrascht. Sie müssen um genauer zu sein über das, was ", der code hat nicht funktioniert für mich, ich habe immer Fehler" bedeutet-und wenn möglich, zeigen Sie uns genau, indem Sie gemeinsam eine reproduzierbare Antwort -- siehe tinyurl.com/reproducible-000
- Ich eigentlich nur am Ende die Schaffung der plot in ggplot2 und es funktionierte. danke aber.
InformationsquelleAutor Ben Bolker
5

Können Sie Ihr Ziel erreichen, indem Sie die gleiche Funktion abline. Verwenden Sie die Funktion abline(a=0, b=1) wo a und b sind das abfangen und slop von der Linie jeweils. Diese zieht man eine gerade Linie von der form Y = a + b*x.

Hoffe, das hilft.

InformationsquelleAutor Alemnew Sheferaw Asrese
3

Benutze ich nur die Segmente-Funktion:
segments(x0=0,y0=0,x1=45,y1=45)

Machen nur die y-Werte gehen nur über die Grenzen von Ihrem Grundstück.

Diese Funktion ist auch nett, da es nicht ändern Sie Ihre x-und y-Grenzen, so kann die Anlage sauber läuft, in Ihrem Diagramm, ohne präzise die x-und y-Werte.
- Dies funktionierte mit plot, nicht sicher über scatterplot
InformationsquelleAutor Denys Kelly

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.