[R] Age-period-cohort model

[SAULEAU Erik-André]

Dear R-list,

I think it's really a newbie question but ...

I try to model age-period-cohort models with polynoms in each effect. I have
for each level of age and period some cases -k- and persons-years -py-.
Models are A_xP_y:
log(k_age,period/py_age,period)=f(polynom(x),age)+f(polynom(y),period). For
exemple I try to adjust an A2P3 :
log(k/py)=a1.age+a2.age*age+p1.period+p2.period*period+p3.period*period*peri
od.  I define factor(effect). My command is:

glm(k~-1+a+a2+p+p2+p3+offset(log(py)), family=poisson).

But I obtain only estimation for a and p and "NA" for a2, p2, p3.
What is the problem?

Thanks for advance.

============================================
Erik-André SAULEAU

SEAIM
Hôpital du Hasenrain
BP 1070
87, Avenue de Altkirch
68051 MULHOUSE Cédex
Tel: 03-89-64-79-95
Mél: sauleauea at ch-mulhouse.fr
============================================




************************************************************************************
Afin d'éviter toute propagation de virus informatique, et en complément des dispositifs en place, ce message (et ses pièces jointes s'il y en a) a été automatiquement analysé par un antivirus de messagerie. 
***********************************************************************************


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._