def dlogpdf_dlink_dtheta(self, f, y, Y_metadata=None):
dlogpdf_dlink_dtheta = np.zeros((self.size, f.shape[<NUM_LIT:0>], f.shape[<NUM_LIT:1>]))<EOL>dlogpdf_dlink_dtheta[<NUM_LIT:0>,:,:] = self.dlogpdf_dlink_dvar(f,y,Y_metadata=Y_metadata)<EOL>return dlogpdf_dlink_dtheta<EOL>
:param link_f: latent variables link(f)
:type link_f: Nx1 array
:param y: data
:type y: Nx1 array
:param Y_metadata: Y_metadata not used in gaussian
:returns: derivative of log likelihood evaluated at points link(f) w.r.t variance parameter
:rtype: Nx1 array