R/fcn.R

In MeasurementError.cor: Measurement Error model estimate for correlation coefficient

"fcn" <-
function(x,exp1,se1,exp2,se2)
{

  r1 <- x[1]
  r2 <- x[2]
  mu1 <- x[3]
  mu2 <- x[4]
  s1 <- x[5]
  s2 <- x[6]
  
  f <- sum(log(2.0)+log(pi)-(exp2^2*se1^2-2.0*se1^2*exp2*mu2+mu2^2*se1^2-2.0*se1*r1*mu1*se2*mu2+2.0*se1*r1*exp2*mu1*se2-2.0*se1*r1*exp1*exp2*se2+2.0*se1*r1*exp1*se2*mu2+2.0*exp1*sqrt(s1*s2)*r2*mu2+2.0*r2*mu1*sqrt(s1*s2)*exp2+exp2^2*s1+exp1^2*se2^2+exp1^2*s2+mu1^2*se2^2+mu1^2*s2-2.0*r2*mu1*mu2*sqrt(s1*s2)-2.0*exp1*r2*sqrt(s1*s2)*exp2+mu2^2*s1-2.0*exp1*s2*mu1-2.0*exp1*se2^2*mu1-2.0*exp2*mu2*s1)/(-se1^2*se2^2-s1*se2^2-se1^2*s2-s1*s2+s1*r2^2*s2+se1^2*r1^2*se2^2+2.0*r1*se1*r2*sqrt(s1*s2)*se2)/2.0+0.5*log(se1^2*se2^2+s1*se2^2+se1^2*s2+s1*s2-s1*r2^2*s2-se1^2*r1^2*se2^2-2.0*r1*se1*r2*sqrt(s1*s2)*se2))

  return(f)
           
}