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