R/addReleve.R

In Rveg: Digitization of Phytosociological Relevés

#'
#' addReleve
#'
#' Digitizing and editing your releves - Rveg database
#' @param DATABASE name of csv files for releve table and header - database
#' @param SAVE name of exporting database
#' @param checklist custom checklist
#' @param extrahead extra rows in header
#' @param start Boolean to start immediately digitizing first releve
#'
#' @returns export two csv files, one for releve and one for header (Rveg database)
#'
#' @examples
#' ## NOT RUN
#' if (interactive()) {
#'   addReleve()
#' }
#'
#' @export
#'
#'

addReleve <- function(DATABASE = "NEW", SAVE = "default", checklist = "default", extrahead = NULL, start = TRUE) {
  # LOAD DATA ------------------------------------------------------------------------------------------

  if (checklist == "default") {
    checklist <- paste0(path.package("Rveg"), "/extdata/DANIHELKA2012rko.txt")
  } # Default checklist
  if (SAVE == "default") {
    SAVE <- file.path(tempdir(), "default") # tempdir
  }

  if (DATABASE == "NEW") {
    DATA <- data.frame(ShortName = character(), stringsAsFactors = FALSE)
    HeaderDATA <- data.frame(
      ShortName = c(
        "ID", "DATE", "SpringDATE", "LOCALITY", "FieldCODE", "Authors",
        "PlotSize", "Latitude", "Longitude", "Accuracy", "CRS", "Slope", "Exposure",
        "E3", "E2", "E1", "Ejuv", "E0", "Note", extrahead
      ),
      stringsAsFactors = FALSE
    ) # create header releve and header template
    write.csv(DATA, paste0(SAVE, "REL.csv"))
    write.csv(HeaderDATA, paste0(SAVE, "HEAD.csv")) # immediate write on disk
  } else {
    DATARELE <- paste0(DATABASE, "REL.csv")
    DATAhead <- paste0(DATABASE, "HEAD.csv")
    DATA <- read.csv(DATARELE, row.names = 1)
    HeaderDATA <- read.csv(DATAhead, row.names = 1) # reading existin db

    write.csv(DATA, paste0(SAVE, "REL.csv"))
    write.csv(HeaderDATA, paste0(SAVE, "HEAD.csv")) # immediate write on disk
  }

  SpLIST1 <- read.delim(checklist, sep = "\t")
  SpLIST <- makeSpLIST(checklist) # creating Species checklist

  # DATA SETUP  ------------------------------------------------------------------------------------------

  while (TRUE) {
    DATA2 <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1) #
    HeaderDATA2 <- read.csv(paste0(SAVE, "HEAD.csv"), row.names = 1) #
    if (DATABASE == "NEW"  & start == TRUE) {
      Header <-
        data.frame(
          ShortName = c(
            "ID", "DATE", "SpringDATE", "LOCALITY", "FieldCODE", "Authors", "PlotSize",
            "Latitude", "Longitude", "Accuracy", "CRS", "Slope", "Exposure", "E3", "E2", "E1", "Ejuv", "E0", "Note", extrahead
          ),
          Value = 0
        )
      RelNew <- data.frame(ShortName = SpLIST[, 2], Value = 0) # new empty rel
      ID <- ncol(DATA2)
      colnames(RelNew)[2] <- (ID - 1)
      Header[1, 2] <- (1)

      ## filling the header
      ab <- readline("DATE?(YYYY/MM/DD) ")
      bb <- readline("SPRINGDATE?(YYYY/MM/DD) ")
      bc <- readline("LOCALITY? ")
      a <- readline("FieldCODE? ")
      b <- readline("Authors? ")
      c <- readline("PlotSize?(m2) ")
      d <- readline("Latitude? ")
      e <- readline("Longitude? ")
      f <- readline("Accuracy? ")
      crs <- readline("Coordinta Reference System? ")
      ff <- readline("Slope?(degrees) ")
      fk <- readline("Exposure? ")
      g <- readline("E3?(%) ")
      h <- readline("E2?(%) ")
      i <- readline("E1?(%) ")
      j <- readline("Ejuv?(%) ")
      k <- readline("E0?(%) ")
      l <- readline("Note? ")

      extraval <- NULL
      if (!is.null(extrahead)) {
        for (val in extrahead) {
          assign(val, readline(paste0(val, "? ")))
          extraval <- c(extraval, eval(as.symbol(val)))
        }
      }
      hh <- c(ab, bb, bc, a, b, c, d, e, f, crs, ff, fk, g, h, i, j, k, l, extraval)
      Header[2:(length(hh) + 1), 2] <- hh
      colnames(Header)[2] <- (1)
      HeaderDATA2 <- data.frame(HeaderDATA2, Header[, 2])
      colnames(HeaderDATA2)[2] <- "X1"
      write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))

      aa <- "NEW" # menu options
      DATABASE <- ""
    } else {
      message(paste0("number of relev\u00e9s: ", ncol(DATA2)-1 , ", number of headers: ", ncol(HeaderDATA2) -1 ))
      aa <-
        toupper(readline(
          "AddReleve?(Y/N/RREL/RHEAD/REMOVEREL/PRINTREL/PRINTHEAD/ADDREL/ADDHEAD/YSP) "
        ))
    }
    if ((aa == "Y" | aa == "YSP") & ncol(DATA2) != ncol(HeaderDATA2)) {
      message("Number of eelev\u00e9s:s and headers must match!!!")
    } ## warning for mismatch

    if (aa == "NEW" | aa == "RREL" | (aa == "Y" & ncol(DATA2) == ncol(HeaderDATA2)) | aa == "ADDREL") {
      if (aa == "Y") {
        # Header <-
        #   data.frame(
        #     ShortName = c(
        #       "ID", "DATE", "SpringDATE", "LOCALITY", "FieldCODE", "Authors", "PlotSize", "Latitude", "Longitude", "Accuracy", "Slope", "Exposure",
        #       "E3", "E2", "E1", "Ejuv", "E0", "Note", extrahead
        #     ),
        #     Value = 0
        #   )

        Header <- data.frame(ShortName = rownames(HeaderDATA2),
                             Value = 0)

        RelNew <- data.frame(ShortName = SpLIST[, 2], Value = 0)
        ID <- ncol(DATA2)
        colnames(RelNew)[2] <- (ID - 1)
        #Header[1, 2] <- (length(colnames(HeaderDATA2)))

        Header[1, 2] <- suppressWarnings(as.numeric(HeaderDATA2[1,ncol(HeaderDATA2)])+1)

        if (is.na(Header[1, 2])) {
          Header[1, 2] <- 1
        }

        ## hinting previous prompts
        hhp <- HeaderDATA2[-1,ncol(HeaderDATA2)]

        message(paste0("RE <- ",hhp[1]))
        ab <- readline("DATE?(YYYY/MM/DD) ")
        message(paste0("RE <- ",hhp[2]))
        bb <- readline("SPRINGDATE?(YYYY/MM/DD) ")
        message(paste0("RE <- ",hhp[3]))
        bc <- readline("LOCALITY? ")
        message(paste0("RE <- ",hhp[4]))
        a <- readline("FieldCODE? ")
        message(paste0("RE <- ",hhp[5]))
        b <- readline("Authors? ")
        message(paste0("RE <- ",hhp[6]))
        c <- readline("PlotSize?(m2) ")
        message(paste0("RE <- ",hhp[7]))
        d <- readline("Latitude? ")
        message(paste0("RE <- ",hhp[8]))
        e <- readline("Longitude? ")
        message(paste0("RE <- ",hhp[9]))
        f <- readline("Accuracy?(m) ")
        message(paste0("RE <- ",hhp[10]))
        crs <- readline("Coordinta Reference System? ")
        message(paste0("RE <- ",hhp[11]))
        ff <- readline("Slope?(degrees) ")
        message(paste0("RE <- ",hhp[12]))
        fk <- readline("Exposure? ")
        message(paste0("RE <- ",hhp[13]))
        g <- readline("E3?(%) ")
        message(paste0("RE <- ",hhp[14]))
        h <- readline("E2?(%) ")
        message(paste0("RE <- ",hhp[15]))
        i <- readline("E1?(%) ")
        message(paste0("RE <- ",hhp[16]))
        j <- readline("Ejuv?(%) ")
        message(paste0("RE <- ",hhp[17]))
        k <- readline("E0?(%) ")
        message(paste0("RE <- ",hhp[18]))
        l <- readline("Note? ")

        if (length(hhp)>18) {
          extrahead <- c()
          for (ext in 1:(length(hhp)-18)) {
            extrahead[ext] <- HeaderDATA2[ext+19,1]
          }
        }

        extraval <- NULL
        if (!is.null(extrahead)) {
          ehcounter <- 1 # extrahead counter
          for (val in extrahead) {
            message(paste0("RE <- ",hhp[18+ehcounter]))
            ehcounter <- ehcounter + 1
            assign(val, readline(paste0(val, "? ")))
            extraval <- c(extraval, eval(as.symbol(val)))
          }
        }



        hh <- c(ab, bb, bc, a, b, c, d, e, f, crs, ff, fk, g, h, i, j, k, l, extraval)

        for (hval in 1:length(hh)) {
          if (toupper(hh[hval]) == "RE") {
            hh[hval] <- hhp[hval]
          }
        }

        Header[2:(length(hh) + 1), 2] <- hh
        colnames(Header)[2] <- (ID - 1)
        HeaderDATA2 <- data.frame(HeaderDATA2, Header[, 2])
        colnames(HeaderDATA2)[2:length(colnames(HeaderDATA2))] <- paste0("X",c(1:(length(colnames(HeaderDATA2)) - 1)))
        write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))
      }
      if (aa == "RREL") {
        DATA2 <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1)
        DATA2 <- DATA2[order(as.numeric(row.names(DATA2))), ]
        while (TRUE) {
          n <- as.numeric(readline("ReleveNumber? ")) # releve to repair
          if (is.na(n) != TRUE) {
            HeaderDATA3 <- read.csv(paste0(SAVE, "HEAD.csv"), row.names = 1)
            print(HeaderDATA3[, n + 1])
            tt <- toupper(readline("CorrectNumber?(Y/N) ")) # double check
            if (tt == "Y") {
              break
            }
          }
        }

        ID <- n + 1
        RelNew <- data.frame(ShortName = SpLIST[, 2], Value = 0)

        # DATA2 <- DATA2[order(as.numeric(row.names(DATA2))),]
        RelNew <- RelNew[order(RelNew[, 1]), ] #
        DATA2 <- DATA2[order(DATA2[, 1]), ] # testing the shortname ordering
        RelNew[RelNew[, 1] %in% DATA2[, 1], ][, 2] <- DATA2[, ID]
        TABLEexp <- RelNew[RelNew[, 2] > 0, ]
        colnames(TABLEexp) <- c("ShortNames", "Cover")
        print(TABLEexp)
      }
      if (aa == "ADDREL") {
        RelNew <- data.frame(ShortName = SpLIST[, 2], Value = 0)
        ID <- ncol(DATA2)
        colnames(RelNew)[2] <- (ID - 1)
        }
      while (TRUE) {
        m <- toupper(readline("AddNewLayer?(Y/N) "))
        if (m == "Y" | m == "N") {
          m <- m
        }

        if (m == "Y") {
          while (TRUE) {
            most <- toupper(readline("Select Layer (3,2,1,J,0) "))
            if (most == 3 | most == 2 | most == 1 | most == "J" | most == 0) {
              break
            }
          }
          while (TRUE) {
            oo <- toupper(readline("P - percentage, BB - Braun B. scale, CS - custom scale "))

            if (oo == "P" | oo == "BB" | oo == "B" | oo == "CS") {
              break
            }
          }
          while (TRUE) {
            m <- toupper(readline("AddSpecies?(GenuSpe/N) "))

            if (nchar(m) == 7) {
              n <- m

              if (oo == "P") {
                while (TRUE) {
                  o <- suppressWarnings(as.numeric(readline("Abundance?(%) ")))
                  if (!is.na(o) && o >= 0 && o <= 100) {break}
                  }
              } else if (oo == "CS") {
                o <- as.character(readline("Abundance?(%) "))
              } else if (oo == "BB" | oo == "B") {
                while (TRUE) {
                  o <- toupper(readline("Abundance?(0,R,+,1,2,M,A,B,3,4,5) "))
                  if (o == "R" | o == "+" | o == "0" | o == "1" | o == "2" | o == "M" | o == "A" | o == "B" |
                    o == "3" | o == "4" | o == "5") {
                    break
                  }
                }
                if (o == "R") {
                  o <- 1
                } else if (o == "+") {
                  o <- 2
                } else if (o == "1") {
                  o <- 3
                } else if (o == "2") {
                  o <- 15
                } else if (o == "M") {
                  o <- 4
                } else if (o == "A") {
                  o <- 8
                } else if (o == "B") {
                  o <- 18
                } else if (o == "3") {
                  o <- 38
                } else if (o == "4") {
                  o <- 63
                } else if (o == "5") {
                  o <- 88
                } else if (o == "0") {
                  o <- 0
                }
              }
              nana <- paste(n, most, sep = "_") # creating the species name
              specie_check <- SpLIST[SpLIST[, 2] == substr(nana, 1, 9), ] # double check
              print(specie_check[3])
              while (TRUE) {
                tt <- toupper(readline("CorrectName?(Y/F(search for name)) ")) # double check
                if (tt == "Y" && !is.na(specie_check[1, 3]) | tt == "F") {
                  break
                }
              }

              if (tt == "F") {
                while (TRUE) {
                  k <- readline("SpeciesFirst3letters?(eg.Che) ") # search correct name in original list
                  if (nchar(k) >= 3) {
                    k <- paste0(toupper(substr(k, 1, 1)), tolower(substr(k, 2, 3)))
                    searchlist <- SpLIST1[grep(paste0("^", k), SpLIST1[, 3]), ]
                    print(searchlist[order(searchlist$FullName), ])
                    while (TRUE) {
                      while (TRUE) {
                        s <- toupper(readline("SpeciesName?(GenuSpe) "))
                        if (nchar(s) == 7) {
                          nana <- paste(s, most, sep = "_")
                          specie_check <- SpLIST[SpLIST[, 2] == substr(nana, 1, 9), ]
                          print(specie_check[3])
                          break
                        }
                      }

                      tt <- toupper(readline("CorrectSpecies?(Y/N)  "))
                      if ((tt == "Y") & !is.na(specie_check[1, 3])) {
                        break
                      }
                    }
                    break
                  }
                }



                nanas <- paste(s, most, sep = "_")

                RelNew[which(RelNew$ShortName == substr(nanas, 1, 9)), ][, 2] <- o
                print(RelNew[(RelNew[, 2] > 0), ])
              } else {
                RelNew[which(RelNew$ShortName == substr(nana, 1, 9)), ][, 2] <- o
                print(RelNew[(RelNew[, 2] > 0), ])
              }
            } else if (m == "N" | m == "n") {
              break
            }
          }
        } else if (m == "N" | m == "n") {
          message("Species_richness")
          print(nrow(RelNew[(RelNew[, 2] > 0), ]))
          if (aa == "NEW" | !isTRUE(start)) {
            DATA2 <- RelNew[(RelNew[, 2] > 0), ] # first releve?
            colnames(DATA2)[2] <- "X1"
            write.csv(DATA2, paste0(SAVE, "REL.csv"))
            start = TRUE
            DATABASE = ""
            break
          }
          if (aa == "Y" | aa == "ADDREL") {
            DATA2 <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1) # New releve
            DATA2 <- createTABLE(SpLIST, RelNew, DATA2)
            colnames(DATA2)[-1] <- paste0("X",1:(length(colnames(DATA2))-1))
            write.csv(DATA2, paste0(SAVE, "REL.csv"))
            break
          }
          if (aa == "RREL") {
            DATA2 <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1) # Repair releve
            DATA2 <- createTABLE(SpLIST, RelNew, DATA2)
            DATA2 <- DATA2[, -c(ID)]
            colnames(DATA2)[length(colnames(DATA2))] <- ID - 1
            DATA2 <- DATA2[, c(colnames(DATA2)[1], 1:(length(colnames(DATA2)) - 1))]
            colnames(DATA2)[-1] <- paste0("X",1:(length(colnames(DATA2)) - 1)) #testX
            write.csv(DATA2, paste0(SAVE, "REL.csv"))
            break
          }

        }
      }
    }

    if (aa == "PRINTHEAD") {
      HeaderDATA3 <- read.csv(paste0(SAVE, "HEAD.csv"), row.names = 1) # actual headers
      print(HeaderDATA3)
    }

    if (aa == "PRINTREL") {
      DATA3 <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1) # actual releves
      print(DATA3)
    }

    if (aa == "RHEAD") {
      HeaderDATA2 <- read.csv(paste0(SAVE, "HEAD.csv"), row.names = 1) #

      while (TRUE) {
        n <- as.numeric(readline("ReleveNumber? "))
        if (is.numeric(n) & is.na(n) == FALSE) {
          if (n <= (ncol(HeaderDATA2) - 1) & n > 0) {
            print(HeaderDATA2[, colnames(HeaderDATA2) == paste0("X", n)])

            tt <- toupper(readline("CorrectColumn?(Y/F) ")) # double check
            if (tt == "Y") {
              break
            }
          }
        }
      }
      while (TRUE) {
        m <- toupper(readline("RepairHeader?(Y/N) "))
        if (m == "Y") {
          print(HeaderDATA2[, 1]) # print header
          while (TRUE) {
            l <- readline("HeaderCharacteristic? ") # choosing the characterstic

            if (l %in% HeaderDATA2[, 1]) {
              break # if exists
            }
          }
          o <- readline("NewValue? ")
          HeaderDATA2[which(HeaderDATA2[, 1] == l), ][, n + 1] <- o # replace
        } else if (m == "N" | m == "n") {
          break
        }
      }
      write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))
    }

    if (aa == "REMOVEREL") {
      DATA2 <-
        read.csv(paste0(SAVE, "REL.csv"), row.names = 1) # rownames
      while (TRUE) {
        n <- readline("ReleveNumber? ")
        n <- as.numeric(n)
        print(HeaderDATA2[, colnames(HeaderDATA2) == paste0("X", n)]) # selection or releve
        tt <- toupper(readline("CorrectNumber?(Y/N) ")) # double check
        if (tt == "Y") {
          DATA2 <- DATA2[, colnames(DATA2) != paste0("X", n),drop = FALSE] # removal of rel
          #HeaderDATA2 <- HeaderDATA2[, HeaderDATA2[1, ] != n] # removal of header
          HeaderDATA2 <- HeaderDATA2[, colnames(HeaderDATA2) != paste0("X", n),drop = FALSE] # removal of header
          colnames(DATA2)[-1] <- paste0("X",1:(length(colnames(DATA2)) - 1)) # reindexing
          #colnames(HeaderDATA2)[2:length(colnames(HeaderDATA2))] <- paste0("X",c(1:(length(colnames(HeaderDATA2)) - 1))) # header rein
          colnames(HeaderDATA2)[-1] <- paste0("X",1:(length(colnames(HeaderDATA2)) - 1))

          write.csv(DATA2, paste0(SAVE, "REL.csv"))
          write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))
          break
        }
      }
    }

    if (aa == "YSP" & ncol(DATA2) == ncol(HeaderDATA2)) {
      while (TRUE) {
        m <- suppressWarnings(as.numeric(readline("How many releves? ")))
        if (!is.na(m) & m>0) {

          Rels <- list()
          DATAtemp <- read.csv(paste0(SAVE, "REL.csv"), row.names = 1) # New releve
          for (i in 1:m) {
            Rels[[paste0("r",i)]] <- data.frame(ShortName = SpLIST[, 2],value = 0)
          }

          break
        }

      } ## how many releves
      while (TRUE) {
        oo <- toupper(readline("P - percentage, BB - Braun B. scale, CS - custom scale "))

        if (oo == "P" | oo == "BB" | oo == "B" | oo == "CS") {
          break
        }
      }
      while (TRUE) {
        n <-toupper(readline("Add new species (GenuSpe)/N?  "))
        if (nchar(n) == 7) {
          while (TRUE){
            l <- toupper(readline("Select Layer (3,2,1,J,0) " ))
            if (l == "1" | l == "2" | l == "3" | l == "J" | l == "0" ) {
              break
            }

          } ## what layer
          nana <- paste(n, l, sep = "_") # creating the species name
          specie_check <- SpLIST[SpLIST[, 2] == substr(nana, 1, 9), ] # double check
          print(specie_check[3])
          while (TRUE) {
            tt <- toupper(readline("CorrectName?(Y/F(search for name)) ")) # double check
            if (tt == "Y" && !is.na(specie_check[1, 3]) | tt == "FFF") {
              break
            } else if (tt == "F") {
              while (TRUE) {
                k <- readline("SpeciesFirst3letters?(eg.Che) ") # search correct name in original list
                if (nchar(k) >= 3) {
                  k <- paste0(toupper(substr(k, 1, 1)), tolower(substr(k, 2, 3)))
                  searchlist <- SpLIST1[grep(paste0("^", k), SpLIST1[, 3]), ]
                  print(searchlist[order(searchlist$FullName), ])
                  while (TRUE) {
                    while (TRUE) {
                      s <- toupper(readline("SpeciesName?(GenuSpe) "))
                      if (nchar(s) == 7) {
                        nana <- paste(s, l, sep = "_")
                        specie_check <- SpLIST[SpLIST[, 2] == substr(nana, 1, 9), ]
                        print(specie_check[3])
                        break
                      }
                    }

                    tt <- toupper(readline("CorrectSpecies?(Y/N)  "))
                    if ((tt == "Y") & !is.na(specie_check[1, 3])) {
                      break
                    }
                  }
                  break
                }
              }

              # nanas <- paste(s, most, sep = "_")
              #
              # RelNew[which(RelNew$ShortName == substr(nanas, 1, 9)), ][, 2] <- o
              # print(RelNew[(RelNew[, 2] > 0), ])

            }
          } # correct name

          RelFake <- data.frame(ShortName = SpLIST[, 2])
          for (i in 1:m) {
            RelFake[[paste0("r",i)]] <- 0
          }

          for (i in 1:m) {
            message(paste0("Relev\u00e9 ",i))

            if (oo == "P") {
              while (TRUE) {
                o <- suppressWarnings(as.numeric(readline("Abundance?(%) ")))
                if (!is.na(o) && o >= 0 && o <= 100) {break}
              }
            } else if (oo == "CS") {
              o <- as.character(readline("Abundance?(%) "))
            } else if (oo == "BB" | oo == "B") {
              while (TRUE) {
                o <- toupper(readline("Abundance?(0,R,+,1,2,M,A,B,3,4,5) "))
                if (o == "R" | o == "+" | o == "0" | o == "1" | o == "2" | o == "M" | o == "A" | o == "B" |
                    o == "3" | o == "4" | o == "5") {
                  break
                }
              }
              if (o == "R") {
                o <- 1
              } else if (o == "+") {
                o <- 2
              } else if (o == "1") {
                o <- 3
              } else if (o == "2") {
                o <- 15
              } else if (o == "M") {
                o <- 4
              } else if (o == "A") {
                o <- 8
              } else if (o == "B") {
                o <- 18
              } else if (o == "3") {
                o <- 38
              } else if (o == "4") {
                o <- 63
              } else if (o == "5") {
                o <- 88
              } else if (o == "0") {
                o <- 0
              }
            }

            Rels[[paste0("r",i)]][which(Rels[[paste0("r",i)]]$ShortName == substr(nana, 1, 9)), ][, 2] <- o
            RelFake[which(RelFake$ShortName == substr(nana, 1, 9)), ][, i+1] <- o
            print(RelFake[(RelFake[,i+1] > 0), ])

          }
          DATA2 <- DATAtemp
          DATA2 <- createTABLE(SpLIST, Rels, DATA2,variation = 2)
          print(DATA2)
          colnames(DATA2)[-1] <- paste0("X",1:(length(colnames(DATA2))-1))
          write.csv(DATA2, paste0(SAVE, "REL.csv"))
          start = TRUE
          DATABASE = ""


        } else if (n == "N") {
          while (TRUE) {
            rs <- toupper(readline("Add headers?(Y/N) "))
            if (rs == "Y") {
              for (i in 1:m) {
                message(paste0("Relev\u00e9 ",i))
                header <- createHEADER(HeaderDATA2)
                HeaderDATA2 <- data.frame(HeaderDATA2, header[, 2])
                colnames(HeaderDATA2)[2:length(colnames(HeaderDATA2))] <- paste0("X",c(1:(length(colnames(HeaderDATA2)) - 1)))
                write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))
                }


              break
            } else if (rs == "N") {

              break
            }

          }
          break

        }

      } ## add new species
    }

    if (aa == "ADDHEAD") {
      header <- createHEADER(HeaderDATA2)
      HeaderDATA2 <- data.frame(HeaderDATA2, header[, 2])
      colnames(HeaderDATA2)[2:length(colnames(HeaderDATA2))] <- paste0("X",c(1:(length(colnames(HeaderDATA2)) - 1)))
      write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv"))
    }

    if (aa == "N") {
      write.csv(DATA2, paste0(SAVE, "REL.csv"))
      write.csv(HeaderDATA2, paste0(SAVE, "HEAD.csv")) # the script ends, database is saved
      break
    }

  }
}