all_ss 里的R2就是你要的值。

您好，还是要麻烦您，直接用nlin命令都不能直接输出R^2，您说复制代码，怎么就有了呢？能不能说具体点

nlin 不直接给出R2, 但是给出了 SS(sum of squares)。 我们先用 ods output 输出 nlin 计算出来的 SS，然后再自己计算出 R2 就行了。

您好，这几天，我一直在试，可是总也用不好ods output。您帮我改改吧！我就想输出您说的数据在计算。
ods output “summary ststisics”=\'e：\\ods\\case\'；
proc nlin data=sim；
var d r；（sim中只有d剂量，r效应）
run；
ods output close；
可是不行啊，怎么试都不对，您能帮我改改吗？拜谢

大神，我把内容输出来啦，可是怎么描述让它做除法呢？能不能给几句命令，提示一下。谢谢

ods output 生成 dataset, 不是生成文件。ods output ANOVA=simu_output 就是把 nlin 生成的 ANOVA表输出到 simu_output 这个 dataset 里去。
你的 proc nlin 里需要指定你的模型和参数。没有这些信息，SAS 就不知道你要什么样的曲线，哪些是曲线的参数。

谢谢大神，我终于把表格输出来啦！太感谢您了。可是，我要用model sum of squares 除以 corrected total 的sum of
squares这个命令怎么编辑呀！再谢

大神，我转置了生成的数据集，又用where命令挑出了ss行，有增加了一列，计算col1/col3。好麻烦，您有没有好办法？能少几步骤

没有，就是这么麻烦。
如果不想用转置可以看七楼的代码，先把两个 SS 分别作为两个 datasets 输出，然后merge 到一起再计算。复杂程度差不多，省不了多少事。

大神，很感谢您的帮助。
在使用proc nlin过程中，simu临时数据集的格式如下图所示可以吗？

要拟合曲线，先要有模型。你先要想想你的自变量和因变量的关系，用个数学公式表达出来，就像 y=f(x1,x2,x3....)。 有了这个，才能做 nlin。 我不知道你的模型是怎么样的，哪是自变量，哪是因变量，所以没法回答你的问题。

大神，您帮我看看哪里不对？怎么a5里的数都一样。拜谢
options nodate nonotes nosource nosource2;
%macro sim (e,f,name);
/*模拟100次*/
%do j=1 %to 100;
/*计算d1、d2、d3模拟效应值,存放在临时数据集simu,计算simu的拟合优度*/
DATA SIMU ;
DO n = 1 to &e ;
%do i= 0 %to 4 ;
d&i=0.1+&i*0.2;
R=1.002+9.023/(1+exp(10.47*(0.5-d&i)));
v&i=R+&f*R*rannor(&j);
%end ;
output ;
end;

/*合并simu集的数据，使变量为v,d,保存在数据集new4         */
data first1;
set simu ;keep d0 v0;
data second1;
set first1;
rename d0=d v0=v;


data first2;
set simu ;keep d1 v1;
data second2;
set first2;
rename d1=d v1=v;

data first3;
set simu ;keep d2 v2;
data second3;
set first3;
rename d2=d v2=v;


data first4;
set simu ;keep d3 v3;
data second4;
set first4;
rename d3=d v3=v;


data first5;
set simu ;keep d4 v4;
data second5;
set first5;
rename d4=d v4=v;



data new1;
set second1 second2 ;
data new2;
set new1 second3;
data new3;
set new2 second4;
data new4;
set new3 second5;


ods output anova=one;
proc nlin data=new4;
model v=a + b/(1 + exp((0.5-d)*c))  ;
parms a=1 b=9 c=10.5 ;
run;
ods output close;

proc transpose data=one out=two;
run;
data three;
set two;
where _NAME_='SS';
run;
DATA four;
SET three;
COL4=COL1/COL3;
RUN;
quit;
/*计算new4集拟合优度，追加在a5中         */
proc append data=four base=a5 ; run

%end ;
%mend sim ;
%sim (40,0.5,a5);

问题在这一行：  
proc append data=four base=a5 ; run

要在 run 的后面加上 ;

嗯嗯，您真厉害！真的解决问题了。大神，真心感谢您啊！可以称呼您高老师么？谢谢！

高老师，您好，帮我在看看，为什么提示“WARNING: PROC NLIN failed to converge”模拟2次就没有，10次就有。

%macro sim (e,k,f,name);
/*模拟100次*/
%do j=1 %to 100;
/*计算d1~~~dn模拟效应值,存放在临时数据集simu,计算simu的拟合优度*/
DATA SIMU ;
DO n=1 to &e ;
%do i= 0 %to &k ;
d&i=0.1+&i*0.8/&k;
R=1.002+9.023/(1+exp(10.47*(0.5-d&i)));
v&i=R+&f*R*rannor(&j);
%end ;
output ;
end;

%do n=0 %to &k;
data first&n;
set simu; keep d&n v&n;
data second&n;
set first&n;rename d&n=d v&n=v;
proc append data=second&n  base=new; %end;
run;quit;

ods output anova=one;
proc nlin data=new;
model v=a + b/(1 + exp((0.5-d)*c))  ;
parms a=1 b=9 c=10.5 ;
run;
ods output close;

proc transpose data=one out=two;
run;
data three;
set two;
where _NAME_='SS';
run;
DATA four;
SET three;
COL4=COL1/COL3;
RUN;
proc append base=aa data=four  ;
proc delete data=new;run;quit;
%end ;
proc export data=work.aa outfile="f:\sasdata\&name" dbms=excel2000 replace;
run;
%mend sim ;
%sim (25,7,1,aa38);


data a;
set aa;
proc means mean;var col4;
run;quit;

你的模拟数据引入的随机误差太大了，造成有时候没法把数据拟合到你的模型上去。把随机误差减一半试试， 就是把 v&i=R+&f*R*rannor(&j);   换成  v&i=R+&f*R*rannor(&j)*0.5;
还不行就把 0.5 换成 0.4, 0.3 等等。

高老师，我希望随机误差能为效应值的一半，比如：效应值是是10，模拟数据围绕在10+-5左右，这样误差也会太大吗？那这个rannonr是不是可以v=r+rannor*r*0.5？还是这个rannor本身就不对？

“随机误差能为效应值的一半”，是指随机误差的平均值，标准差还是最大值？我估计你要的是标准差。v=r+rannor*r*0.5 就是把你的随机误差的标准差控制在 r*0.5， 如果这是你想要的， 用这个公式计算就可以了。

谢谢，高老师，问题已解决，感谢您的耐心帮助