Description Usage Arguments Value Examples

This function identifies patterns of expression across differentially expressed genes and clusters them into groups based on similar patterns of expression across multiple series, such as time series data (I.E. day1, day2, day3). The resulting group data are plotted to show how expression of groups change over the course of the series based on a generalized linear model with a displayed 95 percent confidence interval, or based on the mean expression for each point in the series. Similarity is calculated by clustering and merging genes until all differentially expressed genes have been clustered into the desired number of groups.

1 2 3 |

`res_list` |
A list of DESeq result sets. Results can be calculated individually using DESeq's results() function. Lists of results can be created by creating a list(result1, result2 ... result_N). |

`filename` |
Filename for output plot. Valid extensions are ".pdf" and ".png". File generation for plot file can be turned off using set_output_mode("screen"). Output will be written to the /DE/series_plots/ directory. |

`designVar` |
The design field used when performing results() extractions. Must correspond to a metadata column. String. |

`groupBy` |
The variable by which to group data in the plot. I.E. "Condition_Time". String. |

`numGroups` |
The number of groups to split data into. Default = 5 |

`theme` |
Theme for the layout and color scheme for the plot. Valid selections are integers between 1-6. |

`method` |
Method for drawing. "glm" uses generalized linear model to show overall trends over time. "mean" uses the mean expression value for each group as intersection points. |

`customLabels` |
If customLabels is set to TRUE, the user will be prompted to provide a custom label for each label. |

`returnData` |
If this value is true, this function will return a list of dataframes corresponding to groupwise splits. Default: FALSE |

`writeData` |
If this value is true, this function will write data to file in the DE series folder for each group of genes. Default: FALSE |

`fetchGroup` |
Specifically fetch gene information for a specific group. This will override standard returnData and instead return data for the specified group number. |

If returnData is true, this function will return a list of dataframes corresponding to groupwise splits.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | ```
## Not run:
#This example assumes an experimenal design of ~Condition_Time.
#Prepare a result list.
res.day1 <- results(dds, contrast=c("Condition_Time", "day1_disease", "day1_control"))
res.day2 <- results(dds, contrast=c("Condition_Time", "day2_disease", "day2_control"))
res.day3 <- results(dds, contrast=c("Condition_Time", "day3_disease", "day3_control"))
myResList <- list(res.day1, res.day2, res.day3)
/*
* Cluster genes by similarity into 5 groups, then visualize their expression over the
* course of the series using a generalized linear model.
*/
de_series(res_list=myResList, filename="DE_series_pattern.pdf",
designVar="Condition_Time",
groupBy="Time", numGroups=5, theme=1, method="glm",
returnData=FALSE, writeData=FALSE)
/*
* Cluster genes by similarity into 3 groups, then visualize their expression over the
* course of the series using based on mean group expression values.
*/
de_series(res_list=myResList, filename="DE_series_pattern.pdf",
designVar="Condition_Time",
groupBy="Time", numGroups=3, theme=2, method="mean",
returnData=FALSE, writeData=FALSE)
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.