R 里面怎么写出下面的function？

各位大大，小弟以前习惯用stata编程，现在改用R，所以很多不太熟悉，特此请教各位，在如下一个finite mixture model 当中，是两个OLS混合而来的，所以有如下的概率密度函数：



在stata中，可以如下书写：
----------------------- copy starting from the next line -----------------------

1. * Example generated by -dataex-. To install: ssc install dataex
   
2. clear
   
3.    gen double `f1'=normalden($ML_y1,`xb1',exp(`lns1'))
   
4.    gen double `f2'=normalden($ML_y1,`xb2',exp(`lns2'))
   
5.    tempvar p
   
6.    gen double `p'=exp(`lp')/(1+exp(`lp'))
   
7.    replace `f'=`p'*`f1'+(1-`p')*`f2'
   
8. 
   
9. end

复制代码

------------------ copy up to and including the previous line ------------------

求助各位大大怎么把这段函数转成 R 里面的一个function，不胜感激