Difference between revisions of "Team:Peking/Modeling/Analysis"

Line 32: Line 32:
 
         .individual-picture{max-width:200px;margin-top:15px; margin-bottom:10px; }
 
         .individual-picture{max-width:200px;margin-top:15px; margin-bottom:10px; }
 
         .formula-inline{max-height:30px;}
 
         .formula-inline{max-height:30px;}
         .formula-line{max-height:50px;margin-top:10px; margin-bottom:10px;}
+
         .formula-line{max-height:50px;margin-top:0px; margin-bottom:0px;}
 
         ul.disc {list-style-type:disc;}
 
         ul.disc {list-style-type:disc;}
 
         ul.circle {list-style-type:circle;}
 
         ul.circle {list-style-type:circle;}

Revision as of 11:50, 18 September 2015

Modeling

Reliablity!!!!

Background

To increase the accuracy and specificity of the detection, we developed an assay over our Paired dCas9 Reporter (PC Reporter) System to get more sequence information from the target genome in the purpose of a more reliable result. We designed m pairs of gRNA specific target sites as m markers in the MTB genome. To make sure if the idea mention above actually work, here we used the target gene and the mismatched gene to have a test, respectively. In experimental group, the gRNAs were used to detect the target gene, while in control group, the gRNA were used to detect the mismatched gene. And to reduce the random error, both the experimental and the control group were repeated n times, the result would be shown as the optical power signals, which is generated by our Paired dCas9 Reporter System. Then by comparing the intensity of the optical power signal corresponding to the target gene and mismatched gene, the difference can be seen directly.

Description

  • There is no recognition site for gRNA in the mismatched gene.
  • The measure values Peking-Analysis-X_iY_i.gif of control and experimental group is independent of each other.
  • The measure values from n times repeated test compose the Peking-Analysis-X_iY_i.gif sample set, respectively, and the sample sets are both small.

Model

Model: Wilcoxon Rank Sum Test of Block Design

In view of the unknown distributions and different variances of the signals by our Paired dCas9 Reporter System, we chose a non-parametric statistics method called Wilcoxon Rank Sum Test of Block Design with the data Rank instead of ANOVA.
In the Block Design, we regarded the same gRNA detection of two treatment, i.e. target and mismatch DNA, as a block. To test the difference between two treatments, we test the null hypothesis that two treatment have no difference. The Wilcoxon Rank Sum statistics Peking-Analysis-W_j.gifof each block is calculated first by

Peking-Analysis-Wj%3DsumR_i Peking-Analysis-W_j_range

where Ri indicates the serial number of Xi in the population of both Xj and Yj. Note that Wilcoxon Rank Sum statistics Wj are distribution free and its distribution is known as long as the sample number is known.

For example, if n=3, {x1,x2,x3}={3,3,5}, {y1,y2,y3}={1,4,2}, so {x1,x2,x3,y1,y2,y3}={1,2,3,3,4,5}, which implies that {R1,R2,R3}={2,2,3}
Under the null hypothesis, after calculate all the possible order of two sample sets, the distributions of the statistics are shown as below:

W 6789101112131415
f(W) 0.050.050.100.150.150.150.150.100.050.05

Due to the small sample size, the minimal significance level is 0.05, which means only if Wj=15 leads to a rejection of the null hypothesis, in other words only when the minimum value of Xj was greater than the maximum value of Yj to accept the alternative hypothesis instead of the null hypothesis, the two sets of data is significantly different. So the Wilcoxon Rank Sum Test may face challenge in single block test when the experimental and control group are slightly different.However, by using Block Design, we can integrate data from m blocks similar to the idea of ANOVA. We calculated the sum of Wj(1<=j<=m) as the statistics.

Modeling_Fm5

The Wilcoxon Rank Sum Wj(1<=j<=m) from m blocks are independent and identically distributed (i.i.d), according to the central limit theorem (CLT), as m approaches infinity, the random variable Modeling_Fm6 converges in distribution to a standard normal distribution N(0,1)

Modeling_Fm7

So actually we use the statistics Modeling_Fm6, also we can calculate the p-value Modeling_Fm8, where Modeling_Fm9 is the distribution function of the standard normal distribution. If p-value is less than 0.01 or Modeling_Fm10, then we accept the alternative hypothesis that the two treatment, i.e. target and mismatch DNA, is highly statistic significantly.

Result

Modeling_Analysis_Figure1 Modeling_Analysis_Figure2

Fig. 1 Heatmaps of target and mismatch target. (a) Heatmap of target DNA assay. (b) Heatmap of mismatch target DNA assay.

In our experiment, n=3, so

Modeling_Fm11

(p-value = 7.868593384510819e-23) so target and mismatch DNA, are highly significantly different in signal.

Reference