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