sonad=read.table("http://www.tlu.ee/~jaagup/andmed/keel/korpus/doksonaliigid.txt",
  header=TRUE, sep=",")
s2=scale(sonad[,2:17])
k=prcomp(s2)
#summary(k)
#k$rotation
#plot(k$sdev^2)
factanal(s2, factors=2)