Lognormal distribution in the digg online social network

Size: px
Start display at page:

Download "Lognormal distribution in the digg online social network"

Transcription

1 Eur. Phys. J. B 83, (2011) DOI: /epjb/e Lognorml distribution in the digg online socil network P. Vn Mieghem, N. Blenn nd C. Doerr

2 Eur. Phys. J. B 83, (2011) DOI: /epjb/e Regulr Article THE EUROPEAN PHYSICAL JOURNAL B Lognorml distribution in the digg online socil network P. Vn Mieghem,N.Blenn,ndC.Doerr Delft University of Technology, Fculty of EEMCS, 2628 CD Delft, Netherlnds Received 16 Februry 2011 / Received in finl form 2 June 2011 Published online 15 September 2011 c EDP Sciences, Società Itlin di Fisic, Springer-Verlg 2011 Abstrct. We nlyse the number of votes, clled the digg vlue, which mesures the impct or populrity of submitted informtion in the Online Socil Network Digg. Experiments over five yers indicte tht the digg vlue of story on the first frontpge follows closely lognorml distribution. While the lw of proportionte effect explins lognorml behvior, the proportionlity fctor in tht lw is ssumed to hve constnt men, wheres experiments show tht decreses linerly with time. Our hypothesis, the probbility tht user diggs (votes) on story given tht he observes certin digg vlue m equls m, cn explin observtions, provided tht the popultion of users tht cn digg on tht story is close to Gussin. 1 Introduction In recent yers, online socil networks (OSNs) hve experienced n explosive growth. With hundreds of such services vilble nd subscribed user bse of severl hundred million people, they hve significntly ltered the wy how people spend their time [1] nd how they serch for content [2]. One specific kind of OSNs ddressed in this pper re socil bookmrking services nd, in prticulr, socil news ggregtors such s digg.com, delicious.com or reddit.com. In this type of OSN, users re shring nd commenting on informtion (such s bookmrks, opinions, news, etc.), further clled stories. The community will vote (n ctivity referred to s to digg ) onthesubmittedstories,where the sum of ll votes on story, which is clled the digg vlue of story, is publicly displyed s some rnking informtion nd therefore reflects the impct of story. The life-time of story is lso studied nd used in [3] to predict the populrity of online content. Collbortive tgging in socil medi, s in digg.com, hs been overviewed in [4] from sttisticl physics point of view. Empiricl dt (see e.g. [5,6]) illustrte tht the digg vlue of n rbitrry story follows lognorml distribution s shown in Figure 1. For over hundred yers, lognorml distributions hve been observed in mny different res [7 9], from economy to biology nd now in OSNs. Recent work [10] further demonstrted tht lognorml distribution my further universlly chrcterize dtsets which hve previously been thought of s typicl instnces of power-lw. The fscinting underlying process tht symptoticlly genertes lognorml distribution is the lw of proportionte effect, which is briefly reviewed in Appendix A. Erlier, Wu nd Hubermn [11] hve rgued tht the digg vlue is, indeed, generted by the lw of e-mil: p.f..vnmieghem@tudelft.nl proportionte effect when geing effects re tken into ccount. However, Section 6 questions their rguments. In this rticle, we propose nother process in Section 3 tht leds to lognorml distribution: our counting process voids symptotic limits, but in return requires tht the user popultion is pproximtely Gussin. The new insight is the probbilistic reltion (12) in Section 3.2 of whether user will digg on story, given tht he observes the digg vlue. The reminder of the pper tries to relte both governing processes with experiments. The mechnics of Digg nd the dt extrction method from the Digg OSN, observed over period of five yers, re described in [6,12]. Section 7 concludes our explortion: the lw of proportionte effect lone seems insufficient to explin the experimentl evidence, while our counting process is ble to explin most findings, provided the user popultion is normlly distributed. The key quntity is the proportionlity fctor. Experimentsexhibitroughly linerly decresing with time. The verge of the proportionlity fctor in the lw of proportionte effect (Appendix A) is, however, constnt. Our finite counting process nd the hypothesis (12) gree with experiments: decreses linerly in the number of users tht could digg on story nd pproximtely linerly with time. Our results my hve broder pplicbility, s [13] for exmple lso find log-norml behvior in the populrity of movies t vrious levels, nd hypothesize selfreinforcing process from explicit socil recommendtions s potentil cuse. 2Thediggvlueofstory Let D s denote the digg vlue, the number of diggs, on story s. Ifwessumethtusercnonlydiggonceon story, then N s D s = 1 {user j diggs on s} (1) j=1

