[R] Hi,

Chong Gu chong at stat.purdue.edu
Wed Jan 30 05:15:13 CET 2002


   X-Mailer: Novell GroupWise Internet Agent 6.0.1
   Date: Tue, 29 Jan 2002 21:18:40 -0500
   From: "Michael Roberts" <mroberts at ers.usda.gov>
   Mime-Version: 1.0
   Content-Type: text/plain; charset=US-ASCII
   Content-Transfer-Encoding: 7bit
   Content-Disposition: inline
   Sender: owner-r-help at stat.math.ethz.ch
   Precedence: SfS-bulk

   Hi,

   Sorry for the confusion.

   I would like to estimate a model wherein
   the marginals of z with respect to w1 and w2 
   are smooth functions of  x and y.  I have data
   on z, x, y, w1 and w2. 

   so E[dz/dw1] = f(x,y) and E[dz/dw2] = g(x,y)

   and I would like to estimate f(x,y) and g(x,y)

   I suppose I could try to fit something more general
   using projection pursuit, but the nature of the problem
   suggests the above structure.

   For some reason I thought 

   x:y:z 

   would fit just the interaction term 

   xyz 

   and not expend to 

   x + y + z + xy +xz + yz + xyz

   like x*y*z,  which is why I wrote it the way I did.

   So maybe it should have bern written

   y ~ I(f(x,y)*w1) + I(g(x,y)*w2) + e

   e is a symmetric random error.

   This seems identifiable to me, but am I missing something?



   Michael J. Roberts

   Resource Economics Division
   Production, Management, and Technology
   USDA-ERS
   (202) 694-5557 (phone)
   (202) 694-5775 (fax)


gss would in principal fit y~x*y*w1+x*z*w2 for you, but the current
algorithm is too slow for the sample size needed for fitting such a
model successfully.

As for identifiability, the current form MIGHT be okey, but are you
sure y=e for w1=w2=0?

Chong Gu
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list