Hi R-Help,

I am using R to do functional outlier detection (using PCA to reduce to 2 
dimensions - the functional boxplot methodology used in the Rainbow package), 
and using Hscv.diag function to calculate the bandwidth matrix where this line 
of code is run:
result <- optim(diag(Hstart), scv.mat.temp, method = "Nelder-Mead", control = 
list(trace = as.numeric(verbose)))

Within the optim function, there is a call to an external C function:
.External2(C_optim, par, fn1, gr1, method, con, lower, upper)
Where
Par = (0.339, 0.339),
fn1 = function (diagH)
{
  H <- diag(diagH) %*% diag(diagH)
  if (default.bflag(d = d, n = n))
    scvm <- scv.mat(x.star, H, Gs, binned = binned, verbose = FALSE,
      bin.par = bin.par.star, deriv.order = r)
  else scvm <- scv.mat(x.star, H, Gs, binned = binned, verbose = FALSE,
    deriv.order = r)
  return(scvm)
},
gr1 is Null,
method = "Nelder-Mead"
con is a list of 18 items (trace=0, fnscale=1, parscale=(1,1), ...),
lower and upper are -Inf and Inf respectively.

When this is run, it returns an error
Error in chol.default(H) :
  the leading minor of order 1 is not positive definite

Has anyone any experience with this error, or know where I can locate the code 
for the C_optim function to continue trying to find the issue and a work around?

Thanks for any help,

Kieran

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to