3 252 The Europen Physicl Journl B Fig. 1. (Color online) The probbility density function (pdf) of the digg vlue X d of stories on the first five frontpges in Digg. Ech pdf is fitted by lognorml (2) nd the corresponding prmeters µ nd σ re shown in the legend. where N s is the totl number of users tht hve the opportunity to digg on story s. Theindictorfunction 1 x equls one if the event x is true, else it is zero. 2.1 Experiments Empiricl dt illustrted in Figure 1 show tht only the digg vlue of stories on the first frontpge is well fitted by lognormldensityfunction[14],p.57 f lognorml (u) = 1 uσ 2π exp [ (log u µ)2 2σ 2 ]. (2) The probbility density function (pdf) of the other frontpges decys in the left- nd right til fster thn lognorml distribution. The decrese fster thn lognorml is likely due to the geing of the story (explined in Sect. 2.2), the decresing number of users tht visit subsequentil frontpges nd the convolution effect (tht forces ny distribution towrds Gussin distribution). The ltter is consequence of the definition D s = D s0 + D s1 + D s2 +...,whered sk is the sum of the diggs on story s while on frontpge k =1, 2,...nd D s0 is the digg vlue of story s just before it is promoted from the upcoming section to the first frontpge. Although D sk nd D sm re not independent, we infer tht they re only wekly dependent. 2.2 Users digg dependently It is conceivble tht the probbility tht the j-th user, further clled user j, diggs on story s is dependent on the digg vlue tht he observes, when his eye ctches the story s. Thediggvluethtuserj sees when first encountering the story s is, for j>1, the sum of ll digg vlues of users tht dugg on the story before him/her X s (j) = 1 {user k diggs on s}. (3) The eventul digg vlue of story s is D s = X s (N s +1), while X s (1) = 0 nd E [X s (2)] = p 1,theprobbilitytht the first user diggs on story s. Wemyrguethtthe digging probbility of user is influenced by the story s digg vlue: if the number of diggs is low, user hs n priorifeelingthtthestoryisnotsottrctivendhis motivtion to digg is lowered, while the opposite occurs when high number of diggs is observed. Another phenomenon is geing: fter some time, the novelty of the story my diminish, especilly if the story contins news or temporry informtion. In tht cse, the overll motivtion to digg on the story, independent of the digg vlue, decreses. Wu nd Hubermn[11]explintht the story s growing ttrction, mesured vi its incresing digg vlue, is counterblnced by the geing of the story. As discussed in Section 6, the nlysis nd ssumptions of geing effects in [11] re debtble. Moreover, the user popultion tht visits the k-th frontpge for k>4isdecresing so quickly tht dt is scrce.therefore,inthe sequel, geing is not considered, nor time-dependent effects tht ffect the digg vlue. We lso ignore primcy nd recency effects [15] tht would predict tht the top most nd bottom most elements would ttrct lrger ttention, simply becuse of their position in long list. For exmple, the Digg feture, Top Stories in All Topics,

4 P. Vn Mieghem et l.: Lognorml distribution in the digg online socil network 253 displys 10 stories in n bbrevited formt on ech frontpge, so tht the list cn be seen in one scn by ny visitor (15 items on VGA screen require one scrolling ction). Due to their specil plcement lone (irrespective of their content), Top Stories re likely to receive higher thn norml number of diggs. In the sequel, we do not distinguish between stories. Our min interest is to figure out whether the motivtion to digg on story is linerly dependent on tht story s diggvlue. If liner dependence holds, lognorml distribution is the nturl consequence by the lw of proportionte effect (Appendix A). In order to keep externl fctors s constnt s possible nd guided by the experiments in Figure 1, we minly concentrte on the number of digg vlues of story s long s it is on the first frontpge. 3Agenerldescriptionforthevergedigg vlue We tke dvntge of the fct tht the digg vlue, both X s (j) s D s, is sum of indictors: Bernoulli rndom vribles tht re either zero or one. Since E [ 1 {user k diggs on s} ] = Pr[user k diggs on s] nd since the expecttion opertor is liner (irrespective of dependencies between rndom vribles in the sum), tking the expecttion of (3) yields E [X s (j)] = Pr [user k diggs on s]. By the lw of totl probbility [14], we hve tht m=0 Pr [user k diggs on s] = Pr [user k diggs on s X s (k) =m]pr[x s (k) =m]. When writing the conditionl probbility s Pr [ user k diggs on s X s (k) =m] =g k (m) (4) where g k (x) isnon-negtivefunctionthtmpsx to the intervl [0, 1], we obtin Pr [user k diggs on s] = m=0 g k (m)pr[x s (k) =m] = E [g k (X s (k))] (5) where we hve used the definition of the expecttion of function of rndom vrible [14], p. 17. Thus, we rrive t the generl reltion for the verge number of diggs tht the jth user sees E [X s (j)] = E [g k (X s (k))] (6) from which we obtin the difference eqution E [X s (j)] E [X s (j 1)] = E [g j 1 (X s (j 1))]. (7) The verge totl number D s = X s (N s +1) of diggs on story s by popultion of N s users is N s E [D s ]= E [g k (X s (k))]. We illustrte in Appendix B tht, for positively correlted users s in Digg, the vrince Vr[X s (j)] cn be smll. In fct, the stronger the correltion in the digging behvior between users, the smller the vrince. A smll vrince Vr[X s (j)] implies tht the men E [X s (j)] is good pproximtion for the rndom vrible X s (j). So fr, we hve implicitly ssumed tht N s is constnt. However, N s is, in fct, lso rndom vrible, denoting the number of users tht hs discovered story s within certin time intervl. If N s is rndom vrible, then the bove computtion is vlid for the conditionl expecttion Y s = E [D s N s ], which is the rndom vrible equl to the verge number of diggs on story s given tht the totl number of users equls N s. We believe tht the r.v. Y s is pproximtely mesured. Ech story s hs number D s of diggs tht re recorded, while the totl popultion of potentil diggers is N s.hence,wehvesequenceofstorieswiththeircorresponding diggs nd popultion {(D s, N s )} s 1.Clerly, it follows from the definition (1) tht D s N s.ifthe informtion bout the popultion is omitted, intuitively one feels tht D 1 cnnot be compred to D 2,becuseitis obvious tht if n 2 = n1 2,thenthenumberofdiggsD 2 cnnot be higher thn D1 2.Yet,lldiggvluesofthedifferent stories, the sequence {D s } s 1,replcedin1histogrm. In fct, by doing so, we compre different processes, while histogrmshouldonlyrecordtheoutcomesofsme stochstic process nd ech outcome should be independent of ll others. Therefore, we believe tht the histogrm pproximtes the conditionl rndom vrible D s N s by its best guess, the estimted vlue E [D s N s ]=Y s. The reminder consists of choosing the function g k (x), defined in (4). 3.1 A simple cse: proportionlity A simple mthemticl choice is the liner function g k (x) = k x + b k nd the generl difference (7) becomes E [X s (j)] = (1 + j 1 ) E [X s (j 1)] + b j 1. After m itertions of this difference eqution, we obtin E [X s (j)] = E [X s (j m)] + m l=1 m (1 + j k ) l 1 b j l (1 + j k ) (8)

