Exercise2: Write the following matlab code in R

yDataMatrix = zeros(n,m);

sStar_Hat = zeros(m,1);

for e = 1:m

uVec = rand(n,1);

yData = arrayfun(@(x)QuantileYLN(x,c,t,T,theta,sigma),uVec);

sStar_Hat(e) = sum(yData < c.*(T-t))/n;

yDataMatrix(:,e) = yData;

%waitbar(e/(3*m));

end