Vamos a mirar ADNI, a ver que sale. En principio habria que ajustar tanot por edad como por volumen craneal. A primera vista pareceria una buena opcion, pero una vez que ajusto por edad y volumen total,

library(ADNIMERGE)
pop <- adnimerge[adnimerge$VISCODE=="bl" & (adnimerge$DX=="Dementia" | adnimerge$DX=="CN"),];
m1 <- lm(pop$Ventricles ~ pop$ICV)
pop$adj_Ventricles = pop$Ventricles - m1$coefficients[[2]]*(pop$ICV- mean(pop$ICV,na.rm=TRUE))
m2 <- lm(pop$adj_Ventricles ~ pop$AGE)
pop$vc_Ventricles = pop$adj_Ventricles - m2$coefficients[[2]]*(pop$AGE-mean(pop$AGE,na.rm=TRUE))

Y queda entonces,

uncorrected	corrected

Los valores estan demasiado juntos y la correccion los agrupa mas.

al escoger el hipocampo estoy presuponiendo que hay AD, pero vamos a ver que queda,

uncorrected	corrected

Pues no esta mal.

   UCSFVOLs <- c("FLDSTRENG", "FSVERSION", "IMAGEUID", "Ventricles", 
        "Hippocampus", "WholeBrain", "Entorhinal", "Fusiform", 
        "MidTemp", "ICV")

        ucsfvol$ICV <- ucsfvol$ST10CV
        ucsfvol$Hippocampus <- (ucsfvol$ST29SV + ucsfvol$ST88SV)
        ucsfvol$Entorhinal <- (ucsfvol$ST24CV + ucsfvol$ST83CV)
        ucsfvol$Fusiform <- (ucsfvol$ST26CV + ucsfvol$ST85CV)
        ucsfvol$MidTemp <- (ucsfvol$ST40CV + ucsfvol$ST99CV)
        ucsfvol$Ventricles <- (ucsfvol$ST30SV + ucsfvol$ST37SV + 
            ucsfvol$ST89SV + ucsfvol$ST96SV)
        ucsfvol$WholeBrain <- apply(ucsfvol[, c("ST128SV", "ST17SV", 
            "ST18SV", "ST61SV", "ST16SV", "ST53SV", "ST42SV", 
            "ST29SV", "ST12SV", "ST11SV", "ST65SV", "ST76SV", 
            "ST77SV", "ST120SV", "ST75SV", "ST112SV", "ST101SV", 
            "ST88SV", "ST71SV", "ST70SV", "ST124SV", "ST147SV", 
            "ST148SV", "ST150SV", "ST151SV")], 1, sum)

Usando el diccionario en ADNI podemos hacer algo como,

Fusiform =  rh.fusiform.GrayVol + lh.fusiform.GrayVol
MidTemp = rh.middletemporal.GrayVol + lh.middletemporal.GrayVol

y entonces,

pop$ND = ifelse(pop$DX == "Dementia", 1, 0)
xt <- pop[, c("Hippocampus", "Entorhinal", "Ventricles", "Fusiform", "MidTemp","AGE","ICV","DX")]

Esta es la cosa. Tengo un monton de informacion que podria clasificarme a los sujetos entre Neurodegeneracion SI y Neurodegeneracion NO. Y si tres variables juntas me pueden dar la info mas o menos adecuadad lo suyo es que meta todo lo relevante e intente usarlo al mismo tiempo. Lo primero deberia ser construir el vector de training a partir de los sujetos de ADNI.

> library("e1071")
> library("caret")
> library("caTools")
> library("ADNIMERGE")
> pop <- adnimerge[adnimerge$VISCODE=="bl" & (adnimerge$DX=="Dementia" | adnimerge$DX=="CN"),];
> pop$ND = ifelse(pop$DX == "Dementia", 1, 0)
> xt <- pop[, c("Hippocampus", "Entorhinal", "Ventricles","AGE","ICV","ND")]

Este es el vector de training pero voy a probar el algoritmo tambien asi que lo que hare es dividirlo y ver como funciona. Voy a tomar el 70% para entrenar el algoritmo y el restante 30_% para ver los resultados.

> split <- sample.split(xt, SplitRatio = 0.7)
> train_cl <- subset(xt, split == "TRUE")
> test_cl <- subset(xt, split == "FALSE")
> set.seed(120)
> classifier_cl <- naiveBayes(ND ~ ., data = train_cl)
> y_pred <- predict(classifier_cl, newdata = test_cl)
> cm <- table(test_cl$ND, y_pred)
> confusionMatrix(cm)
Confusion Matrix and Statistics

   y_pred
      0   1
  0 248  29
  1  46  84
                                          
               Accuracy : 0.8157          
                 95% CI : (0.7746, 0.8522)
    No Information Rate : 0.7224          
    P-Value [Acc > NIR] : 7.984e-06       
                                          
                  Kappa : 0.5609          
                                          
 Mcnemar's Test P-Value : 0.06467         
                                          
            Sensitivity : 0.8435          
            Specificity : 0.7434          
         Pos Pred Value : 0.8953          
         Neg Pred Value : 0.6462          
             Prevalence : 0.7224          
         Detection Rate : 0.6093          
   Detection Prevalence : 0.6806          
      Balanced Accuracy : 0.7935          
                                          
       'Positive' Class : 0 

