亲，谁能帮我将下面的R 代码转成SAS？ 谢谢了
pown=function(alpha,sigma,c1,c2,power0){
  n=0
  power=0
  npow=c(0)
  while(power<power0){
    n=n+1
    Z_alpha=qt(1-alpha,n)
    Z1=Z_alpha-(c1-c2)*sqrt(n)/sigma
    Z2=Z_alpha-(c1+c2)*sqrt(n)/sigma
    power=1-(pt(Z1,n)+pt(Z2,n))
  }
  result=c(exp(-c1), exp(-c2), sigma,alpha, power, n)
  result
}
c1=c(-log(0.7),-log(0.8))
c2=c(-log(0.9),-log(0.95),-log(1))
alpha=0.05
sigma=c(0.34,0.29)
power0=c(0.9,0.8)
Res=matrix(0,nrow=24,ncol=6)
colnames(Res)=c("Upp Bound", "TrueRatio", "sigma","alpha", "power", "n")
i=0
for(i1 in 1:length(power0)){
  for(i2 in 1:length(sigma)){
    for(i3 in 1:length(c1)){
      for(i4 in 1:length(c2)){
        i=i+1
        Res[i,]=pown(alpha,sigma[i2],c1[i3],c2[i4],power0[i1])
      }
    }
  }
}