==Sanity Check: Number of genes significantly change==
 
==Sanity Check: Number of genes significantly change==
 
*To determine how many genes fall under each pvalue category, I added a filter to all rows of data and filtered the data according to the values indicated.
 
*To determine how many genes fall under each pvalue category, I added a filter to all rows of data and filtered the data according to the values indicated.
**948 (18.2%) genes have a p value < 0.05
+
**'''948 (18.2%) genes have a p value < 0.05'''
**235 (4.5%) genes have a p value < 0.01
+
**'''235 (4.5%) genes have a p value < 0.01'''
**21 (0.4%) genes have a p value < 0.001
+
**'''21 (0.4%) genes have a p value < 0.001'''
**0 (0%) genes have a p value < 0.0001
+
**'''0 (0%) genes have a p value < 0.0001'''
 
*Two methods were applies in order to correct the p value and adjust it for the number of Ttests I ran. To determine how many genes fell under p value < 0.05, I filtered the data according to the values indicated in the columns indicated.
 
*Two methods were applies in order to correct the p value and adjust it for the number of Ttests I ran. To determine how many genes fell under p value < 0.05, I filtered the data according to the values indicated in the columns indicated.
**0 (0%) genes have a Bonferroni-corrected p value < 0.05
+
**'''0 (0%) genes have a Bonferroni-corrected p value < 0.05'''
**0 (0%) genes have a Benjamini and Hochberg-corrected p value < 0.05
+
**'''0 (0%) genes have a Benjamini and Hochberg-corrected p value < 0.05'''
 
*To determine which genes had a positive or negative fold change, I filtered all of the genes from the "Avg_LogFC_all" column with the indicated conditions.
 
*To determine which genes had a positive or negative fold change, I filtered all of the genes from the "Avg_LogFC_all" column with the indicated conditions.
**352 (6.7%) genes have a positive average log fold change in the "Avg_LogFC_all" column and a p value < 0.05
+
**'''352 (6.7%) genes have a positive average log fold change in the "Avg_LogFC_all" column and a p value < 0.05'''
**596 (11.4%) genes have a negative average log fold change in the "Avg_LogFC_all" column and a p value < 0.05
+
**'''596 (11.4%) genes have a negative average log fold change in the "Avg_LogFC_all" column and a p value < 0.05'''
**339 (6.4%) genes have an significant positive average log fold change (> 0.25) in the "Avg_LogFC_all" column and a p value < 0.05
+
**'''339 (6.4%) genes have an significant positive average log fold change (> 0.25) in the "Avg_LogFC_all" column and a p value < 0.05'''
**578 (11.1%) genes have an significant negative average log fold change (> 0.25) in the "Avg_LogFC_all" column and a p value < 0.05
+
**'''578 (11.1%) genes have an significant negative average log fold change (> 0.25) in the "Avg_LogFC_all" column and a p value < 0.05'''
 
*Question: What criteria did Merrell et al. (2002) use to determine a significant gene expression change? How does it compare to our method?
 
*Question: What criteria did Merrell et al. (2002) use to determine a significant gene expression change? How does it compare to our method?
**Merrell et al. determined statistically significant differences in gene expression using the Statistical Analysis for Microarrays (SAM) program. According to the article, genes that were determined to have statistically significant change were genes that had at least a two fold change in expression. Our method used p values to test for significant differences between experimental and control data.  
+
**'''Merrell et al. determined statistically significant differences in gene expression using the Statistical Analysis for Microarrays (SAM) program. According to the article, genes that were determined to have statistically significant change were genes that had at least a two fold change in expression. Our method used p values to test for significant differences between experimental and control data.'''
    
==Sanity Check: Compare individual genes with known data==
 
==Sanity Check: Compare individual genes with known data==
 
Merrell et al. (2002) report that genes with IDs: VC0028, VC0941, VC0869, VC0051, VC0647, VC0468, VC2350, and VCA0583 were all significantly changed in their data. However, in our data...
 
Merrell et al. (2002) report that genes with IDs: VC0028, VC0941, VC0869, VC0051, VC0647, VC0468, VC2350, and VCA0583 were all significantly changed in their data. However, in our data...
 
*(All genes were checked by searching the gene name in the "ID" column of the Excel spreadsheet.)
 
*(All genes were checked by searching the gene name in the "ID" column of the Excel spreadsheet.)
Exception encountered, of type "Error"
[d4109449] /biodb/fall2015/index.php?diff=cur&oldid=3337&title=Anuvarsh_Week_8 Error from line 434 of /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php: Call to undefined function each()
Backtrace:
#0 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(544): DiffEngine->diag()
#1 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(344): DiffEngine->compareSeq()
#2 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(227): DiffEngine->diffLocal()
#3 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(721): DiffEngine->diff()
#4 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(859): Diff->__construct()
#5 /apps/xmlpipedb/biodb/fall2015/includes/diff/DairikiDiff.php(980): MappedDiff->__construct()
#6 /apps/xmlpipedb/biodb/fall2015/includes/diff/TableDiffFormatter.php(194): WordLevelDiff->__construct()
#7 /apps/xmlpipedb/biodb/fall2015/includes/diff/DiffFormatter.php(140): TableDiffFormatter->changed()
#8 /apps/xmlpipedb/biodb/fall2015/includes/diff/DiffFormatter.php(82): DiffFormatter->block()
#9 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(888): DiffFormatter->format()
#10 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(802): DifferenceEngine->generateTextDiffBody()
#11 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(733): DifferenceEngine->generateContentDiffBody()
#12 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(662): DifferenceEngine->getDiffBody()
#13 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(632): DifferenceEngine->getDiff()
#14 /apps/xmlpipedb/biodb/fall2015/includes/diff/DifferenceEngine.php(453): DifferenceEngine->showDiff()
#15 /apps/xmlpipedb/biodb/fall2015/includes/page/Article.php(795): DifferenceEngine->showDiffPage()
#16 /apps/xmlpipedb/biodb/fall2015/includes/page/Article.php(506): Article->showDiffPage()
#17 /apps/xmlpipedb/biodb/fall2015/includes/actions/ViewAction.php(44): Article->view()
#18 /apps/xmlpipedb/biodb/fall2015/includes/MediaWiki.php(395): ViewAction->show()
#19 /apps/xmlpipedb/biodb/fall2015/includes/MediaWiki.php(273): MediaWiki->performAction()
#20 /apps/xmlpipedb/biodb/fall2015/includes/MediaWiki.php(566): MediaWiki->performRequest()
#21 /apps/xmlpipedb/biodb/fall2015/includes/MediaWiki.php(414): MediaWiki->main()
#22 /apps/xmlpipedb/biodb/fall2015/index.php(44): MediaWiki->run()
#23 {main}