Хотите знать, может ли кто-нибудь помочь мне с использованием ggplot для построения графика результатов линейной смешанной модели? Есть пациенты, перенесшие операцию по повреждению головного мозга, и я хочу показать, как плотность тракта меняется от «до» лечения к «после», также с учетом возраста, пола и объема поражения.
library(lme4)
library(car)
data <- data.frame("Subject" = c(1,1,2,2,3,3,4,4,5,5,6,6),"Time" = c("Pre","Post","Pre","Post","Pre","Post","Pre","Post","Pre","Post","Pre","Post"),"Lesion_Volume" = c(10,10,11,11,10,10,7,7,15,15,12,12),"Tract_Density" = c(100,80,102,79,105,80,92,82,99,51,120,102), "Age" = c(50,50,51,51,55,55,56,56,50,50,53,53), "Sex"=c(1,1,0,0,1,1,1,1,0,0,0,0))
tract_density <- lmer( Tract_Density ~ Lesion_Volume + Age +Sex+Time + (1 | Subject), data=data)
Anova(tract_density)
I считаю, что лучшим способом будет линейный график, где ось x показывает до и после, а ось y показывает плотность тракта, с линией, представляющей каждую тему. но тогда генеральная линия, также представляющая общую тенденцию? Также не уверены, есть ли рекомендуемый способ отображения уровня значимости или тот факт, что возраст / пол / объем поражения также учитывались моделью?
Заранее благодарим за любую помощь или совет !!! Бен