Standing Ovation Problem - v1

globals [
  ;; these are defined in the interface sliders
  ;; error-max    ;;  "errors" ( diversity or E ) will be on [-Error-max, Error-max],
                  ;;  with the largest range used being [-1,1]
  ;; population   ;;  number of people in the theater
  ;; quality      ;;  absolute quality of the show, ie Q
  ;; error-size   ;;
  ;; threshold    ;; threshold signal that causes agent to stand up if randomization is not turned on
                  ;; where signal = Quality + Error
  ;;
  ;;
  pass              ;; which pass, 1 or 2
  T                 ;; threshold for standing on pass 2, same for everyone
  ;; quality        ;; actual quality of the perforance ( ie Q )

  sitting-color     ;;
  standing-color    ;; stand up on pass one
  persuaded-color   ;; stand up on pass two
  percent-standing       ;; count of patrons standing at end of pass 1
  percent-standing-final ;; percent standing at end of pass 2
  percent-persuaded      ;; percent persuaded by peer-pressure in pass 2
  standing-ovation?      ;; true or false
  verbose?               ;; set to true for debugging
]

turtles-own [
  my-signal-threshold   ;;        this is an individual threshold for comparing to perceived quality ( S = Q + error )
  my-peer-pressure-threshold ;;    % of audience standing that will cause  me to stand on pass two, if I am not already standing from pass 1
  my-error       ;;  chosen from a random distribution from -E to E
  my-signal      ;;  after a play, compute this as actual-quality plus my-error, so ranges from -E to (100 + E)
  standing?      ;;  initially false
]

to setup
  clear-all
  set verbose? false
  set standing-ovation? false
  set sitting-color grey
  set standing-color green
  set persuaded-color yellow

  ;; we have a performance, set the quality
  if random-Q? [ set quality random 101 ]

  create-turtles population [ setxy random-xcor random-ycor  set size 3 set shape "person" set color blue]
  ask turtles [

    ;; errors are always randomized, but the error-max can be slid to zero
    set my-error uniform error-max

    set my-signal-threshold threshold ;;
    if random-T?  [set my-signal-threshold random 101 ]  ;; hmmm, question - does it make sense for both  S and T bo be randomized, and if so,
                                                         ;; is random-T chosen on U[0,100] or some other distribution centered on global threshold?

    set my-signal -999  ;;  Set after the show ends to Q + E

    set my-peer-pressure-threshold x-percent
    if random-X? [ set my-peer-pressure-threshold random 101 ]




    ;; they are sitting so
    set standing? false
    set color sitting-color
  ]
  print ( word population " patrons were seated.")
  reset-ticks ;;  ticks starts at 0 not 1
end 

to-report uniform [eee]
     let ranger (2 * eee )
     let xc random (ranger + 1)  ;; x will be an integer between 0 and 2E
     set xc ( xc - eee )         ;; now x is between -E and E inclusive
     report xc
end 

