wo =c(2.2,2.2,2.4,2.4,2.5,2.7,2.8,2.9,3.03,3.03,3.1,3.37,3.4,3.4,3.4,3.5,3.6,3.7,3.7,3.7,3.7,3.77,5.28,28.95)

woo=c(2.2,2.2,2.4,2.4,2.5,2.7,2.8,2.9,3.03,3.03,3.1,3.37,3.4,3.4,3.4,3.5,3.6,3.7,3.7,3.7,3.7,3.77,5.28)



library(MASS)
wo =sort(wo) 
n =length(wo) 
wo = wo[-n] 
xmu = seq(-10,10,length=100) 
mm = md = qr = rep(0,100) 
 for(i in 1:100)
{ 
dc=c(wo,xmu[i])
mm[i] =n*( sd(dc)- sd(wo))
md[i]= n*( mad(dc)- mad(wo))
qr[i]=n*( IQR(dc)- IQR(wo))

}  
plot(xmu,mm, main=" Sensitivity Analysis",ylim=c(-4,8),xlim=c(-10,10),type="l") 
lines(xmu,md,lty=2)
lines(xmu,qr, lty=3)
legend(-10,8,c("STD","MADN","InterQR"),lty =c(1,2,3))