Plotten empirische kumulative Wahrscheinlichkeitsfunktion und deren inverse

Question

Ich habe Daten cdecn:

set.seed(0) 
cdecn <- sample(1:10,570,replace=TRUE) 
a <- rnorm(cdecn,mean(cdecn),sd(cdecn)) 

Ich habe ein Diagramm erstellt, das die kumulativen Wahrscheinlichkeiten anzeigt.

aprob <- ecdf(a) 
plot(aprob) 

Ich frage mich, wie ich die x-Achse und y-Achse umschalten kann ein neues Grundstück zu erhalten, das heißt, die Umkehrung der ECDF.

Gibt es auch für die neue Grafik eine Möglichkeit, eine vertikale Linie hinzuzufügen, wo die meine Kurve 0 schneidet?

Answer 1

Wir können Folgendes tun. Meine Kommentare entlang des Codes sind sehr erklärend.

## reproducible example 
set.seed(0) 
cdecn <- sample(1:10,570,replace=TRUE) 
a <- rnorm(cdecn,mean(cdecn),sd(cdecn)) ## random samples 

a <- sort(a) ## sort samples in ascending order 
e_cdf <- ecdf(a) ## ecdf function 
e_cdf_val <- 1:length(a)/length(a) ## the same as: e_cdf_val <- e_cdf(a) 

par(mfrow = c(1,2)) 

## ordinary ecdf plot 
plot(a, e_cdf_val, type = "s", xlab = "ordered samples", ylab = "ECDF", 
    main = "ECDF") 

## switch axises to get 'inverse' ECDF 
plot(e_cdf_val, a, type = "s", xlab = "ECDF", ylab = "ordered sample", 
    main = "'inverse' ECDF") 

## where the curve intersects 0 
p <- e_cdf(0) 
## [1] 0.01578947 

## highlight the intersection point 
points(p, 0, pch = 20, col = "red") 

## add a dotted red vertical line through intersection 
abline(v = p, lty = 3, col = "red") 

## display value p to the right of the intersection point 
## round up to 4 digits 
text(p, 0, pos = 4, labels = round(p, 4), col = "red") 

Answer 2

cdecn <- sample(1:10,570,replace=TRUE) 
a <- rnorm(cdecn,mean(cdecn),sd(cdecn)) 

aprob <- ecdf(a) 
plot(aprob) 

# Switch the x and y axes 
x <- seq(0,1,0.001754386) 
plot(y=knots(aprob), x=x, ylab = "Fn(y)") 

# Add a 45 degree straight line at 0, 0 
my_line <- function(x,y,...){ 
    points(x,y,...) 
    segments(min(x), y==0, max(x), max(y),...) 
} 
lines(my_line(x=x, y = knots(aprob))) 

Answer 3

Die "Gerade bei x == 0" Bit macht mich vermuten, dass Sie eine qqplot wollen:

qqnorm(a) 
qqline(a)