Saturday, July 5, 2014

Check new rls data



setwd("~/projects/0.network.aging.prj/9.ken")
tb = read.csv('conditionsWeibRedo.csv')
tb$genotype = toupper( as.character(tb$genotype))

str(tb)

tb2 = read.csv("SceORF_name.csv", header=F, colClass=c("character", "character"))

table( tb$genotype %in% tb2$V2 )




tb$flag = tb$genotype %in% tb2$V2
head(tb)

sub = tb[tb$flag, ]
head(sub)
str(sub)
length(unique(sub$genotype))

Some genotype are probably measured in different conditions.




No comments:

Post a Comment