Internal function used compute a sum in FPCA-based covariance updates
f_sum2(y, fixef, mu.q.c, kt, theta)
outcome matrix
current estimate of fixed effects
current value of mu.q.c
number of basis functions
spline basis