data rv1;

do i=1
to
3;

n=ranbin(0,10,0.6);

do j =1
to n-1;

x=100+1*rannor(0);

if i^=lag(i) then y=x;

else y+x;

output;

end;

output;

end;

run;

proc
print
data=rv1;

run;

Obs
i
n
j
x
y

1
1
5
1
100.383
100.383

2
1
5
2

100.172
200.555

3
1
5
3
100.211
300.766

4
1
5
4
101.065
401.831

5
1
5
5
101.065
401.831

6
2
7
1
99.876
99.876

7
2
7
2
100.306
200.182

8
2
7
3
99.236
299.419

9
2
7
4
98.282
397.700

10
2
7
5
100.696
498.396

11
2
7
6
100.921
599.317

12
2
7
7
100.921
599.317

13
3
5
1
99.924
99.924

14
3
5
2
101.336
201.261

15
3
5
3
99.292
300.552

16
3
5
4
100.994
401.546

17
3
5
5
100.994
401.546