attr1 attr2 ...

PC1 PC2 PC3 ...
... ... ...

getS3method("predict", "prcomp")

pca<-prcomp(data)
scale(newdata,with=FALSE]),pca$center,pca$scale)%*%pca$rotation