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)

trimn=function(x)
		{
		trm=(quantile(x,0.25,type=6)+2*quantile(x,0.50,type=6)+quantile(x,0.75,type=6))/4
		return(trm)
		}

library(MASS)
wo =sort(wo) 
n =length(wo) 
wo = wo[-n] 
xmu = seq(-10,10,length=100) 
mm = md = tr = rep(0,100) 
 for(i in 1:100)
{ 
dc=c(wo,xmu[i])
mm[i] =n*( mean(dc)- mean(wo))
md[i]= n*( median(dc)- median(wo))
tr[i]=n*(trimn(dc)-trimn(wo))
}  
plot(xmu,mm, main=" Sensitivity Curve",ylim=c(-6,6),xlim=c(-10,10),type="l") 
lines(xmu,md,lty=2)
lines(xmu,tr, lty=3) 
 
legend(-10,4,c("mean","median","TriMean"),lty =c(1,2,3))