Question

In: Statistics and Probability

The attached code in Chili's Example HW.R models the Chili's restaurant example as seen in the...

The attached code in Chili's Example HW.R models the Chili's restaurant example as seen in the course materials. A version of the example is repeated here for convenience: Let us assume we have two competing models for what effect a proposed Chili’s exterior remodel will have on its stores. One manager is very excited about the effect of exterior renovations and argues that 70% of all Chili’s stores will see at least a 5% increase in sales. Let’s define some notation; the probability of seeing an increase of 5% in sales at any one store is \thetaθ and for this optimistic manager \theta=70\%θ=70% . Another manager argues that the old façade was not really dated and felt that only 20% of remodeled stores will see at least a 5% increase in sales, i.e. \theta = 20\%θ=20% . As a compromise, the managers agree to roll the renovation out to 20 randomly selected stores and then decide whether to remodel all the stores or wait to remodel the other stores. They elect you to make the decision as to who is right. RUN THE CODE LINE-BY-LINE AND DIGEST WHAT HAPPENS WITH EACH LINE. CONFIRM THAT THE FINAL PLOT ECHOES THE MATH DONE FOR THE EXAMPLE AS PRESENTED IN THE COURSE MATERIALS. THEN, MODIFY THE CODE TO ACCOMMODATE THE FOLLOWING CHANGES TO THE BUSINESS ASSUMPTIONS: Assume that you are interested in comparing more than just the two candidate \theta'sθ′s. Modify the code to compare 8 potential \thetaθ values: \theta\in\left\{0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8\right\}θ∈{0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8}. Assume a prior that gives UNEQUAL weight to all managers' models. Assume the pessimistic manager, \theta=10\%θ=10%, wants the prior on his belief to be 65% (i.e.P\left(\theta=10\%\right)=65\%P(θ=10%)=65%) and all the other priors to be 5%. Also, assume that you collect data from FIVE stores. The first, second, and third stores' sales did increase by over 5% and the fourth and fifth stores' sales did not. For coding purposes, assume that a coin flip of heads represents success (i.e. store sales increased by 5% or more) and tails represents a failure (i.e. store sales did not increase). From extracting the proper numbers from the vector of posteriors, find the posterior probability P\left(\theta=10\%\mid D=3\:Successes,\:2\:Failures\right)P(θ=10%∣D=3Successes,2Failures) Enter as a decimal with three significant digits.

#install.packages(c("gridExtra"))
library("ggplot2")
library("gridExtra")
# theta is a parameter representing the probability that a remodeled store sees at least a 5% increase in sales
#Two bosses: Boss1 Theta = 20% and Boss2 Theta = 70%
nThetaVals = 2 ## we have two competing models (i.e. the Boss1 and Boss2) of the world
# Now make the vector of theta values that represent the models of the world that we have:
thetaVals = c(0.2,0.7)

#CALCULATE PRIOR
# pTheta is the vector of prior probabilities on the theta values.
pTheta = c(0.5,0.5) # Makes a uniform belief distribution.

#SPECIFY OBSERVED DATA
# Specify the data. The follwoing are 3 heads and 9 tails and 1 head and 11 tails.
#Data = c(1,1,1,0,0,0,0,0,0,0,0,0)
#Data = c(1,0,0,0,0,0,0,0,0,0,0,0)
Data = c(1) ## one store is a success
nSuccess = sum( Data == 1 ) ##count # of successful stores
nFail = sum( Data == 0 ) ##count # of failure stores

#CALCULATE LIKELIHOOD
# Compute the likelihood of the data for each value of theta:
pDataGivenTheta = thetaVals^nSuccess * (1-thetaVals)^nFail

#CALCULATE POSTERIOR PROBABILITIES
# Compute the posterior:
pData = sum( pDataGivenTheta * pTheta )
pThetaGivenData = pDataGivenTheta * pTheta / pData # This is Bayes' rule!

#make data frame for plotting
plotDF = data.frame(thetaVals, prior = pTheta,likelihood = pDataGivenTheta, posterior = pThetaGivenData)

###code below for plotting the dataframe and showing probabilities
baseLayer = ggplot(data = plotDF, aes(x = thetaVals)) + coord_cartesian(xlim = c(0,1))
priorLayer = geom_bar(aes(y = prior),stat = "identity",width = 0.01, fill = "blue")
likelihoodLayer = geom_bar(aes(y = likelihood),stat = "identity",width = 0.01, fill = "red4")
posteriorLayer = geom_bar(aes(y = posterior),stat = "identity",width = 0.01, fill = "purple4")

p1 = baseLayer + priorLayer + geom_text(aes(y = prior + 0.1, label = signif(prior,3)))
p2 = baseLayer + likelihoodLayer + geom_text(aes(y = 1.2*likelihood, label = signif(likelihood,2)))
p3 = baseLayer + posteriorLayer + geom_text(aes(y = posterior + 0.1, label = signif(posterior,3)))

grid.arrange(p1,p2,p3) ##see line 1 to install gridExtra package for this functionality

Solutions

Expert Solution

Please see the results of the R code below

library("ggplot2")
library("gridExtra")
# theta is a parameter representing the probability that a remodeled store sees at least a 5% increase in sales
#Two bosses: Boss1 Theta = 20% and Boss2 Theta = 70%
nThetaVals = 2 ## we have two competing models (i.e. the Boss1 and Boss2) of the world
# Now make the vector of theta values that represent the models of the world that we have:
thetaVals = c(0.2,0.7)

#CALCULATE PRIOR
# pTheta is the vector of prior probabilities on the theta values.
pTheta = c(0.5,0.5) # Makes a uniform belief distribution.

