I would like to add a factor N/P to the following data frame:
   time   subjects  value
1     0     sub1    10
2     1     sub1    11
3     2     sub1    12
...
10    9     sub1    19
11   10     sub1    20
12    0     sub2    20
13    1     sub2    21
...
21    9     sub2    29
22   10     sub2    30
23    0     sub3    30
24    1     sub3    31
25    2     sub3    32
...
33   10     sub3    40
so that it reads as follows
   time   subjects  value  status
1     0     sub1    10     N
2     1     sub1    11     N
3     2     sub1    12     N
...
10    9     sub1    19     N
11   10     sub1    20     N
12    0     sub2    20     P
13    1     sub2    21     P
...
21    9     sub2    29     P
22   10     sub2    30     P
23    0     sub3    30     N
24    1     sub3    31     N
25    2     sub3    32     N
...
33   10     sub3    40     N
The 'status' factors are stored so far in an array c("N","P","N",....,"P","P") for hundreds of subjects. 
 
     
    