R/orth.r

orth<- function(A){
  ## questa function produce una base ortogonale Q per il range di A.
  ## Dunque t(Q) %*% Q =I, quindi le colonne di Q si trovano nello stesso
  ## spazio delle colonne di A. Il numero di colonne di Q equivale al rango di A.
  
  n=nrow(A)
  p=ncol(A)
  
  Asvd=svd(A)
  r=qr(A)$rank
  Q=Asvd$u[,1:r]
  Q
}

Try the clustrd package in your browser

Any scripts or data that you put into this service are public.

clustrd documentation built on May 8, 2019, 5:03 p.m.