-
Notifications
You must be signed in to change notification settings - Fork 20
/
.Rhistory
133 lines (133 loc) · 2.53 KB
/
.Rhistory
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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
binom.test(14,11, 0.5)
binom.test(14,11)
binom.test(11,14)
binom.test(12,14)
binom.test(11,14)
binom.test(61,86)
?binom.test
?binom.test
binom.test(59, 116, "g")
binom.test(59, 116, 0.5, "g")
binom.test(59, 116, 0.5, "g")
package.install
install.packages
install.packages("gplots")
?read.table
source("plotHeat.R")
traceback()
source("plotHeat.R")
rnames
mat_data
mat_data
source("plotHeat.R")
data.mat
mat
mat_data
array
array(1, 10)
?array
?binom.test
binom.test(5,3, 1/19, "greater")
binom.test(3,5, 1/19, "greater")
1.0 - pbinom(3,19,3/19)
1.0 - pbinom(3,19,2/19)
pbinom(3,19,2/19)
1.0 - pbinom(2,3,2/19)
1.0 - pbinom(4,5,1/3)
pbinom(4,5,1/3)
1.0 - pbinom(3,5,1/3)
1.0 - pbinom(1,3,2/19)
1.0 - pbinom(1,3,3/19)
1.0 - pbinom(1,3,3/19)
1-(1/3)^7*(2/3)
1-(1/3)^7*(2/3)-(1/3)^6*(2/3)*7
1-(1/3)^7*(2/3)-(1/3)^6*(2/3)*7+1/3^7
1-((1/3)^7*(2/3))-((1/3)^6*(2/3)*7+(1/3)^7)
((1/3)^7*(2/3))-((1/3)^6*(2/3)*7+(1/3)^7)
((1/3)^7*(2/3))+((1/3)^6*(2/3)*7+(1/3)^8)+
1.0 - pbinom(1,2,2/18)
1.0 - pbinom(1,3,3/18)
?pbinom
?pbinom
1.0 - pbinom(3,3,3/18)
1.0 - pbinom(3,3,3/18)
pbinom(3,3,3/18)
pbinom(2,3,3/18)
pbinom(1,3,3/18)
1.0 - pbinom(2,3,3/18)
1/19*1/18*1/17
3/19*2/18*1/17
3/18*2/17*1/16
?hclust
?hclust
?hclust
?hclust
install.packages("ROCR")
require(ROCR)
data(ROCR.simple)
pred <- prediction( ROCR.simple$predictions, ROCR.simple$labels)
perf <- performance(pred,"tpr","fpr")
plot(perf)
pred
ROCR.simple$predictions
ROCR.simple$labels
perf <- performance(pred,"tpr","fpr")
perf
plot(perf)
perf
df <- read.table("precRecData/chari2015Train-wang.tab")
df
head(df)
df <- read.table("precRecData/chari2015Train-wang.tab", headers=T)
?read.table
df <- read.table("precRecData/chari2015Train-wang.tab", header=T)
head(df)
df$label
df$seq
df$score
df$seq
df$label
df
pred <- prediction( df$score, df$label)
require(ROCR)
pred <- prediction( df$score, df$label)
pred
perf$x
perf
perf$x.values
dataframe(perf)
?dataframe
??dataframe
?frame
summary(perf)
head(perf)
??S4
?summary
pred
pred
hist(perf$tn)
pred@
require(ROCR)
data.frame([email protected], [email protected])
odf <- data.frame([email protected], [email protected])
names(odf) <- c('x', 'y')
write.table(odf, "out.tmp")
write.table(odf, "out.tmp", row.names=FALSE)
names(odf) <- na
names(odf) <- NULL
write.table(odf, "out.tmp", row.names=FALSE)
require(ROCR)
odf <- data.frame([email protected], [email protected])
odf
summary(odf)
odf <- data.frame([email protected], [email protected])
odf
data.frame([email protected], [email protected])
colnames(odf) <- NULL
odf
?write.table
data.frame([email protected], [email protected])
q