如何在kernlab包的ksvm中自定义一个内核函数？答案

【问题标题】：How to customize a kernel function in ksvm of kernlab package?如何在kernlab包的ksvm中自定义一个内核函数？
【发布时间】：2013-06-08 00:13:36
【问题描述】：

我有经纬度，所以需要将RBF内核重新定义为exp(-1/2||sophere distrance||^2)，也就是说我需要自己重写一个内核函数。我写我的内核如下：

round.kernel <- function(x,y){
  sigma <- 1
  #R <- 6371
  R <- 1
  a <- (sin( (x[1]-y[1])/2 ))^2+cos(x[1])*cos(y[1])*(sin((x[2]-y[2])/2))^2
  c <- 2*atan2(sqrt(a),sqrt(1-a))
  d <- R*c
  res <- exp(-d^2/(2*sigma))
  return (res)
}
class(round.kernel) <- "kernel"

我测试了函数，内核应该是正确的。但是使用以下训练命令我得到了错误：

fit <- ksvm(y=train[,2],x=train[,3:4],kernel=round.kernel,type='eps-svr')

Error in .local(x, ...) : 
  List interface supports only the stringdot kernel.

比较棘手的是，我尝试了ksvm文档中的示例代码：

k <- function(x,y) {(sum(x*y) +1)*exp(-0.001*sum((x-y)^2))}
class(k) <- "kernel"

但我遇到了同样的错误。

有人知道如何正确定义核函数吗？

【问题讨论】：

这真的很有帮助。你能解释一下你是如何推导出核函数的吗？我看不出 exp(-1/2||sphere distnace||^2) 与您如何完成 res 之间的联系。
@momeara，我认为它基本上是基于三角形函数的，所以可以有多种方法使用arctan、arcsin，或者像atan2。检查此链接：movable-type.co.uk/scripts/latlong.html
这个例子对我有用。

标签： r kernel svm kernlab

【解决方案1】：

我的问题解决如下：内核代码是正确的，我应该直接定义一个函数（x，y），并将其类声明为“内核”。问题是即使在文档中 ksvm 确实支持 x,y 样式，但它们实际上不起作用。将其更改为公式数据样式最终可以运行：

fit <- ksvm(Freq~lat+lon,data=train[,2:4],kernel=roundrbf,type='eps-svr')

此外，我还阅读了 rbfdot 的源代码，以及在 kernlab 本身中定义的其他内核。注意他们的代码风格是这样的：

function(params){
  val <- function(x,y){
  # True kernel defined here
  }
  return (new ("kernel_name",.Data=val,kpar=list(params)))
}

但是说真的，我试过了，用这种风格制作内核函数是行不通的。工作方式直接就是这种风格：

k <- function(x,y){
  #calculate the result
}
class(k) <- "kernel"

【讨论】：