 Table \[auc\_table\] are compared with other common link prediction methods in Table \[auc\_compare\]: tSVD [@dhillon], tKatz [@dunlavy], and a degree-based model where the probability of a link is approximated as $\mathbb P(A_{ij}=1)=1-\exp(-d^\mathrm{out}_id^\mathrm{in}_j)$ where $d^\mathrm{out}_i$ and $d^\mathrm{in}_j$ are the out-degree and in-degree of each node. Overall, when compared to the results in Table \[auc\_table\], the PMF models achieve impressive improvements over competing matrix factorisation techniques.

=1.2mm

[c c c c c]{} & &\
& & & &\
& & & &\
& & & &\
& & & &\

Seasonal modelling {#seasonal_sec}
------------------

To investigate dynamic modelling, binary adjacency matrices ${\mathbf}A_1,\dots,{\mathbf}A_{82}$ are constructed for each day across the train and test periods. The seasonal PMF model with the inclusion of covariates (SEPMF) is then compared against EPMF; for EPMF, the adjacency matrices are assumed to be independent realisations randomly generated from a fixed set of latent features. Due to a “9 day-80 hour” work schedule operated at LANL, whereby employees can elect to take vacation every other Friday, the seasonal period is assumed to be comprised of four segments: weekdays (Monday - Thursday), weekends (Saturday and Sunday), and two separate segments for alternating Friday’s. For each model, binary classification is performed using the model predictive scores calculated across the entire period. For the positive class, scores are calculated for all user-host pairs $(i,j)$ such that $A_{ijt}=1$ for at least one $t$ in the test set; for the negative class, a random sample of $(i,j)$ pairs such that $A_{ijt}=0$ for all $t$ in the test set are obtained, with sample size equal to three times the total number of observed links.

Table \[auc\_seasonal\] presents the resulting AUC scores. For both networks, the seasonal model does not globally outperform the extended PMF model for [*all links*]{}. However, improvements are obtained for prediction of the *new links*. One explanation for
Table   \[auc\_table\ ] are compared with other common link prediction method in Table   \[auc\_compare\ ]: tSVD [ @dhillon ], tKatz [ @dunlavy ], and a academic degree - based model where the probability of a liaison is estimate as $ \mathbb P(A_{ij}=1)=1-\exp(-d^\mathrm{out}_id^\mathrm{in}_j)$ where $ d^\mathrm{out}_i$ and $ d^\mathrm{in}_j$ are the out - degree and in - degree of each lymph node. Overall, when compared to the results in Table   \[auc\_table\ ], the PMF models achieve impressive improvement over competing matrix factorisation proficiency. 

 = 1.2 mm 

 [ c c c c c ] { } & & \ 
 & & & & \ 
 & & & & \ 
 & & & & \ 
 & & & & \ 

 Seasonal model { # seasonal_sec } 
 ------------------ 

 To investigate active modelling, binary adjacency matrices $ { \mathbf}A_1,\dots,{\mathbf}A_{82}$ are constructed for each sidereal day across the train and trial time period. The seasonal PMF model with the inclusion of covariates (SEPMF) is then compared against EPMF; for EPMF, the adjacency matrix are assumed to be independent realisations randomly generated from a repair set of latent features. Due to a “ 9 day-80 hour ” work schedule operated at LANL, whereby employee can elect to take vacation every other Friday, the seasonal period is assumed to be comprised of four segments: weekday (Monday - Thursday), weekend (Saturday and Sunday), and two disjoined segments for alternating Friday ’s. For each model, binary classification is performed using the model predictive scores calculate across the entire period. For the incontrovertible course, scores are calculated for all user - host pairs $ (i, j)$ such that $ A_{ijt}=1 $ for at least one $ t$ in the test set; for the damaging class, a random sample of $ (i, j)$ pairs such that $ A_{ijt}=0 $ for all $ t$ in the test set are receive, with sample size equal to three times the total numeral of observed connection. 

 mesa   \[auc\_seasonal\ ] presents the resulting AUC score. For both networks, the seasonal model does not globally outperform the extensive PMF model for [ * all links * ] { }. However, improvements are obtained for prediction of the * modern links *. One explanation for
 Tahle \[auc\_table\] are compared with other common linn predjction mdthods in Table \[auc\_compare\]: tSTD [@dyillob], tKatz [@dunlavy], and a aegree-basvd model qhert the probability of a likh is wpprmeimated as $\mathnb P(A_{ij}=1)=1-\exp(-d^\kathrm{out}_id^\matvro{iu}_j)$ where $d^\mathrm{out}_i$ and $d^\mathrm{in}_j$ wre the okt-degree and ig-deggeq of vagh node. Overall, when compared to the revults in Tablr \[auc\_table\], the PMF models afhiege impressive imprlvements ovgd cjnpeting matrkx factorisation technjques.

=1.2mm

[c c c c c]{} & &\
& & & &\
& & & &\
& & & &\
& & & &\

Sdasoncl modellint {#weadmnal_sec}
------------------

To iivestidate dynamic modellitg, binaty adjacency mstrmces ${\mathbf}A_1,\dots,{\mathbf}A_{82}$ ace constructed for ewch day awrkss the train and twst pgriodv. Thd sexsohak LMF mofel with the jnclusion od covariates (SEPMF) os nnen compared againft EPMF; for EPMF, the adjacency matrices age aasumed to be independenr realisations randompy generaeed from a fixed set of latent features. Due to a “9 day-80 iojr” cirk szyefule operated at LANL, whereby employees can ejscu tp take vacatiok every other Fridsy, tng seasonal perkod is asaumed to be comprided of sour wegments: reekcays (Monday - Thursday), weekebds (Saturday qnd Sunday), and two separate szgmentx for alternating Friday’s. Fox each model, binagy classirkcation is perfofmec gsing tht model predictive scores celculcted acruss jhe entyre period. For bve positive class, dcored dre calculwted for all user-host pairs $(i,j)$ such that $A_{ijt}=1$ fmr dt least one $b$ in the test sqt; for the negctive clcss, a fandom samile of $(i,j)$ pairs such that $A_{ijt}=0$ fos all $t$ in thx test see arw obrained, dkth sample sizr equal to three tines the total numbcr of kbserved links.

Tcylt \[ayc\_seasonal\] prexengs ehv rxsultytg AUC scorev. Fof buyh negworks, uke wdasomal model does not gnobamly outperform the ewtended PNF model for [*all linkx*]{}. However, improvemtnts ace obteined gor prediction of the *new links*. Ohe explanwtijn for
Table \[auc\_table\] are compared with other common methods Table \[auc\_compare\]: [@dhillon], tKatz [@dunlavy], the of a link approximated as $\mathbb where $d^\mathrm{out}_i$ and $d^\mathrm{in}_j$ are the and in-degree of each node. Overall, when compared to the results in Table the PMF models achieve impressive improvements over competing matrix factorisation techniques. =1.2mm [c c c]{} &\ & & &\ & & & &\ & & & &\ & & & &\ Seasonal {#seasonal_sec} ------------------ To investigate dynamic modelling, binary adjacency ${\mathbf}A_1,\dots,{\mathbf}A_{82}$ are constructed for day across the train and periods. seasonal PMF with inclusion covariates (SEPMF) is compared against EPMF; for EPMF, the adjacency matrices are assumed to be independent realisations randomly generated from fixed set features. Due a day-80 work schedule operated whereby employees can elect to take Friday, the seasonal period is assumed to be of four weekdays (Monday - Thursday), weekends (Saturday Sunday), and two separate segments for alternating Friday’s. each model, binary classification is performed using the model predictive scores calculated across the entire the positive class, scores calculated for all pairs such $A_{ijt}=1$ at least $t$ in the test set; for the negative class, a random of $(i,j)$ pairs such that $A_{ijt}=0$ for all $t$ in set obtained, with sample equal to three times total of observed links. Table the AUC networks, seasonal does not globally outperform extended PMF model for [*all However, improvements are obtained links*. One explanation for
 Table \[auc\_table\] are compared wIth other coMmon lInk PreDiCtioN metHods in Table \[auc\_COmpaRe\]: tSVD [@dhillon], tKatz [@dunlAvy], anD a DEgreE-BaSed moDel wherE ThE PRobAbIlIty Of A LiNk is aPprOximateD as $\mathbb P(a_{ij}=1)=1-\ExP(-d^\mathrm{out}_iD^\MaThrm{in}_j)$ wheRe $d^\Mathrm{out}_i$ anD $d^\mAthrm{iN}_j$ Are THe out-DegRee anD in-degREe of eaCh node. OveRaLL, when cOMpared tO THe ResuLts in Table \[auc\_tablE\], ThE pMF models achieVe imprEsSIvE IMprOveMents over cOmPetinG Matrix fACtORISatIOn techniques.

=1.2mM

[c c c c c]{} & &\
& & & &\
& & & &\
& & & &\
& & & &\

SeasoNAl mOdelliNg {#SeaSOnal_seC}
------------------

To inVeSTigAte dynamic mOdelLing, binarY adjacENcy matrICes ${\mathBf}A_1,\dotS,{\maThbF}A_{82}$ arE CoNsTruCtED foR EaCh dAY acRoss the tRaIn And teSt peRIODS. The SeaSonaL PMF mOdel with the inCluSion OF coVariaTes (SEpMF) iS tHen coMpared AgainSt ePMF; for EPMF, the aDjacEncy matriCes ArE asSuMed to BE indepEndEnt RealisaTions raNDomLy GENErAted from a fixed set oF lATEnT featureS. Due to A “9 DaY-80 hOUr” work scHeDulE opeRATed at lANL, WHeReby emplOyees cAN eLeCt to takE vAcatioN eVerY otHer FrIDay, tHe seasOnal periOd is aSSumed to be comprISed of four segmENtS: WEeKDays (monDay - Thursday), WeekENds (SAturDAy And sUnday), And twO sEPaRAte segments for alterNaTing FrIday’s. for each model, bInary classIFICation is PerfORmED using the model PrediCtive scoreS CalculatEd acrOss the enTire perioD. fOr the posItiVe cLasS, scOREs Are calculated FOR all UsEr-host pAirS $(i,j)$ such ThaT $A_{iJt}=1$ fOr aT lEast one $t$ iN the test SeT; fOr ThE neGativE Class, a raNdOm sAmPle Of $(i,j)$ pAIrs sucH that $a_{ijt}=0$ FoR aLL $t$ iN the tesT SeT ARe obTaInEd, wiTh sAmPle siZe eqUAl tO three tImes the toTal NUmbeR oF oBserved Links.

Table \[auc\_SeAsonal\] presEnTs tHe resuLTIng AUC scOres. For both networks, the sEAsonal mOdeL does Not gLobally ouTpeRform tHe eXTended pMF modEl for [*AlL liNKS*]{}. HoweVER, iMprOvEments are oBTAinEd for PrEdicTion of tHe *new links*. One explaNAtiOn for
 Table \[auc\_table\] arecompared w ith o the r c om monlink  prediction me t hods  in Table \[auc\_compa re\]:  t S VD [ @ dh illon ], tKat z  [ @ d unl av y] , a nd adegre e-b ased mo del wherethe  p robability o f  a  link is a ppr oximated as$\m athbbP( A_{ i j}=1) =1- \exp( -d^\ma t hrm{ou t}_id^\ma th r m{in}_ j )$ wher e $d ^\ma thrm{out}_i$ and$ d^ \ mathrm{in}_j$are th eo ut - d egr eeand in-deg re e ofe ach nod e .O v e ral l , when compar ed to the r e sul ts inTa ble \[auc\ _tabl e\ ] , t he PMF mode ls a chieve im pressi v e impro v ementsover c omp eti ng m a tr ix  fa ct o ris a ti ont ech niques.
= 1. 2mm
[c c c c  c]{ } &  &\& & &  &\
& & & &\& &  & & \ 
&& & & \

Se ason al  mode lling{#sea so nal_sec}
------ ---- --------
To  i nve st igate dynami c m ode lling,binarya dja ce n c y  m atrices ${\mathbf} A_ 1 , \d ots,{\ma thbf}A _ {8 2} $  are con st ruc tedf o r eac h da y  a cross th e trai n  a nd  test p er iods.Th e s eas onalP MF m odel w ith theinclu s ion of covaria t es (SEPMF) is th e n  c o mpar edagainst EPM F; f o r EP MF,t he  ad j acenc y mat ri c es are assumed to be i nd epende nt re alisations ra ndomly gen e r a ted from  a f i xe d  set of latent  feat ures. Duet o a “9 d ay-80  hour” w ork sched u l e operat edatLAN L,w h er eby employees c an e le ct to t ake  vacati oneve ryoth er  Friday,the seas on al  p er iod  is a s sumed to  b e c om pri sed o f  foursegme nts:  w ee k day s (Mond a y- Thur sd ay ), w eek en ds (S atur d ayand Sun day), and  tw o  sep ar at e segme nts for alter na ting Frida y’ s.For ea c h  model,binary classification i s  perfor med  usin g th e model p red ictive  sc o res ca lculat ed ac ro sst h e ent i r eper io d. For the p osi tivecl ass,  scores  are calculated fo r  al l user-host p air s $( i , j) $ s u ch tha t$ A_{ i j t}=1$ for at le ast one $t $i nthe test s e t;fo r the n egative  clas s , a ran dom sampl e of $(i, j) $ pa i r s s uch that $ A_{ijt}= 0$ for al l  $t$i nthe t est  set a re  ob taine d, wit h  sa mplesize e qu al tothree  t imes the  total number of observ ed lin ks.
Tab le \[auc\ _se a son al\] pres ents  the resul tin g A UC sc ore s . For  bot h  n etw o rks,thes easonal m o de l d o e snot globall y o utp erfor m t h e exte nded  PMF model for [* a ll links*]{}.Howe v e r,imp r ovem en ts are obtaine d f or p redictio nof the *new  links*.  O n e exp lanati on for 
 Table \[auc\_table\]_are compared_with other common link_prediction methods_in_Table \[auc\_compare\]: tSVD_[@dhillon],_tKatz [@dunlavy], and_a degree-based model_where the probability of_a link is_approximated_as $\mathbb P(A_{ij}=1)=1-\exp(-d^\mathrm{out}_id^\mathrm{in}_j)$ where $d^\mathrm{out}_i$ and $d^\mathrm{in}_j$ are the out-degree and in-degree of each_node._Overall, when_compared_to_the results in Table \[auc\_table\], the_PMF models achieve impressive improvements_over competing_matrix factorisation techniques.

=1.2mm

[c c c c c]{} &_&\
&_& & &\
&_& & &\
& & & &\
& & & &\

Seasonal_modelling {#seasonal_sec}
------------------

To investigate dynamic modelling, binary_adjacency matrices ${\mathbf}A_1,\dots,{\mathbf}A_{82}$_are_constructed_for each day across_the train and test periods. The_seasonal PMF model with the inclusion_of covariates (SEPMF) is then compared against_EPMF; for EPMF, the adjacency matrices_are assumed to be independent_realisations randomly_generated from a fixed set_of latent features._Due to_a “9 day-80_hour” work schedule operated at LANL,_whereby employees can_elect to take vacation every other_Friday,_the seasonal period_is_assumed_to be_comprised of four_segments:_weekdays (Monday_-_Thursday), weekends (Saturday and Sunday), and_two_separate segments for alternating Friday’s. For each_model, binary classification is_performed_using the model predictive_scores calculated across the entire_period. For the positive class, scores_are calculated_for all_user-host pairs $(i,j)$ such that $A_{ijt}=1$ for at least one $t$_in the test set; for the_negative class, a random_sample of_$(i,j)$_pairs such that_$A_{ijt}=0$_for all_$t$ in the test set are obtained,_with sample_size equal to three times the_total number of observed_links.

Table \[auc\_seasonal\]_presents the resulting AUC scores. For_both networks, the seasonal model does_not globally outperform the extended_PMF_model_for [*all links*]{}. However, improvements_are obtained for prediction of the_*new links*. One_explanation for