FYI,
your more likely to get any response to this if you first check if the
issue has been solved in a more recent version of R. Your version is
1.5 years old.
/Henrik
On Thu, Mar 25, 2010 at 2:03 PM, bernardo lagos alvarez
blacert...@gmail.com wrote:
Hello R Community,
I've been run the following codes. However, I've been getting an
unusual segfault
that I'm unable to trace its origin. Please give me a light to
decipher the caught segfault
Thanks for you attention.
Bernardo.
options(STERM='iESS', editor='emacsclient')
rm(list = ls())
source(fgenIGLD.R) #RNG of IGLD(mu,b,sigma)
source(fsYUTH.R) #the fucntions Y,U,Tbbarran and H necessary to
calcualte the hat{b}
source(fMLEPFP.R) #the maximun likelihood profile estimation
of IGLD(mu,sigma)
source(f04P.R) #it to compute the mle, jacknife and bootstrap
source(fprofileloglik06.R) #the log-likelihood profile of IGLD(mu,sigma)
source(fmgfIGLD.R) #read the fmgfIGLD.R file defining the
fmgfIGLD function. It is useful for the prox.
source(f02.R) #here the function kappas.and.deriv(n,
source(fsimbiasIGLDP.R) #it to compute de bias of jacknife and
bootstrap
source(f01P.R) #in this is fmainsimbiasP
ptm-proc.time()
sink(out22031001.txt)
fmainsimbiasP(1,100)
proc.time()-ptm
*** caught segfault ***
address 0x18, cause 'memory not mapped'
Traceback:
1: nlm(f = fprof_deriv, x = x, p = parHInt, b = parFIX, hessian = T,
iterlim = 2000, check.analyticals = F, stepmax = 10)
2: doTryCatch(return(expr), name, parentenv, handler)
3: tryCatchOne(expr, names, parentenv, handlers[[1]])
4: tryCatchList(expr, classes, parentenv, handlers)
5: tryCatch(expr, error = function(e) { call - conditionCall(e)
if (!is.null(call)) { if (identical(call[[1]],
quote(doTryCatch))) call - sys.call(-4) dcall -
deparse(call)[1] prefix - paste(Error in, dcall, : )
LONG - 75 msg - conditionMessage(e) sm -
strsplit(msg, \n)[[1]] if (14 + nchar(dcall, type = w) +
nchar(sm[1], type = w) LONG) prefix -
paste(prefix, \n , sep = ) } else prefix - Error :
msg - paste(prefix, conditionMessage(e), \n, sep = )
.Internal(seterrmessage(msg[1])) if (!silent
identical(getOption(show.error.messages), TRUE)) {
cat(msg, file = stderr()) .Internal(printDeferredWarnings())
} invisible(structure(msg, class = try-error))})
6: try(nlm(f = fprof_deriv, x = x, p = parHInt, b = parFIX, hessian =
T, iterlim = 2000, check.analyticals = F, stepmax = 10), silent =
T)
7: fMLEPFP(datanpboot, FALSE)
8: fboot.jack.mleP(data, B)
9: fsimbiasIGLDP(N[i], parH[, j], R, B)
10: fmainsimbiasP(1, 100)
Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
Selection: Selection: Selection: 2
Save workspace image? [y/n/c]: n
Process R:2 finished at Thu Mar 25 08:06:26 2010
sessionInfo()
R version 2.7.2 beta (2008-08-16 r46368)
i486-pc-linux-gnu
locale:
LC_CTYPE=es_CL.UTF-8;LC_NUMERIC=C;LC_TIME=es_CL.UTF-8;LC_COLLATE=es_CL.UTF-8;LC_MONETARY=C;LC_MESSAGES=es_CL.UTF-8;LC_PAPER=es_CL.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=es_CL.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
version
_
platform i486-pc-linux-gnu
arch i486
os linux-gnu
system i486, linux-gnu
status beta
major 2
minor 7.2
year 2008
month 08
day 16
svn rev 46368
language R
version.string R version 2.7.2 beta (2008-08-16 r46368)
__
R-help@r-project.org mailing list
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.
__
R-help@r-project.org mailing list
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.