Clarisse Burtin a écrit :d'accord je crois comprendre, mais sommes nous obligés de procéder ainsi pour calculer le nombre nb, est-ce fiable?
Personnellement je ne suis pas convaincu depuis le début, mais ce n'est pas une forum de statistiques, donc mes remarques sur la pertinence sont déjà un peu hors-sujet !
Clarisse Burtin a écrit :sinon comment faire directement par groupe de code articles avec tout mon fichier?
On crée une fonction avec nos instructions précédentes...
Code : Tout sélectionner
nb_autour_q75 <- function(p, probs = c(.7, .8)) {
q <- quantile(p, probs)
sum(p >= q[1] & p <= q[2])
}
... qu'on applique à chaque groupe
Code : Tout sélectionner
res1 <- tapply(donnees$prix, donnees$Article.Code, nb_autour_q75)
# pour changer l'intervalle par défaut
res2 <- tapply(donnees$prix, donnees$Article.Code, nb_autour_q75, probs = c(.65, .85))