Noi/NOTES/2012-3-21: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Noi
No edit summary
>Noi
No edit summary
Line 14: Line 14:
GK0344-004 --> GK0050-005a<br>
GK0344-004 --> GK0050-005a<br>
* Generating the list of MFASMAnalysisQValues.GK0XXXXXX.cpg.txt in UCLA_MFASM_Qcpg_list.txt: [[Media:UCLA_MFASM_Qcpg_list.txt]]
* Generating the list of MFASMAnalysisQValues.GK0XXXXXX.cpg.txt in UCLA_MFASM_Qcpg_list.txt: [[Media:UCLA_MFASM_Qcpg_list.txt]]
''' Binomial test '''
* Performed binomial test using asmBinomialTest.pl script written by Dinh: [[Media:asmBinomialTest.txt]]
* Performed binomial test using asmBinomialTest.pl script written by Dinh: [[Media:asmBinomialTest.txt]]
  ./asmBinomialTest.pl UCLA_MFASM_Qcpg_list.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt
  ./asmBinomialTest.pl UCLA_MFASM_Qcpg_list.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt
Line 19: Line 20:
  awk '{if ($98 <=0.01) print $0;}' UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.01.txt
  awk '{if ($98 <=0.01) print $0;}' UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.01.txt
* After lowering p-val cutoff I got total 4,708 SNP:CpG sites.
* After lowering p-val cutoff I got total 4,708 SNP:CpG sites.
''' T-test '''
* Perform t-test using asmTTest.pl script written by Dinh: [[Media:asmTTest.txt]]

Revision as of 01:39, 22 March 2012

ASM analysis of UCLA SZ data set

  • After discussing with Dr. Zhang and Dinh about how to identify sequence-dependent ASM, Dr. Zhang suggested to identify the consistency of sequence preference DNA methylation across individual
  • After having more discussion with Dinh, she helped me out to write scripts for sequence-dependent ASM using binomial test and t-test by testing on ASM HAPMAP data (total 12 samples): [[1]]
on genome-miner: /home/nplongth/Noi_scratch/ASM_HAPMAP1362-1454_2012_02_03/DD_SequenceDependentTest
  • I did the same analysis on UCLA SZ data set
on genome-miner: /home/nplongth/Noi_scratch/ASM_UCLA/combined_ASM_UCLA-4batches_2012_01_24
  • Since there were mislabeling of some samples at the beginning from UCLA and after mapping and ASM analysis, and I generated many files. I still kept all labeling the same as original, but made the correction when generated the methylation matrix or other sample_list files for any analysis. Here are the list of samples used in analysis and correction.

Correct data mislabeling in UCLA.SZ data set
Note: --> = change the labeling to
GK0210-001 --> GK0210-002
GK0210-002 --> GK0210-001
GK0050-005a --> GK0277-002
GK0277-002 --> GK0344-004
GK0344-004 --> GK0050-005a

Binomial test

./asmBinomialTest.pl UCLA_MFASM_Qcpg_list.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt
  • I got total 8,510 SNP:CpG sites showing significant ASM using p-value cutoff 0.05. However, when I look closer to the fraction of methylation level, it showed that some significant sites by this test and p-value cutoff have very close methylation level or had the methylation fraction 0 on both alleles. I may try to increasing stringency by lowering the cufoff p-value and observe the resulting methylation fraction.
awk '{if ($98 <=0.01) print $0;}' UCLA_SZ_SeqDepTest_Binomial_PVal0.05.txt > UCLA_SZ_SeqDepTest_Binomial_PVal0.01.txt
  • After lowering p-val cutoff I got total 4,708 SNP:CpG sites.

T-test