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.

Reply via email to