search for: xscore

...------------------------ cancor.index<-function(object,x,y,center=T,scale=F){ x<-scale(x,center=center,scale=scale) y<-scale(y,center=center,scale=scale) ncor<-length(object$cor) #number of canonical variables nx<-dim(object$xco)[1] #number of X ny<-dim(object$yco)[1] #number of Y xscore<-x%*%(object$xcoe[,1:ncor]) colnames(xscore)<-paste("con",1:ncor,"x",sep=".") yscore<-y%*%(object$ycoe[,1:ncor]) colnames(yscore)<-paste("con",1:ncor,"y",sep=".") #canonical score eigenvalue<-object$cor^2/(1-object$cor^2) #e...

...ed within the validat slot are the predictions of each observation of Y when leave-one-out cross validation is applied. My code is as follows: > mypls <- mvr(Xtrain, Ytrain, method="SIMPLS", validation="CV", ncomp=1, niter=nrow(Ytrain)) > Xhat <- mypls$training$Xscores %*% t(mypls$training$Xload) > R2 <- 1-(sum((Xhat-Xtrain)^2)/sum(Xtrain^2)) > Q2 <- 1-(sum((Ytrain-mypls$validat$Ypred[,,1])^2)/sum(Ytrain^2) Is this right? Thank you Ana ------------------------------------------- Ana Conesa, PhD Bioinformatics Department Centro de Investigac...

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: https://stat.ethz.ch/pipermail/r-help/attachments/20070513/6557c355/attachment.pl