23 May 2014

thermal reactionome

Next post Previous post

May 20 - 23 notes



Is it necessary to parse “Low” transcripts into those that are up-regulated at low temperatures, versus those that are shut-down at high temperatures. To examine this, plotted this for some examples using following code:

A22.low.Ar.int.TPM.dt.sub <- resp.TPM.dt.sub.pred[J(A22.low, Ar.int), "A22")]

for(i in unique(A22.low.Ar.int.TPM.dt.sub$Transcript)[1:10]) {
    s <- A22.low.Ar.int.TPM.dt.sub[i]
    plot(s$val, s$pTPM, ylim = c(0, max(s$pTPM)+.1))
    points(s$val, s$TPM, pch = 8, col = "red")

Some that clearly have activated at high temperatures, while others steadily decline from cold to hot temperatures. Decided to split this into “On low” and “Off high” by taking advantage of quadratic fit; if the expression level at any temperature above 30°C is the minimum, designate the gene as ‘off’ at high temperatures (shallow curve, minimum at or beyond 30°C). Otherwise, designate as ‘on’ at high temperatures.

for(i in unique(A22.low.Ar.int.TPM.dt.sub$Transcript)[1:10]) {
    s <- A22.low.Ar.int.TPM.dt.sub[i]
    if(min(s[which(s$val < 30) , pTPM]) < s[which(s$val == 38.5), pTPM]) print("On low") else {
        print("Off high")


Starting to use DAVID throught the rDAVIDWebService package from Bioconductor. Detailed installation explained here.



Working on DAVID analysis

Writing…sent out draft!


Creative Commons Licence
This work is licensed under a Creative Commons Attribution 4.0 International License.