Are inflammation markers associated with overall arthritis severity?
- Draw on real-world data to inform our study
- Focus on common laboratory measures of inflammation
- Use visualization and simple statistical tools to explore relationships
\[
\text{JDSI}_i =
0.4 \cdot \text{ESR}_i
+ 0.3 \cdot \text{CRP}_i
+ 0.3 \cdot \text{RA}_i
+ \varepsilon_i
\]
Construct JDSI from the imported dataset variables:
arthritis$JDSI <- 0.4*arthritis$ESR +
0.3*arthritis$CRP +
0.3*arthritis$RA +
rnorm(nrow(arthritis), 0, 5)
head(arthritis[, c("ESR","CRP","RA","JDSI")])
## ESR CRP RA JDSI ## 1 23.0 7.37 11.97 14.700838 ## 2 12.0 1.50 10.00 8.535674 ## 5 100.0 92.65 10.00 75.487648 ## 6 32.5 0.88 10.00 10.478171 ## 7 47.5 10.13 10.00 19.743492 ## 8 12.0 2.24 12.59 15.480834
\[
\hat{\beta}_{\text{ESR}}
\;\pm\;
t_{n-p,\,0.975}
\cdot
\text{SE}(\hat{\beta}_{\text{ESR}})
\]
lm(JDSI ~ ESR + CRP + RA, data = arthritis)
## ## Call: ## lm(formula = JDSI ~ ESR + CRP + RA, data = arthritis) ## ## Coefficients: ## (Intercept) ESR CRP RA ## -2.2592 0.4447 0.3223 0.3792