Generování radiálně separovatelného datasetu

Support Vector Machines v R

Kailash Awati

Instructor

Generování 2D rovnoměrně rozložené sady bodů

Generování datasetu s 200 body
- 2 prediktory x1 a x2, rovnoměrně rozdělené mezi -1 a 1.

# Set required number of datapoints
n <- 200
# Set seed to ensure reproducibility
set.seed(42)

# Generate dataframe with 2 predictors x1 and x2 in (-1, 1)
df <- data.frame(x1 = runif(n, min = -1, max = 1),
                 x2 = runif(n, min = -1, max = 1))

Vytvoření kruhové hranice

Vytvoření kruhové rozhodovací hranice s poloměrem 0,7 jednotky.
Kategorická proměnná y je +1 nebo -1 podle toho, zda bod leží vně nebo uvnitř hranice.

radius <- 0.7
radius_squared <- radius ^ 2

#categorize data points depending on location wrt boundary
df$y <- factor(ifelse(df$x1 ^ 2 + df$x2 ^ 2 < radius_squared, -1, 1),
               levels = c(-1, 1))

Vykreslení datasetu

Vizualizace pomocí ggplot.

library(ggplot2)

prediktory vyneseny na 2 osách; třídy rozlišeny barvou.

# Build plot
p <- ggplot(data = df, aes(x = x1, y = x2, color = y)) + 
     geom_point() + 
     scale_color_manual(values = c("-1" = "red", "1" = "blue")) 
# Display plot 
p

Kapitola 3.1 – radiálně separovatelný dataset

Přidání kruhové hranice – část 1

Vytvoříme funkci pro generování kružnice

# Function generates dataframe with points
# lying on a circle of radius r
circle <- 
  function(x1_center, x2_center, r, npoint = 100) {

  # Angular spacing of 2*pi/npoint between points
  theta <- seq(0, 2 * pi, length.out = npoint)
  x1_circ <- x1_center + r * cos(theta)
  x2_circ <- x2_center + r * sin(theta)

  data.frame(x1c = x1_circ, x2c = x2_circ)
}

Přidání kruhové hranice – část 2

Přidání hranice do grafu:
- vygenerovat hranici pomocí funkce circle().
- přidat hranici do grafu pomocí geom_path()

# Generate boundary
boundary <- circle(x1_center = 0,
                   x2_center = 0,
                   r = radius)
# Add boundary to previous plot
p <- p + 
     geom_path(data = boundary,
               aes(x = x1c, y = x2c),
               inherit.aes = FALSE)
# Display plot
p

Kapitola 3.1 – radiálně separovatelný dataset s rozhodovací hranicí

Čas na procvičení!

Support Vector Machines v R