powered by
bt.scalar function
bt.scalar(t, theta)
change in time, see Brix and Diggle (2001)
parameter beta in Brix and Diggle
...