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(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. > ______________________________________________ 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.