An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120731/fd61788e/attachment.pl>
kernlab kpca predict
2 messages · Jessica Streicher
Hm.. seems like its a problem with loading it in the profile.. If i load it again in the console it works fine. Must have something to do with the masking. ---------- Changed load-order with the package which required "stats" that did the masking, and it works now, so, nevermind.
On 31.07.2012, at 16:25, Jessica Streicher wrote:
Hi!
The kernlab function kpca() mentions that new observations can be transformed by using predict. Theres also an example in the documentation, but as you can see i am getting an error there (As i do with my own data). I'm not sure whats wrong at the moment. I haven't any predict functions written by myself in the workspace either. I've tested it with using the matrix version and the dataframe/formula version on kpca.
data(iris)
test <- sample(1:150,20)
kpc <- kpca(~.,data=iris[-test,-5],kernel="rbfdot",kpar=list(sigma=0.2),features=2)
emb <- predict(kpc,iris[test,-5])
Error in UseMethod("predict") :
no applicable method for 'predict' applied to an object of class "c('kpca', 'prc')"
str(kpc)
Formal class 'kpca' [package "kernlab"] with 9 slots
..@ rotated : num [1:130, 1:2] -9.43 -9.17 -9.09 -9.43 -8.53 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:130] "1" "2" "4" "5" ...
.. .. ..$ : NULL
..@ pcv : num [1:130, 1:2] -0.224 -0.218 -0.216 -0.225 -0.203 ...
..@ eig : Named num [1:2] 0.323 0.119
.. ..- attr(*, "names")= chr [1:2] "Comp.1" "Comp.2"
..@ kernelf :Formal class 'rbfkernel' [package "kernlab"] with 2 slots
.. .. ..@ .Data:function (x, y = NULL)
.. .. ..@ kpar :List of 1
.. .. .. ..$ sigma: num 0.2
..@ kpar : list()
..@ xmatrix : num [1:130, 1:4] 5.1 4.9 4.6 5 5.4 4.6 5 4.4 4.9 4.8 ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:130] "1" "2" "4" "5" ...
.. .. ..$ : chr [1:4] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
.. ..- attr(*, "assign")= int [1:4] 1 2 3 4
..@ kcall : language kpca(x = x, data = ..1, kernel = "rbfdot", kpar = ..3, features = 2)
..@ terms :Classes 'terms', 'formula' length 2 ~Sepal.Length + Sepal.Width + Petal.Length + Petal.Width
.. .. ..- attr(*, "variables")= language list(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width)
.. .. ..- attr(*, "factors")= int [1:4, 1:4] 1 0 0 0 0 1 0 0 0 0 ...
.. .. .. ..- attr(*, "dimnames")=List of 2
.. .. .. .. ..$ : chr [1:4] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
.. .. .. .. ..$ : chr [1:4] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
.. .. ..- attr(*, "term.labels")= chr [1:4] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
.. .. ..- attr(*, "order")= int [1:4] 1 1 1 1
.. .. ..- attr(*, "intercept")= num 0
.. .. ..- attr(*, "response")= int 0
.. .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv>
.. .. ..- attr(*, "predvars")= language list(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width)
.. .. ..- attr(*, "dataClasses")= Named chr [1:4] "numeric" "numeric" "numeric" "numeric"
.. .. .. ..- attr(*, "names")= chr [1:4] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
..@ n.action: NULL
[[alternative HTML version deleted]]
______________________________________________ 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.