Displaying 1 result from an estimated 1 matches for "s2eps".
2001 Nov 22
2
optim ???
...165,166,157,147,145,168,180,147,136,151,168,179,129,131)
dantes<-c(130,122,124,104,112,101,121,124,115,102,98,119,106,107,100)
ddespues<-c(125,121,121,106,101,85,98,105,103,98,90,98,110,103,82)
presion<-data.frame( santes, sdespues, dantes, ddespues)
logvero <- function(mu,delta,s2ksi,s2eps,beta,s2eta,x,z)
{
n <- length(x)
sigma <- matrix(c(s2ksi+s2eps,beta*s2ksi-s2eps,beta*s2ksi-s2eps,
beta^2*s2ksi+s2eta+2*s2eps),2,2)
xz <- sum((x-mu)*(z-delta))
A <- matrix(c(sum((x-mu)^2),xz,xz,sum((z-delta)^2)),2,2)
A <- (1/2)*solve(sigma) %*...