Noi/NOTES/2012-6-14

From ZhangLabWiki
Revision as of 20:23, 14 June 2012 by >Noi
Jump to navigation Jump to search

Regression analysis of UCLA SZ data (disease status)

  • From the comment on previous page
    • LDA: random choose samples for training and (have UCLA-SZ and UCSD-GA on training and testing data sets)
  • Training and testing data set (I used the same group of training and testing data as assigned for age analysis)
    • list 1-96 =>UCSD-SZ, 97-169 => UCSD-GA
    • trainAge <- c(1:40, 91:132)
    • testAge <- c(41:90, 133:169)
  • From the previous analysis, I used PC1-PC3 scores as input for LDA analysis. However when I plot scree plot of the variances contributed on each PC, the change in slope in the scree plot occurs at component 4. In this analysis, I'll put PC1-PC4 scores as input.
A=read.table("UCLA_regression_status_pval0.015_MM",header=TRUE,row.names=1)
B=na.omit(A)
B$STDEV=NULL
B$min_RD=NULL
B$mean_RD=NULL
pca1 <- prcomp(t(B), scale=TRUE)
screeplot(pca1, type="line") 
File:Screeplot status pval0.015.png
Script
 lda_status_PC1-4_v2.R
File:Plot LD1 status PC1-4 p.015 v2 training.png File:Plot LD1 status PC1-4 p.015 v2 test.png