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.