#SPECIFY OBSERVED DATA
# Specify the data. The follwoing are 3 heads and 9 tails and 1 head and 11 tails.
#Data = c(1,1,1,0,0,0,0,0,0,0,0,0)
#Data = c(1,0,0,0,0,0,0,0,0,0,0,0)
Data = c(1) ## one store is a success
nSuccess = sum( Data == 1 ) ##count # of successful stores
nFail = sum( Data == 0 ) ##count # of failure stores

#CALCULATE LIKELIHOOD
# Compute the likelihood of the data for each value of theta:
pDataGivenTheta = thetaVals^nSuccess * (1-thetaVals)^nFail

#CALCULATE POSTERIOR PROBABILITIES
# Compute the posterior:
pData = sum( pDataGivenTheta * pTheta )
pThetaGivenData = pDataGivenTheta * pTheta / pData # This is Bayes' rule!

#make data frame for plotting
plotDF = data.frame(thetaVals, prior = pTheta,likelihood = pDataGivenTheta, posterior = pThetaGivenData)

###code below for plotting the dataframe and showing probabilities
baseLayer = ggplot(data = plotDF, aes(x = thetaVals)) + coord_cartesian(xlim = c(0,1))
priorLayer = geom_bar(aes(y = prior),stat = "identity",width = 0.01, fill = "blue")
likelihoodLayer = geom_bar(aes(y = likelihood),stat = "identity",width = 0.01, fill = "red4")
posteriorLayer = geom_bar(aes(y = posterior),stat = "identity",width = 0.01, fill = "purple4")

p1 = baseLayer + priorLayer + geom_text(aes(y = prior + 0.1, label = signif(prior,3)))
p2 = baseLayer + likelihoodLayer + geom_text(aes(y = 1.2*likelihood, label = signif(likelihood,2)))
p3 = baseLayer + posteriorLayer + geom_text(aes(y = posterior + 0.1, label = signif(posterior,3)))

grid.arrange(p1,p2,p3) ##see line 1 to install gridExtra package for this functionality

The results are


Related Solutions

The attached code in Chili's Example HW.R models the Chili's restaurant example as seen in the...
The attached code in Chili's Example HW.R models the Chili's restaurant example as seen in the course materials. A version of the example is repeated here for convenience: Let us assume we have two competing models for what effect a proposed Chili’s exterior remodel will have on its stores. One manager is very excited about the effect of exterior renovations and argues that 70% of all Chili’s stores will see at least a 5% increase in sales. Let’s define some...
Sample Code and Models Run each of the models below and explain the code function and...
Sample Code and Models Run each of the models below and explain the code function and your findings for each system, do they agree/disagree with what you understand and why ?? Matlab Code % Winter 2018 Control Engineering % Lab No.3 - Root Locus problems % Mark Clarke clear s = tf('s') K = 1150; %Proportional Controller Gain, May need to be altered? % Enter Model 1 % This is a model of a simple 2nd order with no zeros...
The code file is attached. Portion of the code is given (1, 2, 3, 4). You...
The code file is attached. Portion of the code is given (1, 2, 3, 4). You have to write the code for the remaining items (5, 6, 7). #include <iostream> using namespace std; struct node { int data; node* next; }; node* head=NULL; void appendNode(int value) { node *newNode, *curr; newNode = new node(); newNode->data=value; newNode->next=NULL; if (!head) head=newNode; else { curr=head; while (curr->next) curr = curr->next; curr->next = newNode; } } void insertNode(int value) { node *newNode, *curr, *previous;...
code in. c++ void seen ); If there is already a Word object in the Words...
code in. c++ void seen ); If there is already a Word object in the Words list, then the number of occurrences for this word is incremented. If there is no Word object for this word already, create a new word object with occurrence =1, and insert this object into the list of Word objects. getNextWord(); Returns the next word of the list and sets the currentItem pointer to the next Word. This may return an empty string, “”, if...
A fast-food restaurant determines the cost and revenue models for its hamburgers. A fast-food restaurant determines...
A fast-food restaurant determines the cost and revenue models for its hamburgers. A fast-food restaurant determines the cost and revenue models for its hamburgers. C = 0.8x + 7100,     0 ≤ x ≤ 50,000 R = 1 10,000 (66,000x − x2),     0 ≤ x ≤ 50,000 (a) Write the profit function for this situation. P =   (b) Determine the intervals on which the profit function is increasing and decreasing. (Enter your answers using interval notation.) increasing     decreasing     (c) Determine how many hamburgers...
You have opened a new motel with an attached restaurant and cocktail lounge. What are the...
You have opened a new motel with an attached restaurant and cocktail lounge. What are the risks that your business faces? What insurance coverages should you consider essential in your insurance program?
Apply The knowledge of “Mapping Models to code” Concept to write Chunks of java code which...
Apply The knowledge of “Mapping Models to code” Concept to write Chunks of java code which maps the following relationships: a. Aggregation b. Composition c. Generalization d. Realization e. Dependency
Are decisions in business based on models? Give an example
Are decisions in business based on models? Give an example
Response in Coral Code and code example in Coral Code please!!!! Respond to the following in...
Response in Coral Code and code example in Coral Code please!!!! Respond to the following in a minimum of 175 words: It is important to program your code efficiently. Efficient code does not include duplicate code that performs the same procedure. Functions are helpful in making code modular. How can functions reduce the amount of written code and make code easier to read/follow? Provide a code example that supports your comments.
George Jetson owns and operates a restaurant and catering business which has seen the business show...
George Jetson owns and operates a restaurant and catering business which has seen the business show an increase in revenue over the last 6 months. With an increase in revenues George realized the expenses would also increase, but expected to see an increase in the cash position of the business. In reviewing the cash position and bank account, he has realized there has been a decrease in the cash position. The business receives cash and credit cards at the restaurant...
ADVERTISEMENT
ADVERTISEMENT
ADVERTISEMENT