[R] error in function AMER

[Marilia Gmail]

Dear all,

I am trying to reproduce the example in the vignette "Using lme4 to fit 
Generalized Additive
Mixed Models" with my dataset.
But...

 > exemplo <- read.table(exemplo.dat,header=T)
 > mod <- amer(pasvig ~ -1 + harvf + tp(dias,by=harvf) + (1 | pac), 
data=exemplo)
Erro em if (from == to) rep.int(from, length.out) else as.vector(c(from,  :
 valor ausente onde TRUE/FALSE necessário

A subset of the dataset is below.
No clue about what is wrong.

Thanks a lot!

Marilia
"pasvig" "dias" "harvf" "pac"
124 0 "HRJB" 1
125 473 "HRJB" 1
125 888 "HRJB" 1
136 1272 "HRJB" 1
161 1614 "HRJB" 1
133 0 "HRV" 2
151 716 "HRV" 2
127 850 "HRV" 2
141 1137 "HRV" 2
135 1382 "HRV" 2
168 0 "HRV" 3
140 169 "HRV" 3
140 555 "HRV" 3
137 1259 "HRV" 3
147 1524 "HRV" 3
150 1842 "HRV" 3
133 1982 "HRV" 3
160 0 "HRV" 4
159 417 "HRV" 4
137 790 "HRV" 4
173 1146 "HRV" 4
149 1761 "HRV" 4
165 2135 "HRV" 4
142 2415 "HRV" 4
137 0 "HRV" 5
144 420 "HRV" 5
145 842 "HRV" 5
127 0 "HRJB" 6
124 365 "HRJB" 6
122 762 "HRJB" 6
138 1134 "HRJB" 6
124 1561 "HRJB" 6
133 0 "HRV" 7
116 502 "HRV" 7
142 1474 "HRV" 7
170 0 "HRV" 8
151 400 "HRV" 8
132 577 "HRV" 8
114 1169 "HRV" 8
146 1547 "HRV" 8
110 1904 "HRV" 8
156 0 "HRV" 9
132 1370 "HRV" 9
154 1784 "HRV" 9
148 2318 "HRV" 9
148 0 "HRV" 10
110 520 "HRV" 10
123 924 "HRV" 10
130 1640 "HRV" 10
142 2075 "HRV" 10
158 0 "HRV" 11
166 274 "HRV" 11
159 1059 "HRV" 11
150 0 "HRV" 12
140 333 "HRV" 12
121 704 "HRV" 12
138 1225 "HRV" 12
110 1406 "HRV" 12
151 0 "HRV" 13
179 318 "HRV" 13
167 913 "HRV" 13
181 1261 "HRV" 13
131 0 "HRJB" 14
123 630 "HRJB" 14
136 1038 "HRJB" 14
138 1463 "HRJB" 14
133 1757 "HRJB" 14
119 2158 "HRJB" 14
144 0 "HRV" 15
150 418 "HRV" 15
145 659 "HRV" 15
141 1099 "HRV" 15
133 1336 "HRV" 15
150 1757 "HRV" 15
140 0 "HRV" 16
127 432 "HRV" 16
118 NA "HRV" 16
151 0 "HRV" 17
140 720 "HRV" 17
150 1455 "HRV" 17
156 0 "HRV" 18
144 402 "HRV" 18
152 807 "HRV" 18
139 1163 "HRV" 18
127 1575 "HRV" 18
126 0 "HRV" 19
141 997 "HRV" 19
143 1505 "HRV" 19
153 1843 "HRV" 19
123 2023 "HRV" 19
135 0 "HRV" 20
128 357 "HRV" 20
135 880 "HRV" 20
151 1561 "HRV" 20
153 1952 "HRV" 20
131 2364 "HRV" 20
122 0 "HRJB" 21
141 462 "HRJB" 21
136 1224 "HRJB" 21
132 1757 "HRJB" 21