-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathload-support-scan.R
29 lines (25 loc) · 1.47 KB
/
load-support-scan.R
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
##
## load data from a minSupport scan
##
run1.graph = "HLAB"
run1.totFRs = 975
run1.minSup = c(300,100,9000,8000,7000,6000,5000,4000,3000,2000,1000,200,300,500)
run1.numFRs = c(646,975, 137, 192, 234, 270, 299, 362, 429, 474, 576,766,646,594)
run1.correct = c(.572,.572,.562,.565,.570,.570,.569,.567,.568,.569,.567,.572,.572,.569)
run1.TPR = c(.558,.559,.405,.540,.564,.555,.557,.555,.555,.557,.555,.557,.558,.552)
run1.FPR = c(.413,.414,.282,.410,.425,.415,.420,.421,.419,.420,.421,.412,.413,.414)
run1.MCC = c(.145,.145,.130,.130,.139,.140,.137,.134,.136,.137,.135,.145,.145,.139)
run1.sensitivity = run1.TPR
run1.specificity = 1.0 - run1.FPR
run1.fracFRs = run1.numFRs/run1.totFRs
run2.graph="SCZ6A+SCZ14C"
run2.totFRs = 1173
run2.minSup = c( 100,200,9000,8000,7000,6000,5000,4000,3000,2000,1000,800,600,400,500,495,480,463,1400)
run2.numFRs = c(1173,1173, 84, 98, 120, 145, 159, 172, 183, 196, 216,222,247,1125,316,350,602,893,206)
run2.correct = c(.595,.595,.553,.553,.563,.558,.557,.557,.557,.557,.566,.567,.577,.594,.587,.588,.593,.593,.563)
run2.TPR = c(.570,.570,.416,.505,.541,.532,.537,.533,.533,.529,.543,.543,.553,.564,.544,.544,.546,.548,.535)
run2.FPR = c(.380,.380,.311,.400,.416,.415,.423,.418,.419,.415,.412,.408,.399,.377,.369,.367,.360,.361,.410)
run2.MCC = c(.190,.190,.109,.106,.125,.117,.114,.115,.114,.114,.132,.135,.154,.188,.176,.177,.187,.188,.125)
run2.sensitivity = run2.TPR
run2.specificity = 1.0 - run2.FPR
run2.fracFRs = run2.numFRs/run2.totFRs