|  |  | 
|  | * 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"
[d727c925] /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}