如何将R程序转成SAS？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

亲，谁能帮我将下面的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])
      }
    }
  }
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R程序 function Power Alpha Sigma function sigma power 程序 如何