Assignment 1 Problems for Nonlinear Statistical Models | ST 762, Assignments of Statistics

Material Type: Assignment; Class: Nonlinear Statistical Models for Univariate and Multivariate Response; Subject: Statistics; University: North Carolina State University; Term: Fall 2007;

Typology: Assignments

Pre 2010

Uploaded on 03/10/2009

koofers-user-9dj
koofers-user-9dj 🇺🇸

10 documents

1 / 4

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
ST 762, HOMEWORK 1 EXTRA PROBLEMS, FALL 2007
These problems are from previous years and are for you to work on or not as you choose; they are
not to be turned in. You should be familiar with the concepts covered by these problems for the
midterm test. Solutions will be posted when the homework problems to be turned in are due.
1. Suppose we have independent pairs (Yj,xj), j= 1,...,n.
(a) Suppose that E(Yj|xj) = f(xj,β), j = 1,...,n, and that the conditional distribution of
Yjgiven xjis Poisson. Write down the corresponding loglikelihood (conditional on xj) for β
and derive the resulting estimating equation to be solved to obtain the maximum likelihood
estimator for β.
(b) Now suppose that E(Yj|xj) = kjf(xj,β), j = 1,...,n, and that Yjis the number of
successes in kjindependent trials, so that Yjgiven xjfollows a binomial distribution.
(c) What features do the estimating equations in (a) and (b) under these different distribu-
tional assumptions have in common?
2. Suppose we have independent pairs (Yj,xj), j= 1,...,n, such that
E(Yj|xj) = f(xj,β),var(Yj|xj) = σ2f2(xj,β), j = 1,...,n,
with σknown. For some of the calculations below, you may find it convenient to define
λβ(xj,β) = ∂/∂β log f(xj,β).
(a) Assuming that the distribution of Yjgiven xjis gamma, write down the corresponding
loglikelihood for β, and derive the form of the estimating equation for the maximum likelihood
estimator for β.
(b) Same as in (a), but now assuming the distribution of Yj|xjis normal.
(c) Same as in (a), but now assuming the distribution of Yj|xjis lognormal. (If ZN(m, γ2),
then Y= exp(Z) has a lognormal distribution.)
(d) Compare each estimating equation in (a) (c) to the equations you found in Problem 1(a)
and (b), commenting specifically on similarities or differences in the general forms of these
equations. Do any of these equations share the common features you noted in 1(c)?
(e) Compare the estimating equations you found in (a) and (b) of this problem. State explic-
itly what features they share and what features are different. Give a possible reason for the
differences you observe.
3. More on transformations for nonlinear models. An assumption underlying the “Transform
Both Sides” (TBS) model discussed on pages 36–38 of the class notes is that a single trans-
formation can be used to achieve both constant variance and normality on the transformed
scale. For independent pairs (Yj,xj), j = 1,...,n, a generalization of this model has
E{h(Yj, λ)|xj}=h{f(xj,β), λ},var{h(Yj, λ)|xj}=σ2q2(xj, θ) (1)
for some transformation hdepending on a scalar parameter λ. Note that this model implies
that, after transformation, variance is not constant but depends on a function qof xjthat is
known if θis known.
1
pf3
pf4

Partial preview of the text

Download Assignment 1 Problems for Nonlinear Statistical Models | ST 762 and more Assignments Statistics in PDF only on Docsity!

ST 762, HOMEWORK 1 EXTRA PROBLEMS, FALL 2007

These problems are from previous years and are for you to work on or not as you choose; they are not to be turned in. You should be familiar with the concepts covered by these problems for the midterm test. Solutions will be posted when the homework problems to be turned in are due.

  1. Suppose we have independent pairs (Yj , xj ), j = 1,... , n. (a) Suppose that E(Yj |xj ) = f (xj , β), j = 1,... , n, and that the conditional distribution of Yj given xj is Poisson. Write down the corresponding loglikelihood (conditional on xj ) for β and derive the resulting estimating equation to be solved to obtain the maximum likelihood estimator for β. (b) Now suppose that E(Yj |xj ) = kj f (xj , β), j = 1,... , n, and that Yj is the number of successes in kj independent trials, so that Yj given xj follows a binomial distribution. (c) What features do the estimating equations in (a) and (b) under these different distribu- tional assumptions have in common?
  2. Suppose we have independent pairs (Yj , xj ), j = 1,... , n, such that

E(Yj |xj ) = f (xj , β), var(Yj |xj ) = σ^2 f 2 (xj , β), j = 1,... , n,

with σ known. For some of the calculations below, you may find it convenient to define

λβ (xj , β) = ∂/∂β log f (xj , β).

