Inference of SNP-Gene Regulatory Networks by Integrating Gene Expressions and Genetic Perturbations

Algorithm 1

Optimization for elastic net in Step 1-2.

(1)  procedure ELASTIC(𝑌,𝑋,𝜆₁,𝜆₂,𝑖,𝜀)⊳𝜆₁ and 𝜆₂ are optimal parameters estimated by cross validation

(2)  while  err>𝜀  do

(3)   𝐛^old_𝑖=𝐛_𝑖,𝐟^old_𝑖=𝐟_𝑖

(4)   for  𝑗←1,𝑀_𝑠  do

(5)    Update 𝑓_𝑖𝑗 via (12)

(6)   end for

(7)   Update 𝐛_𝑖 via (5)

(8)   err=𝐛^old_𝑖−𝐛_𝑖2+𝐟^old_𝑖−𝐟_𝑖2

(9)  end while

(10)  return  𝐛_𝑖 and 𝐟_𝑖

(11) end procedure