5

我正在尝试使用用户定义的内核。我知道 kernlab 在 R 中提供用户定义的内核（自定义内核函数）。我使用了数据垃圾邮件，包括包 kernlab。（变量数=57 示例数=4061）

我定义了内核的形式，

kp=function(d,e){

as=v*d
bs=v*e
cs=as-bs
cs=as.matrix(cs)

exp(-(norm(cs,"F")^2)/2)
}

class(kp)="kernel"

它是高斯核的变换核，其中v是与每个变量的标准差向量相反的连续变化值，例如：

v=(0.1666667,........0.1666667)

训练集定义了 60% 的垃圾邮件数据（保留不同类别的比例）。

如果数据的类型是垃圾邮件，那么对于火车 svm，数据的类型 = 1

m=ksvm(xtrain,ytrain,type="C-svc",kernel=kp,C=10)

但这一步不起作用。它总是在等待回应。

所以，我问你这个问题，为什么？是不是因为例子太多了？是否有任何其他 R 包可以为用户定义的内核训练 SVM？

4

1 回答 1

3

首先，你的内核看起来像一个经典的 RBF 内核，带有v = 1/sigma，那么你为什么要使用它呢？您可以使用内置的 RBF 内核并简单地设置sigma参数。特别是 - 您可以在矢量化矩阵上使用经典欧几里得，而不是在矩阵上使用 frobenius 范数。

第二 - 这工作得很好。

> xtrain = as.matrix( c(1,2,3,4) )
> ytrain = as.factor( c(0,0,1,1) )
> v= 0.01
> m=ksvm(xtrain,ytrain,type="C-svc",kernel=kp,C=10)
> m
Support Vector Machine object of class "ksvm" 

SV type: C-svc  (classification) 
 parameter : cost C = 10 


Number of Support Vectors : 4 

Objective Function Value : -39.952 
Training error : 0

至少有两个原因让你还在等待结果：

  • RBF 内核引发了最难优化 SVM（尤其是大型C）的问题
  • 用户定义内核的效率远低于内置内核

由于我不确定，是否ksvm实际上优化了用户定义的内核计算（实际上我很确定它没有），您可以尝试构建内核矩阵（训练向量K[i,j] = K(x_i,x_j)在哪里x_i）并提供它。你可以通过i'thksvm

K <- kernelMatrix(kp,xtrain)
m <- ksvm(K,ytrain,type="C-svc",kernel='matrix',C=10)

预先计算内核矩阵可能是一个相当长的过程，但优化本身会快得多，所以如果你想测试许多不同的C值（你肯定应该这样做），这是一个好方法。不幸的是，这需要O(n^2)内存，因此如果您使用超过 100 000 个向量，您将需要非常大量的 RAM。

于 2013-09-01T05:10:21.293 回答