Skip to content

Instantly share code, notes, and snippets.

@jergosh

jergosh/gist:8802370

Created February 4, 2014 11:58
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jergosh/8802370 to your computer and use it in GitHub Desktop.
Save jergosh/8802370 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
gistfile1.txt
setwd("~/Documents/projects/paper-pat-seq/")
patsy_table <- read.table("Log/TEST_WT1_vs_MUT1_gcf.tab", header=T)
patsy_data <- patsy_table$runs_ma
names(patsy_data) <- as.character(patsy_table$transcripts)
pasta_table <- read.table("ExtData/PASTA_polyAranking.tab", header=T)
pasta_data <- pasta_table$weighted_length
names(pasta_data) <- as.character(pasta_table$transcripts)
palseq_table <- read.table("ExtData/PAL_Cerevisiae_total.tab")
palseq_data <- palseq_table$Mean.TL
names(palseq_data) <- as.character(palseq_table$Transcript.ID)
tr_union <- intersect(names(palseq_data), names(patsy_data))
cor(patsy_data[tr_union], palseq_data[tr_union])
plot(patsy_data[tr_union], palseq_data[tr_union])
# 0.04678423
tr_union <- intersect(names(pasta_data), names(patsy_data))
cor(patsy_data[tr_union], pasta_data[tr_union])
plot(patsy_data[tr_union], pasta_data[tr_union])
# -0.1607349
tr_union <- intersect(names(pasta_data), names(palseq_data))
cor(palseq_data[tr_union], pasta_data[tr_union])
plot(palseq_data[tr_union], pasta_data[tr_union])
# 0.1428611
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment