The algorithm runs properly, and I ran some descriptive statistics on the final clusters. But it seems to converge rather slowly. So I might come back to this at a later point and try to improve it.