[R] loop calculation

[Ben Bolker]

  something like:
   
## set PIJ to be negative
while (any(PIJ<0)) {
   ## calculate PIJ
}