Editing
Noi/NOTES/2012-5-12
(section)
Jump to navigation
Jump to search
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
== Regression analysis of UCLA SZ data, applying LDA function == * Continued from: http://genome-tech.ucsd.edu/LabNotes/index.php/Noi/NOTES/2012-4-30 * Links ** http://stat.ethz.ch/education/semesters/ss2012/ams/Rscript/scriptV6.1.R ** http://little-book-of-r-for-multivariate-analysis.readthedocs.org/en/latest/src/multivariateanalysis.html * I repeated regression analysis again by changing some criteria to select the candidate significant CpG sites, based on FDR or with raw p-value cutoff. For regression on status, I will not apply multiple test correction since I didn't get any significant signals. However with minimum t statistic of 2.5 (p-value ~0.014) from the previous analysis I could separate out the normal and affected samples. Dr. Zhang suggested that this should be fine. In addition, I will change the color of the control samples in PCA analysis plots. LDA will be used then the number of samples in the two groups classified by the line will be counted. === PCA analysis, plot the data in scatter plot === '''Scripts''' [[Media:plot_age_10FDR-2.txt| plot_age_10FDR-2.R]] [[Media:plot_sex_10FDR-2.txt| plot_sex_10FDR-2.R]] [[Media:plot_status_p0.015-2.txt| plot_status_p0.015-2.R]] '''Plots''' [[Media:UCLA-GA_age_10FDR_pcaplot1.pdf | UCLA-GA_age_10FDR_pcaplot1]] [[Media:UCLA-GA_age_10FDR_pcaplot2.pdf | UCLA-GA_age_10FDR_pcaplot2]] [[Media:UCLA-GA_sex_10FDR_pcaplot1.pdf | UCLA-GA_sex_10FDR_pcaplot1]] [[Media:UCLA-GA_sex_10FDR_pcaplot2.pdf | UCLA-GA_sex_10FDR_pcaplot2]] [[Media:UCLA-GA_status_pval0.015_pcaplot1.pdf | UCLA-GA_status_pval0.015_pcaplot1]] [[Media:UCLA-GA_status_pval0.015_pcaplot2.pdf | UCLA-GA_status_pval0.015_pcaplot2]] * I generate PC score (PC1-PC3). The reason is that the three components seem to contribute to largest variance and the slope change in the scree plot occurs at component 3 (age) as the picture below. Actually, I should try to analyze with different number of component and see if the classification change or not. Note the plot of status seem to be different from age because the slope change occur at component4 or 5. [[File:age_pca_var.png| 200px]] A=read.table("10%FDR_UCLA-GA_age_methylMatrix.txt",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) summary(pca1) v<-data.frame(pca1$x[,1],pca1$x[,2],pca1$x[,3]) names(v)<-c("PC1","PC2","PC3") * Then add the labels (age, sex or disease status) to each sample in the last column and used this data as the input for LDA V1 V2 V3 V4 1 0.8641308 3.1103773 -1.2768520 adult 2 -7.3738694 -3.7825118 2.4297650 adult 3 6.3925767 0.5246122 0.9405213 young adult 4 6.0525615 2.0565664 4.4626373 young adult '''Scripts''' [[Media:lda_age_10FDR.txt| lda_age_10FDR.R]] [[Media:lda_sex_10FDR.txt| lda_sex_10FDR.R]] [[Media:lda_status_10FDR.txt| lda_status_10FDR.R]] '''Plots''' [[Media:plot_LD1_age_10FDR.pdf | plot_LD1_age_10FDR]] [[Media:plot_LD1_sex_10FDR.pdf | plot_LD1_sex_10FDR]] [[Media:plot_LD1_status_p.015.pdf | plot_LD1_status_p.015]]
Summary:
Please note that all contributions to ZhangLabWiki may be edited, altered, or removed by other contributors. If you do not want your writing to be edited mercilessly, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource (see
ZhangLabWiki:Copyrights
for details).
Do not submit copyrighted work without permission!
Cancel
Editing help
(opens in new window)
Navigation menu
Personal tools
Not logged in
Talk
Contributions
Create account
Log in
Namespaces
Page
Discussion
English
Views
Read
Edit
View history
More
Search
Navigation
Main Page
Current events
Recent changes
Random page
Investigators
Matt Cai
Song Chen
Eric Chu
Dinh Diep
Elizabeth Duong
Shicheng Guo
Alan Fung
Daniel Jacobsen
Blue Lake
Huy Lam
Alice Li
Andrew Richards
Brandon Sos
Chris Wei
Yan Wu
Kun Zhang
Tools
What links here
Related changes
Special pages
Page information