-
Notifications
You must be signed in to change notification settings - Fork 0
/
joint_model_excluded_analysis_ks.r
87 lines (65 loc) · 9.48 KB
/
joint_model_excluded_analysis_ks.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
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
#!/usr/bin/env Rscriptlibrary(dgof)
# KS Test Code for pilot_data_final
d_Qi <- c(255.502281988, 325.03014496363636, 340.0255063032728, 355.0475437087273, 375.09480188424243, 415.1508959652727, 476.3684698739999, 520.3727085783031, 580.1384137520607, 685.1984424550909)
m_Qi <- c(270.4798520463664, 330.1460056546598, 340.06379819393567, 350.0177170783116, 360.08919136973856, 390.16047297866106, 481.1306915550341, 520.177189689558, 571.2964388391151, 675.2678693374806)
ks.test(d_Qi, m_Qi, alternative='two.sided')
d_Ben <- c(252.284822094, 320.15884852054546, 340.1596114709091, 355.16218160245455, 375.2046823663636, 410.02596441227274, 460.1240865122727, 530.2354699976363, 610.0819194214545, 760.1461231321819)
m_Ben <- c(250.06118918011057, 330.0919903941218, 345.05128209635546, 350.22598084877484, 365.12121841704555, 390.18245453596, 445.2340932972472, 545.2802528204758, 650.2560777796871, 841.2564662291451)
ks.test(d_Ben, m_Ben, alternative='two.sided')
d_Henry <- c(260.113788988, 320.06199930151513, 355.04916563157576, 380.0664899989091, 415.1318858153333, 461.21729048142424, 515.9670440330913, 580.2269543888788, 650.1695738107878, 749.8303148791813)
m_Henry <- c(260.10958144953685, 335.0751788313502, 375.1021533367173, 405.2005173013477, 450.2276358109083, 512.2807751288931, 585.8985338328889, 650.3394362512173, 715.2869041257219, 819.8461603136587)
ks.test(d_Henry, m_Henry, alternative='two.sided')
d_Esme <- c(252.77957371, 320.1280230036363, 330.1030085263636, 340.01811516, 350.0783000536363, 370.17925640545457, 415.2541211536364, 495.10015895454546, 560.2303559954545, 633.4851081009084)
m_Esme <- c(256.17419187366505, 325.0782627161711, 325.2058387670011, 330.1145999078986, 335.01566766337163, 340.07337267722085, 355.18224976243977, 445.22039624326857, 530.1439767217032, 590.2389506981798)
ks.test(d_Esme, m_Esme, alternative='two.sided')
d_Xi <- c(255.125279584, 300.1734665695354, 320.0382280123636, 335.1281125844849, 355.15519710965657, 400.2115100998788, 465.04242657212114, 510.23113283313137, 580.1106773217778, 695.1100247105454)
m_Xi <- c(262.32212166560765, 295.11580438186155, 305.13657796720713, 315.1415215028196, 325.2289756851097, 360.0739678774221, 450.26541639084473, 505.04681555524877, 580.0755119191106, 695.28409495015)
ks.test(d_Xi, m_Xi, alternative='two.sided')
d_Andrew <- c(265.08592903000005, 325.0685480755354, 345.14864531440406, 360.24935173175754, 385.12237089399997, 420.1814198046364, 470.1861727745151, 530.1745165591415, 600.164951542505, 750.421707039)
m_Andrew <- c(284.0533606674497, 340.1643110928304, 355.1101818148281, 365.08038149388597, 385.1023228456451, 420.18029307549165, 481.5875219594302, 560.2021422273781, 650.1864480682073, 849.6917369904597)
ks.test(d_Andrew, m_Andrew, alternative='two.sided')
d_Kaisei <- c(272.50699346199997, 305.03725841245455, 320.0344103919091, 340.12148882545455, 405.0279841280303, 460.03234870330306, 490.10845962863635, 520.1840115732728, 580.1239371068788, 688.6532878518177)
m_Kaisei <- c(280.1205611154334, 300.052354168622, 305.23087626125607, 320.089710111475, 425.28732698799155, 480.24878336984676, 500.2230274209716, 520.1945901786752, 570.2094547170744, 683.6058945672103)
ks.test(d_Kaisei, m_Kaisei, alternative='two.sided')
d_Sophia <- c(267.51631705200003, 320.06729654723233, 335.0266177873737, 345.2617681448485, 365.09681566072726, 405.17456648913134, 485.3333648866667, 525.1714811358788, 580.2366903122222, 700.1889382134545)
m_Sophia <- c(295.05687254746954, 325.0750373379489, 330.252616243536, 335.22777240765697, 350.006162874964, 370.2296138015713, 500.20855648655385, 530.0384556269963, 580.2170951553171, 710.1490911100068)
ks.test(d_Sophia, m_Sophia, alternative='two.sided')
d_random <- c(256.94558703399997, 315.11045005248485, 335.0205187218182, 350.15401671903027, 375.1217602999798, 420.0700376940606, 475.1699060929697, 525.1586529829293, 590.3573593079597, 700.3344677318181)
m_random <- c(271.8910777750025, 320.0464938106733, 330.1753933701325, 340.1573888070685, 360.0225683627771, 405.03017814491824, 480.20033038184033, 530.2801863130269, 605.1546664263561, 720.2583565159725)
ks.test(d_random, m_random, alternative='two.sided')
d_Ryan <- c(242.79266809400002, 365.0798432755656, 415.1158042784142, 474.34680104578797, 505.30514638916156, 533.5076828739394, 565.1829599539394, 595.2798451025758, 655.2037476812828, 768.4622123679088)
m_Ryan <- c(250.1986464179889, 380.27722096748226, 475.25389890456466, 520.366604664392, 540.2882245454687, 563.494540512154, 585.1484669787578, 610.2256249513513, 665.2584207582095, 813.016855778751)
ks.test(d_Ryan, m_Ryan, alternative='two.sided')
d_Riccardo <- c(290.11908416, 390.16717078745455, 455.26521451678786, 505.1987140270909, 530.227892134606, 550.3323376543636, 575.0743187349091, 600.0803639410303, 640.3426531589697, 695.3212617159999)
m_Riccardo <- c(324.7726778028116, 481.89444104418095, 535.1439533211701, 560.3047835430201, 570.3966994291521, 585.2228437634672, 595.2120673530022, 610.2113786500622, 640.2678668708353, 670.3046707321979)
ks.test(d_Riccardo, m_Riccardo, alternative='two.sided')
d_Govind <- c(291.66982916, 355.0721905199192, 415.0984181718384, 470.0381688780909, 490.14463658236366, 510.2270617354041, 534.9895773531516, 560.0422482487272, 597.4121626961717, 665.3693760787273)
m_Govind <- c(326.44450478636105, 370.10312219024263, 470.26081304055555, 500.080406935121, 510.23762237469776, 520.2849092883462, 535.000101970627, 540.2821854469249, 567.3957430176174, 611.5670721479631)
ks.test(d_Govind, m_Govind, alternative='two.sided')
d_Yushan <- c(278.514466336, 345.18741643199996, 390.10918731288893, 444.6189804613333, 490.0715924248889, 515.2639758062222, 540.139213744, 570.3158383182222, 615.1746030106666, 700.275733408)
m_Yushan <- c(305.52470627702064, 345.2453106758585, 415.09864236810887, 484.7777404117889, 515.2704843833006, 530.1605338377412, 550.0446406912841, 565.1858762830899, 592.1876394715747, 660.5069027532647)
ks.test(d_Yushan, m_Yushan, alternative='two.sided')
d_Xiaoyu <- c(277.524942694, 385.032964967899, 430.19777909658586, 465.15246158169697, 497.3884739666061, 525.3347479725657, 550.2251889324848, 580.0318839281213, 635.3472919576363, 730.1565902676363)
m_Xiaoyu <- c(312.52345113424155, 440.2533412065847, 485.19857488500037, 510.1487399101691, 532.3034332822388, 550.1119390234644, 560.1498549703239, 575.380702328141, 630.2566828930866, 717.8654362285481)
ks.test(d_Xiaoyu, m_Xiaoyu, alternative='two.sided')
d_Tiangan <- c(285.04643373799996, 359.3000744952929, 385.0957649091919, 435.17592885875763, 475.07080653387874, 500.1736754078586, 525.2559952856061, 555.2703743342828, 600.0805381907879, 680.884343024)
m_Tiangan <- c(320.3853479738107, 369.3586757994825, 395.23620262780497, 459.0533646060218, 490.14980308800693, 505.26092544911455, 520.0788357304908, 540.0349641443491, 565.3851766087205, 628.1432756086696)
ks.test(d_Tiangan, m_Tiangan, alternative='two.sided')
d_Zoey <- c(251.2507933, 355.10035229494946, 383.2466311612121, 435.8569420509091, 498.2836567866666, 535.2683732127273, 560.239672390303, 600.209904199192, 660.2947125042424, 776.5944497436363)
m_Zoey <- c(258.21888326273086, 360.20247156034196, 388.2721334055191, 476.28329822611147, 540.149905224997, 560.2873478779063, 580.2325739147338, 625.2245410769016, 699.0334741643583, 811.4269373360625)
ks.test(d_Zoey, m_Zoey, alternative='two.sided')
d_Sriram <- c(260.792910924, 333.7365065023838, 354.2306257701414, 404.75206238327274, 470.1757149942828, 515.0279465244445, 535.2394182011515, 570.2260076585454, 640.2403476935556, 821.0291923429093)
m_Sriram <- c(291.7508692576053, 330.11799038297653, 339.29367071491424, 394.6184859775916, 500.1266837169368, 530.0491967697222, 538.7673123961553, 560.3100517548746, 636.1210798757173, 933.1142769006258)
ks.test(d_Sriram, m_Sriram, alternative='two.sided')
d_Monde <- c(288.87370543599997, 330.13620774503033, 375.08692416351516, 445.13989636366665, 485.0803339738889, 515.3244962809596, 555.1523967670303, 600.190647148505, 660.004403178808, 750.9800414024542)
m_Monde <- c(305.5559606999359, 324.39150418712705, 370.12339492516935, 480.25980631451125, 515.0865219719524, 537.1705929300845, 575.1423517161193, 620.1794314171682, 670.1578375569524, 755.3404805720218)
ks.test(d_Monde, m_Monde, alternative='two.sided')
d_Thales <- c(255.95837770600002, 373.8896258031313, 460.17855856579797, 490.1102241592121, 500.28533881094955, 515.0414085794141, 525.0747692219394, 535.312324676909, 560.178602296303, 620.3847723172727)
m_Thales <- c(290.2793827229856, 453.89850366373327, 505.1655550558129, 515.0489790147034, 515.2149978449005, 520.1120816453396, 520.2168221944213, 525.1030671880571, 530.2062423416902, 550.1640065135032)
ks.test(d_Thales, m_Thales, alternative='two.sided')
d_Rushank <- c(272.699893192, 340.28547949555553, 425.31622340333337, 490.20173529557576, 511.0016001475555, 525.2055260751919, 540.2829276343031, 580.4031143295757, 630.1819748167879, 675.2745958321818)
m_Rushank <- c(302.0008362248576, 340.1193482327884, 500.2903653701554, 525.2055814587604, 530.9460568574076, 535.3435690013371, 545.167676519208, 580.3244224422477, 620.9955799025772, 645.1430073084422)
ks.test(d_Rushank, m_Rushank, alternative='two.sided')
d_Bill <- c(285.17701646, 345.1577092633535, 365.1959783543434, 395.05426746993936, 460.14612577727274, 495.3309730801414, 515.1412895860606, 535.2922037140203, 563.3416866037778, 623.9309867403636)
m_Bill <- c(311.9036492796989, 345.16632584986974, 360.083214809399, 375.1719559452432, 472.2244539491782, 500.16184737470786, 510.1034992144179, 515.1846573436538, 525.2412638381262, 555.2489096538833)
ks.test(d_Bill, m_Bill, alternative='two.sided')