请教一个sas计算的问题

抱歉，那个命令错了一个字母。请用如下命令。

title"Computation of min sample size for test of fit";
data one;
rmsea0=.05;*null hyp rmsea;
rmseaa=.08;*alt hyp rmsea;
d=15;*degrees of freedom;
alpha=.05;*alpha level;
powd=.80;*desired power;
*initialize values;
powa=0.0;
n=0;
*begin loop for finding initial level of n;
do until (powa>powd);
n+100;
ncp0=(n-1)*d*rmsea0**2;
ncpa=(n-1)*d*rmseaa**2;
*compute power;
if rmsea0>rmseaa then do;
cval=cinv(alpha,d,ncp0);
powa=probchi(cval,d,ncpa);
end;
if rmsea0<rmseaa then do;
cval=cinv(1-alpha,d,ncp0);
powa=1-probchi(cval,d,ncpa);
end;
end;
*begin loop for interval halving;
dir=-1;
newn=n;
intv=200;
powdiff=powa-powd;
do until(powdiff<.001);
intv=intv*.5;
*compute new power;
ncp0=(newn-1)*d*rmsea0**2;
ncpa=(newn-1)*d*rmseaa**2;
*compute power;
if rmsea0>rmseaa then do;
cval=cinv(alpha,d,ncp0);
powa=probchi(cval,d,ncpa);
end;
if rmsea0<rmseaa then do;
cval=cinv(1-alpha,d,ncp0);
powa=1-probchi(cval,d,ncpa);
end;
powdiff=abs(powa-powd);
if powa<powd then dir=1;else dir=-1;
end;
minn=newn;
output;
proc print data=one;
var rmsea0 rmseaa powd alpha d minn powa; run;

liuliqingdong 发表于 2012-9-28 11:35
抱歉，那个命令错了一个字母。请用如下命令。

title"Computation of min sample size for test of fit"; ...

还是没有run statement

liuliqingdong 发表于 2012-9-28 11:35
抱歉，那个命令错了一个字母。请用如下命令。

title"Computation of min sample size for test of fit"; ...

另外，如果run得位置如下：
title"Computation of min sample size for test of fit";
data one;
rmsea0=.05;*null hyp rmsea;
rmseaa=.08;*alt hyp rmsea;
d=15;*degrees of freedom;
alpha=.05;*alpha level;
powd=.80;*desired power;
*initialize values;
powa=0.0;
n=0;
*begin loop for finding initial level of n;
do until (powa>powd);
n+100;
ncp0=(n-1)*d*rmsea0**2;
ncpa=(n-1)*d*rmseaa**2;
*compute power;
if rmsea0>rmseaa then do;
cval=cinv(alpha,d,ncp0);
powa=probchi(cval,d,ncpa);
end;
if rmsea0<rmseaa then do;
cval=cinv(1-alpha,d,ncp0);
powa=1-probchi(cval,d,ncpa);
end;
end;
*begin loop for interval halving;
dir=-1;
newn=n;
intv=200;
powdiff=powa-powd;
do until(powdiff<.001);
intv=intv*.5;
*compute new power;
ncp0=(newn-1)*d*rmsea0**2;
ncpa=(newn-1)*d*rmseaa**2;
*compute power;
if rmsea0>rmseaa then do;
cval=cinv(alpha,d,ncp0);
powa=probchi(cval,d,ncpa);
end;
if rmsea0<rmseaa then do;
cval=cinv(1-alpha,d,ncp0);
powa=1-probchi(cval,d,ncpa);
end;
powdiff=abs(powa-powd);
if powa<powd then dir=1;else dir=-1;
end;
minn=newn;
output;
run;

不好意思，公司的电脑很久没有给出答复，可能数据计算量太大

sunset1986 发表于 2012-9-28 12:15
另外，如果run得位置如下：
title"Computation of min sample size for test of fit";
data one;

没关系。我刚刚下载了SAS，我再试试。不过还是非常感谢你。

liuliqingdong 发表于 2012-9-28 12:23
没关系。我刚刚下载了SAS，我再试试。不过还是非常感谢你。

别客气，多交流~
您这个是在作哪方面的分析呢？

主要是有关SEM中，power analysis和sample size estimation相关的分析。目前一个project需要用到这个。你是哪个行业或领域的。我是psychology的。

data test;
input acct $ date amount;
cards;
10001 201201 130
10001 201201 140
10001 201301 124
10002 201201 130
10002 201201 120
        ;
run;
proc sql;
create table chk as
select distinct acct , sum(amount) as amount
from test
group by acct
;
quit;

liuliqingdong 发表于 2012-9-28 12:31
主要是有关SEM中，power analysis和sample size estimation相关的分析。目前一个project需要用到这个。你是 ...

我是做金融风险分析的~ 呵呵

用sql, group by语句