* Click on the drop-down arrow on your "Pvalue" column.  Select "Custom".  In the window that appears, set a criterion that will filter your data so that the Pvalue has to be less than 0.05.
 
* Click on the drop-down arrow on your "Pvalue" column.  Select "Custom".  In the window that appears, set a criterion that will filter your data so that the Pvalue has to be less than 0.05.
 
** '''''How many genes have p value < 0.05? and what is the percentage (out of 5221)?'''''
 
** '''''How many genes have p value < 0.05? and what is the percentage (out of 5221)?'''''
 +
*** 948 genes, 18%
 
** '''''What about p < 0.01? and what is the percentage (out of 5221)?'''''
 
** '''''What about p < 0.01? and what is the percentage (out of 5221)?'''''
 +
*** 235 genes, 4.5%
 
** '''''What about p < 0.001? and what is the percentage (out of 5221)?'''''
 
** '''''What about p < 0.001? and what is the percentage (out of 5221)?'''''
 +
*** 24 genes, 0.46%
 
** '''''What about p < 0.0001? and what is the percentage (out of 5221)?'''''
 
** '''''What about p < 0.0001? and what is the percentage (out of 5221)?'''''
 +
*** 2 genes, 0.038%
 
* When we use a p value cut-off of p < 0.05, what we are saying is that you would have seen a gene expression change that deviates this far from zero less than 5% of the time.
 
* When we use a p value cut-off of p < 0.05, what we are saying is that you would have seen a gene expression change that deviates this far from zero less than 5% of the time.
 
* We have just performed 5221 T tests for significance.  Another way to state what we are seeing with p < 0.05 is that we would expect to see this magnitude of a gene expression change in about 5% of our T tests, or 261 times. (Test your understanding: [http://xkcd.com/882/ http://xkcd.com/882/].) Since we have more than 261 genes that pass this cut off, we know that some genes are significantly changed.  However, we don't know ''which'' ones.  To apply a more stringent criterion to our p values, we performed the Bonferroni and Benjamini and Hochberg corrections to these unadjusted p values.  The Bonferroni correction is very stringent.  The Benjamini-Hochberg correction is less stringent.  To see this relationship, filter your data to determine the following:
 
* We have just performed 5221 T tests for significance.  Another way to state what we are seeing with p < 0.05 is that we would expect to see this magnitude of a gene expression change in about 5% of our T tests, or 261 times. (Test your understanding: [http://xkcd.com/882/ http://xkcd.com/882/].) Since we have more than 261 genes that pass this cut off, we know that some genes are significantly changed.  However, we don't know ''which'' ones.  To apply a more stringent criterion to our p values, we performed the Bonferroni and Benjamini and Hochberg corrections to these unadjusted p values.  The Bonferroni correction is very stringent.  The Benjamini-Hochberg correction is less stringent.  To see this relationship, filter your data to determine the following:
 
** '''''How many genes are p < 0.05 for the Bonferroni-corrected p value? and what is the percentage (out of 5221)?'''''
 
** '''''How many genes are p < 0.05 for the Bonferroni-corrected p value? and what is the percentage (out of 5221)?'''''
 +
*** 0 genes, 0%
 
** '''''How many genes are p < 0.05 for the Benjamini and Hochberg-corrected p value? and what is the percentage (out of 5221)?'''''
 
** '''''How many genes are p < 0.05 for the Benjamini and Hochberg-corrected p value? and what is the percentage (out of 5221)?'''''
 +
*** 0 genes, 0%
 
* In summary, the p value cut-off should not be thought of as some magical number at which data becomes "significant".  Instead, it is a moveable confidence level.  If we want to be very confident of our data, use a small p value cut-off.  If we are OK with being less confident about a gene expression change and want to include more genes in our analysis, we can use a larger p value cut-off.   
 
* In summary, the p value cut-off should not be thought of as some magical number at which data becomes "significant".  Instead, it is a moveable confidence level.  If we want to be very confident of our data, use a small p value cut-off.  If we are OK with being less confident about a gene expression change and want to include more genes in our analysis, we can use a larger p value cut-off.   
 
* The "Avg_LogFC_all" tells us the size of the gene expression change and in which direction.  Positive values are increases relative to the control; negative values are decreases relative to the control.
 
* The "Avg_LogFC_all" tells us the size of the gene expression change and in which direction.  Positive values are increases relative to the control; negative values are decreases relative to the control.
 
** Keeping the (unadjusted) "Pvalue" filter at p < 0.05, filter the "Avg_LogFC_all" column to show all genes with an average log fold change greater than zero.  '''''How many are there? (and %)'''''
 
** Keeping the (unadjusted) "Pvalue" filter at p < 0.05, filter the "Avg_LogFC_all" column to show all genes with an average log fold change greater than zero.  '''''How many are there? (and %)'''''
 +
*** 352 genes, 6.7%
 
** Keeping the (unadjusted) "Pvalue" filter at p < 0.05, filter the "Avg_LogFC_all" column to show all genes with an average log fold change less than zero.  '''''How many are there? (and %)'''''
 
** Keeping the (unadjusted) "Pvalue" filter at p < 0.05, filter the "Avg_LogFC_all" column to show all genes with an average log fold change less than zero.  '''''How many are there? (and %)'''''
 +
*** 596, 11.4%
 
** '''''What about an average log fold change of > 0.25 and p < 0.05? (and %)'''''
 
** '''''What about an average log fold change of > 0.25 and p < 0.05? (and %)'''''
Exception encountered, of type "Error"
[2bbb1232] /biodb/fall2015/index.php?diff=3459&oldid=3448&title=Kevin_Wyllie_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}