Einzelner Wert auf der X-Achse anzeigen

[FranziskaBönte]

Liebe Kommunity,

leider bin ich ein R-Anfänger.

Ich male ein Bildchen mit plot(), bei dem ich nur wenige Werte anzeigen möchte. So auch beta auf der X-Achse (unter der roten Geraden). Aber ich bekomme das nicht passend hin. Könnt Ihr mir helfen? Es geht um die rot eingefärbte Zeile.

Hier mein Code:

plot(function(x)dnorm(x ,mean=1,sd=2) ,-10 ,10 ,main=" Konsistenz von Schätzfunktionen " ,xlab="",
ylab= "" , col="black",xaxt="n", ,yaxt="n" )
plot(function(x)dnorm(x ,mean=-2,sd=4) ,-10 ,10 ,add=TRUE,col=" blue ",xaxt="n",yaxt="n" )
abline(v=4, col="red")
text(-5.5, 0.09, bquote(T==10), col="blue")
text(-2, 0.15, bquote(T==100), col="black")
text(5, 0.12, bquote(T==1.000), col="red")
text(4, 0, bquote(beta), col="red")
title(ylab=bquote("f(" ~hat(beta)~")"), line=0, cex.lab=1.2)
title(xlab=bquote(hat(beta)), line=1, cex.lab=1.2)

Danke für Eure Hilfe, Franziska

[bigben]

Hallo Franziska,

bei mir schreibt die rote Zeile einen roten Buchstaben beta ins Diagramm:

Rplot17.png (5.33 KiB) 844 mal betrachtet

Was ist bei Dir das Problem?

LG,
Bernhard

[FranziskaBönte]

Mein Problem ist, dass ich das rote beta in dem Bildchen nicht mitten in der Linie stehen haben möchte, sondern da, wo quasi der zugehörige x-Achsenwert stehen würde.

[bigben]

Achsenbeschriftungen macht man mit axis(), wobei ich auf Anhieb gerade nicht weiß, wie man das rot kriegt:

Code: Alles auswählen

plot(function(x)dnorm(x ,mean=1,sd=2) ,-10 ,10 ,main=" Konsistenz von Schätzfunktionen " ,xlab="",
     ylab= "" , col="black",xaxt="n", ,yaxt="n" )
plot(function(x)dnorm(x ,mean=-2,sd=4) ,-10 ,10 ,add=TRUE,col=" blue ",xaxt="n",yaxt="n" )
abline(v=4, col="red")
text(-5.5, 0.09, bquote(T==10), col="blue")
text(-2, 0.15, bquote(T==100), col="black")
text(5, 0.12, bquote(T==1.000), col="red")
#text(4, 0, bquote(beta), col="red")
title(ylab=bquote("f(" ~hat(beta)~")"), line=0, cex.lab=1.2)
title(xlab=bquote(hat(beta)), line=1, cex.lab=1.2)
axis(1, at = 4, labels = expression(beta), col = "red")

LG,
Bernhard

[bigben]

Nachtrag: Google hat geholfen:

Code: Alles auswählen

axis(1, at = 4, labels = expression(beta), col = "red", col.axis = "red")

[FranziskaBönte]

super! Das war es, was ich suchte. Perfekt. 1.000 Danke