log.vero {geoComp} | R Documentation |
Usage
log.vero(theta1, dados.comp, print.pars = F)
Arguments
theta1 |
|
dados.comp |
|
print.pars |
|
Examples
##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function(theta1,dados.comp,print.pars=F){
if(theta1[1] <= 0) return(.Machine$double.xmax^0.5)
if(theta1[2] <= 0) return(.Machine$double.xmax^0.5)
if(theta1[3] <= 0) return(.Machine$double.xmax^0.5)
if(theta1[4] < 0) return(.Machine$double.xmax^0.5)
if(abs(theta1[5]) > 1) return(.Machine$double.xmax^0.5)
## Pegando o Y
y <- dados.comp[[2]][[1]]
X <- cbind(rep(1:0,length=length(y)),rep(0:1,length=length(y)))
####
V <- monta.V(theta1, dados.comp = dados.comp)
ldetV <- determinant(V,log=TRUE)$modulus[1]
mu <- drop(solve(crossprod(X,solve(V,X))) %*%crossprod(X,solve(V,y)))
Qe <- drop(crossprod(y,solve(V,y))-2*crossprod(y,solve(V,X%*%mu))+crossprod(mu,crossprod(X,solve(V,X%*%mu))))
if(Qe < 0) return(Qe=.Machine$double.xmax^0.5)
n <- length(y)
s1 <- sqrt(Qe/n)
# drop diminui a dimensao do objeto
ll <- drop(-0.5*(n*log(2*pi)+n*log(s1^2)+ ldetV + (1/(s1^2))*Qe))
if(print.pars) print(c(theta1,ll))
return(-ll)
}
[Package
geoComp version 0.1-0
Index]