gistfile1.txt

#### DATA EXPLORATION ####
ggplot(df,aes(x=date_interval,
                        y=feed_provided, col=as.factor(lcn_id))) + 
  geom_point(alpha=0.1) + 
  theme_bw()+
  facet_wrap(~ogn_id)+
  theme(legend.position = "none")
df%>%
  filter(feed_provided_cum>0)%>%
  ggplot(.,aes(x=date_interval)) + 
  geom_point(aes(y=feed_provided_cum),alpha=0.1) + 
  theme_bw()+
  facet_wrap(~ogn_id)
df%>%
  ggplot(.,aes(x=curve_day)) + 
  geom_point(aes(y=feed_provided_cum),alpha=0.1) + 
  theme_bw()+
  facet_wrap(~ogn_id)
df%>%
  filter(feed_provided_cum>0)%>%
  ggplot(.,aes(x=curve_day)) + 
  geom_point(aes(y=aantal_voerpannen),alpha=0.1) + 
  theme_bw()+
  facet_wrap(~ogn_id)
df%>%
  filter(feed_provided_cum>0)%>%
  ggplot(.,aes(x=curve_day)) + 
  geom_point(aes(y=feed_provided_cum),alpha=0.1) + 
  theme_bw()+
  facet_wrap(~ogn_id)
df%>%
  filter(feed_provided_cum>0)%>%
  ggplot(.,aes(x=date_interval)) + 
  geom_point(aes(y=comp1_provided),alpha=0.1, color="green") + 
  geom_point(aes(y=comp2_provided),alpha=0.1, color="red") + 
  geom_point(aes(y=comp3_provided),alpha=0.1, color="orange") + 
  geom_point(aes(y=comp4_provided),alpha=0.1, color="blue") + 
  theme_bw()+
  facet_wrap(~ogn_id)
df%>%
  filter(!is.na(lcn_id))%>%
  ggplot(.,aes(x=date_interval)) + 
  geom_line(aes(y=feed_provided, 
                 group=lcn_id, 
                 colour=as.factor(lcn_id)),alpha=0.5) + 
  theme_bw()+
  theme(legend.position = "none")