程序有句話不懂，求解？能改編不？

是 否 +2 論壇幣

k人 參與回答

經(jīng)管之家送您一份

應(yīng)屆畢業(yè)生專(zhuān)屬福利!

求職就業(yè)群

趙安豆老師微信：zhaoandou666

經(jīng)管之家聯(lián)合CDA

送您一個(gè)全額獎(jiǎng)學(xué)金名額~ !

立即領(lǐng)取

感謝您參與論壇問(wèn)題回答

經(jīng)管之家送您兩個(gè)論壇幣！

+2 論壇幣

/* Compute the moving average within a BY-group of last n observations.
  For the first n-1 observations within the BY-group, the moving average
  is set to missing. */
data ds1;
do patient='A','B','C';
do month=1 to 7;
  num=int(ranuni(0)*10);
  output;
end;
end;
run;
proc sort;
by patient;
%let n = 4;
data ds2;
set ds1;
by patient;
retain num_sum 0;
if first.patient then do;
  count=0;
  num_sum=0;
end;
count+1;
last&n=lag&n(num);
if count gt &n then num_sum=sum(num_sum,num,-last&n);
else num_sum=sum(num_sum,num);
if count ge &n then mov_aver=num_sum/&n;
else mov_aver=.;
run;

proc print;
run;

不知道紅色那句是什么意思呢？

根據(jù)這個(gè)程序，能不能改編成moving standard deviation of last n observation.覺(jué)得貌似可以，但是好像又沒(méi)那么簡(jiǎn)單。不知大家有何見(jiàn)解？

掃碼加我 拉你入群

請(qǐng)注明：姓名-公司-職位

以便審核進(jìn)群資格，未注明則拒絕

分享0 收藏0 回帖

關(guān)鍵詞：observations observation Deviation observat Standard average moving within count 程序

不懂也頂一下~~~~

其實(shí)見(jiàn)識(shí)一個(gè)SUM函數(shù)  求和罷了

last&n=lag&n(num);   *//取前第n個(gè)數(shù)lastn//;
if count gt &n then num_sum=sum(num_sum,num,-last&n);  *//by分組，如果count>n,則求移動(dòng)累和:加上當(dāng)前觀測(cè)，減去第前n個(gè)數(shù)lastn//;
else num_sum=sum(num_sum,num);    *//如果count<=n，則直接求累和，不必減去前第n個(gè)值lastn//;
if count ge &n then mov_aver=num_sum/&n;   *//移動(dòng)累和轉(zhuǎn)為移動(dòng)平均//;
else mov_aver=.;  //前n個(gè)觀測(cè)置為缺失//;

注意：程序里分組的第一個(gè)觀測(cè)，count＝0，而不是1。（其實(shí)應(yīng)該為1好點(diǎn)）

移動(dòng)標(biāo)準(zhǔn)差的話,把觀測(cè)值的平方進(jìn)行移動(dòng)累和，減去n被移動(dòng)平均值，最后除以n就行了。