Sorry The problem was with my eyes, not with R, predict works as it is supposed to on a subset, I was misled by the labels, which looked like they were giving the wrong thing. Sorry again Peter