[R] Problem in lme longitudinal model

From: Bernardo Rangel Tura <tura_at_centroin.com.br>
Date: Wed 18 May 2005 - 01:11:06 EST


Hi R-masters!

I trying model Heart disease mortality in my country with a lme model like this:

m1.lme<-lme(log(rdeath)~age*year,random=~age|year,data=dados)

where: rdeath is rate of mortality per 100000 person per age and year

            age: age of death (22 27 32 37 42 47 52 57 62 67 72 77 82)
            year: year of death (1980:2002)

I donīt have problem to fit the model, but in residual analysis I have one problem.

If i type acf(m1.lme$residuals) the graph show 4 plot with a sin curve (a example in attach) and I get this variogram:

> Variogram(m1.lme)

       variog dist n.pairs
1  0.3939065    1     276
2  0.6486452    2     253
3  0.9679870    3     230
4  1.2765094    4     207
5  1.4158147    5     184
6  1.4264685    6     161
7  1.4048608    7     138
8  1.1902613    8     115
9  0.9619330    9      92
10 0.7037427   10      69
11 0.8591257   11      46
12 1.2215657   12      23

Well
I need help for solution this problem.... Somebody can help me?

Thanks in advance

Bernardo Rangel Tura, MD, MSc
National Institute of Cardiology Laranjeiras Rio de Janeiro Brazil

No virus found in this outgoing message. Checked by AVG Anti-Virus.



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed May 18 01:17:04 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:48 EST