小弟不才 只会用sas做简单的处理,对于宏循环批量处理尝试了下论坛里的方法却总是报错,还请大大们指导。
以下是单个文件的导入和处理代码,希望能够实现批量导入output1.csv到output100.csv的文件,经过批量处理,导出ca,ka,ks三个参数的拟合结果到E:根目录下的excel格式文件里
proc import
datafile='E:\SOC\yanzheng\output1.csv'
out = test
dbms=csv replace;
getnames=yes;
run;
proc nlin data=test method=marquardt hougaard;
parameters Ka=0.001 to 0.02 by 0.001 Ks=0.0001 to 0.001 by 0.0001 Ca=10 to 1000 by 100;
Model
csoc=ca*exp(-ka*day)+(csoc-cr-ca)*exp(-ks*day)+Cr*exp(-kr*day);
der.ca=ka*exp(-ka*day)-ks*exp(-ks*day);
der.ka=ca*exp(-ka*day)-ca*exp(-ka*day)*ka*day;
der.ks=(csoc-cr-ca)*exp(-ks*day)-(csoc-cr-ca)*exp(-ks*day)*ks*day;
run;