Преобразование кода sas в код R
Привет!
Я хотел бы преобразовать приведенный ниже код SAS в R-код. Это было бы большим подспорьем для меня!
Заранее спасибо!
%пусть NN=1000;
proc iml;
cov={1 .1 .1 .1 .1,
.1 1 .1 .1 .1,
.1 .1 1 .1 .1,
.1 .1 .1 1 .1,
.1 .1 .1 .1 1
};
M=5;
h={[5].2};
Rcvr={[5].2};
T={[&NN]0};
VL={[&NN]0};
VR={[&NN]0};
s={[&NN]0};
Р=0.05;
rv=rannor(repeat(1200,&NN,5));
Г=РВ*(корень(ков));
спред=0;
i=1;
j=1;
do while (i<=&NN);
do while (j<=5);
x=cdf('Normal',Y[i,j]);
Y[i,j]= -(log(1-x))/h[j];
j=j+1;
конец;
T[i]=min(Y[i,]);
j=1;
к=0.5;
делать пока (к&ЛТ;Т[Я] &ампер; к&ЛТ;=м);
VL[i]=VL[i]+EXP(-r*k);
к=к+0.5;
конец;
do while(j<=5);
если (Т[Я]=Г[Я,J] & Т[Я] И Л;М), то ВР[я]=ВР[я]+(1-Rcvr[Дж])*(ЕХР(-р*т[я]));
j=j+1;
конец;
j=1;
если VL[i]^=0, то s[i]=VR[i]/VL[i];
спред=спред+s[i];
i=i+1;
конец;
spread=спред/&NN;
распространение печати;
бежать;
выходить;
Что я уже пробовал:
Я stucked в несколько шагов, в моих усилий, чтобы переводить код
Richard MacCutchan
К сожалению, этот сайт не предоставляет код или преобразование кода на заказ.