gistfile1.txt

## Create summary time-series datasets for quick relationship building
length(levels(arc$customer))
length(levels(arc$locality))
length(levels(arc$region))
length(levels(arc$strain))
length(levels(arc$generation1))
length(levels(arc$`generation(a/s)`))
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(daysinsea, na.rm=TRUE),
                    .type="floor")
## Create summary data by Customer and date
arc_company<-arc%>%group_by(customer)%>%
  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),
                    weightclass=median(weightclasshalv, na.rm=TRUE),
                    dayssea=median(daysinsea, na.rm=TRUE),
                    .type="floor")
plot(arc_month[, c(3:10)])
plot(arc_company[, c(3:10)])