Histogramme und Ausreißer

Kreditrisikomodellierung in R

Lore Dirick

Manager of Data Science Curriculum at Flatiron School

Funktion hist() verwenden

hist(loan_data$int_rate)

Histogram of interest rate

Funktion hist() verwenden

hist(loan_data$int_rate, main = "Histogram of interest rate", xlab = "Interest rate")

Histogram of interest rate

hist() für annual_inc verwenden

hist(loan_data$annual_inc, xlab = "Annual Income", main = "Histogram of Annual Income")

Screen Shot 2020-06-12 at 1.55.03 PM.png

hist() für annual_inc verwenden

hist_income <- hist(loan_data$annual_inc,
                    xlab = "Annual Income",
                    main = "Histogram of Annual Income")
hist_income$breaks

0  500000 1000000 1500000 2000000 2500000 3000000 3500000 4000000 4500000 ...

Das breaks-Argument

n_breaks <- sqrt(nrow(loan_data)) # n_breaks = 170.5638
hist_income_n <- hist(loan_data$annual_inc, breaks = n_breaks, 
                      xlab = "Annual Income", main = "Histogram of Annual Income")

Screen Shot 2020-06-12 at 1.55.58 PM.png

annual_inc

plot(loan_data$annual_inc, ylab = "Annual Income")

Screen Shot 2020-06-12 at 1.56.24 PM.png

annual_inc

plot(loan_data$annual_inc, ylab = "Annual Income")

Screen Shot 2020-06-12 at 1.56.53 PM.png

Ausreißer

Wann ist ein Wert ein Ausreißer?
- Expertenurteil
- Faustregel, z. B.
  - Q1 - 1,5 * IQR
  - Q3 + 1,5 * IQR
- Meist: Kombination aus beidem

Expertenurteil

„Annual salaries > $3 million are outliers"

# Find outlier
index_outlier_expert <- which(loan_data$annual_inc > 3000000)

# Remove outlier from dataset
loan_data_expert <- loan_data[-index_outlier_expert, ]

Faustregel

Ausreißer, wenn größer als Q3 + 1,5 * IQR

# Calculate Q3 + 1.5 * IQR
outlier_cutoff <- quantile(loan_data$annual_inc, 0.75) + 1.5 * IQR(loan_data$annual_inc)

# Identify outliers
index_outlier_ROT <- which(loan_data$annual_inc > outlier_cutoff)

# Remove outliers
loan_data_ROT <- loan_data[-index_outlier_ROT, ]

hist(loan_data_expert$annual_inc,
     sqrt(nrow(loan_data_expert)), 
     xlab = "Annual income")

hist(loan_data_ROT$annual_inc,
     sqrt(nrow(loan_data_ROT)), 
     xlab = "Annual income")

Bivariates Diagramm

plot(loan_data$emp_length, loan_data$annual_inc, 
     xlab= "Employment length", ylab= "Annual income")

Screen Shot 2020-06-12 at 1.58.14 PM.png

Bivariates Diagramm

plot(loan_data$emp_length, loan_data$annual_inc, 
     xlab= "Employment length", ylab= "Annual income")

Screen Shot 2020-06-12 at 1.58.34 PM.png

Lass uns üben!

Kreditrisikomodellierung in R