liud4 / rVMAP

Data Management code for VMAC-MAP study
Other
3 stars 0 forks source link

CMRO2 and OEF derivation #36

Closed omair-a-khan closed 4 years ago

omair-a-khan commented 4 years ago

From Kim on 20191127 on Asana > VMAC Biostatistics Collaboration > MAP Derivation Requests:

CMRO2 = asl_rest_grey_matter_hct * (((asl_3t_trust_spo2 + asl_3t_trust_spo2_1+asl_3t_trust_spo2_2)/3)-oef_yv_hct) * (55.6 * bld_c_hgb)

oef_oef_hct = (((asl_3t_trust_spo2 + asl_3t_trust_spo2_1+asl_3t_trust_spo2_2)/3)-oef_yv_hct))/(oef_yv_hct*100)

oef_oef = (((asl_3t_trust_spo2 + asl_3t_trust_spo2_1+asl_3t_trust_spo2_2)/3)-oef_yv)/(oef_yv*100)
omair-a-khan commented 4 years ago

Kim had an error in the equations she specified for these variables. The code should be:

oef.ya <- (asl.3t.trust.spo2 + asl.3t.trust.spo2.1 + asl.3t.trust.spo2.2) / 3
cmro2 <- asl.rest.grey.matter.hct * (oef.ya - oef.yv.hct) * (55.6 * bld.c.hgb)
oef.oef <- 100 * ((oef.ya - oef.yv) / oef.ya)
oef.oef.hct <- 100 * ((oef.ya - oef.yv.hct) / oef.ya)