#######################
### affectIndiv
### Some trajectories
traj <- gald3d()["traj"]
### 4 clusters centers
center <- traj[runif(4,1,nrow(traj)),,]
### Affectation of each individual
part <- affectIndiv3d(traj,center)
#################
### K-means simulation (3 steps)
plot(clusterLongData3d(traj),parTraj=parTRAJ(col=part+1))
for (i in 1:3){
center <- calculTrajMean3d(traj,part)
part <- affectIndiv3d(traj,center)
plot(clusterLongData3d(traj),parTraj=parTRAJ(col=part+1))
}
Run the code above in your browser using DataLab