Frequently Asked Questions


 
Contents

Practical questions

  • How much costs Biomapper?
  • Does Biomapper work with Idrisi32?
  • I'm working with Arcview. How can I use Biomapper?
  • On which platform does Biomapper run?
  • How to convert an ESRI grid into an Idrisi/Biomapper raster map?
  • In which programming language is Biomapper developped?
  • Could you kindly send me the user manual.
  • How to quote Biomapper / ENFA?
  • Why is it not possible to use most of the modules of biomapper with Idrisi16 (.img) files?
  • Statistics

  • I don't understand what mean all the validation statistics
  • What is the "modified Kappa coefficient"?
  • Do you have an example data set to test Biomapper?
  • I have very good absence data. Can I use them with Biomapper?
  • How to use abundance data with Biomapper?
  • What is a Box-Cox transformation and why is it needed?
  • What is the "broken-stick advice"?
  • The number of categories per factor for the making of the HS map changes by steps of two. Is this normal?
  • I have computed an ENFA model and I would like to extrapolate it on a wider / other area. Is there some equation I can use to do it?
  • How is computed the score matrix?
  • Whate are global marginality, specialisation and tolerance? 
  • How do you compute the HS value for each cell from the scores matrix?
  • How to interpret factor biological meaning?
  • Box-Cox normalisation fails with some EGV maps. Should I discard them?
  • How are the ROC curves and kappa calculated if Biomapper does not use absence data?
  • How can I compare results from GLM and ENFA?
  • What is the difference between unidimensional and multidimensional histogram algorithms for the computation of HS maps and how strong is the impact on the resulting HS maps?
  • Can I compare the global marginality and specialisation coefficients of different species if I use the same area but different set of ecogeographical maps, in particular if have had to discard different correlated maps in the process?
  • You say that sometimes the marginality factor takes also a part of the specialisation into account. Where can I find how much?
  • What's the difference between "explained variance" and "explained information"?
  • I know that the habitat suitability of my species is linearly related to some variable. Nevertheless, in the HS map, the optimum for this variable seems not to be at an extremum. Why?
  • Procedures

    Importing files

    Problems

    Miscellanies

    Your Question
    If you don't find here an answer to your question, please feel free to ask me.
    Practical questions
    How much costs Biomapper?
    Biomapper is free. You just have to download it.
    Does Biomapper work with Idrisi32?
    Yes, Biomapper can read both Idrisi16 and Idrisi32 file format.
    I'm working with Arcview. How can I use Biomapper?
    Several of Biomapper's users are using Arcview. You will have to prepare your maps in Arcview before converting them to Idrisi/Biomapper format. Once the whole process has been completed, you can reimport the resulting maps into Arcview for further analysis or display.

    If you are using Arcview 3.0, you can use the Biomapper module Grid Convertor.

    How to convert an ESRI grid into an Idrisi/Biomapper raster map and conversely?
    There are several possibilities:

    1° ArcView 3.x extension
      There is an extension doing the job on the ESRI site, made by  Holger Schäuble. Look for Grid Converter  (av2idrisi.zip) on http://gis.esri.com/arcscripts/index.cfm. It works for ArcView 3.x.

    2° Biomapper's GridConvertor
      If you have Arcview 3.0 or 3.1 and Spatial Analyst, you can use the Biomapper module Grid Convertor. It allows to convert several files in only one operation.

    3° Manual conversion
      If you have another version of ArcView, GridConvertor may not work (always problems due to ESRI proprietary policy).You will therefore have to convert your grids manually. In the Biomapper help file you will find a description of the Idrisi/Biomapper file format that should help you to transfer your data.You can try this: in Arcview, go into the file menu and select Export Data Source and export your image as a binary raster (has an flt extension). Exit Arcview and then change the file extension from .flt to .rst. You should then be able to create a document file (*.rdc) for this image with the information supplied in the help file. The flt and rst format are identical. You will probably have to choose "real" data type.

    4° Otherwise...
      You can also find useful programs at http://www.pierssen.com/idrisi/grid.htm

    In which programming language is Biomapper developped?
    I'm using Borland Delphi for all my programming work. It allow me to quickly conceive fast running procedures cloaked in a user-friendly interface. The source code is not open but you will find the crucial procedures in the annexes of my PhD thesis.
    Could you kindly send me the user manual.
    There is no user manual for Biomapper. There is the help file though, and on the web, you will find the FAQ ( www.unil.ch/biomapper/faq.html ) and more general information.
    How can I use satellite imagery pictures (or other file format) with Biomapper?
    Biomapper has no importing capabilities. It can build a map from raw data but cannot convert alien file formats. It works only with the Idrisi (16 and 32) file format but Idrisi itself has an extensive set of conversion tools.
    How to quote Biomapper / ENFA?
    You can quote the main paper:

    Hirzel, A.H., Hausser, J., Chessel, D., Perrin, N., 2002. Ecological-niche factor analysis: How to compute habitat- suitability maps without absence data? Ecology 83, 2027-2036.

    and the Biomapper software itself:

    Hirzel, A., Hausser, J., Perrin, N., 2002. Biomapper 2.0. Lausanne, Lab. for Conservation Biology. URL: http://www.unil.ch/biomapper.

    Why is it not possible to use most of the modules of Biomapper with Idrisi16 (.img) files?
      Only the main Biomapper program has a button to switch from one format to the other. The other module use one of them according to the folder in which they are placed (it is not very elegant, I agree, and I will have to modify that some time in the future).

      When the module is launched, it checks if there is an "idrisi.env" file in the same folder as the module executable. If there is one, it switches to the Idrisi16 mode, else it switch to Idrisi32 mode.

      You can therefore enforce one mode or the other by placing a fake "idrisi.env" in your Biomapper folder.

    On which platform does Biomapper run?
      Biomapper is developped for Windows32 platforms, i.e. on all Windows since Win95. I have personally tested it on Win95, Win98, Win NT4 and Win2000. I have been told that it worked well on a Mac with a PC-card, but there is no version specially developed for Mac, Linux, Unix or any other platforms.
    Statistics
    I don't understand what mean all the validation statistics
      I agree that this validation part is not yet very well documented... It is the fruit of hard work to find a way to evaluate a model without absence data. I tried thus many methods that are still present in the output, without any further explanation... I shall try here to unveil a few details on this subject. This will finally be incorporated into the help file.

      You must have parted your sample into two sets. For validation purpose, you must use as "validation map" the set you did NOT used for model calibration. This observed map is therefore a boolean map indicating species presence. You can then evaluate the habitat suitability map produced with the ENFA model.

      When no reliable absence data are available, evaluation consist to compare various statistics computed on the "predicted map" (the habitat suitability map): 1° On the whole study area 2° Only on the validation points.

      The best way to understand this is to look at the box-plot displayed after the validation process. A good model should produce high species HS-value (80-100). The global box-plot  gives you insight on how marginal the species is in the studied area and thus how much these results could have been obtained by chance only.

      The results window gives a few statistics to resume this box-plot. It is composed of three parts: 1° species, 2° global statistics and 3° Comparisons.

      The two first part give identical information: first, common distribution statistics are computed (mean, median, quartiles, etc.). Then a few more specifics statistics: the one I find the most useful is the "proportion of presence cells >50": this is the proportion of validation points that have a predicted HS-value over 50. The higher this value, the better your model.The "Prob. to be above this value by chance" statistic use a bootstrap procedure to assess how much this value could have been obtained by chance. Practically, it is not very interesting as I have allways got here a 0.000 probability (good news)... Then you can see the 90th percentile and 95th percentile. You can compare all the previous values between global and validation sets, but only the validation set gives a really objective idea of the quality of your model.

      The last part give three comparison values between the two sets: The "Kappa" coefficient is a modified kappa statistics that integrate both how good is the model and how far from random it is. "Prop.of pix being significantly above 50" is the difference between the two "Prop.>50" values. And the "Probability to be over 50%" is the probability to be over 50% by chance, computed on the global set distribution histogram.

      These three last values are interesting to assess how the model is different from what could be achieved by a random model but it says nothing about the absolute quality of your model. They are highly related to the global HS of the study area and thus, if your species is not very marginal nor very specialised, the model could be very good but get a very bad "far from random" score.

    What is the "modified Kappa coefficient"?
    I gave this name to a home-made statistic whose behaviour was to be similar to the Cohen's Kappa coefficient. It is computed as follows: 

    Kmod = ( MS-MG ) / (1- MG )

    Do you have an example data set to test Biomapper?
    Alas, none of our data are in public domain and we cannot give them away.
    I have very good absence data. Can I use them with Biomapper?
    If you are really sure that your absence data are good and that no historical or spatial factors could have biased them, you will probably get a better model by using a presence/absence-based method, like Generalized Linear Model (GLM) or Generalized Additive Model (GAM). But you can ever put your absence data aside and apply ENFA on presence data only.
    How to use abundance data with Biomapper?
    You can use it to weight the presence data (weighted boolean map). Just use a map with integere weights in place of 0 and 1 as species map. 
    What is a Box-Cox transformation and why is it needed?
    Box and Cox (1964) developed a procedure for estimating the best transformation to normality within the family of power transformations:

    Y' = (YL -1)/L  (for L<>0)
    Y' = ln Y           (for L = 0)

    See "Biometry", Sokal & Rohlf, 1995, pp.417-419 for further explanations.

    Biomapper uses the Box-Cox algorithm to normalise as well as possible the ecogeographic variables. Empirically, we have found that normality was not a crucial factor and this step could as well be ignored.

    What is the "broken-stick advice"?
      The distribution of the eigenvalue is compared to the distribution of Mac Arthur's broken-stick. It is the expected distribution when breaking a stick randomly. Therefore, the eigenvalues that are larger than what would have been obtained randomly may be considered "significant". You can also keep only the factors with an eigenvalue larger than 1. These are objective means to choose how many factors to keep for HS map computing. These are just indications designed as a support when selecting the factors.
    The number of categories per factor for the making of the HS map changes by steps of two. I mean, you can only chose 2, 4, 6 and so on classes per factor. Is this normal?
      Yes. It is because the HS computation is based on the median of the factor distributions and the median must fall between two classes and so  the number of classes must be even.
    I have computed an ENFA model and I would like to extrapolate it on a wider / other area. Is there some equation I can use to do it?
      Presence-only models are difficult to extrapolate to other areas. Indeed they are based on the comparison between the locations where the species has been observed settling down and the available habitat. Although it can accurately make prediction on the study area, exporting the model to another place can be very tricky. In particular if you are comparing areas very far from each other.

      Even extrapolation in a closer area can be tricky if the environmental layers have not been collected in the same way. An environmental variable as simple as mean water temperature can be done in very different ways (time of the day, depth, season, etc.) which could prevent the model to make good predictions. Or it is sufficient to move the study window a few kilometers to make it cover a different habitat distribution, which will bring unpredictable disturbances.

      Finally, Biomapper has not been done for that purpose and there is presently no easy way to conduct such an extrapolation.

    How is computed the scores matrix?
      The full mathematical details of this operation are described in the main paper in Ecology. You can get an intuitive understanding by looking at the help file or on this site at http://www.unil.ch/biomapper/enfa.html . Here is a short view of this process:

      The eigenvalues and eigenvectors are extracted as follows: Compute the matrix W=Rs-1 * Rg where Rg is the global correlation matrix and Rs the species covariance matrix. From W, extract the marginality factor (I don't give here the mathematical procedure), which gives us the matrix W*. The specialisation factors are computed by extracting the eigensystem from W*.

    What are global marginality, specialisation and tolerance?
        Global Marginality = M = Sqrt(Si=1,v[Mi2] )/1.96
        Global Specialisation = S = Sqrt(Si=1,v(li)/V)
        Global tolerance = T = 1/S

    where Mi are the coefficients of the marginality factor, Sqrt() is the square root function, V is the number of variables and l are the eigenvalues.

    How do you compute the HS value for each cell from the scores matrix?
      This is a rather complex procedure. The full mathematical details of this operation are described in the main paper in Ecology.

      Shortly said, for each retained factor, a frequency histogram is computed over all the cells of the map. The median of this distribution is computed. The further a cell is from this median, the lower is its suitability for this factor. The global suitability is then obtained by computing a weighted mean on these "partial suitabilities". Marginality has a weight of 1, the sum of specialisation factors has a weight of 1, proportionally to their eigenvalue.

    How to interpret factor biological meaning?
      Look at the scores matrix. The first column of this matrix is the marginality factor. The other columns are the V-1 specialisation factors. (V is the number of variables). The rows are the EGV contributions to each factor.
    Box-Cox normalisation fails with some EGV maps. Should I discard them?
      For myself, when the Box-Cox fails, I keep the original map. A "Box-Coxised" map gives better results than a "brute" map, but a "brute" map is still better than no map at all. 

      Anyway, you may include it at the beginning (to compute the big time-consuming covariance map. Once it is computed, you can remove easily variables from it, but if you add new variables, the whole matrix will have to be recomputed) and then try to remove it to see how it affects the result.

    How are the ROC curves and kappa calculated if Biomapper does not use absence data?
      The kappa and the ROC curves that you find in the validation dialog box assume that blank values are true absences. They are therefore not suitable for most data sets where absences are unreliable. I put them here for the cases where you can rely upon absences.
    How can I compare results from GLM and ENFA?
      Comparing ENFA and GLM is a tricky stuff. In my recent paper (Hirzel et al, 2001, Ecological Modelling 145), I was able to compare them because I was using virtual data and so I had access to the "reality", the "truth". But in the general case, we do not know it and so we are constrained to use the standard statistics (Kappa, ROC, etc.). There are three main problems when comparing presence/absence to presence-only methods:

    1° If absences are thought to be unreliable to build a model, there is no sense in using them to validate it afterwards. So, the standard statistics are not useful. I tried to develop a few statistics to replace them (see the FAQ on the Biomapper site) but the perfect statistic has still to be invented.

    2° As it is based on presence data only, the ENFA is more efficient to model the areas with average to high suitability; its predictions for low-suitabitility regions should be taken with prudence. By contrast, presence/absence methods and in particular GLMs, will tend to model good versus bad areas causing such a kind of "stepped" response which is different from the linear one of the ENFA.

    3° Without absences - i.e. bad habitat points - to "fix the floor", ENFA must scale its suitability index to the ceiling. That means that, on an ENFA HS map, you will always, by construction, have at least one cell with a HS of 1 (or 100) ; with GLM, it is generally not the case as it is computing "probability of presence" and so the maximum values are generally quite lower than 1.

      Thus, when comparing visually GLM and ENFA maps (computed on well-known species at equilibrium, i.e., when absence data are largely reliable), the results are obviously quite similar. But if you try to compare them statistically, you get strange results biased either for one or for the other depending on which base hypothesis you use. To compare them you must correct both results to make them comparable:

    1° Remove the "ceiling effect" by stretching the GLM results between 0 and 1.

    2° Synchronise the "step effects" by transforming both GLM and ENFA results into boolean maps (by choosing a threshold).

    3° Then you can compare these results with standard presence/absence statistics.

    What is the difference between unidimensional multidimensional histogram algorithm for the computation of HS maps and how strong is the impact on the resulting HS maps?
      The ful detail of the unidimensional algorithm are in the main paper in Ecology. The multidimensional is by now obsolete but I let here the explanation for history's sake. I give here a "feeling" of what they do and how they differ:

    Unidimensional algorithm:

      Once the ENFA factors have been computed, it is possible to compute for every cell of the map a value along each of them (in fact, usually, one computes it only for a few of the first factors). The distribution histogram of each factor is then computed taking into account only the presence points. Each histogram will be used to attribute a "partial suitability" value to every cell (the more the cell factor value depart from the median of the histogram, the lower its "partial suitability". Then, a weighted sum of these partial suitability values is made for each map cell, and finally these sum are stretched in order to have their maximum at 100.
    Multidimensional algorithm:
      Here we do not address the selected factors independently. By crossing all factors together, the factor space is divided in small units (hypercubes). Then we count how many cells belong to each hypercube : this is the multidimensional histogram. This multidimensional distribution is computed both for all cells (global distribution) and the presence cells (species distribution). Finally, the HS value is computed for each hypercube by dividing the species hypercube by the global hypercube (and multiplied by 100).
      The problem with the multi algorithm is that it need very huge amount of presence data to be accurate, in particular if you want to include more than two factors, which is generally the case. It is also very memory-consuming. So far, I never got good results with this algorithm and it is why it will not be described in your paper. We strongly advise Biomapper users to use only the unidimensional algorithm. In fact, I could well have removed it from the interface...

      I am currently working on new kinds of algorithms, but it is another story... Stay tuned! Biomapperians will be informed of all the new developments.

    Can I compare the global marginality and specialisation coefficients of different species if I use the same area but different set of ecogeographical maps, in particular if have had to discard different correlated maps in the process?
      Strictly speaking, you cannot. Practically, the main biasing effect is the study area. When you have to remove a variable, it is because it does not contain more information than is already included into the model, so removing it should not alter significantly the global marginality and specialisation coefficients. Thus, provided the map sets do not differ drastically, you can still comparer them by these statistics. Anyway, do not assume too much significance to small differences in marginality or specialisation between species

      To be true, I never tested this. You could test it by building a common minimal EGV set and apply it to all your species. You will then see how the marginality et specialisation differ between the common data-set and the species-optimal one. Tell what you get, should you decide to try this.

    You say that sometimes the marginality factor takes also a part of the specialisation into account. Where can I find how much?
      The amount of variance explained by the first factor is in fact the amount of specialisation. It generally ranges from 10% to 70%. This value is given in the eigenvalue table.

      To summarise, the marginality factor explains always 100% of the marginality and some part of the specialisation. It is why it has always a great weight (minimum 0.5) for HS computation.

    What's the difference between "explained variance" and "explained information"?
      The Ecological Niche factors are conveying two kinds of information: marginality and specialisation. In fact, the first factor explains always 100% of the marginality and some varying part of specialisation; the subsequent factors explain no marginality and the rest of the specialisation.

      Historically, I was using the concept of "explained variance", which in fact was related only to "explained specialisation". Marginality was not used in this value. Therefore, when using this index for deciding how many factors had to be kept for the HS analysis, the user could be misled by this value. Accordingly, I introduced the concept of "explained information" to cope with this by giving marginality and specialisation the same "information power".

      Mathematically, if L1,L2,...,Lf are the eigenvalues of the f retained factors, and SL is the sum of all n eigenvalues (= L1+L2+...+Lf+...+Ln), we have :

    Explained variance = Explained specialisation = Se = (L1+L2+...+Lf)/SL

      In the explained information index, the marginality gets the same weight as the specialisation and we have therefore:

    Explained information = Ie = (Se+1)/2

      Ie is therefor always greater than Se and can never be smaller than 0.5. If the user chooses to keep only the first factor, he will explain at least half of the information, the one included into marginality. This Ie value is to be seen as a decision support value to choose how many factors are to be included into the HS analysis.

    I know that the habitat suitability of my species is linearly related to some variable. Nevertheless, in the HS map, the optimum for this variable seems not to be at an extremum. Why?
      The HS computing algorithm is not linear but bases itself on the observed distribution. Your problem arises generally on the marginality factor. Let's imagine a species linearly related to the frequency of forest and that this variable is strongly correlated to the marginality factor (but the reasoning hold also for specialisation factors) ; it means that, the more forest there is around a given cell, the more suitable it is for the species, the maximum being a frequency of 100%. This optimum is what we know from our  knowledge of the species, field studies, etc. Now, let's see the Biomapper's point of view: To it, the optimal frequency is the one where the species is the most frequent (More precisely, the median of the species distribution along the marginality factor. As the distribution of this factor is generally unimodal and more or less symmetrical, the median corresponds also to the most frequent.) Therefore, if large forests are rare in the study area, points with 100% forest freq. will be rare too, and the optimum for the species will not be 100% but lower (say 80%). Then, when computing the HS index, 80% freq. will provide the highest partial suitability value and this will decrease when freq. is either increasing or decreasing. The rarer the large forests, the steeper the rate of decrease.

      Sometimes this effect is welcome (when dealing with median optima) and sometimes counterproductive (with extreme optima). I am currently working at new algorithms which will hopefully address this problem.

    Procedures
    Do you have a step-by-step explanation of how doing a whole habitat suitability analysis?
      Yes. You will find it in the menu Help/Modus operandi. The whole process is outlined there and by clicking on the main titles you access a more detailed step-by-step description of the process.
    How can I create the species-presence map?
      There are several possibilities depending on what kind of data you have at hand:
      • List of observations coordinates: Put them in an ASCII file, using a structure x-coordinate tabulation y-coordinate (You can do this with Microsoft Excel) and the use the "Convertor" module to create a boolean presence map from this file.
      • Observation map in a point-vector-file: Simply rasterise it (with Idrisi function "PointRas"), using the same resolution and window as your ecogeographical maps.
      • Population map in a polygon-vector-file: First, rasterise it (with Idrisi function "PolyRas"), using the same resolution and window as your ecogeographical maps. Then make this map boolean (1=inside populations). Finally, use the module "Sampler" to divide this map into a calibration and a validation data sets.
    How to convert a point vector map into a species map?
    In Idrisi32:
    • Menu Reformat/Raster-vector conversion/POINTRAS
    • Select your vector map, e.g. "Species.vct"
    • Choose a name for the "image file to be updated". It must be a new name (not an existing raster) (e.g. "Species_bl.rst")
    • Choose "Change cells to record the presence of 1 or more points"
    Idrisi asks you if you want to bring up INITIAL: Answer YES.
    • In the INITIAL dialog box:
    • Select "Copy spatial parmeters..."
    • Select one of your EGV maps in "Image to copy parameters from"
    • In "Output data type", select "Byte"
    • Click OK
      And this should do the job. Species_bl can now be used as species map. You may want to partition it into calibration and validation data sets. You can do this with the Biomapper module Sampler.

      Rememember that all the EGV maps and the species map must be in the same directory.

    How to convert a polygon vector map into a species map?
    In Idrisi32:
  • Menu Reformat/Raster-vector conversion/POLYRAS
  • Select your vector map, e.g. "Species.vct"
  • Choose a name for the "image file to be updated". It must be a new name (not an existing raster) (e.g. "Species.rst") 
  • Idrisi asks you if you want to bring up INITIAL: Answer YES.
    • In the INITIAL dialog box:
    • Select "Copy spatial parmeters..."
    • Select one of your EGV maps in "Image to copy parameters from"
    • In "Output data type", select "Byte"
    • Click OK
    Then, you must booleanise the rasterised polygons:
    • Menu Analysis/Database query/image calculator
    • Select "Logical expression"
    • Type Species_BL = [Species]>0
    • Click OK
      Species_BL can now be used as species map. You may want to partition it into calibration and validation data sets. You can do this with the Biomapper module Sampler.

      Rememember that all the EGV maps and the species map must be in the same directory.

    How to insert the species-presence map?
        Once the species-presence map has been created, you must insert it in Biomapper in order to use it in the analyses. The species-map must be inserted in the Work maps list (NOT the EGV maps). This can be done in the  Files/Work maps/Add maps menu. Once inserted in this list, you must declare it as the current species map by selecting it, right-clicking on it and selecting "Mark as species map". The current species map is then marked by a red circle in front of its name.
    Why have the EGV maps to be quantitative?
      The ENFA is based on quantitative computations. For example, the marginality factor is based on means. Therefore, if you use a purely categorical map these computations will be misleading. Example: say you have a vegetation map with 1.grassland, 2.forest, 3.agriculture field, 4.bushes, but which is mainly constitute of types 1 and 4. The ENFA will computes that the global mean for this map is about 2.5 ( 0.5*(4+1) ), which at least doesn't mean anything and at worst is strongly misleading (could be taken as forest or agriculture land).
    I have a categorical map (say a vegetation type map). How would you handle preparing this map?
      I would use BOOLEANISATOR to get a boolean map of each relevant category and then feed these boolean maps into DISTAN and CIRCAN to get distance and frequency maps of them. You could even choose several radius for the moving window in order to take into account several influence distances.CIRCAN is useful for all resources (food, shelter, etc.) variables, which the species could need in its home range. DISTAN is useful to makes disturbance (mainly human impact like tourism, noise, pollution,etc.) variables. Generally, I compute both distance and frequency maps for all my boolean variables, compute the ENFA and look at the score to see if there is any difference between how the species is sensible to both aspects.

      You could also consider to use some fragmentation index according to the species you are mapping, as "Border length" in CIRCAN (It works well for animals/plants living or feeding near forest boundaries, like Ungulates)
     

    When I verify the EGV maps, I get a warning message telling "Map "xxxxx" is not continuous enough. What does it mean and how to fix this problem?
        This indicates that the EGV maps is either boolean or nearly boolean, i.e. the map contains of almost only two values. This may happen if you took too small a radius in Circan with too sparse an original boolean map; sometimes, the Box-Cox algorithm may have this effect too; in this case, just take the map withou Box-Coxing it. Feeding "nearly boolean" maps into the ENFA could cause "negative eigenvalue" or "very large eigenvalue" problems. This warning is just to make you aware of a possible cause of problems. In fact, visually, the EGV-maps should be beautiful rainbows of colors and not almost black and white to get the best results.

        To fix this problem you have four options:
       

      1. Don't change anything and try to perform the ENFA. If you get negative or very large eigenvalues, Biomapper will protest and you will know this is probably due to these maps. But perhaps it will work anyway. It's worth a try.
      2. Remove the faulty map(s) and perform the ENFA, but, alas! with less information.
      3. Try to increase the CircAn radius for the faulty map(s).
      4. Use DistAn for these maps.
    Is it possible to obtain a better model by reducing the amount of explained variance (i.e. the number of factors used) and, consquently, increasing the number of classes per factor?
        Yes. You must find the best trade-off between explained variance and smoothness of the HS model. Note that it is generally not useful to select more than 10 classes.
    How to compare models obtained by various factor number/class number trade-off ?
        You can use the validation module of Biomapper. Look at the box-plot it generates. Focus on the species box; it must be as high as possible. The higher and the narrower the better. (The global-box is not useful to assess model quality; it is here to see how different from randomness is your model. If you built your model on an area globally good for your species, you will get a good model simply because the species can live everywhere.)
    How can I print a map from biomapper? I guess I can do it from Idrisi but I´d like to use the rainbow palette.
        Display the map and save it as a BMP file. This file you can insert in a word document or print with any picture software. Alternatively, you will find the rainbow palette in the Biomapper package: Biomapper.smp.
    What are these "Biomapper extensions" used for? Can I ignore them?
      These extensions are made to simplify the browsing and help the user to select among all the maps, those having the right data type. But these extensions are not used by Biomapper to verify the maps; it uses the raster documentation file (*.rdc). Thus you can as well ignore the biomapper extensions.
    In the Options dialog box, it is possible to switch between correlation and covariance matrix, and to change the norming of the eigenvectors. However, this doesn't seem to affecte ENFA outputs.
    These options are intended for the Principal Components Analysis. They do not affect ENFA indeed.
    Problems
    I get an error message "... is not a valid floating point value"
      It happens because your computer is not set to use "." as decimal separator. It is mandatory if you want to use Biomapper and greatly advised if you want to use Idrisi. You can change this setting in the Start Menu/parameters/Control panel/Regional settings/Numbers.
    When I try to open a project, I get a message "Documentation file not found".
        This message means that Biomapper could not find a map. This happen generally because the map is not in the same directory as the project file (*.bio). All EGV maps must be in the same directory. This constraint was introduced in autumn 2000 in order to make the transfer of project between different computers transparent. By now, when you create a new project, Biomapper will verify that all your files are in the same directory. However, you could have problem with project created with previous versions.
    I get error messages when I add EGV maps, or try to normalise them, or mask them, etc.
      This kind of problems generally happen when the metadata file is incorrect.
    • If you are working with ArcView grids that you have imported into Idrisi format, your maps probably suffer from the infamous -9999 background bug. Check in the metadata file (*.rdc) what's the flag value: if it is -9999, this bug is probably  at work. There is a tool in the MapManager module that will help you to fix this problem.
    • If this does not solve your problmen, check the metadata file (*.rdc. You can open it in the NotePad, for example). Check for inconsistancies between your EGV maps and species maps, in particular in the following fields: columns, rows, Max X, Min X, Max Y, Min Y and ref. system.
    •  If this doesn't help, please contact me.
    Miscellanies
    Do you have a published paper about ENFA?
      There are already several papers. More are coming. You can get a list of related publications on the Bibliography page. There, you can also download a PDF version of my PhD thesis. Registered users will be informed when new papers are coming out. 
    Is Biomapper better suitable for plants or animals modelling?
      So far, Biomapper has been applied on a wide range of the life tree: Mammals, Birds, Insects, tropical Trees and Ferns. Most of them are animals but it's mainly because I use to meet more zoologists than botanists...

      In fact, absence data tend to be poorer for animals and so Biomapper could be particularly adapted to this case.
     

    I would like to look for any word in the Biomapper help file. How can I do it?
      The Windows help-file system allows you to generate a word search from any help file. Open the help file and click on the "Find" tab (rightmost tab). The first time you do this, you can choose the size of the database you want to use; generally I use the smallest size, but if you want more words to be included in the database, you can select a larger size. Then, the "Find" tab will provide you with the kind of tool you are looking for. Just type a word and the you will see a list of the "chapters" containing this word.
    Some checkboxes, buttons and options are not documented in the help file.
      The help file is actually written along the thread of the modus operandi (or step-by-step guide) to guide the user through Biomapper to allow her/him to compute and validate a HS map easily and understand the main outputs. I agree there are a lot of undocumented small controls (checkboxes, buttons, etc.) I added for the comfort of the user. Whilst only a few are described in the help file or in this FAQ on the website, all of them have a built-in help message that appears when you hover the mouse cursor on the control; the help message - or "mouse hint" - appears both in the bottom bar of Biomapper main window and in a floating temporary panel; they are present for most controls in all modules. These hints try to explain what's the use of the control. When a more complete explanation was needed, I generally put it either in the help file or in the FAQ. If I did not, please send me an e-mail and I shall clarify the situation.
    I found inconsistencies and undocumented options between the software and the help file.
      There are some minor inconsistencies between the help file and the program as  Biomapper is evolving continuously and I don't modify the help file each time. I hope you understand that I am a scientific researcher and not a professional programmer. My main interest is in solving ecological questions by developing more efficient algorithms and applying them to real cases. Biomapper is free and I'm alone to program it, therefore, I cannot devote - and am not interested in devoting - too much time and energy in help-file writing and updating work. So far, Biomapper development is a pleasure; I try to keep it so. Should it become a pain, I would probably have to find a way to get some money-compensation.
        However, I am always happy to answer the questions of the Biomapperians and do my best to fulfill their needs. All their questions and my answers are filed in this FAQ, which is therefore a good knowledge-base. This FAQ is updated regularly, each time I answer a new question, in fact. The help-file is more painful and time-consuming to update, and so I update it only about once a year, if there is enough matter to add and modify.
    Why is it not possible to use EGV maps having different scales, or not overlayable?
      The ENFA requires to have for every cell of the map a value for each EGV. If the maps are not overlayable, they would have to be made such internally. That would mean to include into Biomapper a complete interpolation algorithm. Although it would be possible, interpolation is a whole world in itself and I could not implement it exhaustively into Biomapper nor keep up-to-date with all the new development in the domain. My strategy is therefore to leave this aspect to the softwares dedicated to this problem. Idrisi for example has a whole geostatistic module allowing to do all kind of interpolations, from the simpler  to the most complex ones.
      The function of Biomapper is to go where other GIS softwares do not go, not to repeat what they are already doing. There is of course some overlap but I keep it minimal.

    Contact the webmaster
    [Top] [Home] [LBC] [UNIL]
    Updated 02.12.2002