AnalysisScripts/classifydataframe.R at master · TaliaferroLab/AnalysisScripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
myvector <- rep('none',nrow(test))
myvector2 <- rep('no', nrow(test))
myvector[which(test$kmer == 'GAGGG')] = 'HNRNPA2B1'
myvector[which(test$kmer == 'CAGGG')] = 'HNRNPA2B1'
myvector[which(test$kmer == 'AGGGC')] = 'HNRNPA2B1'
myvector[which(test$kmer == 'GGGAG')] = 'SRSF'
myvector[which(test$kmer == 'GGAGG')] = 'SRSF'
myvector[which(test$kmer == 'CAGGA')] = 'SRSF'
#myvector[which(test$kmer == 'UGCAU')] = 'RBFOX'
#myvector[which(test$kmer == 'GCAUG')] = 'RBFOX'
#myvector[which(test$kmer == 'UGCUU')] = 'MBNL'
#myvector[which(test$kmer == 'AUGCU')] = 'MBNL'
#myvector[which(test$kmer == 'UGCUG')] = 'MBNL'
#myvector[which(test$kmer == 'ACUAA')] = 'QKI'
#myvector[which(test$kmer == 'CUAAC')] = 'QKI'
#myvector[which(test$kmer == 'GUAUG')] = 'HNRNPU'
myvector2[which(test$corrected_pvalue < 0.1)] = 'yes'


test2 <- data.frame(test, protein = myvector, sigenriched = myvector2)