Science Forums Biology Forum Molecular Biology Forum Physics Chemistry Forum

Science Forums Biology Forum Molecular Biology Forum Physics Chemistry Forum (http://www.molecularstation.com/forum/)
-   BioStatistics Forum (http://www.molecularstation.com/forum/biostatistics-forum/)
-   -   Chi-Squared doesn't work!! (http://www.molecularstation.com/forum/biostatistics-forum/83318-chi-squared-doesnt-work.html)

kyr 08-09-2011 07:04 AM

Chi-Squared doesn't work!!
 
Hey guys

I need to do a Chi-Squared between two sequence and their reference.

I have a table of 2 columns and 2 rows.
column 1 : cell type 1
column 2 : cell type 2
row 1 : reference nucleotide
row 2 : mutation (SNP)

I wanna know if the difference between the reference and mutation is significant. So we counted the number of read (from sequencer) aligned to the reference and the mutation.

In this objective, I decided to do a Chi-squared with a 2x2 Contingency table with the following formula:
(AD-BC)^2 . N
(A+B)(C+D)(A+C)(B+D)

The problem is everytime that we have a de novo SNP with no read aligned (value 0) to the reference the chi-squared and p-value are returned as null (coz with this formulas you cannot time by 0).
For example if
A=0
B=0
C=22
D=239
==> FOR US this means that we have a de novo SNPs so it is significant. BUT X^2 RETURNS a value of 0...

Does anyone has a better idea to improve my method of finding de novo and significant mutation ?

Thanks

Abbreviations used in X^2 table :
ROW 1 / CELL 1 : A
ROW 1 / CELL 2 : B
ROW 2 / CELL 1 : C
ROW 2 / CELL 2 : D

A + B + C + D = N


All times are GMT. The time now is 10:35 PM.

Powered by vBulletin® Version 3.8.4
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.
Copyright 2005 - 2012 Molecular Station | All Rights Reserved

Page generated in 0.08913 seconds with 11 queries