How can I write the code of the following formula?

y i,t = (1−θ) y i,t−1 + θαi + θβx i,t + θγz t + ε i,t

i have to use panel data, fixed effect model. I tried to use plm package but I'm having some difficulties because I don't know how to indicate the parameter θ and how to insert it in the formula. That's the code I used:

femodel = plm(Y ~ lag(Y) + lag(Y)-1 + X1 + X2 + X3 + Z1 + Z2, index = c("ID","ANNO"), data = pdata, model = "within", effect = "twoways").

Could you help me to insert the other parameter θ?