Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

GeneralizedLinearModelFit and offsert for poisson regression

  • To: mathgroup at smc.vnet.net
  • Subject: [mg112975] GeneralizedLinearModelFit and offsert for poisson regression
  • From: Ester Vilaprinyó <evilaprinyo at cmb.udl.cat>
  • Date: Sat, 9 Oct 2010 06:33:13 -0400 (EDT)

Hello, I would like to solve a poisson regression with the 
GeneralizedLinearModelFit function.
I have a matrix with the following form:

{{x11,x12, x13, y1},
{x21,x22, x23, y2},
...
}

The GeneralizedLinearModelFit for an exponential family Poisson will return
e^(b0+b1 f1+ b2 f2+ b3 f3)
Can I fix that b3=1?

In my case y refers to death counts and x3 to population. In other 
programs such as Stata this would be solved by offset(x3). I tried to 
solve this using LinearOffsetFunction but I did not succeed with the 
following steps: 1) created a new matrix for data = {{x11,x12, 
y1},{x21,x22, y2},...}. 2) Create a function that will return the offset 
for each input f[x1_,x2_]:=x3 3) 
GeneralizedLinearModelFit[data,{x1,x2},{x1,x2},ExponentialFaily->"Poisson", 
LinearOffsetFunction->(f[#1,#2]&)]. Why it does not work?

Thanks!




  • Prev by Date: Further expanding d[xy]/dx in the output of a Dt calculation
  • Next by Date: Re: A list based table join?
  • Previous by thread: Re: Further expanding d[xy]/dx in the output of a Dt calculation
  • Next by thread: Select, from table data