Nothing
context("Kernelized Implicitly Constrained Least Squares Classifier")
data(testdata)
modelform <- testdata$modelform
classname<-all.vars(modelform)[1]
D <- testdata$D
D_test <- testdata$D_test
X <- testdata$X
X_u <- testdata$X_u
y <- testdata$y
X_test <- testdata$X_test
y_test <- testdata$y_test
test_that("Kernel and Linear give same result: 2 class", {
g_kernel <- KernelICLeastSquaresClassifier(X,y,X_u,kernel=kernlab::vanilladot(),lambda=0.0000001, scale = TRUE,projection="supervised",y_scale=TRUE,x_center=TRUE)
g_linear <- ICLeastSquaresClassifier(X,y,X_u,intercept=FALSE,scale=TRUE,x_center=TRUE,y_scale=TRUE)
sum(loss(g_linear,X_test,y_test))
sum(loss(g_kernel,X_test,y_test))
expect_equal(predict(g_kernel,X_test), predict(g_linear,X_test))
expect_equal(loss(g_kernel,X_test,y_test), loss(g_linear,X_test,y_test),tolerance =10e-5)
})
test_that("Kernel and Linear give the same result for supervised projection settings", {
# Same for supervised projection
g_linear<-ICLeastSquaresClassifier(X,y,X_u,
projection="supervised", intercept=FALSE,
x_center=TRUE, scale=TRUE,y_scale=TRUE)
g_kernel<-KernelICLeastSquaresClassifier(X,y,X_u,
projection="supervised",
kernel=kernlab::vanilladot(), lambda=0.0000001,
x_center=TRUE, scale = TRUE, y_scale=TRUE)
expect_equal(mean(loss(g_linear,X_test,y_test)),mean(loss(g_kernel,X_test,y_test)), tolerance=10e-6)
expect_equal(as.numeric(g_linear@theta), as.numeric(t(g_kernel@Xtrain) %*% g_kernel@theta),tolerance=10e-5)
expect_equal(decisionvalues(g_kernel,X_test), decisionvalues(g_linear,X_test),tolerance=10e-5)
expect_equal(g_linear@scaling, g_kernel@scaling)
expect_equal(g_linear@y_scale, g_kernel@y_scale)
})
test_that("Kernel and Linear give the same result for semi-supervised projection", {
g_linear <- ICLeastSquaresClassifier(X,y,X_u,
projection="semisupervised", intercept=FALSE,
x_center=TRUE, scale=TRUE, y_scale=TRUE)
g_kernel <- KernelICLeastSquaresClassifier(X,y,X_u,
projection="semisupervised",
kernel=kernlab::vanilladot(), lambda=10e-10,
x_center=TRUE, scale = TRUE, y_scale=TRUE)
expect_equal(mean(loss(g_linear,X_test,y_test)), mean(loss(g_kernel,X_test,y_test)), tolerance=10e-5)
expect_equal(mean(loss(g_linear,X,y)), mean(loss(g_kernel,X,y)), tolerance=10e-6)
expect_equal(as.numeric(g_linear@theta), as.numeric(t(g_kernel@Xtrain) %*% g_kernel@theta), tolerance=10e-5)
expect_equal(decisionvalues(g_kernel,X_test), decisionvalues(g_linear,X_test), tolerance=10e-5)
expect_equal(g_linear@scaling, g_kernel@scaling)
expect_equal(g_linear@y_scale, g_kernel@y_scale)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.