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