| Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision |
| neuroimagen:bioface_atn [2021/09/10 09:20] – [random Forest] osotolongo | neuroimagen:bioface_atn [2021/09/12 10:23] – [Métodos no lineales] osotolongo |
|---|
| </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?}} | |
| |
| rf_cl <- randomForest(ND ~ ., data = xt, na.action = na.omit) | rf_cl <- randomForest(ND ~ ., data = xt, na.action = na.omit) |
| base$ND <- predict(rf_cl, newdata = base) | 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")) | 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> | </code> |
| |
| {{ :neuroimagen:rf_classifier_bioface.png?600 |}} | ^ Clasificación ^ Probabilidades de N+ ^ |
| | | {{ :neuroimagen:rf_classifier_bioface.png?600 |}} | {{ :neuroimagen:rf_prob_profile_bioface.png?600 |}} | |
| |
| **FACEHBI** | **FACEHBI** |
| <code> | <code> |
| face$ND <- predict(rf_cl, newdata = face) | 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")) | 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> | </code> |
| |
| {{ :neuroimagen:rf_classifier_facehbi.png?600 |}} | ^ Clasificación ^ Probabilidades de N+ ^ |
| | | {{ :neuroimagen:rf_classifier_facehbi.png?600 |}} | {{ :neuroimagen:rf_prob_profile_facehbi.png?600 |}} | |
| |
| === support-vector machine === | === support-vector machine === |
