Noi/NOTES/2012-4-30: Difference between revisions
Jump to navigation
Jump to search
>Noi No edit summary |
>Noi No edit summary |
||
Line 1: | Line 1: | ||
* [[http://genome-tech.ucsd.edu/LabNotes/index.php/noi:DMR220k_LabNotes '''Link to calendar''']] | * [[http://genome-tech.ucsd.edu/LabNotes/index.php/noi:DMR220k_LabNotes '''Link to calendar''']] | ||
* Good tutorial for LDA by R: https://onlinecourses.science.psu.edu/stat857/book/export/html/17 | * Good tutorial for LDA by R: https://onlinecourses.science.psu.edu/stat857/book/export/html/17, http://miner.chem.purdue.edu/Lectures/Lecture10.pdf | ||
== Regression analysis of UCLA SZ data, applying LDA function == | == Regression analysis of UCLA SZ data, applying LDA function == | ||
* Continued from: http://genome-tech.ucsd.edu/LabNotes/index.php/Noi/NOTES/2012-4-17 | * Continued from: http://genome-tech.ucsd.edu/LabNotes/index.php/Noi/NOTES/2012-4-17 |
Revision as of 16:51, 3 May 2012
- [Link to calendar]
- Good tutorial for LDA by R: https://onlinecourses.science.psu.edu/stat857/book/export/html/17, http://miner.chem.purdue.edu/Lectures/Lecture10.pdf
Regression analysis of UCLA SZ data, applying LDA function
- Continued from: http://genome-tech.ucsd.edu/LabNotes/index.php/Noi/NOTES/2012-4-17
- I repeated regression analysis again by changing some criteria to select the candidate significant CpG sites, based on FDR or with minimum raw p-value. For status regression analysis, I will not apply multiple test correction since I didn't get any signal. However with minimum t statistic of 2.5 from the previous analysis I could separate out the normal and affected samples. Dr. 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.
Color code for PCA plot
- Age
mycol=c(1,1,4,4,4,1,1,4,1,1,4,4,1,4,1,1,4,1,1,4,4,4,4,1,1,4,4,1,1,4,4,4,1,1,4,4,4,1,1,4,1,1,4,4,1,1,4,4,1,1,4,4,4,4,1,1,4,4,1,4,4,4,4,1,1,4,4,1,1,4,4,1,1,4,4,4,1,1,4,4,1,1,4,4,1,1,4,4,1,1,4,4,1,1,4,4,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2) mycol[mycol==1]<- "black" #1=UCLA parents mycol[mycol==4]<- "#32CD32" #4=UCLA children(limegreen) mycol[mycol==2]<- "#3A5FCD" #2=UCSD-GA:all old adults(royalblue3)
- sex
mycol=c(1,4,1,1,1,1,4,1,1,4,4,1,4,1,1,4,1,1,4,4,1,1,1,1,4,1,1,1,4,1,4,1,1,4,1,4,1,1,4,1,1,4,1,1,4,1,1,1,1,4,1,4,1,1,1,4,4,1,1,1,1,4,1,1,4,1,1,1,4,1,1,1,4,1,1,1,1,4,1,1,1,4,4,1,1,4,1,1,1,4,1,1,1,4,1,4,2,8,2,2,8,2,2,8,2,2,2,2,8,8,2,8,2,2,2,2,2,8,2,8,2,8,8,8,8,8,8,2,2,2,8,8,2,8,8,8,2,2,2,2,2,2,8,8,8,2,8,8,8,8,2,8,8,2,8,8,2,8,8,8,2,2,8,8,8,8,2,2,2) mycol[mycol==1]<- "#000080" mycol[mycol==4]<- "#CD1076" mycol[mycol==2]<- "#1874CD" mycol[mycol==8]<- "#CD69C9"
- status
mycol=c(4,4,1,4,1,4,4,4,4,4,4,4,4,1,4,4,4,4,4,4,1,4,1,4,4,4,1,4,4,4,4,1,4,4,4,4,1,4,4,4,4,4,4,1,4,4,1,4,4,4,1,4,4,1,4,4,4,4,4,4,1,4,4,4,4,1,4,4,4,4,1,4,4,1,4,1,4,4,1,1,4,4,4,4,4,4,1,4,4,4,1,4,4,4,1,4,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2) mycol[mycol==1]<- "#EE0000" #1=UCLA, affected (red2) mycol[mycol==4]<- "#0000CD" #4=UCLA, unaffected (mediumblue) mycol[mycol==2]<- "#00BFFF" #2=UCSD-GA(deepskyblue)
- Directory: /home/nplongth/Noi_scratch/Data_analysis/UCLA_SZ_data_analysis_2012_04_05/UCLA_regression_analysis/2012_05_01_regression/
- I used the files generated from previous analyses. Age, sex and status data with the calculated p-value and q-value were printed out in separated files.
PCA analysis, required p-value <= 0.01
- Print out the CpG sites with the R square value > 0.1 and p-value <=0.01
awk '{if ($2>0.1 && $8<=0.01) print $0}' sorted_UCLA_regression_age_qval.txt>UCLA_regression_age_pval0.01 1038 UCLA_regression_age_pval0.01 awk '{if ($2>0.1 && $8<=0.01) print $0}' sorted_UCLA_regression_sex_qval.txt > UCLA_regression_sex_pval0.01 678 UCLA_regression_sex_pval0.01 awk '{if ($2>0.1 && $8<=0.01) print $0}' sorted_UCLA_regression_status.txt > UCLA_regression_status_pval0.01 110 UCLA_regression_status_pval0.01
- Extract methylation value into the matrix
./extract_methyFreqMatrix.pl ../UCLA-GA_methyl_min70_minSTD0.05 ./UCLA_regression_age_pval0.01 > UCLA_regression_age_pval0.01_MM 1038 UCLA_regression_age_pval0.01_MM ./extract_methyFreqMatrix.pl ./UCLA-GA_methyl_min70_minSTD0.05 ./UCLA_regression_sex_pval0.01 > UCLA_regression_sex_pval0.01_MM 678 UCLA_regression_sex_pval0.01_MM ./extract_methyFreqMatrix.pl ./UCLA-GA_methyl_min70_minSTD0.05 ./UCLA_regression_status_pval0.01 > UCLA_regression_status_pval0.01_MM 110 UCLA_regression_status_pval0.01_MM