Is this along the lines of what you are trying to do?
sim.data-data.frame(matrix(rnorm(350*10),350,10))
day-seq(1:350)
sim.data-data.frame(day,sim.data)
pc1.load.all-NULL
pc2.load.all-NULL
for (i in seq(0,300,by=50)){
sim.data.i-subset(sim.data,sim.data$dayisim.data$day(i+50))
.
Ignacio
-Original Message-
From: Roger Bivand [mailto:[EMAIL PROTECTED]
Sent: Tuesday, April 26, 2005 3:00 AM
To: Ignacio Colonna
Cc: R-help@stat.math.ethz.ch
Subject: Re: [R] multiple autocorrelation coefficients in spdep?
On Mon, 25 Apr 2005, Ignacio Colonna wrote:
Hello
The code below gives almost identical results for a split-block analysis in
lme and SAS proc mixed, in terms of variance components and F statistics. It
just extends the example in Pinheiro Bates (p.162) to a split block
design.
I am including below the SAS code and the data in case you want to
Hello,
Has anyone modified the errorsarlm in the R package spdep to
allow for more than a single spatial autocorrelation coefficient (i.e.
'lambda')?
Or, if not, any initial suggestions on how to make that
modification? I have looked at the source code for the function
Maybe aggregate() is what you are looking for?
e.g. say your data frame is called 'mydata'
sum.by.CAT-aggregate(mydata,list(CAT),sum)
this will give you sums by CAT for all the variables in the data set and
will yield 'NA' for any character variables you may have.
Ignacio
-Original
José,
Notice that the values to the left of the | in your stem plot are
all even. Odd numbers are included in the same line.
Try
stem(time,scale=2)
The decimal point is 1 digit(s) to the right of the |
3 | 2
4 | 17
5 | 09
6 | 4667789
7 | 38
8 | 3
9 | 0335
10 |
Emma,
I am not an expert, but I have been trying to fit similar models.
Adding to Keith's reply to your question, I can suggest what I concluded was
the most reasonable model for my case. Based on Keith's Model1, you might
also want to allow for a correlation among years within each
Owen,
I think this gives the plot you are looking for. There may be other
better ways to do it, this is just the one I know. Inside 'panel' you would
need to use 'ltext()' instead of 'text()', as in the example you provided.
xyplot(V1~V2, data=a, groups=V3,
panel = function(x, y, groups)