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(10000,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(10000, 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.
