Maxent model


This page contains some analysis of the Maxent model result, created Sun Sep 02 20:03:16 CEST 2018 using 'dismo' version 1.1-4 & Maxent version 3.4.1. If you would like to do further analyses, the raw data used here is linked to at the end of this page.


Analysis of omission/commission

The following picture shows the omission rate and predicted area as a function of the cumulative threshold. The omission rate is is calculated both on the training presence records, and (if test data are used) on the test records. The omission rate should be close to the predicted omission, because of the definition of the cumulative threshold.


The next picture is the receiver operating characteristic (ROC) curve for the same data. Note that the specificity is defined using predicted area, rather than true commission (see the paper by Phillips, Anderson and Schapire cited on the help page for discussion of what this means). This implies that the maximum achievable AUC is less than 1. If test data is drawn from the Maxent distribution itself, then the maximum possible test AUC would be 0.978 rather than 1; in practice the test AUC may exceed this bound.



Some common thresholds and corresponding omission rates are as follows. If test data are available, binomial probabilities are calculated exactly if the number of test samples is at most 25, otherwise using a normal approximation to the binomial. These are 1-sided p-values for the null hypothesis that test points are predicted no better than by a random prediction with the same fractional predicted area. The "Balance" threshold minimizes 6 * training omission rate + .04 * cumulative threshold + 1.6 * fractional predicted area.

Cumulative thresholdCloglog thresholdDescriptionFractional predicted areaTraining omission rate
1.0000.017Fixed cumulative value 10.0950.003
5.0000.129Fixed cumulative value 50.0560.019
10.0000.302Fixed cumulative value 100.0430.043
0.0060.000Minimum training presence0.3910.000
16.3500.47310 percentile training presence0.0350.099
9.7940.293Equal training sensitivity and specificity0.0430.043
3.0620.077Maximum training sensitivity plus specificity0.0650.006
1.1980.022Balance training omission, predicted area and threshold value0.0900.003
4.8850.126Equate entropy of thresholded and original distributions0.0560.019


(A link to the Explain tool was not made for this model. The model uses product features, while the Explain tool can only be used for additive models.)



Analysis of variable contributions


The following table gives estimates of relative contributions of the environmental variables to the Maxent model. To determine the first estimate, in each iteration of the training algorithm, the increase in regularized gain is added to the contribution of the corresponding variable, or subtracted from it if the change to the absolute value of lambda is negative. For the second estimate, for each environmental variable in turn, the values of that variable on training presence and background data are randomly permuted. The model is reevaluated on the permuted data, and the resulting drop in training AUC is shown in the table, normalized to percentages. As with the variable jackknife, variable contributions should be interpreted with caution when the predictor variables are correlated.

VariablePercent contributionPermutation importance
GebCo_TB_ext42.734.1
bio1425.932.5
bio415.527.2
GebCo_TB_ext_slope4.30.1
iwahashi40
bio63.81.1
KG1.30.3
bio1510.5
bio1313.3
bio120.30
GebCo_TB_ext_aspect0.10
bio10.10.8
bio500.2




Raw data outputs and control parameters


The data used in the above analysis is contained in the next links. Please see the Help button for more information on these.
The model applied to the training environmental layers
The coefficients of the model
The omission and predicted area for varying cumulative and raw thresholds
The prediction strength at the training and (optionally) test presence sites
Results for all species modeled in the same Maxent run, with summary statistics and (optionally) jackknife results


Regularized training gain is 2.884, training AUC is 0.980, unregularized training gain is 2.993.
Algorithm terminated after 500 iterations (47 seconds).

The follow settings were used during the run:
323 presence records used for training.
10036 points used to determine the Maxent distribution (background points and presence points).
Environmental layers used: GebCo_TB_ext GebCo_TB_ext_aspect GebCo_TB_ext_slope KG(categorical) bio1 bio12 bio13 bio14 bio15 bio4 bio5 bio6 iwahashi(categorical)
Regularization values: linear/quadratic/product: 0.050, categorical: 0.250, threshold: 1.000, hinge: 0.500
Feature types used: hinge product linear quadratic
outputdirectory: /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent
samplesfile: /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/presence
environmentallayers: /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/absence
autorun: true
visible: false
prefixes: false
Command line used: autorun -e /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/absence -o /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent -s /Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/presence -z

Command line to repeat this species model: java density.MaxEnt nowarnings noprefixes -E "" -E species outputdirectory=/Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent samplesfile=/Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/presence environmentallayers=/Users/micha/Documents/Masterarbeit/Results/maps/models/Rupicapra_pyrenaica_maxent/absence autorun novisible noprefixes -t KG -t iwahashi