to go
  set pass ( ticks + 1 )     ;; because ticks start at 0
  if ( pass > 2 ) [ stop ]  ;; two passes is all that is needed to run the model

  if verbose? [ print (word "this is pass #" pass) ]
  if (ticks = 0) [
    print (word  "The show quality was " quality )
   ]  ;; the show happens

  if ( pass = 1 )
 [ ask turtles   [
    if verbose? [print ( word "Patron " who " is doing the first pass with show Q = " quality )]


       ;; compute my signal, that is, how I perceive the show
        ;;;;;;;;;;;;;;;;;;;;;;;;;set my-error uniform error-max

       set my-signal (quality + my-error) ;; the play's quality as I see it


    if verbose? [    print ( word ".... my-error = " my-error ", my-signal = " my-signal ", and my-signal-threshold = " my-signal-threshold)]

      ;; then I stand up if my-signal is over my-signal-threshold

       if-else ( my-signal > my-signal-threshold ) [
         if verbose? [  print "....standing up" ]
         set standing? true
         set color standing-color
        ]
        [
      if verbose? [ print ".... remaining sitting."]
    ]
    if verbose? [ print " " ]
  ]
    set percent-standing  ( 100 * count turtles with [ standing? = true ]) / ( count turtles )
    print (word "At the end of the first pass,  percent standing = " precision percent-standing 0)
    if verbose? [print " "]
  ]  ;; end of pass 1

   if (pass = 2 )
   [


    if verbose? [print (word "Starting second pass, with percent standing after passs 1 = " percent-standing)]
   ;; if ( percent-standing > X-percent )[ ask turtles with [standing? = false]  ;;  in version 1, this is wrong, tests global setting

   ask turtles with [standing? = false]  ;;  in version 2, this is better
      [
          if ( percent-standing > my-peer-pressure-threshold )[                      ;; this tests each individual person's threshhold
               set standing? true
               if verbose? [print (word "Patron " who " now stands up")]
               set color persuaded-color ]
      ]

    ;; ok now count again
    set percent-standing-final 100 * ( count turtles with [ standing? = true ]) / ( count turtles )
    set percent-persuaded 100 * count turtles with [color = persuaded-color] / (count turtles )
    if-else  percent-standing-final = 100 [ set standing-ovation? true ][ set standing-ovation? false]

    wrap-up-show
  ] ;; end of pass 2

  tick
end 

to wrap-up-show
    print " ========================================="
    print ( word " The show quality was " quality )
    print ( word " At the end of the first pass,  percent standing  : " precision percent-standing 0)
    print ( word " Then, this percent was also persuaded to stand up: " precision percent-persuaded 0)
    print ( word " So, the final total percent standing             : " precision percent-standing-final 0)
    if-else ( standing-ovation?  )
          [ print "   Standing ovation!"]
          [ print "   Did not get a full standing ovation"]

  ;;  print " ========================================="
  ;;  print " update graphs and output here "

  let show-text (word "For Q=" quality ", T=" threshold ", E="  error-max ", X=" X-percent ", %pass-1=" precision percent-standing 0 ", %pass-2="  precision percent-persuaded 0 " %total=" precision percent-standing-final 0)                                   ;;////////////////////////////////////////////////////////////////////////////
  print show-text
  ;;  print " ========================================="

  if population = 1 [ display-visual ]
end 

to display-visual ;; for only one turtle


  let my-T 0
  let my-E 0
  let my-S 0
  let my-Q quality

  ask turtle 0 [ set hidden? true
  set my-T my-signal-threshold
  set my-E my-error
  set my-S my-signal
  ]

  analyze-one-Emax error-max my-Q   my-T my-S  ;; this is the one to come home to
end 

to analyze-one-Emax [ amax my-Q   my-T my-S ]

  ask patches [ set pcolor gray]
  ask patches with [ pxcor > -10 and pxcor < 10 and pycor < 3 and pycor > -3 ] [ set pcolor white ] ;; 21 patches wide, each about 5 points
  let white-count count patches with [ pcolor = white ]

  ;; make a vertical black bar at T
  let plot-T ( my-T / 5)  - 10
  ask patches with [ pxcor = round plot-T and pycor < 3 and pycor > -3 ] [ set pcolor black ]

  ;; make a blue short vertical bar at Q
  let plot-Q ( my-Q / 5) - 10 ;
  ask patches with [pxcor = round plot-Q and pycor > -1 and pycor < 6 ] [ set pcolor blue ]

  ;; make a red short vertical bar at S ( what the dice actually rolled )
  let plot-S ( my-S / 5) - 10 ;
  ask patches with [ pxcor = round plot-S and pycor < 3 and pycor > 0 ] [ set pcolor red ]


  ;; clear out old blue bar
  ;; ask patches with [  pycor = 3 ]  [ set pcolor white]

  ;; make a new blue horizontal bar showing what the dice MIGHT have rolled ( ie, Q - E to Q + E)
  let lefter  (my-Q - amax)
  let righter (my-Q + amax)
  set lefter  round ( ( lefter  / 5 ) - 10);
  set righter round ( ( righter / 5 ) - 10);
  ask patches with [pxcor >= lefter and pxcor <= righter and pycor = 3 ]  [ set pcolor blue]

  ;; mark possible places the red bar could land that would result in standing up on the first pass
  ask patches with [pxcor > plot-T and pxcor <= righter   and pycor < 3 and pycor > -3 ][ set pcolor green]

  ;; compute the odds of standing with these settings ( Q and T fixed,  S has uniform distribution around Q of half-width E )
  let green-count count patches with [ pcolor = green and pycor = 0 ]  ;; count values of S that result in standing
  let blue-count  count patches with [ pcolor = blue  and pycor = 3 ]  ;; count possible values S could take on
  let stand-pct ( 100 * green-count / blue-count)

  ;; BUG -- truncated pct at 100 max instead of locating the bug in my code!  ////////////////////////////////////////////////////

  if (stand-pct > 100 ) [ set stand-pct 100 ]
  print (word "Given E = " amax  " ,  Q = " my-Q " , T= " my-T ",  the odds of standing in pass 1 are " round stand-pct " percent" )

  ;; and clean up the display out of order to put the red on top of the green
  ask patches with [ pxcor = round plot-S and pycor < 3 and pycor > 0 ] [ set pcolor red ]
end