(在延續這個post)問:KNN中的R - 奇怪的行爲
有誰知道爲什麼下面的KNN R代碼裏面給出了不同種不同的預測?這很奇怪,因爲K <-5,因此大多數都是明確的。另外,浮點數很大 - 因此不會出現數據問題的精確度+數據被縮放和居中。
library(class)
from = -(2^30)
to = -(from)
seed <- -229881389
set.seed(seed)
K <- 5
m = as.integer(runif(1, K, 20))
n = as.integer(runif(1, 5, 1000))
train = matrix(runif(m*n, from, to), nrow=m, ncol=n)
trainLabels = sample.int(2, size = m, replace=T)-1
test = matrix(runif(n, from, to), nrow=1)
sc<-function(x){(x-mean(x))/sd(x)}
train<-apply(train,2,sc)
test<-t(apply(test,1,sc))
seed <- as.integer(runif(1, from, to))
set.seed(seed)
pred_1 <- knn(train=train, test=test, cl = trainLabels, k=K)
message("predicted: ", pred_1, ", seed: ", seed)
seed <- as.integer(runif(1, from, to))
set.seed(seed)
pred_2 <- knn(train=train, test=test, cl = trainLabels, k=K)
message("predicted: ", pred_2, ", seed: ", seed)
手動檢查:
euc.dist <- function(x1, x2) sqrt(sum((x1 - x2)^2))
result = vector(mode="numeric", length=nrow(train))
for(i in 1:nrow(train)) {
result[i] <- euc.dist(train[i,], test)
}
a <- data.frame(result, trainLabels)
names(a) = c("RSSE", "labels")
b <- a[with(a, order(result, decreasing =T)), ]
headK <- head(b, K)
message("Manual predicted K: ", paste(K," class:", names(which.max(table(headK[,2])))))
print(b)
將使預測0,對於前k個(= 5)。
可能重複的[問:KNN在R-奇怪的行爲](http://stackoverflow.com/questions/38932289/q-knn-in-r-strange-behavior) – nrussell
這不是一個重複的,因爲這個問題是做縮放和居中(雖然[這一個](http://stackoverflow.com/questions/38932289/q-knn-in-r-strange-behavior)不是) – erans
順便說一句,對於一個可重複的例子,你應該不要使用隨機本身的種子......也就是說,在應用knn時存在一些隨機性,這是處理關係所需要的。所以,使用不同的種子你會得到不同的結果。 – DatamineR