Reference no: EM13989984
Detailed Question:
The homework is for computational statistics course (Optimization, R coding, regression).
This is an advanced course. THE R CODES ARE IMPORTANT AND SHOULD BE INCLUDED IN THE ANSWER.
Instructions:
A scientist studied the effects of a very strong storm in trees in a forest in Minnesota as a function of the tree diameters and the severity of winds in the areas where the trees were located. The data for n = 659 trees are given at tree.txt file
The data set consists of the following variables:
D = Diameter of the tree.
S = Severity of the storm at the location of the tree in the scale of 0 to 1, with 1 being most severe.
y = 0 if the tree survived and 1 if the tree is blown down.
Here, y1, ... , yn is a set of n observations where yi ∼ Bernouli(Πi). We are to estimate Πi = P(yi = 1), the probability that a tree is blow out, as a function of D and S. In particular we parameterize this probability by the logit function
Πi(β) = exp(xiTβ)/1 + exp(xiTβ) for i = 1, ......n (1)
where xiT = (xi1, xi2, xi3), where xi1 = 1 is the intercept, xi2 = log(D), and xi3 = S, and β = (β1, β2, β3)T is a set of p parameters to be estimated.
(a) Write the likelihood function l(β) as a function of yi's and Πi(β).
(b) Obtain the first differential dΠi(dβ) and the second differential ddΠi(dβ, dβ).
(c) Obtain the first differential dl(dβ) and use it to derive formulas for ∂l(β)/∂β)i for B = 1, 2, 3.
(d) Write an R code to implement steepest ascent method with step-halving to obtain the maximum likelihood estimates for β.
(e) Obtain the second differential ddl(dβ, dβ).
(f) Obtain the expected value dd[dl(dβ, dβ)] and use it derive formulas for the elements of the Fisher information matrix.
(g) Write an R code to implement the Fisher-scoring algorithm to obtain the maximum likelihood estimates for β. [Do not use iteratively reweighted least squares to run Fisher scoring.]
(h) Use the glm function in R to obtain maximum likelihood estimates for β.
(i) Plot the three-dimensional graph of Π(β), using your maximum likelihood estimates and as a function of X1 = S and X2 = log(D), For your plot, X1 should range in the interval (0,1) and X2 should range in the range of observed values for log(D). Moreover, draw the constant value contours and give an interpretation of the values.
(j) What is the prediction of your model for the probability that a tree with diameter 10 would blow down as a function of the severity of wind, as it ranges from 0 to 1. Show the probability using a graph.
(h) What is the prediction of your model for the probability that a tree would be blown down with wind severity 0.4 and diameters ranging from 5 to 30.
Attachment:- tree.txt
|
Knowledge of conflict-negotiations and power in organization
: How can knowledge of conflict, negotiations, and power in organizations help you be more successful in your career? Please be specific in your response for each of the three topics.
|
|
Explain the action of indicators
: How Ostwald's theory and quininoid theory explain the action of indicators
|
|
Discuss the interests of facebooks major stakeholder groups
: Define and discuss the interests of Facebook’s major stakeholder groups. What measures does the company take to manage its different stakeholders? How does it prioritize or balance their different needs?
|
|
Expected rates of return for different products if the paint
: Gregory's company has ordered a large supply of paint, but the manufacturer says there is a possibility that the batch they received was defective. However, the company has no choice but to start filling orders, and they are trying to decide where th..
|
|
What is the prediction of your model for the probability
: What is the prediction of your model for the probability that a tree with diameter 10 would blow down as a function of the severity of wind, as it ranges from 0 to 1 - Write the likelihood function l(β) as a function of yi's and Πi(β).
|
|
Does an elective abortion lead to nagetive psychology effect
: Article- Does an Elective Abortion Lead to Nagetive Psychology Effects? Read the article, please read it and write 3 pages
|
|
Imagine that youre working as a video game developer at the
: Imagine that you're working as a video game developer at the Chicago office of Electronic Arts (EA), and you have been offered a transfer to the EA development studio in Tokyo, Japan.
|
|
Definition of linear transformation and write it
: explain mathematically what the plots show, namely, the given transformation y=f(x) is a linear one. For 25 points, show the plots without the convincing argument.find a definition of linear transformation and write it (with citation).
|
|
What is the charge q on the capacitor time
: An uncharged capacitor with C = 14.9 μF, a resistor with R = 24.3 k Omega, and a V 25.7V are connected in series, What is the charge q (in C)on the capacitor Time t = 2 time constants, after the switch is dosed?
|