Question: Fwd: feature weights in 'superpc' package
0
6.4 years ago by
Yukti Choudhury10 wrote:
I am running 'superpc' to model a supervised principal component predictor to predict survival. My dataset has 18 features and 218 samples. After running superpc.predict.red to form the reduced model, I am trying to extract feature weights that are used to construct the reduced predictor. According to documentation for 'superpc', "wt" corresponding to "Weight for each feature, in constructing the reduced predictor" is one of the values of the output list from superpc.predict.red. However, after running superpc.predict.red, I do not find this value in the output. The same is true when I run the example script provided in the documentation, as below: set.seed(332) #generate some data x<-matrix(rnorm(1000*40),ncol=40) y<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) ytest<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) censoring.status<- sample(c(rep(1,30),rep(0,10))) censoring.status.test<- sample(c(rep(1,30),rep(0,10))) featurenames <- paste("feature",as.character(1:1000),sep="") data<-list(x=x,y=y, censoring.status=censoring.status, featurenames=featurenames) data.test<-list(x=x,y=ytest, censoring.status=censoring.status.test, featurenames= featurenames) a<- superpc.train(data, type="survival") fit<- superpc.predict(a, data, data.test, threshold=1.0, n.components=1, prediction.type="continuous") fit.red<- superpc.predict.red(a,data, data.test, threshold=.6) fit.red does not include a value called "wt". I am trying to derive a formula based on the weights of selected features which will to assign a supervised principal components score, with which survival outcome can be determined. Any advice on the use of this function will be much appreciated. Yukti Choudhury Postdoctoral Fellow Institute of Bioengineering and Nanotechnology, A*STAR, Singapore [[alternative HTML version deleted]]
survival assign • 900 views
modified 6.4 years ago by Tim Triche4.2k • written 6.4 years ago by Yukti Choudhury10
Answer: Fwd: feature weights in 'superpc' package
0
6.4 years ago by
Tim Triche4.2k
United States
Tim Triche4.2k wrote:
why did you not at least cc: the maintainer of the package (Robert Tibshirani) on your email? The odds are quite good that he will understand your query much better than those of us who are merely users of the methodology, as he and Eric Bair developed it. Cc:'ing Dr. Tibshirani re: weights in superpc. Best, --t On Mon, May 6, 2013 at 2:28 AM, Yukti Choudhury <yukti.c@gmail.com> wrote: > I am running 'superpc' to model a supervised principal component predictor > to predict survival. > > My dataset has 18 features and 218 samples. After running > superpc.predict.red to form the reduced model, I am trying to extract > feature weights that are used to construct the reduced predictor. > > According to documentation for 'superpc', "wt" corresponding to "Weight for > each feature, in constructing the reduced predictor" is one of the values > of the output list from superpc.predict.red. However, after running > superpc.predict.red, I do not find this value in the output. > > The same is true when I run the example script provided in the > documentation, as below: > > set.seed(332) > #generate some data > > x<-matrix(rnorm(1000*40),ncol=40) > y<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) > ytest<-10+svd(x[1:60,])$v[,1]+ .1*rnorm(40) > censoring.status<- sample(c(rep(1,30),rep(0,10))) > censoring.status.test<- sample(c(rep(1,30),rep(0,10))) > > featurenames <- paste("feature",as.character(1:1000),sep="") > data<-list(x=x,y=y, censoring.status=censoring.status, > featurenames=featurenames) > data.test<-list(x=x,y=ytest, censoring.status=censoring.status.test, > featurenames= featurenames) > > a<- superpc.train(data, type="survival") > > fit<- superpc.predict(a, data, data.test, threshold=1.0, n.components=1, > prediction.type="continuous") > > fit.red<- superpc.predict.red(a,data, data.test, threshold=.6) > > > > fit.red does not include a value called "wt". > > > I am trying to derive a formula based on the weights of selected features > which will to assign a supervised principal components score, with which > survival outcome can be determined. > > Any advice on the use of this function will be much appreciated. > > > > Yukti Choudhury > Postdoctoral Fellow > > Institute of Bioengineering and Nanotechnology, A*STAR, > Singapore > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > -- *A model is a lie that helps you see the truth.* * * Howard Skipper<http: cancerres.aacrjournals.org="" content="" 31="" 9="" 1173.full.pdf=""> [[alternative HTML version deleted]]