search for: obs6

...ects in lme/lmer ? should I try to define a "new" correlation structure corMultiAR1 ? If so, where can I find help on how to write such a piece of code ( nlme:::corAR1 is not clear to me) ? Or is there a way to define a general parametrised covariance matrix in gls ? Olivier > obs6 <- matrix( c(1,2,3,4,5,6, 1,1,2,2,3,3), byrow=F, nc=2) > dimnames(obs6) <- list(NULL, c("y","time")) > obs6 <- data.frame(obs6) > obs6 y time 1 1 1 2 2 1 3 3 2 4 4 2 5 5 3 6 6 3 > gls (y~1, correl=corAR1(0.0,~time), data=obs6) Error in Init...

...1 15/02/2000 6 n 2 1 23/02/2000 I want to select the observations whose difference of "var4" is less than 7 for the cases with the same values of var1,var2 andvar3. The obervations have the same var1, var2 and var3 are, part1 (obs2 and obs3) and part2 (obs4,obs5, and obs6). For obs2 and obs3, their date difference is less than 7, so we donot need to delete any of them. For obs4,obs5, and obs6,we can see that obs6 should be deleted becuase its date is over 7 dyas longer than obs4. So the final dataset should obs1,obs2,obs3,obs4, and obs5. I have a lot of obse...

...DA DID6 DID7 DI 0580: 44 38 20 45 42 41 53 20 43 50 53 50 20 45 45 43 D8 EBAS CPSP EEC 0590: 50 20 45 56 43 50 20 58 42 41 53 20 4f 42 53 31 P EVCP XBAS OBS1 05a0: 20 4f 42 53 32 20 4f 42 53 33 20 4f 42 53 34 20 OBS2 OBS3 OBS4 05b0: 4f 42 53 35 20 4f 42 53 36 20 4f 42 53 37 20 4f OBS5 OBS6 OBS7 O 05c0: 42 53 38 20 00 48 05 41 54 4d 43 08 50 54 4d 43 BS8 .H.ATMC.PTMC 05d0: 08 41 54 52 41 08 50 54 52 41 08 50 4e 48 4d 20 .ATRA.PTRA.PNHM 05e0: 54 42 41 42 20 54 42 41 48 20 52 54 49 50 08 54 TBAB TBAH RTIP.T 05f0: 53 4f 44 08 41 54 50 43 08 50 54 50 43 08 50 46 SOD.ATPC.PTPC...