求助关于sas中的循环问题

布莱特

如：现有

No cost

1 3767.11

2 4683.89

3 3111.59

4 1111.03

5 4030.38

6 3791.35

7 5361.34

8 1380.84

9 7691.02

10 2079.44

11 3890.88

12 623.45

13 5287.07

14 653

15 4426.4

16 11651.7

17 1569.17

18 4751.24

19 1850.73

20 3267.78

第一个问题：我如果要生成一个变量m，它等于连续两个cost值得平均值。如m2=（cost1+cost2）/2 然后m3=(cost2+cost3)/2；这样一直循环下去.

第二个问题，生成一变量r；当cost得值按顺序如果成同方向变化的时候，等于0，反方向变化的时候等于1。如这里数据中cost2>cost1，如cost3>cost2，则等于0，如cost3<cost2则等于１，这个时候形成新的方向；如cost4<cost3则等于0，如cost4>cost3则等于1，依此类推循环。

waterlwh

<br />
/*****************************************Q1**********************************/<br />
data a;<br />
input cost;<br />
cards;<br />
3767.11<br />
4683.89<br />
3111.59<br />
1111.03<br />
4030.38<br />
3791.35<br />
5361.34<br />
1380.84<br />
7691.02<br />
2079.44<br />
3890.88<br />
623.45<br />
5287.07<br />
653<br />
4426.4<br />
11651.7<br />
1569.17<br />
4751.24<br />
1850.73<br />
3267.78<br />
;<br />
run;<br />
proc transpose data=a out=b(drop=_name_) prefix=c;<br />
run;<br />
data ab(keep=m);<br />
set b;<br />
array aa(*) _numeric_;<br />
do i=2 to dim(aa);<br />
m=mean(aa(i-1),aa(i));<br />
output;<br />
end;<br />
run;<br />
data test1;<br />
merge a ab;<br />
run;<br />
<br />
/*********************************************Q2************************************/<br />
data temp(keep=r);<br />
set b;<br />
array bb(*) _numeric_;<br />
do i=3 to dim(bb);<br />
if (bb(i-2)-bb(i-1))*(bb(i-1)-bb(i))>0 then do; r=0 ;output;end;<br />
else do; r=1;output;end;<br />
end;<br />
run;<br />
<br />
data test2;<br />
merge a temp;<br />
run;<br />

fairy

这个问题，其实用lagn（x）函数很容易就解决了：

data a1;

input no $ cost;

cards;

1 3767.11

2 4683.89

3 3111.59

4 1111.03

5 4030.38

6 3791.35

7 5361.34

8 1380.84

9 7691.02

10 2079.44

11 3890.88

12 623.45

13 5287.07

14 653

15 4426.4

16 11651.7

17 1569.17

18 4751.24

19 1850.73

20 3267.78

;

run;

data a2;

set a1;

m=(cost+lag1(cost))/2;

if cost ge lag1(cost) then r=0;

else r=1;

run;

waterlwh

[quote]引用第2楼fairy于2007-10-25 17:54发表的“”:

这个问题，其实用lagn（x）函数很容易就解决了：

data a1;

input no $ cost;

cards;

1 3767.11

.......[/quote]

第一个问题可以这样做,第二个问题楼上的需要再看看

fairy

[quote]引用第3楼waterlwh于2007-10-25 21:29发表的“”:

第一个问题可以这样做,第二个问题楼上的需要再看看[/quote]

嗯，没看后面的变化了，那就结合2阶函数一起使用好了：

data a2;

set a1;

m=(cost+lag1(cost))/2;

if (cost-lag1(cost))*(lag1(cost)-lag2(cost)) gt 0 then r=0;

else if (cost-lag1(cost))*(lag1(cost)-lag2(cost)) lt 0 then r=1;

else r=2;

if no in ('1','2') then r=.;

run;

hongtianli

data a;

input No cost;

costnew=lag(cost);

m=(costnew+cost)/2;

if dif(cost)>0 then r=0;

else r=1;

cards;

1 3767.11

2 4683.89

3 3111.59

4 1111.03

5 4030.38

6 3791.35

7 5361.34

8 1380.84

9 7691.02

10 2079.44

11 3890.88

12 623.45

13 5287.07

14 653

15 4426.4

16 11651.7

17 1569.17

18 4751.24

19 1850.73

20 3267.78

;

proc print;run;

试试这个！

王笑权

data a;

input No cost;

mean=lag(cost)+dif(cost)/2;

if mean= . then mean= cost;

if mean^= . then mean=mean;

if dif(cost)>0 then r1=0;

else r1=1;

if dif2(cost)> dif(cost) then r2=0;

else r2=1;

if r1=1 & r2=0|r1=0 &r2=1 then r=0;

if r1=1 & r2=1|r1=0 & r2=0 then r=1;

drop no y r1 r2;

cards;

1 3767.11

2 4683.89

3 3111.59

4 1111.03

5 4030.38

6 3791.35

7 5361.34

8 1380.84

9 7691.02

10 2079.44

11 3890.88

12 623.45

13 5287.07

14 653

15 4426.4

16 11651.7

17 1569.17

18 4751.24

19 1850.73

20 3267.78

;

proc print;

run;

结果

Obs cost mean r

1 3767.11 3767.11 1

2 4683.89 4225.50 0

3 3111.59 3897.74 0

4 1111.03 2111.31 1

5 4030.38 2570.71 0

6 3791.35 3910.87 0

7 5361.34 4576.35 0

8 1380.84 3371.09 0

9 7691.02 4535.93 0

10 2079.44 4885.23 0

11 3890.88 2985.16 0

12 623.45 2257.17 0

13 5287.07 2955.26 0

14 653.00 2970.04 0

15 4426.40 2539.70 0

16 11651.70 8039.05 1

17 1569.17 6610.44 0

18 4751.24 3160.21 0

19 1850.73 3300.99 0

20 3267.78 2559.26 0

http://www.pinggu.org/bbs/dispbbs.asp?boardid=68&id=271259