library(datasets)
euclid <- function(x1,x2) #x musi byc postaci x1 - x2
{
sqrt(rowSums((x1-x2)^2))
}
wymiar <- 2
k_means <- function(dane, k, max_iter, wymiar)
#k - liczba klastrow, max_iter - maksymalna liczba iteracji
{
#
# losowanie punktow poczatkowych
#
means <- (sample(1:dim(dane)[1],k))
clusters <- dane[means,]
print('klastry')
print(clusters)
#do kazdego klastra zostana przypisane najblizej lezace punkty
#dla kazdego punktu obliczane sa odleglosci od wszystkich klastrow
#przypisywany jest do najblizszego
clustlist = list() #inicjalizacja listy na macierze klastrow
w = list()
for (i in 1:max_iter)
{
for (j in 1:k) #petla po wszystkich klastrach
{
# macierz klastrow
clustlist[[j]] <- t(matrix(clusters[j,],wymiar,dim(dane)[1]))
w[[j]]<-euclid(dane,unlist(clustlist[[j]]))
#print(w)
}
}
return(w)
}
a<-k_means(iris[, 3:4], 2, 1,2)
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}