Hi all, I asked this before the holiday, didn't get any response. So would like to resend the message, hope to get any fresh attention. Since this is not purely lme technical question, so I also cc-ed R general mailing list, hope to get some suggestions from there as well.
I asked some questions on how to analyze longitudinal study with only 2 time points (baseline and a follow-up) previously. I appreciate many useful comments from some members, especially Dennis Murphy and Marc Schwartz who refered the following paper addressing specifically this type of study with only 2 time points: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1121605/ Basically, with only 2 time points (baseline and one follow-up), ANCOVA with follow-up as dependent variable and baseline as covariate should be used: follow-up = a + b*baseline + treatment Now I have a regular longitudinal study with 6 time points, 7 treatments (vehicle, A, B, C, D, F, G), measuring a response variable "y". The dataset is attached. I have some questions, and appreciate any suggestions on how to analyze the dataset. dat<-read.table("dat.txt",sep='\t',header=T,row.names=NULL) library(MASS) dat$trt<-relevel(dat$trt,'vehicle') xyplot(y~time, groups=trt, data=dat, ylim=c(3,10),col=c(1:6,8),lwd=2,type=c('g','a'),xlab='Days',ylab="response", key = list(lines=list(col=c(1:6,8),lty=1,lwd=2), text = list(lab = levels(dat$trt)), columns = 3, title = "Treatment")) So as you can see that there is some curvature between glucose level and time, so a quadratic fit might be needed. dat$time2<-dat$time*dat$time A straight fit like below seems reasonable: fit<-lmer(y~trt*time+trt*time2+(time|id),dat) Checking on random effects, it appears that variance component for random slope is very small, so a simpler model with random intercept only may be sufficient: fit<-lmer(y~trt*time+trt*time2+(1|id),dat) Now, I want to incorporate baseline response into the model in order to account for any baseline imbalance. I need to generate a new variable "baseline" based on glucose levels at time=0: dat<-merge(dat, dat[dat$time==0,c('id','y')], by.x='id',by.y='id',all.x=T) colnames(dat)[c(4,6)]<-c('y','baseline') so the new fit adding baseline into the mixed model is: fit<-lmer(y~baseline+trt*time+trt*time2+(1|id),dat) Now my question is 1). Is the above model a reasonable thing to do? 2) when baseline is included as a covariate, should I remove the data points at baseline from the dataset? I am kind of unsure if it's reasonable to use the baseline both as a covariate and as part of the dependent variable values. Next thing I want to do with this dataset is to do multiple comparisons between each treatment (A, B, C, D, F, G) vs. vehicle at a given time point, say time=56 (the last time points) after adjusting the baseline imbalance. This seems to be done using Dunnet test. When I say "after adjusting baseline imbalance", I mean the comparisons should be done based on the difference between time=56 and time=0 (baseline), i.e. is there any difference in the change from baseline for treatment A (or B, C, D, F, G) vs. vehicle?. How can we test this? Will glht() in multcomp work for a lmer fit? If yes, how can I specify the syntax? Finally, with the above model, how to estimate the difference (and the standard error) between time=56 and time=0 (baseline) for each treatment groups? Thank you all for your attention. John
id trt time y 1 A 14 8 1 A 28 4 1 A 7 9 1 A 42 7 1 A 0 2 1 A 56 4 2 vehicle 14 3 2 vehicle 28 2 2 vehicle 7 6 2 vehicle 42 2 2 vehicle 0 3 2 vehicle 56 3 3 A 0 3 3 A 56 4 3 A 14 5 3 A 28 4 3 A 7 8 3 A 42 7 4 A 14 4 4 A 28 8 4 A 7 7 4 A 42 5 4 A 0 5 4 A 56 8 5 A 14 7 5 A 28 3 5 A 7 7 5 A 42 6 5 A 0 3 5 A 56 7 6 F 0 6 6 F 56 30 6 F 14 17 6 F 28 14 6 F 7 15 6 F 42 30 7 B 14 5 7 B 28 8 7 B 7 5 7 B 42 1 7 B 0 4 7 B 56 2 8 B 14 7 8 B 28 5 8 B 7 8 8 B 42 1 8 B 0 7 8 B 56 5 9 B 0 7 9 B 56 5 9 B 14 10 9 B 28 7 9 B 7 15 9 B 42 9 10 B 0 2 10 B 56 3 10 B 14 4 10 B 28 5 10 B 7 1 10 B 42 2 11 vehicle 14 11 11 vehicle 28 5 11 vehicle 7 9 11 vehicle 42 7 11 vehicle 0 7 11 vehicle 56 3 12 B 14 3 12 B 28 2 12 B 7 3 12 B 42 1 12 B 0 3 12 B 56 12 13 B 14 7 13 B 28 20 13 B 7 7 13 B 42 9 13 B 0 7 13 B 56 8 14 F 14 4 14 F 28 5 14 F 7 4 14 F 42 5 14 F 0 5 14 F 56 5 15 vehicle 14 26 15 vehicle 28 4 15 vehicle 7 11 15 vehicle 42 21 15 vehicle 0 7 15 vehicle 56 4 16 C 14 3 16 C 28 3 16 C 7 3 16 C 42 5 16 C 0 3 16 C 56 2 17 C 14 6 17 C 28 5 17 C 7 4 17 C 42 3 17 C 0 3 17 C 56 5 18 vehicle 0 3 18 vehicle 56 5 18 vehicle 14 5 18 vehicle 28 3 18 vehicle 7 7 18 vehicle 42 5 19 C 14 1 19 C 28 3 19 C 7 1 19 C 42 3 19 C 0 1 19 C 56 5 20 C 0 5 20 C 56 5 20 C 14 5 20 C 28 6 20 C 7 7 20 C 42 3 21 F 0 2 21 F 56 2 21 F 14 3 21 F 28 2 21 F 7 2 21 F 42 2 22 F 14 17 22 F 28 7 22 F 7 18 22 F 42 6 22 F 0 7 22 F 56 15 23 C 14 3 23 C 28 3 23 C 7 3 23 C 42 5 23 C 0 4 23 C 56 5 24 D 0 4 24 D 56 17 24 D 14 13 24 D 28 2 24 D 7 3 24 D 42 3 25 D 0 6 25 D 56 0 25 D 14 7 25 D 28 3 25 D 7 7 25 D 42 12 26 F 14 5 26 F 28 6 26 F 7 3 26 F 42 3 26 F 0 5 26 F 56 5 27 D 14 1 27 D 28 7 27 D 7 5 27 D 42 4 27 D 0 1 27 D 56 3 28 vehicle 0 2 28 vehicle 56 3 28 vehicle 14 3 28 vehicle 28 2 28 vehicle 7 4 28 vehicle 42 5 29 D 0 6 29 D 56 3 29 D 14 5 29 D 28 7 29 D 7 5 29 D 42 4 30 D 0 0 30 D 28 6 30 D 56 6 31 D 14 8 31 D 28 2 31 D 7 9 31 D 42 5 31 D 0 11 31 D 56 3 32 A 0 3 32 A 14 3 32 A 28 6 32 A 7 3 32 A 42 2 32 A 56 15 33 B 0 5 33 B 14 13 33 B 28 3 33 B 7 4 33 B 42 6 33 B 56 5 34 F 0 7 34 F 14 12 34 F 28 9 34 F 7 10 34 F 42 8 34 F 56 18 35 C 0 3 35 C 14 20 35 C 28 4 35 C 7 7 35 C 42 6 35 C 56 7 36 G 0 0 36 G 14 7 36 G 28 5 36 G 7 4 36 G 42 5 36 G 56 3 37 A 0 7 37 A 14 3 37 A 28 7 37 A 7 3 37 A 42 3 37 A 56 7 38 A 0 3 38 A 14 5 38 A 28 4 38 A 7 3 38 A 42 4 38 A 56 3 39 D 0 7 39 D 14 7 39 D 28 16 39 D 7 7 39 D 42 5 39 D 56 14 40 F 0 3 40 F 14 2 40 F 28 3 40 F 7 5 40 F 42 3 40 F 56 3 41 B 0 5 41 B 14 12 41 B 28 6 41 B 7 7 41 B 42 11 41 B 56 10 42 D 0 3 42 D 14 3 42 D 28 4 42 D 7 4 42 D 42 3 42 D 56 3 43 D 0 7 43 D 56 21 44 B 0 7 44 B 14 7 44 B 28 25 44 B 7 7 44 B 42 6 44 B 56 8 45 C 0 5 45 C 14 6 45 C 28 5 45 C 7 10 45 C 42 8 45 C 56 7 46 G 0 2 46 G 14 5 46 G 28 4 46 G 7 3 46 G 42 5 46 G 56 8 47 A 0 3 47 A 14 8 47 A 28 5 47 A 7 10 47 A 42 8 47 A 56 4 48 A 0 7 48 A 14 5 48 A 28 20 48 A 7 9 48 A 42 7 48 A 56 12 49 A 0 5 49 A 14 13 49 A 28 5 49 A 7 5 49 A 42 13 49 A 56 9 50 B 0 2 50 B 14 9 50 B 28 7 50 B 7 10 50 B 42 5 50 B 56 6 51 G 0 6 51 G 14 12 51 G 28 9 51 G 7 4 51 G 42 20 51 G 56 12 52 A 0 5 52 A 14 7 52 A 28 7 52 A 7 2 52 A 42 2 52 A 56 5 53 G 0 7 53 G 14 8 53 G 28 19 53 G 7 17 53 G 42 14 53 G 56 12 54 vehicle 0 7 54 vehicle 14 11 54 vehicle 28 8 54 vehicle 7 14 54 vehicle 42 16 54 vehicle 56 35 55 B 0 6 55 B 14 5 55 B 28 6 55 B 7 8 55 B 42 5 55 B 56 7 56 C 0 5 56 C 14 9 56 C 28 11 56 C 7 7 56 C 42 5 56 C 56 9 57 A 0 5 57 A 14 10 57 A 28 9 57 A 7 7 57 A 42 5 57 A 56 5 58 B 0 3 58 B 14 5 58 B 28 7 58 B 7 5 58 B 42 5 58 B 56 3 59 A 0 3 59 A 14 2 59 A 28 1 59 A 7 3 59 A 42 3 59 A 56 3 60 C 0 5 60 C 14 26 60 C 28 9 60 C 7 7 60 C 42 7 60 C 56 12 61 C 0 7 61 C 14 15 61 C 28 2 61 C 7 5 61 C 42 3 61 C 56 15 62 D 0 5 62 D 14 5 62 D 28 4 62 D 7 5 62 D 42 5 62 D 56 3 63 B 0 1 63 B 14 6 63 B 28 2 63 B 7 4 63 B 42 8 63 B 56 5 64 F 0 7 64 F 14 7 64 F 28 9 64 F 7 9 64 F 42 26 64 F 56 10 65 B 0 3 65 B 14 1 65 B 28 3 65 B 7 1 65 B 42 3 65 B 56 3 66 D 0 5 66 D 14 20 66 D 28 18 66 D 7 10 66 D 42 25 66 D 56 6 67 vehicle 0 7 67 vehicle 14 10 67 vehicle 28 7 67 vehicle 7 1 67 vehicle 42 7 68 vehicle 0 3 68 vehicle 14 4 68 vehicle 7 5 68 vehicle 56 5 69 A 0 3 69 A 14 5 69 A 28 7 69 A 7 7 69 A 42 6 69 A 56 9 70 vehicle 0 4 70 vehicle 14 19 70 vehicle 28 11 70 vehicle 7 16 70 vehicle 42 6 70 vehicle 56 6 71 A 0 7 71 A 14 15 71 A 28 11 71 A 7 30 71 A 42 15 71 A 56 8 72 vehicle 0 3 72 vehicle 14 3 72 vehicle 28 3 72 vehicle 7 3 72 vehicle 42 3 72 vehicle 56 5 73 A 0 3 73 A 14 12 73 A 28 15 73 A 7 5 73 A 42 3 73 A 56 14 74 F 0 7 74 F 14 17 74 F 28 10 74 F 7 20 74 F 42 15 74 F 56 26 75 F 0 3 75 F 14 9 75 F 28 9 75 F 7 5 75 F 42 10 75 F 56 8 76 D 0 6 76 D 14 10 76 D 7 7 77 D 0 5 77 D 14 6 77 D 28 3 77 D 7 5 77 D 42 6 77 D 56 7 78 C 0 3 78 C 14 6 78 C 28 3 78 C 7 5 78 C 42 2 78 C 56 2 79 C 0 3 79 C 14 9 79 C 28 4 79 C 7 12 79 C 42 5 79 C 56 6 80 B 0 3.5 80 B 14 4 80 B 28 5 80 B 7 4.5 80 B 42 3 80 B 56 9.5 81 F 0 6 81 F 14 6 81 F 28 15.5 81 F 7 24 81 F 42 9.5 81 F 56 8 82 vehicle 0 2.5 82 vehicle 14 6.5 82 vehicle 28 10 82 vehicle 7 9 82 vehicle 42 9 82 vehicle 56 2.5 83 C 0 3.5 83 C 14 6 83 C 28 4 83 C 7 7 83 C 42 5 83 C 56 4 84 G 0 7 84 G 14 5 84 G 28 6 84 G 7 4.5 84 G 42 5 84 G 56 12 85 C 0 3 85 C 14 2 85 C 28 1 85 C 7 2 85 C 42 2.5 85 C 56 1 86 D 0 6 86 D 14 2 86 D 28 2 86 D 7 3 86 D 42 2.5 86 D 56 3.5 87 B 0 3 87 B 14 6 87 B 28 4 87 B 7 3.5 87 B 42 4.5 87 B 56 9 88 D 0 2 88 D 14 5.5 88 D 28 3 88 D 7 2 88 D 42 5 88 D 56 5 89 C 0 3 89 C 14 3 89 C 28 2 89 C 7 6 89 C 42 4 89 C 56 6 90 F 0 6.5 90 F 14 9 90 F 28 4 90 F 7 17 90 F 42 5.5 90 F 56 3 91 G 0 5.5 91 G 14 5.5 91 G 28 2.5 91 G 7 6 91 G 42 3 91 G 56 7.5 92 D 0 2 92 D 14 2 92 D 28 2 92 D 7 2 92 D 42 2.5 92 D 56 2 93 G 0 2 93 G 14 3.5 93 G 28 2 93 G 7 4 93 G 42 2.5 93 G 56 2 94 F 0 6 94 F 14 12 94 F 28 8 94 F 7 8.5 94 F 42 7 94 F 56 7.5 95 vehicle 0 6 95 vehicle 14 2.5 95 vehicle 28 7 95 vehicle 7 3.5 95 vehicle 42 6.5 95 vehicle 56 5.5 96 D 0 3 96 D 14 6 96 D 28 9 96 D 7 6 96 D 42 5 96 D 56 5 97 G 0 2 97 G 7 3.5 98 B 0 4 98 B 14 6 98 B 28 1 98 B 7 2.5 98 B 42 3.5 98 B 56 3 99 G 0 2 99 G 14 3 99 G 28 5 99 G 7 5 99 G 42 2 99 G 56 3.5 100 D 0 5 100 D 14 5 100 D 28 4 100 D 7 5.5 100 D 42 4 100 D 56 6 101 B 0 6 101 B 14 2 101 B 28 10 101 B 7 5 101 B 42 7 101 B 56 10 102 B 0 2 102 B 14 3 102 B 28 2 102 B 7 1 102 B 42 2 102 B 56 2 103 G 0 2 103 G 14 3 103 G 28 3 103 G 7 3 103 G 42 4 103 G 56 4 104 B 0 5 104 B 14 6 104 B 28 4 104 B 7 5 104 B 42 2 104 B 56 2 105 F 0 4 105 F 14 5 105 F 28 5 105 F 7 5 105 F 42 2 105 F 56 3 106 B 0 3 106 B 14 5 106 B 28 8 106 B 7 6 106 B 42 3 106 B 56 3 107 vehicle 0 5 107 vehicle 14 2 107 vehicle 28 13 107 vehicle 7 3 107 vehicle 42 1 107 vehicle 56 5 108 G 0 1 108 G 14 1 108 G 28 2 108 G 7 4 108 G 42 2 108 G 56 2 109 G 0 3 109 G 14 2 109 G 28 2 109 G 7 3 109 G 42 2 109 G 56 1 110 C 0 2 110 C 14 2 110 C 28 1 110 C 7 1 110 C 42 1 110 C 56 1 111 F 0 3 111 F 14 2 111 F 28 2 111 F 7 8 111 F 42 2 111 F 56 8 112 B 0 2 112 B 14 3 112 B 28 7 112 B 7 2 112 B 42 5 112 B 56 4 113 B 0 5 113 B 14 1 113 B 28 2 113 B 7 3 113 B 42 2 113 B 56 1 114 B 0 5 114 B 14 5 114 B 28 8 114 B 7 2 114 B 42 3 114 B 56 2 115 G 0 6 115 G 14 17 115 G 28 22 115 G 7 16 115 G 42 10 115 G 56 29 116 D 0 5 116 D 14 6 116 D 28 4 116 D 7 2 116 D 42 6 116 D 56 6 117 F 0 4 117 F 14 3 117 F 28 1 117 F 7 2 117 F 42 2 117 F 56 2 118 vehicle 0 1 118 vehicle 14 2 118 vehicle 28 7 118 vehicle 7 4 118 vehicle 42 3 118 vehicle 56 2 119 C 0 3 119 C 14 5 119 C 28 3 119 C 7 3 119 C 42 4 119 C 56 2 120 A 0 5 120 A 14 7 120 A 28 7 120 A 7 7 120 A 42 5 120 A 56 6 121 vehicle 0 2 121 vehicle 14 8 121 vehicle 28 3 121 vehicle 7 2 121 vehicle 42 1 121 vehicle 56 5 122 C 0 1 122 C 14 3 122 C 28 2 122 C 7 1 122 C 42 4 122 C 56 3 123 vehicle 0 2 123 vehicle 14 5 123 vehicle 28 4 123 vehicle 7 2 123 vehicle 42 4 123 vehicle 56 4 124 G 0 6 124 G 14 3 124 G 28 4 124 G 7 6 124 G 42 5 124 G 56 5 125 F 0 9 125 F 14 2 125 F 7 3 125 F 56 8 126 A 0 14 126 A 14 4 126 A 7 7 126 A 56 4 127 vehicle 0 3 127 vehicle 14 2 127 vehicle 28 4 127 vehicle 7 2 127 vehicle 42 2 127 vehicle 56 3 128 F 0 7 128 F 56 5 129 C 0 3 129 C 14 10 129 C 28 4 129 C 7 8 129 C 42 2 129 C 56 2 130 D 0 7 130 D 14 10 130 D 28 8 130 D 7 9 130 D 42 11 130 D 56 6 131 C 0 7 131 C 14 15 131 C 28 9 131 C 7 15 131 C 42 14 131 C 56 25 132 F 0 3 132 F 14 10 132 F 28 5 132 F 7 3 132 F 42 3 132 F 56 13 133 F 0 2 133 F 14 6 133 F 28 3 133 F 7 3 133 F 42 5.5 133 F 56 5.5 134 B 0 2 134 B 14 2 134 B 28 2 134 B 7 2 134 B 42 2 134 B 56 2 135 A 0 2 135 A 14 2 135 A 28 1.5 135 A 7 2 135 A 42 2 135 A 56 2 136 D 0 6 136 D 14 4 136 D 28 7 137 B 0 3 137 B 14 4 137 B 28 2 137 B 7 3 137 B 42 4 137 B 56 18 138 D 0 6 138 D 14 9 138 D 28 6 138 D 7 6 138 D 42 10.5 138 D 56 11 139 G 0 2 139 G 14 2 139 G 28 7 139 G 7 3 139 G 42 6 139 G 56 7 140 F 0 4 140 F 14 9 140 F 28 3 140 F 7 7 140 F 42 4.5 140 F 56 2.5 141 C 0 7 141 C 14 9 141 C 28 17.5 141 C 7 6 141 C 42 2 141 C 56 9 142 A 0 1.5 142 A 14 1 142 A 28 2.5 142 A 7 2 142 A 42 2 142 A 56 2 143 B 0 3 143 B 14 6.5 143 B 28 7 143 B 7 7 143 B 42 4 143 B 56 3.5 144 A 0 1.5 144 A 14 2 144 A 28 1 144 A 7 6 144 A 42 3 144 A 56 2 145 C 0 2 145 C 14 2 145 C 28 3 145 C 7 2 145 C 42 1 145 C 56 3 146 C 0 6 146 C 14 6 146 C 28 6 146 C 7 4 146 C 42 6 146 C 56 6 147 B 0 7 147 B 14 6 147 B 28 5 147 B 7 2.5 147 B 42 3 147 B 56 2.5 148 vehicle 0 3.5 148 vehicle 14 8 148 vehicle 28 1 148 vehicle 7 9 148 vehicle 42 2 148 vehicle 56 3 149 F 0 1 149 F 14 3 149 F 28 1 149 F 7 1 149 F 42 3 149 F 56 2 150 D 0 5 150 D 14 12 150 D 28 9 150 D 7 15 150 D 42 9 150 D 56 10 151 G 0 5 151 G 14 7 151 G 28 30 151 G 7 5 151 G 42 10 151 G 56 7 152 C 0 4 152 C 14 9 152 C 28 18 152 C 7 8 152 C 42 18 152 C 56 7 153 vehicle 0 5 153 vehicle 14 10 153 vehicle 28 30 153 vehicle 7 3 153 vehicle 42 10 153 vehicle 56 15 154 G 0 1 154 G 14 4 154 G 28 5 154 G 7 2 154 G 42 7 154 G 56 2 155 G 0 7 155 G 14 6 155 G 7 4 156 B 0 4 156 B 14 3 156 B 28 2 156 B 7 2 156 B 42 7 156 B 56 6 157 B 0 2 157 B 14 2 157 B 28 4 157 B 7 1 157 B 42 1 157 B 56 1 158 vehicle 0 4 158 vehicle 14 9 158 vehicle 28 5 158 vehicle 7 6 158 vehicle 42 5 158 vehicle 56 6 159 A 0 5 159 A 14 8 159 A 28 7 159 A 7 7 159 A 42 7 159 A 56 10 160 B 0 7 160 B 14 5 160 B 28 6 160 B 7 5 160 B 42 13 160 B 56 4 161 F 0 7 161 F 14 20 161 F 28 15 161 F 7 10 161 F 42 16 161 F 56 10 162 D 0 2 162 D 14 6 162 D 28 4 162 D 7 5 162 D 42 1 162 D 56 2 163 D 0 1 163 D 14 5 163 D 28 5 163 D 7 5 163 D 42 5 163 D 56 1 164 F 0 5 164 F 14 20 164 F 28 12 164 F 7 8 164 F 42 13 164 F 56 17 165 C 0 6 165 C 14 16 165 C 28 16 165 C 7 12 165 C 42 21 165 C 56 16 166 F 0 4 166 F 14 8 166 F 28 9 166 F 7 8 166 F 42 4 166 F 56 6 167 G 0 6 167 G 14 25 167 G 28 15 167 G 7 6 167 G 42 32 167 G 56 29 168 A 0 3 168 A 14 6 168 A 28 6 168 A 7 5 168 A 42 5 168 A 56 4 169 vehicle 0 3 170 C 0 3 170 C 14 3 170 C 28 3 170 C 7 3 170 C 42 4 170 C 56 3 171 B 0 3 171 B 14 4 171 B 28 3 171 B 7 4 171 B 42 4 171 B 56 3 172 A 0 4 172 A 14 5 172 A 28 2 172 A 7 3 172 A 42 6 172 A 56 7 173 G 0 6 173 G 14 6 173 G 28 8 173 G 7 4 173 G 42 3 173 G 56 5 174 G 0 2 174 G 14 3 174 G 28 3 174 G 7 1 174 G 42 3 174 G 56 3 175 G 0 4 176 D 0 1 176 D 14 1 176 D 28 1 176 D 7 1 176 D 42 1 176 D 56 5 177 G 0 6 177 G 14 10 177 G 28 1 177 G 7 28 177 G 42 4 177 G 56 3 178 D 0 3 178 D 14 8 178 D 28 7 178 D 7 9 178 D 42 11 178 D 56 8 179 C 0 4 179 C 14 9 179 C 28 5 179 C 7 8 179 C 42 8 179 C 56 15 180 A 0 2 180 A 14 0 180 A 28 1 180 A 7 1 180 A 42 1 180 A 56 1 181 vehicle 0 5 181 vehicle 14 7 181 vehicle 28 8 181 vehicle 7 5 181 vehicle 42 9 181 vehicle 56 8 182 G 0 7 182 G 14 9 182 G 28 7 182 G 7 5 182 G 42 6 182 G 56 8 183 vehicle 0 2 183 vehicle 14 1 183 vehicle 28 4 183 vehicle 7 12 183 vehicle 42 1 183 vehicle 56 11 184 vehicle 0 7 184 vehicle 14 2 184 vehicle 28 1 184 vehicle 7 2 184 vehicle 42 13 184 vehicle 56 2 185 C 0 5 185 C 14 5 185 C 28 8 185 C 7 4 185 C 42 8 185 C 56 4 186 G 0 2 186 G 14 11 186 G 28 7 186 G 7 5 186 G 42 9 186 G 56 5 187 D 0 2 187 D 14 10 187 D 28 0 187 D 7 5 187 D 42 6 187 D 56 10 188 B 0 3 188 B 14 5 188 B 28 2 188 B 7 2 188 B 42 5 188 B 56 1 189 G 0 2 189 G 14 5 189 G 28 9 189 G 7 1 189 G 42 7 189 G 56 8 190 A 0 3 190 A 14 4 190 A 28 5 190 A 7 8 190 A 42 6 190 A 56 7 191 G 0 6 191 G 14 10 191 G 28 8 191 G 7 13 191 G 42 6 191 G 56 10 192 vehicle 0 2 192 vehicle 14 7 192 vehicle 28 16 192 vehicle 7 5 192 vehicle 42 5 192 vehicle 56 5 193 A 0 3 193 A 14 5 193 A 28 2 193 A 7 1 193 A 42 5 193 A 56 5 194 D 0 1 194 D 14 9 194 D 28 8 194 D 7 9 194 D 42 1 194 D 56 5 195 C 0 3 195 C 14 6 195 C 28 5 195 C 7 5 195 C 42 10 195 C 56 8 196 vehicle 0 1 196 vehicle 14 6 196 vehicle 28 3 196 vehicle 7 2 196 vehicle 42 4 196 vehicle 56 6 197 A 0 6 197 A 7 8 197 A 56 15 198 B 0 7 198 B 14 10 198 B 28 7 198 B 7 2 198 B 42 7 198 B 56 11 199 B 0 2 199 B 14 11 199 B 28 13 199 B 7 9 199 B 42 15 199 B 56 10 200 vehicle 0 7 200 vehicle 14 7 200 vehicle 28 14 200 vehicle 7 11 200 vehicle 42 13 200 vehicle 56 10 201 D 0 7 201 D 14 2 201 D 28 3 201 D 7 6 201 D 42 7 201 D 56 8 202 A 0 7 202 A 14 6 202 A 28 3 202 A 7 7 202 A 42 12 202 A 56 1 203 D 0 2 203 D 14 8 203 D 28 1 203 D 7 3 203 D 42 18 203 D 56 5 204 vehicle 0 7 204 vehicle 14 20 204 vehicle 28 14 204 vehicle 7 7 204 vehicle 42 18 204 vehicle 56 10 205 F 0 2 205 F 14 5 205 F 28 3 205 F 7 5 205 F 42 8 205 F 56 4 206 vehicle 0 3 206 vehicle 14 1 206 vehicle 28 1 206 vehicle 7 0 206 vehicle 42 1 206 vehicle 56 0 207 vehicle 0 6 207 vehicle 14 5 207 vehicle 28 1 207 vehicle 7 5 207 vehicle 42 5 207 vehicle 56 2 208 F 0 5 208 F 14 0 208 F 28 7 208 F 7 0 208 F 42 7 208 F 56 0 209 G 0 7 209 G 14 1 209 G 28 20 209 G 7 12 209 G 42 12 209 G 56 5 210 vehicle 0 1 210 vehicle 14 1 210 vehicle 28 1 210 vehicle 7 1 210 vehicle 42 1 210 vehicle 56 1 211 A 0 2 211 A 14 13 211 A 28 5 211 A 7 6 211 A 42 5 211 A 56 6 212 A 0 4 212 A 14 2 212 A 28 5 212 A 7 2 212 A 42 4 212 A 56 4 213 vehicle 0 1 213 vehicle 14 1 213 vehicle 28 0 213 vehicle 7 1 213 vehicle 42 0 213 vehicle 56 2 214 D 0 1 214 D 14 3 214 D 28 20 214 D 7 18 214 D 42 14 214 D 56 12 215 F 0 2 215 F 14 2 215 F 28 0 215 F 7 4 215 F 42 0 215 F 56 5 216 C 0 2 216 C 14 6 216 C 28 5 216 C 42 3 216 C 56 3 217 F 0 2 217 F 14 2 217 F 28 1 217 F 7 1 217 F 42 5 217 F 56 5 218 A 0 5 218 A 14 7 218 A 28 4 218 A 7 5 218 A 42 8 218 A 56 5 219 C 0 1 219 C 14 8 219 C 28 8 219 C 7 5 219 C 42 3 219 C 56 8 220 F 0 2 220 F 14 1 220 F 28 1 220 F 7 2 220 F 42 2 220 F 56 1 221 F 0 4 221 F 14 1 221 F 28 5 221 F 7 5 221 F 42 12 221 F 56 2 222 vehicle 0 7 222 vehicle 14 8 222 vehicle 28 9 222 vehicle 7 10 222 vehicle 42 10 222 vehicle 56 10 223 B 0 1 223 B 14 1 223 B 7 1 223 B 56 5 224 F 0 7 224 F 56 10 225 B 0 2 225 B 14 5 225 B 28 10 225 B 7 5 225 B 42 12 225 B 56 9 226 A 0 7 226 A 14 9 226 A 28 28 226 A 7 1 226 A 42 3 226 A 56 12 227 B 0 1 227 B 14 20 227 B 28 14 227 B 7 4 227 B 42 7 227 B 56 9 228 D 0 5 228 D 14 2 228 D 28 2 228 D 7 5 228 D 42 10 228 D 56 9 229 F 0 2 229 F 14 8 229 F 28 8 229 F 7 3 229 F 42 5 229 F 56 5 230 vehicle 0 2 231 D 0 5 231 D 14 0 231 D 28 4 231 D 7 2 231 D 42 3 231 D 56 6 232 G 0 4 232 G 14 12 232 G 28 3 232 G 7 5 232 G 42 12 232 G 56 4 233 F 0 6 234 G 0 2 234 G 14 6 234 G 28 1 234 G 7 4 234 G 42 2 234 G 56 6 235 A 0 7 235 A 14 10 235 A 28 6 235 A 7 16 235 A 42 1 235 A 56 1 236 D 0 3 236 D 14 2 236 D 28 0 236 D 7 3 236 D 42 18 236 D 56 0 237 vehicle 0 1 237 vehicle 14 0 237 vehicle 28 0 237 vehicle 7 10 237 vehicle 42 6 237 vehicle 56 6 238 B 0 7 238 B 14 16 238 B 28 7 238 B 7 16 238 B 42 12 238 B 56 8 239 vehicle 0 4 239 vehicle 14 4 239 vehicle 28 6 239 vehicle 7 8 239 vehicle 42 2 239 vehicle 56 0 240 B 0 0 240 B 14 6 240 B 28 4 240 B 7 6 240 B 42 1 240 B 56 0 241 C 0 4 241 C 14 0 241 C 28 1 241 C 7 24 241 C 42 0 241 C 56 2 242 B 0 0 242 B 56 0 243 F 0 6 243 F 14 8 243 F 28 18 243 F 7 8 243 F 42 7 243 F 56 9 244 D 0 4 244 D 14 2 244 D 28 2 244 D 7 3 244 D 42 2 244 D 56 3 245 G 0 1 245 G 14 5 245 G 28 3 245 G 7 1 245 G 42 8 245 G 56 4 246 D 0 7 246 D 14 20 246 D 28 12 246 D 7 7 246 D 42 16 246 D 56 13 247 vehicle 0 5 247 vehicle 14 1 247 vehicle 28 4 247 vehicle 7 6.5 247 vehicle 42 3 247 vehicle 56 4 248 D 0 1 248 D 14 5 248 D 28 2 248 D 7 2 248 D 42 2 248 D 56 5 249 A 0 5 249 A 14 5 249 A 28 15 249 A 7 7 249 A 42 9 249 A 56 18 250 C 0 1 250 C 14 3 250 C 28 3 250 C 7 1 250 C 42 2 250 C 56 1 251 F 0 2 251 F 14 4 251 F 28 2 251 F 7 7 251 F 42 2 251 F 56 5 252 F 0 1 252 F 14 4 252 F 28 3 252 F 7 2 252 F 42 3 252 F 56 3 253 A 0 2 253 A 14 4 253 A 28 7 253 A 7 3 253 A 42 9 253 A 56 4 254 D 0 5 254 D 14 17 254 D 28 4 254 D 7 23 254 D 42 28 254 D 56 18 255 C 0 1 255 C 14 1 255 C 28 4 255 C 7 3 255 C 42 1 255 C 56 2 256 F 0 7 256 F 14 10 256 F 28 17 256 F 7 20 256 F 42 15 256 F 56 26 257 D 0 7 257 D 14 9 257 D 28 10 257 D 7 7 257 D 42 11 257 D 56 7 258 vehicle 0 7 258 vehicle 14 6 258 vehicle 28 5 258 vehicle 7 2 258 vehicle 42 3 258 vehicle 56 6 259 G 0 5 259 G 14 8 259 G 28 10 259 G 7 7 259 G 42 13 259 G 56 7.5 260 vehicle 0 7 260 vehicle 14 3 260 vehicle 28 6 260 vehicle 7 4 260 vehicle 42 4 260 vehicle 56 4 261 G 0 7 261 G 14 9 261 G 28 4 261 G 7 7 261 G 42 6.5 261 G 56 6 262 C 0 3 262 C 14 5 262 C 28 4 262 C 7 5 262 C 42 5 262 C 56 17 263 C 0 6 263 C 14 5.5 263 C 28 4 263 C 7 5 263 C 42 4 263 C 56 3 264 A 0 4.5 264 A 14 1 264 A 28 7 264 A 42 4 264 A 56 6 265 vehicle 0 4 265 vehicle 14 3 265 vehicle 28 4 265 vehicle 7 9 265 vehicle 42 5 265 vehicle 56 10 266 D 0 6 266 D 14 13 266 D 7 5 266 D 56 8 267 F 0 7 267 F 14 5 267 F 28 8 267 F 7 6 267 F 42 5 267 F 56 11 268 B 0 6.5 268 B 14 20 268 B 28 20 268 B 7 8 268 B 42 10 268 B 56 11 269 vehicle 0 4.5 269 vehicle 14 8 269 vehicle 28 2.5 269 vehicle 7 5 269 vehicle 42 6.5 269 vehicle 56 3 270 G 0 4.5 270 G 14 10 270 G 28 5 270 G 7 9 270 G 42 30 270 G 56 10 271 vehicle 0 7 271 vehicle 28 7 271 vehicle 7 15 271 vehicle 42 10 271 vehicle 56 10 272 F 0 4 272 F 14 27 272 F 28 30 272 F 7 21 272 F 42 18 272 F 56 26 273 D 0 4 273 D 14 10 273 D 28 16.5 273 D 7 15 273 D 42 21 273 D 56 16 274 C 0 4 274 C 14 5 274 C 28 13 274 C 7 12 274 C 42 5 274 C 56 7 275 A 0 4 275 A 14 2.5 275 A 28 4 275 A 7 12 275 A 42 3 275 A 56 3 276 F 0 5 276 F 14 11 276 F 28 15 276 F 42 16 276 F 56 5 277 G 0 7 277 G 14 33 277 G 28 10 277 G 7 20 277 G 42 12 277 G 56 10 278 C 0 1 278 C 14 8 278 C 28 1 278 C 7 1 278 C 56 7.5 279 A 0 5 279 A 14 0 279 A 28 11 279 A 7 0 279 A 42 0 279 A 56 6 280 F 0 5.8 280 F 14 15 280 F 28 8.8 280 F 7 5 280 F 42 10 280 F 56 11 281 F 0 6 281 F 14 11.8 281 F 28 35 281 F 7 16 281 F 42 12 281 F 56 12 282 C 0 5.8 282 C 14 10 282 C 28 5 282 C 7 5 282 C 42 18 282 C 56 35 283 C 0 6 283 C 14 15 283 C 28 5 283 C 7 11 283 C 42 8 283 C 56 21 284 G 0 7 284 G 14 35 284 G 28 35 284 G 7 35 284 G 42 18 284 G 56 35 285 D 0 6 285 D 14 16 285 D 28 11 285 D 7 6 285 D 42 13 286 vehicle 0 5.5 286 vehicle 14 19 286 vehicle 28 35 286 vehicle 7 22.8 286 vehicle 42 12 286 vehicle 56 8 287 C 0 0 287 C 14 6 287 C 28 0 287 C 7 6 287 C 42 5 287 C 56 15 288 vehicle 0 6 288 vehicle 14 8 288 vehicle 28 6 288 vehicle 7 7 288 vehicle 42 8 288 vehicle 56 7 289 B 0 5 289 B 14 5 289 B 28 6 289 B 7 6 289 B 42 8 289 B 56 6 290 C 0 3 290 C 14 3 290 C 28 2 290 C 7 3 290 C 42 2 290 C 56 3 291 A 0 6 291 A 14 7 291 A 28 15 291 A 7 5 291 A 42 5 291 A 56 24 292 G 0 6 292 G 14 5 292 G 28 8 292 G 7 5 292 G 42 8 292 G 56 6 293 B 0 6 293 B 14 12 293 B 28 7 293 B 7 7 293 B 42 12 293 B 56 16 294 A 0 3 294 A 14 1 294 A 28 3 294 A 7 4 294 A 42 2 294 A 56 2 295 D 0 5 295 D 14 7 295 D 28 5 295 D 7 6 295 D 42 7 295 D 56 6 296 C 0 6 296 C 14 6 296 C 28 7 296 C 7 7 296 C 42 6 296 C 56 5 297 A 0 6 297 A 14 8 297 A 28 10 297 A 7 7 297 A 42 7 297 A 56 8 298 F 0 6 298 F 14 7 298 F 28 5 298 F 7 7 298 F 42 5 298 F 56 7 299 B 0 5 299 B 14 10 299 B 28 7 299 B 7 6 299 B 42 8 299 B 56 11 300 D 0 6 300 D 14 4.5 300 D 28 7.5 300 D 7 5 300 D 42 5 300 D 56 4 301 G 0 5.5 301 G 14 5 301 G 28 5.5 301 G 7 5 301 G 42 10 301 G 56 10 302 vehicle 0 4 302 vehicle 14 9 302 vehicle 28 9 302 vehicle 7 2 302 vehicle 56 14 303 vehicle 0 3 303 vehicle 14 6 303 vehicle 28 3 303 vehicle 7 4 303 vehicle 42 2 303 vehicle 56 5 304 F 0 5.5 304 F 14 11 304 F 28 14 304 F 7 14 304 F 42 13 304 F 56 20 305 C 0 5.5 305 C 14 9 305 C 28 4.5 305 C 7 5 305 C 42 7 305 C 56 5 306 B 0 6.5 306 B 14 17 306 B 28 12.5 306 B 7 5.5 306 B 42 10 306 B 56 14 307 A 0 5 307 A 14 5 307 A 28 4 307 A 7 7.5 307 A 42 5 307 A 56 7.5 308 G 0 5 308 G 14 7.5 308 G 28 10 308 G 7 7 308 G 42 14 308 G 56 8.5 309 A 0 3 309 A 14 12 309 A 28 4 309 A 7 6 309 A 42 8.5 309 A 56 6 310 G 0 6.5 310 G 14 5.5 310 G 28 11 310 G 7 6 310 G 42 12.5 310 G 56 9.5 311 D 0 6.5 311 D 14 7 311 D 28 7.5 311 D 7 5 311 D 42 6 311 D 56 5 312 B 0 5.5 312 B 14 6.5 312 B 28 5.5 312 B 7 3 312 B 42 6.5 312 B 56 10.5 313 A 0 4 313 A 14 4 313 A 28 4 313 A 7 3 313 A 42 5 313 A 56 5 314 A 0 6 314 A 14 10 314 A 28 11 314 A 7 7 314 A 42 8 315 C 0 7 315 C 14 4 315 C 28 8 315 C 7 5 315 C 42 3 315 C 56 6 316 C 0 6 316 C 14 5 316 C 28 7 316 C 7 6 316 C 42 7 316 C 56 4 317 D 0 6 317 D 14 7 317 D 28 5 317 D 7 4 317 D 42 4 317 D 56 5 318 G 0 5 318 G 14 3 318 G 28 6 318 G 7 5 318 G 42 8 318 G 56 5 319 F 0 5 319 F 14 10 319 F 28 10 319 F 7 10 319 F 42 14 319 F 56 17 320 G 0 4 320 G 14 6 320 G 28 18 320 G 7 4 320 G 42 15.5 320 G 56 6 321 vehicle 0 3 321 vehicle 14 7 321 vehicle 28 5 321 vehicle 7 7 321 vehicle 42 6 321 vehicle 56 5 322 F 0 7 322 F 14 21 322 F 28 20 322 F 7 18 322 F 42 16 322 F 56 21 323 G 0 6 323 G 14 5 323 G 28 12 323 G 7 6 323 G 42 12 323 G 56 15 324 C 0 3 324 C 14 6.5 324 C 28 8 324 C 7 3 324 C 42 4 324 C 56 7 325 C 0 3 325 C 14 5 325 C 28 8.5 325 C 7 3 325 C 42 10 325 C 56 6.5 326 G 0 6.5 326 G 14 8 326 G 28 9 326 G 7 8.5 326 G 42 8 326 G 56 10 327 D 0 3 327 D 14 7 327 D 7 5
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.