invitro = read.csv("data/invitro/plast_vs_dist_sjostrom_hausser_2006.csv")

invitro["epsp_ratio"] = invitro["epsp_ratio"] / 100
invitro = invitro[, c("epsp_ratio", "risetime")]

invitro
insilico = read.csv("data/insilico/plast_vs_dist.csv")

insilico
library(fasano.franceschini.test)

fasano.franceschini.test(S1=invitro, S2=insilico)

    Fasano-Francheschini Test

data:  invitro and insilico
D-stat = 0.39226, p-value = 5.578e-05
sample estimates:
    dff,1     dff,2 
0.3925676 0.3919595 
library(ks)

kde.test(x1=invitro, x2=insilico)
$Tstat
[1] 0.1573882

$zstat
[1] 3.20971

$pvalue
[1] 0.0006643442

$mean
[1] 0.03861666

$var
[1] 0.001369283

$var.fhat1
[1] 0.01707485

$var.fhat2
[1] 0.02114026

$n1
[1] 74

$n2
[1] 100

$H1
            [,1]        [,2]
[1,]  0.02318962 -0.08411697
[2,] -0.08411697  0.73436593

$H2
           [,1]       [,2]
[1,] 0.06318114 0.01550384
[2,] 0.01550384 0.14156182

$psi1
[1] 0.1627154

$psi12
[1] 0.08741366

$psi21
[1] 0.1022355

$psi2
[1] 0.184322
LS0tCnRpdGxlOiAiRGlzdGFuY2UtZGVwZW5kZW50IHBsYXN0aWNpdHkgaW4gTDIzUEMgdG8gTDVUVFBDIGNvbm5lY3Rpb25zIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCmBgYHtyfQppbnZpdHJvID0gcmVhZC5jc3YoImRhdGEvaW52aXRyby9wbGFzdF92c19kaXN0X3Nqb3N0cm9tX2hhdXNzZXJfMjAwNi5jc3YiKQoKaW52aXRyb1siZXBzcF9yYXRpbyJdID0gaW52aXRyb1siZXBzcF9yYXRpbyJdIC8gMTAwCmludml0cm8gPSBpbnZpdHJvWywgYygiZXBzcF9yYXRpbyIsICJyaXNldGltZSIpXQoKaW52aXRybwpgYGAKYGBge3J9Cmluc2lsaWNvID0gcmVhZC5jc3YoImRhdGEvaW5zaWxpY28vcGxhc3RfdnNfZGlzdC5jc3YiKQoKaW5zaWxpY28KYGBgCmBgYHtyfQpsaWJyYXJ5KGZhc2Fuby5mcmFuY2VzY2hpbmkudGVzdCkKCmZhc2Fuby5mcmFuY2VzY2hpbmkudGVzdChTMT1pbnZpdHJvLCBTMj1pbnNpbGljbykKYGBgCgoKYGBge3J9CmxpYnJhcnkoa3MpCgprZGUudGVzdCh4MT1pbnZpdHJvLCB4Mj1pbnNpbGljbykKYGBg