data(iris, package="datasets")
dat = scale(iris[,-5])
res_pca = i_pca(dat[1:50,-5])
aa = seq(50,150,50)
for (k in c(1:2)) {
res_pca = update(res_pca, dat[c((aa[k]+1):aa[k+1]),])
}
plot(res_pca, animation = FALSE)
Run the code above in your browser using DataLab