R/findrow.R

Defines functions findrow

Documented in findrow

findrow <-
function(autocorr = 0.2, DIFFMAT=data$DIFF50, VERBOSE=FALSE) {

  #--------------------------------------------------------------
  # 
  # TITLE:     findrow()
  # AUTHOR:    TARMO REMMEL
  # DATE:      26 October 2016
  # CALLS:     N/A
  # CALLED BY: singlemap()
  # NEEDS:     LOOKUP TALBE DIFF50 FOR BIAS CORRECTION AND
  #            A RHO VALUE TO LOOK UP
  # NOTES:     USED TO FIND THE CLOSEST VALUE IN THE DIFF50
  #            LOOKUP TABLE MATRIX FOR RHO BIAS CORRECTION
  #            USED IN CONJUNCTION WITH findcol()
  #--------------------------------------------------------------

  # RUN ONLY IF THE AUTOCORRELATION IS > 0 and < 0.2499999
  if(autocorr > 0 & autocorr < 0.2499999) {

    minrow <- max(which(rownames(DIFFMAT) < autocorr))
    maxrow <- min(which(rownames(DIFFMAT) > autocorr))

    val11 <- abs(as.numeric(rownames(DIFFMAT)[minrow]) - autocorr)
    val22 <- abs(as.numeric(rownames(DIFFMAT)[maxrow]) - autocorr)

    if(VERBOSE) {
      cat("minrow: ", minrow, "\n", sep="")    
      cat("maxrow: ", maxrow, "\n", sep="")
      cat("val12: ", val11, "\n", sep="")    
      cat("val22: ", val22, "\n", sep="")
    } # END IF

    if(val11 < val22) row <- as.numeric(minrow) else row <- as.numeric(maxrow)
  } # END IF
  else {
    # A SPECIAL CASE, INDICATING THAT THE CORRECTION IS NEGLIGIBLE AT THIS LEVEL
    row <- 99
  } # END ELSE
  
  return(row)

} # END FUNCTION: findrow

Try the PatternClass package in your browser

Any scripts or data that you put into this service are public.

PatternClass documentation built on March 14, 2020, 1:07 a.m.