Y no es perfecto pero esta bien!

Asi que para evaluar los sujetos habria que construir una matriz de datos identica a la de ADNI. Esto no es tan dificil, asi que vamos,

[osotolongo@brick03 bioface]$ xnat_pullfs.pl -s aseg -x bioface19 -o bf_base_aseg.csv
[osotolongo@brick03 bioface]$ xnat_pullfs.pl -s aparc -x bioface19 -o bf_base_aparc.csv
[osotolongo@brick03 bioface]$ join -t, bf_base_aseg.csv bf_base_aparc.csv > bf_base.csv
[osotolongo@brick03 bioface]$ awk -F"," '{print $3","$5}' bioface_np.csv | sed 's/Subjecte/Subject_ID/; s/Edad/AGE/' > bf_age.csv
[osotolongo@brick03 bioface]$ join -t, bf_base.csv bf_age.csv > bf_data.csv

Lo demas voy a hacerlo en R,

> base <- read.csv("bf_data.csv")
> base$Hippocampus = base$Left.Hippocampus + base$Right.Hippocampus
> base$Entorhinal = base$lh.entorhinal.GrayVol + base$rh.entorhinal.GrayVol
> base$Ventricles <- base$Left.Inf.Lat.Vent + base$Right.Inf.Lat.Vent + base$Left.Lateral.Ventricle + base$Right.Lateral.Ventricle
> base$ICV = base$eTIV

Ahora el procedimiento,

> classifier_cl <- naiveBayes(ND ~ ., data = xt)
> base$ND <- predict(classifier_cl, newdata = base)

A ver como queda,

> plot(base$AGE, base$Hippocampus, main = "Hippocampus Volume versus Age", xlab="Age", ylab="HV", pch=19, col=ifelse(base$ND==1,"red","blue"))

Nice!

A partir de la base de datos de ADNI, tomo los sujetos diagnosticados como Dementia y CN en el baseline, así como su edad y ciertos valores calculados a partir de la segmentacion. Tomo el diagnostico de Dementia como positivo en neurodegeneracion y negativo en caso contrario.

library("e1071")
library("caret")
library("caTools")
library("ADNIMERGE")
pop <- adnimerge[adnimerge$VISCODE=="bl" & (adnimerge$DX=="Dementia" | adnimerge$DX=="CN"),];
pop$ND = as.factor(ifelse(pop$DX == "Dementia", 1, 0))
xt <- pop[, c("Hippocampus", "Entorhinal", "Ventricles", "MidTemp", "AGE", "ICV", "ND")]

Con estos datos, entreno el clasificador,

classifier_cl <- naiveBayes(ND ~ ., data = xt)

Ahora, tomo un proyecto y reuno los datos de segmentacion y demograficos. No es neceario todo pero la edad y la segmentacion basica si. Esto puede cambiar de proyeto en proyecto pero es mas o menos asi, (lo repito por tenerlo todo junto)

xnat_pullfs.pl -s aseg -x bioface19 -o bf_base_aseg.csv
xnat_pullfs.pl -s aparc -x bioface19 -o bf_base_aparc.csv
join -t, bf_base_aseg.csv bf_base_aparc.csv > bf_base.csv
awk -F"," '{print $3","$5}' bioface_np.csv | sed 's/Subjecte/Subject_ID/; s/Edad/AGE/' > bf_age.csv
join -t, bf_base.csv bf_age.csv > bf_data.csv

Cargo estos datos y construyo las variables identicas a la de los datos de ADNIMERGE.

setwd("/home/data/bioface")
base <- read.csv("bf_data.csv")
base$Hippocampus = base$Left.Hippocampus + base$Right.Hippocampus
base$Entorhinal = base$lh.entorhinal.GrayVol + base$rh.entorhinal.GrayVol
base$Ventricles <- base$Left.Inf.Lat.Vent + base$Right.Inf.Lat.Vent + base$Left.Lateral.Ventricle + base$Right.Lateral.Ventricle
base$MidTemp = base$lh.middletemporal.GrayVol + base$rh.middletemporal.GrayVol
base$ICV = base$eTIV

Ahora puedo calcular el valor de (N),

base$ND <- predict(classifier_cl, newdata = base)

o la probabilidad de que exista neurodegeneración

base$post <- predict(classifier_cl, newdata = base, type="raw")

Clasificación	Probabilidades de N+

Y esto voy a hacerlo para FACEHBI también para compararlo,

setwd("/home/data/facehbi")
read.csv("face_data.csv") -> face
face$Hippocampus = face$Left.Hippocampus + face$Right.Hippocampus
face$Entorhinal = face$lh.entorhinal.GrayVol + face$rh.entorhinal.GrayVol
face$Ventricles <- face$Left.Inf.Lat.Vent + face$Right.Inf.Lat.Vent + face$Left.Lateral.Ventricle + face$Right.Lateral.Ventricle
face$MidTemp = face$lh.middletemporal.GrayVol + face$rh.middletemporal.GrayVol
face$ICV = face$eTIV
face$ND <- predict(classifier_cl, newdata = face)
face$post <- predict(classifier_cl, newdata = face, type="raw")