• Chinese Physics B
  • Vol. 29, Issue 8, (2020)
Hong-Li Zeng1,2,† and Erik Aurell3,4
Author Affiliations
  • 1School of Science, New Energy Technology Engineering Laboratory of Jiangsu Province, Nanjing University of Posts and Telecommunications, Nanjing 20023, China
  • 2Nordita, Royal Institute of Technology, and Stockholm University, SE-10691 Stockholm, Sweden
  • 3KTH – Royal Institute of Technology, AlbaNova University Center, SE-10691 Stockholm, Sweden
  • 4Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, 30-38 Kraków, Poland
  • show less
    DOI: 10.1088/1674-1056/ab8da6 Cite this Article
    Hong-Li Zeng, Erik Aurell. Inverse Ising techniques to infer underlying mechanisms from data[J]. Chinese Physics B, 2020, 29(8): Copy Citation Text show less

    Abstract

    As a problem in data science the inverse Ising (or Potts) problem is to infer the parameters of a Gibbs–Boltzmann distributions of an Ising (or Potts) model from samples drawn from that distribution. The algorithmic and computational interest stems from the fact that this inference task cannot be carried out efficiently by the maximum likelihood criterion, since the normalizing constant of the distribution (the partition function) cannot be calculated exactly and efficiently. The practical interest on the other hand flows from several outstanding applications, of which the most well known has been predicting spatial contacts in protein structures from tables of homologous protein sequences. Most applications to date have been to data that has been produced by a dynamical process which, as far as it is known, cannot be expected to satisfy detailed balance. There is therefore no a priori reason to expect the distribution to be of the Gibbs–Boltzmann type, and no a priori reason to expect that inverse Ising (or Potts) techniques should yield useful information. In this review we discuss two types of problems where progress nevertheless can be made. We find that depending on model parameters there are phases where, in fact, the distribution is close to Gibbs–Boltzmann distribution, a non-equilibrium nature of the under-lying dynamics notwithstanding. We also discuss the relation between inferred Ising model parameters and parameters of the underlying dynamics.
    P(\boldsymbols)=exp(β(iθisi+i<jJijsisj))Z,(1)

    View in Article

    Z=\boldsymbolsexp(β(iθisi+i<jJijsisj)).(2)

    View in Article

    {θi,Jij}ML=argmax[iθisii<jJijsisj1βlogZ],(3)

    View in Article

    P(\boldsymbols)W\boldsymbols,\boldsymbols=P(\boldsymbols)W\boldsymbols,\boldsymbols.(4)

    View in Article

    ri(,\boldsymbolsi)=ri(+,\boldsymbolsi)e2βθi2βjJijsj,(5)

    View in Article

    ri(\boldsymbols)=γi(\boldsymbolsi)eβΔiE(\boldsymbols),(6)

    View in Article

    mi=tanh(hi+jJijmj).(7)

    View in Article

    cij=sisjsisj=mihj(8)

    View in Article

    Jij,nMF=(c1)ij.(9)

    View in Article

    hi,nMF=tanh1mijJij,nMFmj.(10)

    View in Article

    mi=tanh(hi+jJijmjmijJij2(1mj2)).(11)

    View in Article

    Jij,TAP+2mimj(Jij,TAP)2=(c1)ij.(12)

    View in Article

    P(si|\boldsymbolsi)=exp(β(θisi+jJijsisj))s=±1exp(β(θis+jJijssj)).(13)

    View in Article

    (θi,Jij)PML,i=argmax[θisijJijsisj1βlogζi],(14)

    View in Article

    ddtp(s1,,sL;t)=iωi(si)p(s1,,si,,sL;t)iωi(si)p(\boldsymbols;t),(15)

    View in Article

    ωi(\boldsymbols)=γ1+exp[2si(θi+jJijsj)]=γ2[1sitanh(θi+jJijsj)].(16)

    View in Article

    Pe(\boldsymbols,t)=1Ns=1N1\boldsymbols,\boldsymbolss(t)(17)

    View in Article

    Hi(t)=jJijsj(t)+θi,(18)

    View in Article

    ωi(\boldsymbols,t)=γ2[1si(t)tanh(Hi(t))].(19)

    View in Article

    ωTOT(\boldsymbols,t)=iωi(\boldsymbols,t).(20)

    View in Article

    si(t+δt)={+1withprobability1/{1+exp[2βHi(t)]},1withprobability1/{1+exp[2βHi(t)]}.()

    View in Article

    {γ/{1+exp[2βHi(t)]},whensi(t)=1,γ/{1+exp[2βHi(t)]},whensi(t)=1.()

    View in Article

    Jij=Jijs+kJijas,k0,(21)

    View in Article

    Jijs=Jjis,Jijas=Jjias.(22)

    View in Article

    Jijs2=Jijas2=g2N11+k2.(23)

    View in Article

    p(Jij)exp((Jijμ)22σ2)(24)

    View in Article

    mi(t)=si(t),cij(t0+τ,t0)=si(τ+t0)sj(t0)mi(τ+t0)mj(t0).(25)

    View in Article

    dmi(t)dt=mi(t)+tanh[Hi(t)],(26)

    View in Article

    dsi(t)sj(t0)dt=si(t)sj(t0)+tanh[Hi(t)sj(t0)].(27)

    View in Article

    bi=θi+jJijmj(28)

    View in Article

    Hibi+jJijδsj(t),(29)

    View in Article

    dmi(t)dt=mi(t)+tanhbi(t)(lowestorderclosure).(30)

    View in Article

    si(t)sj(t0)+dsi(t)sj(t0)dt=mimj+(1mi2)(jJikδsk(t)δsj(t0)),(31)

    View in Article

    Cij(t,t0)=δsi(t)δsj(t0)=si(t)sj(t0)mimj,(32)

    View in Article

    Dij(t,t0)=Cij(t,t0)+dCij(t,t0)dt,(33)

    View in Article

    Dij(t,t0)=(1mi2)kJikCkj(t,t0).(34)

    View in Article

    J,asynnMF=A1DC1,(35)

    View in Article

    Hi(t)=bimikiJik2(1mk2)+kJikδsk(t),(36)

    View in Article

    J,asynTAP=A1DC1,(37)

    View in Article

    Aij=δij(1mi2)[1(1mi2)jJij2(1mj2)].(38)

    View in Article

    Jt+1=A(Jt)1DC1.(39)

    View in Article

    Fi=(1mi2)jJij2(1mj2)(40)

    View in Article

    JijasynTAP=Vij(1mi2)(1Fi),(41)

    View in Article

    Fi(1Fi)2jVij2(1mj2)1mi2=0.(42)

    View in Article

    JijasynTAP=JijasynnMF1Fi.(43)

    View in Article

    μi1(1)dt1μin(n)dtn×exp(μ(1)t1μ(2)(t2t1)μ(n+1)(ttn)),()

    View in Article

    L=i,tlog[(1γδt)δsi(t+δt),si(t)+γδtesi(t+δt)Hi(t)2coshHi(t)].(44)

    View in Article

    δJijLJij=flips[si(t+δt)tanh(Hi(t))]sj(t)+γδt2noflipsqi(t)si(t+δt)sj(t),(45)

    View in Article

    δJijγ1C˙ij(0)+Cij(0)tanh(Hi(t))sj(t),(46)

    View in Article

    ε=ij(JijJijTrue)2N(N1),(47)

    View in Article

    si(t)={1,ift[tif,min(tif+1,tin+X))withXexp(γ1),1,otherwise,(48)

    View in Article

    δθi=η(siDatasiModel),δJij=η(sisjDatasisjModel).()

    View in Article

    si(t)={1,ifVi(t,Δt)Vthi,1,ifVi(t,Δt)<Vthi.(49)

    View in Article

    Jij=C(0)ij1.()

    View in Article

    Jij=11mi2(dC(τ)dτ|τ=0C(0)1)ij.()

    View in Article

    QJ,J=i,jJijJiji,jmax(Jij,Jij)2.(50)

    View in Article

    F(\boldsymbols)=F0+ifisi+ijfijsi,sj.(51)

    View in Article

    P(\boldsymbols)t|sel=P(\boldsymbols)(F(\boldsymbols)F(\boldsymbols)),(52)

    View in Article

    P(\boldsymbols)t|rec=rξ,\boldsymbolsC(ξ)[Q(\boldsymbols1,\boldsymbols2)P2(\boldsymbols1,\boldsymbols2)Q(\boldsymbols,\boldsymbols)P2(\boldsymbols,\boldsymbols)],(53)

    View in Article

    P2(\boldsymbols1,\boldsymbols2)=P(\boldsymbols1)P(\boldsymbols2).(54)

    View in Article

    PAB(xA,xB)=1Ns=1N1xA(s),xA1xB(s),xB,(55)

    View in Article

    Jij=fijrcij,(56)

    View in Article

    cij=ξC(ξ)(ξi(1ξj)+(1ξi)ξj).(57)

    View in Article

    fij=Jijrcij,(58)

    View in Article

    Hong-Li Zeng, Erik Aurell. Inverse Ising techniques to infer underlying mechanisms from data[J]. Chinese Physics B, 2020, 29(8):
    Download Citation