IL6 activity

Author

Gibran Hemani

Published

May 7, 2025

Background

n <- 50000

g_il6 <- rbinom(n, 2, 0.5)
g_il6r <- rbinom(n, 2, 0.5)
il6 <- rnorm(n, 2.19, 3.13) / (1000 * 23.7)
il6r <- rnorm(n, 102.46, 53.59) / 50
sgp130 <- rnorm(n, 307.13, 63.47) / 100


sd(il6)
sd(il6r)
sd(sgp130)^2


dimer <- 0.5 * (il6 + il6r + 0.5) - sqrt(
    ((il6 + il6r + 0.5)/2)^2 - (il6 * il6r)
)


trimer <- 0.5 * (dimer + sgp130 + 0.05) - sqrt(
    ((dimer + sgp130 + 0.05)/2)^2 - (dimer * sgp130)
)


var(dimer, na.rm = TRUE)
var(sgp130)


var(trimer, na.rm = TRUE)


hist(dimer / trimer)

var(sqrt(
    ((dimer + sgp130 + 0.05)/2)^2 - (dimer * sgp130)
), na.rm = TRUE)


range(sqrt(
    ((dimer + sgp130 + 0.05)/2)^2 - (dimer * sgp130)
), na.rm = TRUE)


var(0.5 * (dimer + sgp130 + 0.05), na.rm = TRUE)

cor(dimer, trimer, use="pair")

dimer2 <- dimer - trimer

cor(dimer2, trimer, use="pair")

hist(dimer2/trimer)

hist(trimer/dimer2)


cor(trimer / dimer2, dimer / trimer, use="pair")


cor(dimer, dimer2, use="pair")

cor(trimer/dimer2, il6, use="pair")
cor(trimer/dimer2, il6r, use="pair")
cor(trimer/dimer2, sgp130, use="pair")

cor(trimer/dimer2, sgp130, use="pair")

cor(dimer2/trimer, sgp130, use="pair")

hist(dimer2/trimer, breaks=100)


dimer / trimer - contaminated activity score
dimer2 / trimer - activity score
trimer / dimer2 - inactivity score
sessionInfo()
R version 4.5.0 (2025-04-11)
Platform: aarch64-apple-darwin20
Running under: macOS Sonoma 14.6.1

Matrix products: default
BLAS:   /Library/Frameworks/R.framework/Versions/4.5-arm64/Resources/lib/libRblas.0.dylib 
LAPACK: /Library/Frameworks/R.framework/Versions/4.5-arm64/Resources/lib/libRlapack.dylib;  LAPACK version 3.12.1

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

time zone: Europe/London
tzcode source: internal

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
 [1] htmlwidgets_1.6.4 compiler_4.5.0    fastmap_1.2.0     cli_3.6.5        
 [5] tools_4.5.0       htmltools_0.5.8.1 yaml_2.3.10       rmarkdown_2.29   
 [9] knitr_1.50        jsonlite_2.0.0    xfun_0.52         digest_0.6.37    
[13] rlang_1.1.6       evaluate_1.0.3