gistfile1.txt

arc_month<-arc%>%
  summarise_by_time(sampledate,
                    .by="month",
                    Weight=median(wt.s,na.rm=TRUE),
                    rochef=median(rochefant, na.rm=TRUE),
                    nirala=sum(`nir-ala`, na.rm=TRUE),
                    nirdha=median(`nir-dha`, na.rm=TRUE),
                    nirdpa=sum(`nir-dpa`, na.rm=TRUE),
                    nirepa=median(`nir-epa`,na.rm=TRUE),
                    nireta=sum(`nir-eta`,na.rm=TRUE),
                    nirfasummon=median(`nir-fasummon`,na.rm=TRUE),
                    nirfasumo3=median(`nir-fasumo3`,na.rm=TRUE),
                    nirfasumo6=median(`nir-fasumo6`,na.rm=TRUE),
                    nirfasumsatfa=median(`nir-fasumsatfa`,na.rm=TRUE),
                    nirastax=median(nirastax., na.rm=TRUE),
                    nirfat=median(nirfat, na.rm=TRUE),
                    nirtotpig=median(nirtotpig, na.rm=TRUE),
                    nirlino=median(`nir-linoleicacid`, na.rm=TRUE),
                    cfactor=median(c.factor, na.rm=TRUE),
                    slaughtloss=median(slaug.loss, na.rm=TRUE),
                    length=median(length, na.rm=TRUE),
                    weightslaughter=median(weig.slaug, na.rm=TRUE),
                    weightwhole=median(weig.whole, na.rm=TRUE),
                    dayssea=median(log(daysinsea), na.rm=TRUE),
                    .type="floor")
arc_month%>%filter(!is.na(sampledate))%>%
                     plot_time_series(.data=.,
                                      .date_var=sampledate,
                                      .value=rochef,
                                      .interactive=FALSE)