On 2011-02-18 11:16, Juergen Rose wrote:> If do:
>> library("e1071")
>> example(svm)
>
> I get:
>
>
> svm> data(iris)
>
> svm> attach(iris)
>
> svm> ## classification mode
> svm> # default with factor response:
> svm> model<- svm(Species ~ ., data = iris)
>
> svm> # alternatively the traditional interface:
> svm> x<- subset(iris, select = -Species)
>
> svm> y<- Species
>
> svm> model<- svm(x, y)
>
> svm> print(model)
>
> Call:
> svm.default(x = x, y = y)
>
>
> Parameters:
> SVM-Type: C-classification
> SVM-Kernel: radial
> cost: 1
> gamma: 0.25
>
> Number of Support Vectors: 51
>
>
> svm> summary(model)
>
> Call:
> svm.default(x = x, y = y)
>
>
> Parameters:
> SVM-Type: C-classification
> SVM-Kernel: radial
> cost: 1
> gamma: 0.25
>
> Number of Support Vectors: 51
>
> ( 8 22 21 )
>
>
> Number of Classes: 3
>
> Levels:
> setosa versicolor virginica
>
>
>
>
> svm> # test with train data
> svm> pred<- predict(model, x)
>
> svm> # (same as:)
> svm> pred<- fitted(model)
>
> svm> # Check accuracy:
> svm> table(pred, y)
> y
> pred setosa versicolor virginica
> setosa 50 0 0
> versicolor 0 48 2
> virginica 0 2 48
>
> svm> # compute decision values and probabilities:
> svm> pred<- predict(model, x, decision.values = TRUE)
>
> svm> attr(pred, "decision.values")[1:4,]
> setosa/versicolor setosa/virginica versicolor/virginica
> 1 1.196152 1.091460 0.6705626
> 2 1.064621 1.056332 0.8479934
> 3 1.180842 1.074534 0.6436474
> 4 1.110699 1.053143 0.6778595
>
> svm> # visualize (classes by color, SV by crosses):
> svm> plot(cmdscale(dist(iris[,-5])),
> svm+ col = as.integer(iris[,5]),
> svm+ pch = c("o","+")[1:150 %in% model$index + 1])
>
> *** caught segfault ***
> address (nil), cause 'unknown'
>
> Traceback:
> 1: .Call("La_rs", x, only.values, PACKAGE = "base")
> 2: eigen(-x/2, symmetric = TRUE)
> 3: cmdscale(dist(iris[, -5]))
> 4: plot(cmdscale(dist(iris[, -5])), col = as.integer(iris[, 5]),
> pch = c("o", "+")[1:150 %in% model$index + 1])
> 5: eval.with.vis(expr, envir, enclos)
> 6: eval.with.vis(ei, envir)
> 7: source(tf, local, echo = echo, prompt.echo = paste(prompt.prefix,
> getOption("prompt"), sep = ""), continue.echo =
paste(prompt.prefix,
> getOption("continue"), sep = ""), verbose = verbose,
max.deparse.length
> = Inf, encoding = "UTF-8", skip.echo = skips, keep.source =
TRUE)
> 8: example(svm)
>
> Possible actions:
> 1: abort (with core dump, if enabled)
> ..
>
> I did already "update.packages(), what can I still do.
Works just fine for me. What's your sessionInfo()?
Here's mine:
> sessionInfo()
R version 2.12.1 Patched (2010-12-27 r53883)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=English_Canada.1252 LC_CTYPE=English_Canada.1252
[3] LC_MONETARY=English_Canada.1252 LC_NUMERIC=C
[5] LC_TIME=English_Canada.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] e1071_1.5-24 class_7.3-3
loaded via a namespace (and not attached):
[1] tools_2.12.1
Peter Ehlers
>>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.