Revision a300fe2b8c94aa973f39d050e5d9ddf04f609e1a authored by M. Helena Gon\xe7alves on 17 February 2013, 00:00:00 UTC, committed by Gabor Csardi on 17 February 2013, 00:00:00 UTC
1 parent e9f9bbd
f2var.f
C CCCC
C function to calculate the integral for omega for the second order
C dependence model with random effects
C dimension of x and vectors assume the maximum values
CCCCC
DOUBLE PRECISION FUNCTION f2var(v,i)
DOUBLE PRECISION logL,prob,beta1,bt1,lpsi1,x1,
*theta1,work1,omega1,v
INTEGER m, mpar, y1,i
DIMENSION x1(4500,10),theta1(4500),
*work1(4500),y1(4500),prob(4500),lpsi1(2),
*beta1(10),bt1(10)
COMMON/param/x1,theta1,work1,
*y1,lpsi1,beta1,bt1,m,mpar,omega1
beta1(1) = v +bt1(1)
CALL mlik2i(logL,prob,mpar,m)
f2var = dexp(logL-(v**2)/(2*dexp(omega1)))*
*((v**2-dexp(omega1))/(2*dexp(2*omega1)))
return
end
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...