我正在嘗試使用R型提供程序來擬合和預測支持向量機模型。我能夠擬合模型,但是當我嘗試預測返回的矢量與訓練矢量具有相同的長度時,它不應該有這個長度。爲什麼R.predict.svm返回錯誤大小的列表?
我在R中直接嘗試了等效代碼,並且返回的列表具有正確的長度。
這是怎麼發生的?
下面是一個例子:
open System
open RDotNet
open RProvider
open RProvider.stats
open RProvider.e1071
// Random number generator
let rng = Random()
let rand() = rng.NextDouble()
// Generate fake X1 and X2
let X1s = [ for i in 0 .. 9 -> 10. * rand() ] // length = 10
let X2s = [ for i in 0 .. 9 -> 5. * rand() ] // length = 10
let Z1s = [ for i in 0 .. 5 -> 10. * rand() ] // length = 6
let Z2s = [ for i in 0 .. 5 -> 5. * rand() ] // length = 6
// Build Ys
let Ys = [0;1;0;1;0;1;0;1;0;1]
let XMat =
["X1", box X1s; "X2", box X2s]
|> namedParams
|> R.cbind
let ZMat =
["Z1", box Z1s; "Z2", box Z2s]
|> namedParams
|> R.cbind
let svm_model =
["x", box XMat; "y", box Ys ; "type", box "C"; "gamma", box 1.0]
|> namedParams
|> R.svm
let svm_predict = R.predict(svm_model, ZMat)
let res =
if svm_predict.Type = RDotNet.Internals.SymbolicExpressionType.IntegerVector then
svm_predict.AsInteger()
|> List.ofSeq
else failwithf "Expecting a Numeric but got a %A" svm_predict.Type
printfn "The predicted values are: %A" res
// The predicted values are: [1; 2; 1; 2; 1; 2; 1; 1; 1; 2]
這裏是最初的R代碼:
library(stats)
library(e1071)
// Random number generator
x1 <- 10 * rnorm(10)
x2 <- 5 * rnorm(10)
x = cbind(x1, x2)
z1 <- 10 * rnorm(5)
z2 <- 5 * rnorm(5)
z = cbind(z1, z2)
zs <- c(0,1,0,1,0,1,0,1,0,1)
svm_fit = svm(x=x,y=zs,type="C",gamma=1.0)
svm_pred = predict(svm_fit, z)
print(svm_pred)
1 2 3 4 5
1 0 1 1 1
Levels: 0 1
你的原始R代碼是什麼樣的? –
我添加了R代碼。感謝您的關注。 – Soldalma
我發現'let svm_predict = R.predict(svm_model,ZMat)'中的第二個參數是不相關的。無論我輸入什麼,結果都是一樣的。即使我忽略了這個論點。 – Soldalma