(a) Assuming that the distribution of Yj given xj is gamma, write down the corresponding loglikelihood for β, and derive the form of the estimating equation for the maximum likelihood estimator for β. (b) Same as in (a), but now assuming the distribution of Yj |xj is normal. (c) Same as in (a), but now assuming the distribution of Yj |xj is lognormal. (If Z ∼ N (m, γ^2 ), then Y = exp(Z) has a lognormal distribution.) (d) Compare each estimating equation in (a) – (c) to the equations you found in Problem 1(a) and (b), commenting specifically on similarities or differences in the general forms of these equations. Do any of these equations share the common features you noted in 1(c)? (e) Compare the estimating equations you found in (a) and (b) of this problem. State explic- itly what features they share and what features are different. Give a possible reason for the differences you observe.

  1. More on transformations for nonlinear models. An assumption underlying the “Transform Both Sides” (TBS) model discussed on pages 36–38 of the class notes is that a single trans- formation can be used to achieve both constant variance and normality on the transformed scale. For independent pairs (Yj , xj ), j = 1,... , n, a generalization of this model has

E{h(Yj , λ)|xj } = h{f (xj , β), λ}, var{h(Yj , λ)|xj } = σ^2 q^2 (xj , θ) (1)

for some transformation h depending on a scalar parameter λ. Note that this model implies that, after transformation, variance is not constant but depends on a function q of xj that is known if θ is known.

(a) By an argument similar to that on page 37 of the notes, show that (1) is roughly equivalent to a certain mean-variance model for E(Yj |xj ) and var(Yj |xj ). Give the form of E(Yj |xj ) and var(Yj |xj ). Hint: Define ej = [h(Yj , λ) − h{f (xj , β), λ}]/{q(xj , θ)}. (b) Suppose we have data for which we believe var(Yj |xj ) follows the model in (2.14) on page 40 of the class notes. Would it be possible to use model (1) with a suitable choice of q and the Box-Cox transformation on page 35 to arrive at a model on the transformed scale for which the variance is constant on the transformed scale? Explain. (c) Consider the usual TBS model on page 36 of the notes with the Box-Cox transformation, for which g(xj , θ) ≡ 1. Suppose the response is such that Y takes on only positive values. Is it always possible to transform such a response to normality using (2.11) for any value of λ? Give an argument justifying your answer.

  1. Still more on transformations for nonlinear models. The Michaelis-Menten (MM) model

f (x, β) = V x K + x

= {β 0 + β 1 /x}−^1 , β 1 = 1/V, β 2 = K/V, β = (β 1 , β 2 )T

is widely used to model data in biological, biochemical, and other situations. In fisheries research, it is called the Beverton-Holt spawner-recruit model. Often, data that are well represented by this model also exhibit nonconstant variance. A number of methods have been proposed to estimate β when it is assumed that E(Y |x) = f (x, β). Here, we explore how these may all be viewed as special cases of (1) in the previous problem. Specifically, consider the general TBS model for independent pairs (Yj , xj ), j = 1 ,... , n given in (1) for some transformation h depending on a scalar parameter λ with q(xj , θ) = xθj. Throughout this problem, take h to be the Box-Cox transformation on page 35 of the notes. Note that x is a scalar here. (a) Rather than fitting the (nonlinear) MM model to data on the original scale, perhaps modeling nonconstant variance explicitly, it has been traditional to try to “linearize” the model by using some sort of transformation. Two common such approaches are to consider the following models, which are usually written in the “classical” response = mean + error form:

(i) Lineweaver-Burk: 1/Y = β 0 + β 1 /x + e, where e is taken to be mean zero with variance σ^2. (ii) Woolf: 1/Y = β 0 + β 1 /x + e/x, where e is taken to be mean zero with variance σ^2.

Note that in both cases, the model for conditional mean is linear in β 0 and β 1 , which is the appeal of these approaches. Show that each of (i), (ii) is a special case of the general TBS model described above by finding the values of λ and θ to which they correspond. Using your result to Problem 3, give the approximate mean-variance relationship on the original scale that is being assumed for “small” error. (b) Another popular approach has been to fit the so-called Scatchard model, which is usually written as Y /x = (1/β 1 ) − (β 0 /β 1 )Y + e for “errors” e with mean 0 and constant variance σ^2. Verify that the Scatchard model corresponds to assuming that Y |x has constant coefficient of variation, and find this coefficient

the “test” to the 5th-from-last observation, and so on until an observation is less than 10% likely. Now hold these fixed and obtain start values for β 1 , β 2 by considering the form of (2).

(b) Write a program that implements your method. To test it, apply it to the data in the file biexp.dat available on the class web page. Using your IRWLS program, assuming var(Yj |xj ) = σ^2 f 2 θ(xj , β) with θ = 1.0, use your starting values to fit (2) to these data.