# NOT RUN {
activity.gm <- 10
activity.wm <- activity.gm * 0.4
percent.gm <- matrix(seq(0.1, 1, by=0.1), nrow=2)
percent.wm <- -percent.gm + 1
activity.obs <- percent.gm * rnorm(n=length(percent.gm), mean=activity.gm, sd=5) +
rnorm(n=length(percent.wm), mean=activity.wm, sd=5)
activity.corrected <- partialVolumeCorrection(activity.obs, percent.gm, percent.wm)
# }
Run the code above in your browser using DataLab