“…Note: Here, pm is used for proposed model. library(AdequacyModel) data=c(14, 9,8,5,20,1,3333,9,4,178,240,3993,69,72,93,6,7,1,4,19,65,3,2,60,182,6,41,6,200,7,1,2,5,30,812,12,3) sort(data) data ################# PDF of the proposed model pdf_pm <function(par, plot(x,ecdf,l-ty=1,lwd=4,type="s",xlab="x",ylab="G(x; 0.3147927, 0.2818076)", ylim=c(0,1),xlim=c(min(x),max(x)),col="black") par(new=TRUE) plot(x,proposedcdf,xlab="x",ylab="G(x; 0.3147927, 0.2818076)", ylim=c(0,1),xlim=c(min(x),max(x)),col="red",lty =5, lwd=4,type="l") par(new=TRUE) legend(1500, 0.4,c("Real Data","Proposed Model"),-col=c(1,2), lty =c(1,5), bty="n", cex=1.2) B.3. R Code for Plotting the Fitted Survival Function.…”