5 254 The Europen Physicl Journl B nd for j = N s +1sinceD s = X s (N s +1), E [D s ]=E [X s (N s +1 m)] + N s l=n s+1 m b l N s k=l+1 N s k=n s+1 m (1 + k ). (1 + k ) With the initil condition E [X s (2)] = p 1,wehve N s E [D s ]=p 1 k=2 N s (1 + k )+ l=2 b l N s k=l+1 (1 + k ). (9) When choosing ll k = nd ll b k = b, thenexpression (8) simplifies to E [X s (j)] = E [X s (j m)] (1 + ) m + b = such tht, for j = N s +1, E [D s ]= m (1 + ) l 1 l=1 ( E [X s (j m)] + b ) (1 + ) m b (10) ( E [X s (N s +1 m)] + b ) (1 + ) m b. (11) Hence, the verge digg vlue E [D s ] = E [D s1 ] of story just before it disppers from the frontpge cn be expressed s function of the verge digg vlue E [X s (N s +1 m)] = E [D s0 ]justfterthestoryppers on the frontpge 1.Duringthetimeonthefrontpge, precisely m users hve hd the opportunity to digg on tht story. Suppose tht the number of users m in (11) is Gussin N ( µ, σ 2), then the rndom vrible Y s1 = E [D s1 m] islognorml,providede [X s (N s +1 m)] = E [D s0 ]isknownconstnt.indeed, [( Pr [Y s1 x] =Pr E [D s0 ]+ b ) (1 + ) m b ] x log =Pr m = 1 σ 2π x+ b E[D s0]+ b log (1 + ) log x+ b E[D s0 ]+ b log(1+) ] (t µ)2 exp [ 2 σ 2 dt. 1 Similrly, we cn express the verge number on the k-th frontpge by considering E [D s] = E [D sk ], E [X s (N s +1 m)] = E [D s, ]ndm is the number of users tht hd the opportunity to digg on the story during its sty on the k-th frontpge. such tht f Ys1 (x) = d Pr [Y s1 x] dx 1 = ( x + b ) σ log (1 + ) 2π [ exp (log ( ( x+ ) b log E [Ds0 ]+ b 2 σ 2 log 2 (1 + ) ) ] µ log (1+)) 2 which is recognized from (2) s lognorml distribution in u = x + b with prmeters µ =log( E [D s0 ]+ b ) + µ log (1 + ) nd σ = σ log (1 + ). Hence, ssuming Gussin popultion, the digg vlue D sk of the story on frontpge k cn be described by lognorml distribution, possibly with different prmeters nd b when geing is tken into ccount. Clerly, the expressions simplify considerbly if b =0. In the cse b =0,ourderivtionisinlinewiththelwof proportionte effect s shown in Section 4. However, the prmeters of the lognorml distribution re different. 3.2 Interprettion The choice g j (x) = j x + b j implies, s follows from (4), tht Pr [user j diggs on s X s (j) =m] = j m + b j. This generl liner form dependent on user j cn be useful to specify the digging behvior of friends nd non-friends of the origintor of the story s. Incsethtuserj is friend of the origintor of the story s,he/she my be insensitive to the digg vlue X s (j) =m tht he/she observes when first encountering the story, becuse his/her friendship reltion with the origintor outweighs the judgements of other diggers. Hence, j =0ndb j = b 1, constnt vlue tht expresses the fith or depth of the friendship reltion of friends in the origintor. On the other hnd, non-friends hve lmost no fith in the origintor, but in their peers, such tht j = 1ndb j =0.Consequently, the generl liner reltion(9)becomes E [D s ]=p 1 (1 + ) (N non-friends 2) + b (N friends 2) where N friends is the number of friends of the origintor of the story s nd N s = N non-friends + N friends. Assuming tht both popultions of friends nd nonfriends re Gussin-like distributed (with possibly different men nd vrince), this expression shows tht E [D s N friends, N non-friends ]isthesumofgussinnd lognorml,whichisginhevy-tiled.sincen friends is usully smller thn N non-friends,lognormldistribution is expected to dominte [6]. When returning to the simplest liner cse for g k (x) = j x + b j,where j = nd b j =0,then Pr [ user j diggs on s X s (j) =m] =m (12)

6 P. Vn Mieghem et l.: Lognorml distribution in the digg online socil network 255 implying tht >0 is ssumed to be constnt for ech user j. Thus, 1 N s in order tht ll conditionl probbilities, for ech j N s,resmllerthnorequlto1. Thus, (12) illustrtes tht decreses with N s,whichis in line with Figure 4. Furthermore, reltion (5) reduces to Pr [user k diggs on s] =E [X s (k)]. If the user is the k-th user, counted since the story is on the j-th frontpge, then, with (10), we hve tht Pr [user k diggs on s] =E [D s,j 1 ](1+) k which shows n exponentilly incresing digging probbility in k. 4Applictionofthelwofproportionte effect to the digg vlue The lw of proportionte effect, explined in Appendix A, is more generl thn the nlysis in Section 2.2: it pplies to the rndom vrible X s (n) instedofthemen E [X s (n)]. On the other hnd, it is n symptotic result (n ), which implicitly ssumes rpid convergence towrds Gussin in order to observe for finite n lredy the lognorml distribution. In our cse, where the digg vlue X j = X s (j +1) is sum of indictors s defined in (3), we hve tht X s (j +1) X s (j) =1 {user j diggs on s}. According to the lw of proportionte effect X j = (1 + α j ) X j 1,wededucetht 1 {user j diggs on s} = α j X s (j). (13) Since α j nd X s (j) reindependentsssumedinthe lw of proportionte effect, tking the expecttion leds to Pr [user j diggs on s] =E [α j ] E [X s (j)] nd we conclude from (5) with g k (x) = j x tht E [α j ]= j.thelwofproportionteeffectssumesthtllrndom vribles α j re i.i.d. with men E [α], such tht E [α j ]=E[α] nd j = for ech j. Hence,theconditionl probbility (12) is mnifesttion of the lw of proportionte effect nd provides nother wy to verify proportionl behvior. We will now compute the prmeter µ = E [log (1 + α)] in the lognorml limit lw for X n = X s (n +1) for lrge n (see Appendix A). Tking the m-th power of (13) yields 1 {user j diggs on s} = α m j X m s (j) from which it follows tht Pr [user j diggs on s] =E [ α m ] j E [X m s (j)] = E [α m ] E [X m s (j)]. Since E [X s (j)] > 1fornottoosmllj, wemyconclude tht E [α m ] < 1forllm 1. Hence, the series [ ] E [log (1 + α)] = E ( 1) m 1 α m m = m=1 m=1 ( 1) m 1 E [α m ] m =Pr[userj diggs on s] m=1 ( 1) m 1 me [X m s converges nd the prmeter µ = E [log (1 + α)] equls ( 1) m 1 µ =Pr[userj diggs on s] me [Xs m (j)]. m=1 (j)] The ltter lternting series with decresing terms is bounded by 1 E [X s (j)] 1 2E [Xs 2 (j)] < Invoking the bounds yields nd Pr [user j diggs on s] E [X s (j)] m=1 ( 1) m 1 me [Xs m (j)] < 1 E [X s (j)]. ( 1 E [X ) s (j)] 2E [Xs 2 (j)] < <µ Pr [user j diggs on s] E [X s (j)] µe [X s (j)] < Pr [user j diggs on s] < E [X s (j)] µ ( 1 E[Xs(j)] 2E[X 2 s (j)] ). Hence, we find tht µe [X s (j)] < 1. However, experimentl results (Fig. 1) indicte tht both µ 5.2 nd E [X s (j)] re lrger thn 1, contrdicting µe [X s (j)] < 1. The nlysis indictes tht the lw of proportionte effect only sets in when certin vlue of X s (j), or equivlently j, isreched. This conclusion is supported by the nlysis in Section 3.1 (nd Fig. 3 below). Using (11), ssuming tht E [X s (N s +1 m)] = D s0,thediggvlueofstorys just before it is promoted to the frontpge, is known constnt nd m is normlly distributed s N ( µ, σ 2),we obtin lognorml with men µ =log ( D s0 + b ) + µ log (1 + ) nd 1 m.confiningtothecsewhereb =0ndusing the verge µ s good estimte for m such tht = 1 µ, then ( µ log (1 + ) µ log 1+ 1 µ ) ( 1 = µ µ 1 ( )) 1 2 µ 2 + O µ 3 =1 1 2 µ + O ( 1 µ 2 ) < 1

7 256 The Europen Physicl Journl B Fig. 2. (Color online) The time t of story s on which user hs dug on tht story. Five typicl stories re shown. The horizontl lines indicte the durtion when the story ws on the frontpge nd the time t seems to be pproximtely liner with the the number of digging users while the story is on the frontpge. A story stys, on verge, 2.4 h on the first frontpge, on which most Digg users re ctive. nd µ log (D s0 ) 1. Since µ 5.2 infigure1,wededucethtd s0 e 4.2, which grees with experiments [6]. 5Theincrementlincreseofthediggvlue with time If the lw of proportionte effect is correct, then it follows from (A.2) tht X n X m = X m n (1 + α j ) 1. j=m+1 For not too smll m, (13)showsthtα j 1 X j < 1. Multiplying out nd neglecting terms with products yields X n X m X m n j=m+1 α j X m (n m) E [α] nd we rrive, with = E [α] sshowninsection4,t X n X m (n m) X m. Figure 2 shows the ge t of story s function of the number of users tht dugg on the story s. Inotherwords, ech dt point reflects the time t which user diggs on the story s. Theremrkbleobservtionisthtthe number of digging users on stories, when stories re on the frontpge (time intervl between two horizontl lines in Fig. 2), seem to be liner in time t. Figure2illustrtesfive stories with different eventul digg vlue. Hence, Figure 2 suggests tht the user s increment n m for stories on the frontpge is proportionl with time t, i.e.the n-th user nd his ppernce re relted s n = βt n + δ. In tht cse (nd ssuming tht lso non-digging users rrive ccording to the liner lw), we hve pproximtely X n X m X m β (t n t m ). (14) Figure 3 shows four five minute intervls in which the increment f (x) =X n X m is drwn versus the digg vlue x = X m.whenthetimeintervlisshort,hrdlynyliner correltion s suggested by (14) is observed. When the time intervl is longer, for exmple 1 h insted of 5 min, Figure 4 strts revelingthelwin(14): the increments re proportionl to the digg vlue t the beginning of the time intervl. Hence, Figure 4 supports the clim tht the digg vlues of story re dependent nd obey the lw of proportionte effect, when the digg vlue is sufficiently high. The lter condition is mthemticlly not precisely determined, becuse the lw of proportionte

8 P. Vn Mieghem et l.: Lognorml distribution in the digg online socil network 257 Fig. 3. (Color online) The increments during 5 min of stories on the frontpge versus their diggvlue t the beginning of the 5 min intervl. Four 5 min intervls re shown nd the liner fit of the increments f (x) versusthediggvluex. Noticethtthe liner fit y = x on log-log scle is line with slope 1 (or 45 degrees) nd t x =1equltolog. Fig. 4. (Color online) The increments of story s versus its digg vlue x during 1 h. The only difference with Figure 3 is the durtion of the intervl (1 h versus 5 min).

9 258 The Europen Physicl Journl B Fig. 5. (Color online) The slope, fitted from the increments during one hour versus the digg vlue, versus the sequentil time intervls of 1 h. The fit of the first hours is lso included. effect is n symptotic lw (see Appendix A), yet observble in relity, i.e. when the digg vlues re finite. Figure 4 indictes tht the tendency towrds the symptotic regime is rther slow. In view of the rther slow tendency towrds the symptotic regime, it is surprising tht the distribution Pr [X d = k] infigure1isclosetolognorml.these observtions question tht the underlying process responsible for the nice lognorml on the front pge is generted by the lw of proportionte effect, becuse the time scles do not mtch. Rther, they suggest tht the user popultion must be Gussin-like distributed (s needed for finite digg vlues in Sect. 3.1) nd tht collective dependence, in proportionl fshion ccording to the hypothesis (12) constitutes the generting process. Figure 5 shows the fits β (t n t m ) in (14), where t n t m =1hoursuchthtβ ct γ n.additionllypproximting γ 1ndusingn = βt n + δ, wededuce tht c c βt n n nd this result (with c 0.3 < 1) grees with 1 N s derived in Section 3.2, but disgrees with the lw of proportionte effect. For longer time intervls (exceeding 17 h), we observe in Figure 5 devitions in the fit of the slopes. For these long times, other effects (minly due to geing, the role of the user interfce nd the principles of humn ttention) strt dominting. 6AnlysisofWundHubermnin[11] After n initilly fst increse, the digg vlue of story flttens with time becuse the novelty of the story hs pssed. Wu nd Hubermn [11] hve tken the geing of story into ccount. The effect of geing mens tht the set {α k } k 1 in (A.1) strts decresing fter some time threshold k c. Bsed on fitting experiments, Wu nd Hubermn [11] propose α k = r k b k, where r k = exp ( 0.4k 0.4) is the geing fctor nd the set {b k } k 1 of rndom vribles is i.i.d. nd with finite men nd vrince. The Lindeberg conditions in [16] for the CLT stte tht σk 2 =Vr[α k]shouldbesmllcompred to N s j=1 σ2 k nd the ltter sum should tend to infinity when N s. The decrese in α k is so fst, tht n lim n j=1 σ2 k is finite, in which cse the limiting distribution is not Gussin nd, consequently, lognorml cnnot be explined from the lw of proportionte effect! In tht cse, the function g k (x) in(4)cnnotbeliner function. Nevertheless, Wu nd Hubermn [11] clim convergence to Gussin (lognorml) by referring to Embrechts nd Mejim [17], Theorem 2, who show tht Z = j=1 c j (λ) X j converges to Gussin with rte of the order of O ( λ αβ),whereα>0ndwhere{x j } j 0 is setofi.i.drndomvribleswithmen0ndvrince1, nd where c j (λ) =O ( λ β) with β>0. Yet, it is not cler whether this theorem pplies to demonstrteconvergence to Gussin (lognorml), becuse r j =exp ( 0.4j 0.4) is compred to c j (λ) nd the r j s depend on j, while the convergence in Embrechts nd Mejim ssumes tht c j (λ) =O ( λ β),forllj. While we hve shown in Section 5 tht the lw of proportionte effect cnnot explin the experiments, the rguments of Wu nd Hubermn t lest lck rigor, s sketched bove. Aprt from mthemticl rigor, since geing does not ply dominnt role on the first frontpge,

10 P. Vn Mieghem et l.: Lognorml distribution in the digg online socil network 259 their pproch is essentilly equl to the lw of proportionte effect, which cnnot explin the nice lognorml on the first frontpge. 7Conclusion We hve investigted two different nlyses tht led to lognorml distribution, observed (see Fig. 1) for the digg vlue of story, while on the first few frontpges. Usully, the lognorml distribution is the chrcteristic fingerprint of the lw of proportionte effect. However, we found tht the lw of proportionte effect only seems to hold when certindiggvlueofthestoryisreched,ndnotfrom the beginning of the digg counting. Moreover, the proportionlity fctor is not constnt s required by the lw of proportionte effect, but dependent on the number of users s shown in our nlysisofsection3ndexperimentilly verified in Section 5. Finlly, the governing difference eqution (A.1) of the lw of proportionte effect is not experimentlly observed t times the story is on the first frontpge. Hence, the lw of proportionte effect cnnot explin the fst convergence towrds the lognorml distribution of digg vlues on the first frontpge. The second nlysis just sums dependent Bernoulli rndom vribles. We show tht the dependence mong users is reflected by (12), Pr [user j diggs on s X s (j) =m] =m. This conditionl probbility (12) is mnifesttion of proportionte effect nd provides notherwytoverifyproportionl behvior. The conditionl probbility (12) illustrtes how individul humn behvior is ffected by tht of others, given tht the individul cn observe how the others rect, e.g. vi the digg vlue X s (j) = m. Also, the bove conditionl probbility shows tht 1 m for ech user m. Atlst,ournlysisofSection3demonstrtes tht lognorml distribution of the digg vlue D s is obtined when the popultion N s of potentil diggers is normlly distributed. Gussin user popultions re pproximtely mesured 2. At lst, our nlysis underlines the importnce of the user popultion specifics, which re, by its symptotic nture (N s ), not relevnt in the lw of proportionte effect. We re grteful to Siyu Tng for her useful comments. Nobert Blenn is funded by TRANS ( Appendix A: The lw of proportionte effect As mentioned in [7] nd in [18], Kpteyn considered in 1903 the eqution X j X j 1 = α j f (X j 1 ) 2 See e.g. Alex the Web Informtion Compny, AlexInternet,Inc.,2010. where the set {α j } 1 j n of rndom vribles is mutully independent nd identiclly distributed, equl to the distribution of the rndom vrible α with men E [α] nd vrince Vr[α]. Moreover, the set {α j } 1 j n of rndom vribles is lso independent of the rndom vribles X 1,X 2,...,X n.thespecilcsewheref(x) =x reduces to X j =(1+α j ) X j 1 (A.1) nd the process tht determines the sequence X 1, X 2,...,X n, given X 0, is sid to obey the lw of proportionte effect, which ws first introduced by Gibrt [19]. After iterting the eqution (A.1), we obtin n X n = X 0 (1 + α j ). (A.2) j=1 By the Centrl Limit Theorem [14] nd ssuming tht ny α j > 1, the sum S n = n j=1 log (1 + α j)ofthe i.i.d. rndom vribles {log (1 + α j )} j 1,echwithdistribution identicl to tht of log (1 + α) with(finite)men E [log (1 + α)] = µ nd vrince σ 2 =Vr[log(1+α)], converges to S n nµ σ d N (0, 1) n d which implies tht S n = log( n j=1 (1 + α j)) N ( nµ, nσ 2).Equivlently,e Sn = n j=1 (1 + α j)tends,for lrge n, tolognormldistributionwithprmetersnµ nd nσ 2.Hence,wehveshowntht,forlrgen, X n is symptoticlly lognormlly distributed with prmeters nµ nd nσ 2,thtrelinerinn. There is continuous vrint of the lw of proportionte effect. In biology, the growth in the number n (t) of items of sme species over time t cn be modelled by the following first order differentil eqution dn (t) dt = r (t) n (t). which reltes the growth (chnge in the popultion) s proportionl to the popultion n (t) ndtheproportion- lity fctor r (t) istimedependent.thegenerlsolution is, for t>, log n (t) =logn ()+ t r (u) du. When we dditionlly ssume tht r (t) chnges t some times = t 0 <t 1 <t 2 <...<t m = t, wheret j re rndom time moments, then where R j = tj t r (u) du = m j=1 tj t j 1 r (u) du = m j=1 t j 1 r (u) du = r (ξ j )(t j t j 1 ) nd ξ j [t j 1,t j ] R j

11 260 The Europen Physicl Journl B ρ = {Pr [ user k diggs on s user l diggs on s] Pr [user k diggs on s]} Pr [user l diggs on s] Pr [user k diggs on s](1 Pr [user k diggs on s]) Pr [user l diggs on s](1 Pr [user l diggs on s]). is rndom vrible with men µ j = E [r (ξ j )(t j t j 1 )]. Assuming tht the Centrl Limit Theorem cn be pplied, the set of rndom vribles {R j } 1 j m tends to GussinN ( mµ, mσ 2),ndthelognormldistribution of n (t) forlrget then follows in the usul wy. Agin, the men is liner in the time t becuse t m = E [ t], the verge time-spcing. Appendix B: The vrince of X s (j) In order to compute the vrince Vr[X s (j)] = E [ Xs 2 (j)] (E [X s (j)]) 2,wefirstrewriteXs 2 (j) using the definition (3) s Xs 2 (j) = 1 {user k diggs on s} 1 {user l diggs on s} nd l=1 = X 2 s (j) =X s (j) 1 {user k diggs on s} +2 1 {user k diggs on s} 1 {user l diggs on s} l= {(user k diggs on s) (user l diggs on s)}. l=1 Tking the expecttion gives E [ X 2 s (j) ] = E [X s (j)] +2 Pr [(user k diggs on s) (user l diggs on s)] l=1 where user l proceeds user k. Invoking the conditionl probbility, finlly results in E [ X 2 s (j) ] = E [X s (j)] +2 Pr [user k diggs on s user l diggs on s] l=1 Pr [user l diggs on s]. (B.1) In generl, we cn upper bound E [ X 2 s (j) ] s E [ Xs 2 (j) ] E [X s (j)] + 2 Pr [user l diggs on s] l=1 = E [X s (j)] + 2 E [X s (k)] (j 1) 2 from which we obtin n upper bound of the vrince Vr [X s (j)] E [X s (j)] + 2 E [X s (k)] (E [X s (j)]) 2. If users were to digg independently so tht Pr [user k diggs on s user l diggs on s] = Pr [user k diggs on s] then the conditionl probbility in (B.1) would reduce to E [ X 2 s (j)] = E [X s (j)] Pr [user k diggs on s] E [X s (k)]. j 1 +2 In Digg, users re positively correlted 3 such tht Pr [user k diggs on s user l diggs on s] Pr [user k diggs on s] which leds to higher second moment E [ X 2 s (j)] in (B.1) thn for independent users. However, the men E [X s (j)] j 1inDiggislsolrger.SinceVr[X s (j)] (j 1) 2 (E [X s (j)]) 2,themendecresesthevrince qudrticlly, which my led toreltivelysmllvrince. In other words, the stronger the proportionte effect (mesured vi the proportionlity fctor in (12)), the smller the vrince Vr[X s (j)] nd the better the men E [X s (j)] pproximtes the rndom vrible X s (j). References 1. Rob McGnn, Internet Edges Out Fmily Time More Thn TV Time (ClickZ.com, 2005) 2. C. Nuttll nd D. Gelles, Fcebook becomes bigger hit thn Google (Finncil Times, 2010) 3 Applying the definition [14], p. 30 of the liner correltion coefficient ρ yields see eqution bove.

12 P. Vn Mieghem et l.: Lognorml distribution in the digg online socil network G. Szbo, B.A. Hubermn, Commun. ACM 53, 80 (2010) 4. C. Cstellno, S. Fortunto, V. Loreto, Revi. Mod. Phys. 81, 591(2009) 5. K. Lermn, A. Glstyn, Proceedings of the First Workshop on online Socil Networks (WOSP 08) (ACM, New York, 2008), pp S. Tng, N. Blenn, C. Doerr, P. Vn Mieghem, to pper in IEEE Trnsctions on Multimedi 7. E.L. Crow, K. Shimizu, Lognorml distributions, Theory nd Applictions (Mrcel Dekker, Inc. New York, 1988) 8. E. Limpert, W.A. Sthel, M. Abbt, Bioscience 51, 341 (2001) 9. W. Shockley, Proceedings of the IRE 45, 279(1957) 10. F. Rdicchi, S. Fortunto, C. Cstellno, Proc. Ntl. Acd. Sci. USA 105, 17268(2008) 11. F. Wu, B.A. Hubermn, Proc. Ntl. Acd. Sci. USA 104, (2007) 12. C. Doerr, S. Tng, N. Blenn, P. Vn Mieghem, Are Friends Overrted? A Study of the Socil News Aggregtor Digg.com (IFIP Networking, 2011) 13. R.K. Pn, S. Sinh, New J. Phys. 12, (2010) 14. P. Vn Mieghem. Performnce Anlysis of Communictions Systems nd Networks (Cmbridge University Press, Cmbridge, U.K., 2006) 15. Deese, Kufmn, J. Exp. Psychol. 54, 180(1957) 16. W. Feller, An Introduction to Probbility Theory nd Its Applictions, 2ndedn.(JohnWiley&Sons,NewYork, 1971), Vol P. Embrechts, M. Mejim, Zeitung für Whrscheinlichkeitstheorie und verwndte Gebiete 68, 191(1984) 18. M. Armte, Mthémtiques et sciences humines 129, 5 (1995) 19. R. Gibrt, Bulletin de l Sttistique Générle de l Frnce 19, 469(1930)

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

1 Probability Density Functions

1 Probability Density Functions Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our

More information

Chapter 5 : Continuous Random Variables

Chapter 5 : Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

CS667 Lecture 6: Monte Carlo Integration 02/10/05

CS667 Lecture 6: Monte Carlo Integration 02/10/05 CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of

More information

Lecture 3 Gaussian Probability Distribution

Lecture 3 Gaussian Probability Distribution Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning

Solution for Assignment 1 : Intro to Probability and Statistics, PAC learning Solution for Assignment 1 : Intro to Probbility nd Sttistics, PAC lerning 10-701/15-781: Mchine Lerning (Fll 004) Due: Sept. 30th 004, Thursdy, Strt of clss Question 1. Bsic Probbility ( 18 pts) 1.1 (

More information

Entropy and Ergodic Theory Notes 10: Large Deviations I

Entropy and Ergodic Theory Notes 10: Large Deviations I Entropy nd Ergodic Theory Notes 10: Lrge Devitions I 1 A chnge of convention This is our first lecture on pplictions of entropy in probbility theory. In probbility theory, the convention is tht ll logrithms

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

ECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance

ECO 317 Economics of Uncertainty Fall Term 2007 Notes for lectures 4. Stochastic Dominance Generl structure ECO 37 Economics of Uncertinty Fll Term 007 Notes for lectures 4. Stochstic Dominnce Here we suppose tht the consequences re welth mounts denoted by W, which cn tke on ny vlue between

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

3.4 Numerical integration

3.4 Numerical integration 3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Tests for the Ratio of Two Poisson Rates

Tests for the Ratio of Two Poisson Rates Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17 EECS 70 Discrete Mthemtics nd Proility Theory Spring 2013 Annt Shi Lecture 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion,

More information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

The steps of the hypothesis test

The steps of the hypothesis test ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of

More information

1B40 Practical Skills

1B40 Practical Skills B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need

More information

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

Math 31S. Rumbos Fall Solutions to Assignment #16

Math 31S. Rumbos Fall Solutions to Assignment #16 Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0. STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

Part I: Basic Concepts of Thermodynamics

Part I: Basic Concepts of Thermodynamics Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

Riemann is the Mann! (But Lebesgue may besgue to differ.)

Riemann is the Mann! (But Lebesgue may besgue to differ.) Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >

More information

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17 CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

Operations with Polynomials

Operations with Polynomials 38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils

More information

du = C dy = 1 dy = dy W is invertible with inverse U, so that y = W(t) is exactly the same thing as t = U(y),

du = C dy = 1 dy = dy W is invertible with inverse U, so that y = W(t) is exactly the same thing as t = U(y), 29. Differentil equtions. The conceptul bsis of llometr Did it occur to ou in Lecture 3 wh Fiboncci would even cre how rpidl rbbit popultion grows? Mbe he wnted to et the rbbits. If so, then he would be

More information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Heat flux and total heat

Heat flux and total heat Het flux nd totl het John McCun Mrch 14, 2017 1 Introduction Yesterdy (if I remember correctly) Ms. Prsd sked me question bout the condition of insulted boundry for the 1D het eqution, nd (bsed on glnce

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model: 1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

and that at t = 0 the object is at position 5. Find the position of the object at t = 2. 7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we

More information

A SHORT NOTE ON THE MONOTONICITY OF THE ERLANG C FORMULA IN THE HALFIN-WHITT REGIME. Bernardo D Auria 1

A SHORT NOTE ON THE MONOTONICITY OF THE ERLANG C FORMULA IN THE HALFIN-WHITT REGIME. Bernardo D Auria 1 Working Pper 11-42 (31) Sttistics nd Econometrics Series December, 2011 Deprtmento de Estdístic Universidd Crlos III de Mdrid Clle Mdrid, 126 28903 Getfe (Spin) Fx (34) 91 624-98-49 A SHORT NOTE ON THE

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

Student Activity 3: Single Factor ANOVA

Student Activity 3: Single Factor ANOVA MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether

More information

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior Reversls of Signl-Posterior Monotonicity for Any Bounded Prior Christopher P. Chmbers Pul J. Hely Abstrct Pul Milgrom (The Bell Journl of Economics, 12(2): 380 391) showed tht if the strict monotone likelihood

More information

Math 42 Chapter 7 Practice Problems Set B

Math 42 Chapter 7 Practice Problems Set B Mth 42 Chpter 7 Prctice Problems Set B 1. Which of the following functions is solution of the differentil eqution dy dx = 4xy? () y = e 4x (c) y = e 2x2 (e) y = e 2x (g) y = 4e2x2 (b) y = 4x (d) y = 4x

More information

MATH SS124 Sec 39 Concepts summary with examples

MATH SS124 Sec 39 Concepts summary with examples This note is mde for students in MTH124 Section 39 to review most(not ll) topics I think we covered in this semester, nd there s exmples fter these concepts, go over this note nd try to solve those exmples

More information

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve. Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F

More information

Math 426: Probability Final Exam Practice

Math 426: Probability Final Exam Practice Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by

More information

221A Lecture Notes WKB Method

221A Lecture Notes WKB Method A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e

More information

Non-Linear & Logistic Regression

Non-Linear & Logistic Regression Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

Module 6: LINEAR TRANSFORMATIONS

Module 6: LINEAR TRANSFORMATIONS Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for

More information

2008 Mathematical Methods (CAS) GA 3: Examination 2

2008 Mathematical Methods (CAS) GA 3: Examination 2 Mthemticl Methods (CAS) GA : Exmintion GENERAL COMMENTS There were 406 students who st the Mthemticl Methods (CAS) exmintion in. Mrks rnged from to 79 out of possible score of 80. Student responses showed

More information

IN GAUSSIAN INTEGERS X 3 + Y 3 = Z 3 HAS ONLY TRIVIAL SOLUTIONS A NEW APPROACH

IN GAUSSIAN INTEGERS X 3 + Y 3 = Z 3 HAS ONLY TRIVIAL SOLUTIONS A NEW APPROACH INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 8 (2008), #A2 IN GAUSSIAN INTEGERS X + Y = Z HAS ONLY TRIVIAL SOLUTIONS A NEW APPROACH Elis Lmpkis Lmpropoulou (Term), Kiprissi, T.K: 24500,

More information

A PREY-PREDATOR MODEL WITH COVER FOR THE PREY AND AN ALTERNATIVE FOOD FOR THE PREDATOR AND CONSTANT HARVESTING OF BOTH THE SPECIES *

A PREY-PREDATOR MODEL WITH COVER FOR THE PREY AND AN ALTERNATIVE FOOD FOR THE PREDATOR AND CONSTANT HARVESTING OF BOTH THE SPECIES * Jordn Journl of Mthemtics nd Sttistics (JJMS) (), 009, pp. 43-54 A PREY-PREATOR MOEL WITH COVER FOR THE PREY A A ALTERATIVE FOO FOR THE PREATOR A COSTAT HARVESTIG OF BOTH THE SPECIES * K. LAKSHMI ARAYA.PATTABHIRAMACHARYULU

More information

A New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation

A New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation Xiodong Zhung A ew Sttistic Feture of the Short-Time Amplitude Spectrum Vlues for Humn s Unvoiced Pronuncition IAODOG ZHUAG 1 1. Qingdo University, Electronics & Informtion College, Qingdo, 6671 CHIA Abstrct:

More information

4.4 Areas, Integrals and Antiderivatives

4.4 Areas, Integrals and Antiderivatives . res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

Shear and torsion interaction of hollow core slabs

Shear and torsion interaction of hollow core slabs Competitive nd Sustinble Growth Contrct Nº G6RD-CT--6 Sher nd torsion interction of hollow core slbs HOLCOTORS Technicl Report, Rev. Anlyses of hollow core floors December The content of the present publiction

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

Name Solutions to Test 3 November 8, 2017

Name Solutions to Test 3 November 8, 2017 Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier

More information

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from

More information

Green s functions. f(t) =

Green s functions. f(t) = Consider the 2nd order liner inhomogeneous ODE Green s functions d 2 u 2 + k(t)du + p(t)u(t) = f(t). Of course, in prctice we ll only del with the two prticulr types of 2nd order ODEs we discussed lst

More information

Math 135, Spring 2012: HW 7

Math 135, Spring 2012: HW 7 Mth 3, Spring : HW 7 Problem (p. 34 #). SOLUTION. Let N the number of risins per cookie. If N is Poisson rndom vrible with prmeter λ, then nd for this to be t lest.99, we need P (N ) P (N ) ep( λ) λ ln(.)

More information

Rel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited

More information

Applicable Analysis and Discrete Mathematics available online at

Applicable Analysis and Discrete Mathematics available online at Applicble Anlysis nd Discrete Mthemtics vilble online t http://pefmth.etf.rs Appl. Anl. Discrete Mth. 4 (2010), 23 31. doi:10.2298/aadm100201012k NUMERICAL ANALYSIS MEETS NUMBER THEORY: USING ROOTFINDING

More information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019 ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil

More information

Generalized Fano and non-fano networks

Generalized Fano and non-fano networks Generlized Fno nd non-fno networks Nildri Ds nd Brijesh Kumr Ri Deprtment of Electronics nd Electricl Engineering Indin Institute of Technology Guwhti, Guwhti, Assm, Indi Emil: {d.nildri, bkri}@iitg.ernet.in

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Probability Distributions for Gradient Directions in Uncertain 3D Scalar Fields

Probability Distributions for Gradient Directions in Uncertain 3D Scalar Fields Technicl Report 7.8. Technische Universität München Probbility Distributions for Grdient Directions in Uncertin 3D Sclr Fields Tobis Pfffelmoser, Mihel Mihi, nd Rüdiger Westermnn Computer Grphics nd Visuliztion

More information

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information