####################### Funciones identificadoras de outliers en Boxplots en S-PLUS ###########
idenbox<-
function(y)
{   q1<-quantile(y)[2]
    q3<-quantile(y)[4]
    boxplot(y)
    may<-which(y>q3+(1.5*(q3-q1)))
    men<-which(y<q1-(1.5*(q3-q1)))
    for (i in 1:length(may)) {text(20,y[may[i]],y[may[i]])}
    for (i in 1:length(men)) {text(20,y[men[i]],y[men[i]])}
}

######### Ejemplo: idenbox(obs)

idenbox.numero<-
function(y)
{   q1<-quantile(y)[2]
    q3<-quantile(y)[4]
    boxplot(y)
    may<-which(y>q3+(1.5*(q3-q1)))
    men<-which(y<q1-(1.5*(q3-q1)))
    for (i in 1:length(may)) {text(20,y[may[i]],y[may[i]],label=may[i])}
    for (i in 1:length(men)) {text(20,y[men[i]],y[men[i]],label=men[i])}
}


####################################################################
     Función identificadora de outliers en boxplots en R
####################################################################

 identify(rep(1,length(obs)),obs) #### Identifica el número de observación