07 Problem Solving

Stat405 Problem solving

Hadley Wickham
Wednesday, 16 September 2009

Homework
Great improvement. A few tips:
Headings should be in sentence case, and be a
single sentence summary of the section.
Include ggsave calls in your code.
Don’t use " and " to get “ and ”.
Use ` ` and ' '.
(Also, I think I’ve been grading quite hard. Think
of a 5 as an A+, a 4 as an A, and a 3 as an A-/B+)


Feedback
Speed seems about right.
Frustration is a natural part of the learning
experience. You will get better! In
particularly, it will get easier to ﬁnd
interesting plots (helps to start homework
early!)
Split homework help session? Also, feel
free to email me to set up time to talk.


1. Recap of problem
2. Task
3. Basic strategy
4. Turning ideas into code


Slots
Casino claims that slot machines have
prize payout of 92%. Is this claim true?
mean(slots$payoff)
t.test(slots$prize, mu = 0.92)
qplot(prize, data = slots, binwidth = 1)
Can we do better?


Bootstrapping
Slot machine is basically a random
number generator - three numbers
between 1 and 7. We can generate a new
number by sampling from the data.
Doing that lots of times, and calculating
the payoff should give us a better idea of
the payoff. But ﬁrst we need some way to
calculate the prize from the windows


DD DD DD 800
windows <- c("DD", "C", "C")
7 7 7 80
# How can we calculate the
BBB BBB BBB 40 # payoff?
B B B 10
C C C 10
Any bar Any bar Any bar 5
C C * 5
C * C 5
C * * 2
* C * 2 DD doubles any winning
combination. Two DD
* * C 2 quadruples.

Your turn

There are three basic cases of prizes.
What are they? Take 3 minutes to
brainstorm with a partner.


Cases

1. All windows have same value
2. A bar (B, BB, or BBB) in every window
3. Cherries and diamonds


Same values


Same values

1. Determine that all windows are the
same. How?


Same values

same. How?
2. Given single window, look up prize
value. How?


Same values

same. How?
2. Given single window, look up prize
value. How?

With a partner, brainstorm
for 2 minutes on how to
solve one of these problems

# Same value
same <- length(unique(windows)) == 1

# OR
same <- windows[1] == windows[2] &&
windows[2] == windows[3]

if (same) {
# Lookup value
}


If
if (cond) expression
Cond should be a logical vector of length 1
Use && and || to combine sub-conditions
Not & and | - these work on vectors
&& and || are “short-circuiting”: they
do the minimum amount of work


if (TRUE) {
# This will be run
}

if (FALSE) {
# This won't be run
} else {
# This will be
}

# Single line form: (not recommended)
if (TRUE) print("True!)
if (FALSE) print("True!)


if (TRUE) {
# This will be run
}

if (FALSE) {
# This won't be run
} else {
# This will be
}
Note indenting.
Very important!
# Single line form: (not recommended)
if (TRUE) print("True!)
if (FALSE) print("True!)


x <- 5
if (x < 5) print("x < 5")
if (x == 5) print("x == 5")

x <- 1:5
if (x < 3) print("What should happen here?")

if (x[1] < x[2]) print("x1 < x2")
if (x[1] < x[2] && x[2] < x[3]) print("Asc")
if (x[1] < x[2] || x[2] < x[3]) print("Asc")


if (window[1] == "DD") {
prize <- 800
} else if (windows[1] == "7") {
prize <- 80
} else if (windows[1] == "BBB") ...

# Or use subsetting
c("DD" = 800, "7" = 80, "BBB" = 40)
c("DD" = 800, "7" = 80, "BBB" = 40)["BBB"]
c("DD" = 800, "7" = 80, "BBB" = 40)["0"]
c("DD" = 800, "7" = 80, "BBB" = 40)[window[1]]


Your turn

Complete the previous code so that if all
the values in win are the same, then prize
variable will be set to the correct amount.


All bars

How can we determine if all of the
windows are B, BB, or BBB?
(windows[1] == "B" ||
windows[1] == "BB" ||
windows[1] === "BBB") && ... ?


All bars

How can we determine if all of the
windows are B, BB, or BBB?
(windows[1] == "B" ||
windows[1] == "BB" ||
windows[1] === "BBB") && ... ?

Take 2 minutes to brainstorm
possible solutions

windows[1] %in% c("B", "BB", "BBB")
windows %in% c("B", "BB", "BBB")

allbars <- windows %in% c("B", "BB", "BBB")
allbars[1] & allbars[2] & allbars[3]
all(allbars)

# See also ?any for the complement


Your turn

Complete the previous code so that the
correct value of prize is set if all the
windows are the same, or they are all
bars


payoffs <- c("DD" = 800, "7" = 80, "BBB" = 40,
"BB" = 25, "B" = 10, "C" = 10, "0" = 0)

allbars <- all(windows %in% c("B", "BB", "BBB"))

if (same) {
prize <- payoffs[windows[1]]
} else if (allbars) {
prize <- 5
}


Cherries

Need numbers of cherries, and numbers
of diamonds (hint: use sum)
Then need to look up values (like for the
ﬁrst case) and multiply together


cherries <- sum(windows == "C")
diamonds <- sum(windows == "DD")

c(0, 2, 5)[cherries + 1] *
c(1, 2, 4)[diamonds + 1]


payoffs <- c("DD" = 800, "7" = 80, "BBB" = 40,
"BB" = 25, "B" = 10, "C" = 10, "0" = 0)


if (same) {
prize <- 5
} else {

prize <- c(0, 2, 5)[cherries + 1] *
c(1, 2, 4)[diamonds + 1]
}


Writing a function

Now we need to wrap up this code in to a
reusable fashion. We need a function
Have used functions a lot, and this will be
our ﬁrst time writing one


Basic structure
• Name
• Input arguments
• Names/positions
• Defaults
• Body
• Output value


calculate_prize <- function(windows) {
payoffs <- c("DD" = 800, "7" = 80, "BBB" = 40,
"BB" = 25, "B" = 10, "C" = 10, "0" = 0)


if (same) {
prize <- 5
} else {

prize <- c(0, 2, 5)[cherries + 1] *
c(1, 2, 4)[diamonds + 1]
}
prize
}


Homework

What did we miss? Find out the problem
with our scoring function and ﬁx it.
Also, some readings related to group
work.


07 Problem Solving

More Related Content

Similar to 07 Problem Solving

More from Hadley Wickham

Recently uploaded

07 Problem Solving