User Tools

Site Tools


neuroimagen:bioface_atn

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
neuroimagen:bioface_atn [2021/09/09 09:47]
osotolongo [What we got so far]
neuroimagen:bioface_atn [2021/09/12 10:21]
osotolongo [What we got so far]
Line 206: Line 206:
  
 <code> <code>
 +setwd("/home/data/bioface")
 base <- read.csv("bf_data.csv") base <- read.csv("bf_data.csv")
 base$Hippocampus = base$Left.Hippocampus + base$Right.Hippocampus base$Hippocampus = base$Left.Hippocampus + base$Right.Hippocampus
Line 228: Line 229:
  
 ^ Clasificación ^ Probabilidades de N+ ^ ^ Clasificación ^ Probabilidades de N+ ^
-| {{ :neuroimagen:neurodegeneration_classifier.png?600 | my pretty cool neurodegeneration classifier output }} | {{ :neuroimagen:nplus_probabilities_histogram.png?600 | how the probabilities are distributed? }} |+| {{ :neuroimagen:bioface_nplus_classifier.png?600 | my pretty cool neurodegeneration classifier output }} | {{ :neuroimagen:nplus_probabilities_histogram.png?600 | how the probabilities are distributed? }} |
  
 Y esto voy a hacerlo para FACEHBI también para compararlo, Y esto voy a hacerlo para FACEHBI también para compararlo,
  
 <code> <code>
 +setwd("/home/data/facehbi")
 read.csv("face_data.csv") -> face read.csv("face_data.csv") -> face
 face$Hippocampus = face$Left.Hippocampus + face$Right.Hippocampus face$Hippocampus = face$Left.Hippocampus + face$Right.Hippocampus
Line 243: Line 245:
 </code> </code>
  
 +^ Clasificación ^ Probabilidades de N+ ^
 | {{ :neuroimagen:facehbi_nbayes_classifier.png?600 | my pretty cool neurodegeneration classifier output (FACEHBI)}} | {{ :neuroimagen:facehbi_nplus_prob_histogram.png?600 |how the probabilities are distributed?}} | | {{ :neuroimagen:facehbi_nbayes_classifier.png?600 | my pretty cool neurodegeneration classifier output (FACEHBI)}} | {{ :neuroimagen:facehbi_nplus_prob_histogram.png?600 |how the probabilities are distributed?}} |
 +
 +==== Métodos no lineales ====
 +
 +Vamos a intentar con métodos no lineales. Ojo, estos  métodos no tienen interpretación probabilistica. Son construcciones matemáticas que funcionan.
 +
 +=== random Forest ===
 +
 +** BIOFACE **
 +<code>
 +rf_cl <- randomForest(ND ~ ., data = xt, na.action = na.omit)
 +base$ND <- predict(rf_cl, newdata = base)
 +base$post <- predict(rf_cl, newdata = base, type="prob")
 +plot(base$AGE, base$Hippocampus, main = "Hippocampus Volume versus Age", xlab="Age", ylab="HV", pch=19, col=ifelse(base$ND==1,"red","green"))
 +hist(base$post[,2], col="red", density = 50, probability=TRUE, breaks = 10, main="Probabilities profile", xlab = "Probability of N+")
 +</code>
 +
 +{{ :neuroimagen:rf_classifier_bioface.png?600 |}}
 +
 +**FACEHBI** 
 +
 +<code>
 +face$ND <- predict(rf_cl, newdata = face)
 +face$post <- predict(rf_cl, newdata = face, type = "prob")
 +plot(face$AGE, face$Hippocampus, main = "Hippocampus Volume versus Age", xlab="Age", ylab="HV", pch=19, col=ifelse(face$ND==1,"red","green"))
 +hist(face$post[,2], col="red", density = 50, probability=TRUE, breaks = 10, main="Probabilities profile", xlab = "Probability of N+")
 +</code>
 +
 +{{ :neuroimagen:rf_classifier_facehbi.png?600 |}}
 +
 +=== support-vector machine ===
 +
 +**BIOFACE**
 +
 +<code>
 +svm_cl <- svm(ND ~ ., data = xt, na.action = na.omit)
 +base$ND <- predict(svm_cl, newdata = base)
 +plot(base$AGE, base$Hippocampus, main = "Hippocampus Volume versus Age", xlab="Age", ylab="HV", pch=19, col=ifelse(base$ND==1,"red","green"))
 +</code>
 +
 +{{ :neuroimagen:svm_classifier_bioface.png?600 |}}
 +
 +** FACEHBI **
 +
 +<code>
 +face$ND <- predict(svm_cl, newdata = face)
 +plot(face$AGE, face$Hippocampus, main = "Hippocampus Volume versus Age", xlab="Age", ylab="HV", pch=19, col=ifelse(face$ND==1,"red","green"))
 +</code>
 +
 +{{ :neuroimagen:svm_classifier_facehbi.png?600 |}}
neuroimagen/bioface_atn.txt · Last modified: 2022/06/14 08:22 by osotolongo