[R] Problem with script

[Carlos López]

Hello :-) I am trying to run the next script, it generates "random 
areas" inside a map of the american continent,
and then plot it, it´s suppose that every frame gives you the evolution 
of the program but at some point it stops
with the weirdest of the errors I´ve ever seen in R, I don´t even have a 
line like that :-|

Error in if (random_matrix[x_list[random_sq] - 1, y_list[random_sq]] ==  :
  missing value where TRUE/FALSE needed


the script is this:

# ______________

                                        # ----- Define a function for 
plotting a matrix ----- #
myImagePlot <- function(x, ...){
  min <- min(x)
  max <- max(x)
  yLabels <- rownames(x)
  xLabels <- colnames(x)
  title <-c()
                                        # check for additional function 
arguments
  if( length(list(...)) ){
    Lst <- list(...)
    if( !is.null(Lst$zlim) ){
      min <- Lst$zlim[1]
      max <- Lst$zlim[2]
    }
    if( !is.null(Lst$yLabels) ){
      yLabels <- c(Lst$yLabels)
    }
    if( !is.null(Lst$xLabels) ){
      xLabels <- c(Lst$xLabels)
    }
    if( !is.null(Lst$title) ){
      title <- Lst$title
    }
  }
                                        # check for null values
  if( is.null(xLabels) ){
    xLabels <- c(1:ncol(x))
  }
  if( is.null(yLabels) ){
    yLabels <- c(1:nrow(x))
  }
 
  layout(matrix(data=c(1,2), nrow=1, ncol=2), widths=c(4,1), heights=c(1,1))
 
                                        # Red and green range from 0 to 
1 while Blue ranges from 1 to 0
  ColorRamp <- rgb( seq(0,1,length=256),  # Red
                   seq(0,1,length=256),  # Green
                   seq(1,0,length=256))  # Blue
  ColorLevels <- seq(min, max, length=length(ColorRamp))
 
                                        # Reverse Y axis
  reverse <- nrow(x) : 1
  yLabels <- yLabels[reverse]
  x <- x[reverse,]
 
                                        # Data Map
  par(mar = c(3,5,2.5,2))
  image(1:length(xLabels), 1:length(yLabels), t(x), col=ColorRamp, xlab="",
        ylab="", axes=FALSE, zlim=c(min,max))
  if( !is.null(title) ){
    title(main=title)
  }
  axis(BELOW<-1, at=1:length(xLabels), labels=xLabels, cex.axis=0.7)
  axis(LEFT <-2, at=1:length(yLabels), labels=yLabels, las= HORIZONTAL<-1,
       cex.axis=0.7)
 
                                        # Color Scale
  par(mar = c(3,2.5,2.5,2))
  image(1, ColorLevels,
        matrix(data=ColorLevels, ncol=length(ColorLevels),nrow=1),
        col=ColorRamp,
        xlab="",ylab="",
        xaxt="n")
 
  layout(1)
}
                                        # ----- END plot function ----- #

                                        # Read the table to create the 
matrix needed

setwd("/home/natorro/Batmaps/")
mexico_matrix <- read.table("/home/natorro/Batmaps/mexico_matrix.dat", 
header=TRUE)

                                        # Let's generate a matrix full 
of zeroes
random_matrix <- matrix(0, nrow=175, ncol=175)

                                        # Size of range (area of 
distribution)
range_size <- 100

                                        # Let's choose initial point

initial_point_passed <- FALSE
initial_point <- sample(175, 2)    

while (initial_point_passed == FALSE){
  if (mexico_matrix[initial_point[1], initial_point[2]] == 0)
    initial_point <- sample(175, 2) else {
      random_matrix[initial_point[1], initial_point[2]] <- 1
      initial_point_passed <- TRUE
    }
}

                                        # Define vectors x_list and y_list

x_list <- matrix(0, nrow=range_size, ncol=1)
y_list <- matrix(0, nrow=range_size, ncol=1)

x_list[1] <- initial_point[1]
y_list[1] <- initial_point[2]


                                        # Range counter
range_counter <- 1
flag <- 0

while (range_counter < range_size){
  random_sq <- round((runif(1) * range_counter ) + 1)
  flag <- 0    
  while (flag == 0) {
    
    if(random_matrix[x_list[random_sq]-1, y_list[random_sq]]== 0 ||
       random_matrix[x_list[random_sq] + 1, y_list[random_sq]] == 0 ||
       random_matrix[x_list[random_sq], y_list[random_sq] + 1] == 0 ||
       random_matrix[x_list[random_sq], y_list[random_sq] - 1] == 0) {
      
      aleatorio <- runif(1)
      if (aleatorio < 1) {
        if (random_matrix[x_list[random_sq],y_list[random_sq]-1]== 0 ) {
          random_matrix[x_list[random_sq], y_list[random_sq] - 1] <- 1
          flag <- 1
          range_counter <- range_counter + 1
          x_list[range_counter] <- x_list[random_sq]
          y_list[range_counter] <- y_list[random_sq] - 1
          myImagePlot(mexico_matrix + random_matrix)
        }
      }
      
      if (aleatorio < 0.75){
        if (random_matrix[x_list[random_sq],y_list[random_sq] + 1] == 0) {
          random_matrix[x_list[random_sq], y_list[random_sq] + 1] <- 1
          flag <- 1
          range_counter <- range_counter + 1
          x_list[range_counter] <- x_list[random_sq]
          y_list[range_counter] <- y_list[random_sq] + 1
          myImagePlot(mexico_matrix + random_matrix)
        }
      }
      
      
      if (aleatorio < 0.50){
        if(random_matrix[x_list[random_sq]+1, y_list[random_sq]] == 0) {
          random_matrix[x_list[random_sq] + 1, y_list[random_sq]] <- 1
          flag <- 1
          range_counter <- range_counter + 1
          x_list[range_counter] <- x_list[random_sq] + 1
          y_list[range_counter] <- y_list[random_sq]
          myImagePlot(mexico_matrix + random_matrix)
        }
      }
      
      if (aleatorio < 0.25) {
        if(random_matrix[x_list[random_sq]-1, y_list[random_sq]] == 0 ) {
          random_matrix[x_list[random_sq] - 1, y_list[random_sq]] <- 1
          flag <- 1
          range_counter <- range_counter + 1
          x_list[range_counter] <- x_list[random_sq] - 1
          y_list[range_counter] <- y_list[random_sq]
          myImagePlot(mexico_matrix + random_matrix)
        }
      }
      
      
    }
  }
}

sum(random_matrix)


#-------------------------------


I do the "sum(random_matrix)" just to check out if the number of zeroes 
needed have been set (never has reached 100 so far :-()

I´m not sure what I am doing wrong, I´ve been working with this file in 
windows and mac at the same time in a colaboration.
Any help will be greatly appreciated.

Carlos


-- 
This message has been scanned for viruses and
dangerous content by MailScanner, and is
believed to be clean.