proc model下计算参数

现在有这公式：d=p**a , a是我所要求的参数。
我想在sas proc model  模块下完成对a 进行迭代求解，保证每个plot下每个group下的sum(d)=0.067*num, 并且0<d<1。
数据见附件（有很多个plot，这里我只列除了两三个）。请高手帮个忙，着急做出结果毕业，谢谢了！现在我的代码如下，有挺大问题的。
proc import datafile="c:\data.xls" replace
out=one;
getnames=yes;
data one; set one;
proc means noprint;  by plot group;
    var p;
    output out=two n=m1 ;
data three; merge one two; by plot group;
proc model data=three;by plot group;
sum=0;
do i=1 to m1;
sum=sum+p**a;
end;
eq.fn=0.067*num-sum;
solve a/ out=four ;

高手呢？帮个忙吧！