Scribd
Upload
5 views3 pages

Support Vector Machine R程式練習2

The document provides an R programming exercise focused on Support Vector Machines (SVM) using the e1071 package. It includes code for generating synthetic data, training an SVM model with both radial and linear kernels, and evaluating the model's performance on training and test datasets. The results include summaries of the model parameters and confusion matrices for predictions.

Uploaded by

hubertkuo418
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
5 views3 pages

Support Vector Machine R程式練習2

The document provides an R programming exercise focused on Support Vector Machines (SVM) using the e1071 package. It includes code for generating synthetic data, training an SVM model with both radial and linear kernels, and evaluating the model's performance on training and test datasets. The results include summaries of the model parameters and confusion matrices for predictions.

Uploaded by

hubertkuo418
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Support Vector Machine : R 程式練習 2

411210002 郭玉皓

2024-12-24
library(e1071)

## Warning: 套件 'e1071' 是用 R 版本 4.4.2 來建造的

set.seed(1)
x = matrix(rnorm(50*2), ncol=2)
y = c(rep(0, 25), rep(1, 25))
x[y==0, 2] = x[y==0, 2]+2
dat = data.frame(x=x, y=as.factor(y))
par(mfrow=c(1, 1))
plot(x, col=(y+1))

svmfit = svm(y~., data=dat, kernel="radial", cost=10, gamma=1)


plot(svmfit, dat)
library(ISLR)

names(Khan)

## [1] "xtrain" "xtest" "ytrain" "ytest"

dim(Khan$xtrain)

## [1] 63 2308

dim(Khan$xtest)

## [1] 20 2308

length(Khan$ytrain)

## [1] 63

length(Khan$ytest)

## [1] 20

table(Khan$ytrain)

##
## 1 2 3 4
## 8 23 12 20

table(Khan$ytest)
##
## 1 2 3 4
## 3 6 6 5

dat = data.frame(x=Khan$xtrain, y=as.factor(Khan$ytrain))


out = svm(y~., data=dat, kernel="linear", cost=10)
summary(out)

##
## Call:
## svm(formula = y ~ ., data = dat, kernel = "linear", cost = 10)
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: linear
## cost: 10
##
## Number of Support Vectors: 58
##
## ( 20 20 11 7 )
##
##
## Number of Classes: 4
##
## Levels:
## 1 2 3 4

table(out$fitted, dat$y)

##
## 1 2 3 4
## 1 8 0 0 0
## 2 0 23 0 0
## 3 0 0 12 0
## 4 0 0 0 20

dat.te = data.frame(x=Khan$xtest, y=as.factor(Khan$ytest))


pred.te = predict(out, newdata=dat.te)
table(pred.te, dat.te$y)

##
## pred.te 1 2 3 4
## 1 3 0 0 0
## 2 0 6 2 0
## 3 0 0 4 0
## 4 0 0 0 5

You might also like