Package dplyr et modèle linéaire

[namaa ghassan]

Bonjour Jeunes Gens :p

Je suis un peu perdue là !!! je ne sais pas comment coder ce que j'ai en tête !!!

Alors, J'ai une base contenant 4 colonnes ( Code - Date - Moyenne Qte - Moyenne Prix ) et je veux regrouper les Mean Qte et prix du même mois dans un vecteur, puis appliquer un ML sur chaque vecteur et récupérer le a.
J'ai pensé à mutate mais après pour les mettre en vecteur je ne sais pas du tt comment faire !!! j'ai mis vect puis c puis stack mais rien ne marche hahaha !!

Voilà Un petit aperçu sur ce que j'ai:

Code : Tout sélectionner

> str(rslt)
'data.frame':   1464 obs. of  4 variables:
 $ CodeFam  : Factor w/ 4 levels "351","361","381",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ DateT    : Date, format: "2015-10-01" ...
 $ MeanQ    : num  NA NA NA NA NA NA NA NA NA NA ...
 $ MeanP    : num  NA NA NA NA NA NA NA NA NA NA ...


ce que je veux faire après avoir regrouper mes données

Code : Tout sélectionner

model <- lm(Q~P)
a <- cbind(a,coef(model)[2] )


D'avance merci :)

Cordialement,
Namaâ.

[Logez Maxime]

Bonjour,

Si par CodeFam et par mois et en distinguant les années :

Code : Tout sélectionner

rslt %>% mutate(Mois = format(DateT, "%Y-%m")) %>%
  group_by(CodeFam, Mois) %>%
    do(Mods = lm(MeanQ ~ MeanP, data = .)) %>% mutate(b = coef(Mods)[2])

Si juste les mois (pas CodeFam) :

Code : Tout sélectionner

rslt %>% mutate(Mois = format(DateT, "%Y-%m")) %>%
  group_by(Mois) %>%
    do(Mods = lm(MeanQ ~ MeanP, data = .)) %>% mutate(b = coef(Mods)[2])

Si tu veux pas distinguer les années alors tu remplaces "%Y-%m" par "%m".

Cordialement,
Maxime

[namaa ghassan]

Merciiiiiiiiiiiiiiiii bcp je vais essayer d'adapter ceci à mon problème :)

Cdt,
Namaâ.

[namaa ghassan]

Le code marche trooop bien :)

je voulais aussi faire la même chose mais l'appliquer juste sur deux dates consécutives !! je me suis inspirée d'un truc que t'avais posté

Code : Tout sélectionner

test2 <- tabfam %>% group_by(CodeFam)  %>%
        do(Mods = lm(LogQ ~ LogP, data = cbind(.[(1:2),]) )) %>%
   mutate(b = round(coef(Mods)[2],5))


mais apparemment il ne prend pas le cbind hhhhhhh !!

Cdt,
Namaâ