Dear R community,
I am using the R package `interval` in order to perform some modelling
tests of the
NPMLE convergence in the case of censoring. So all I am doing is drawing a
sample
from exponential distribution, making it a censored sample and computing
the
NPMLE of
its distribution function. But when run on Linux Calculate 10.4 the
program
keeps
crashing and reporting a segmentation fault
after the call to the `icfit` function when the sample size gets to 70.
When run on Windows 7 it seems to be fine.
That is why I am totally confused and have decided to ask for help.
I have attached the code I am running which results in a segmentation
fault
if run on Linux Calculate.
It has the seed set to the value which leads to this error. But it is
important to note
that if the parameters used in the program and the seed are changed it
doesn't
necessarily crash.
Here is the description of my R version and OS:
sessionInfo()
R version 2.10.1 (2009-12-14)
i686-pc-linux-gnu
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
After calling `icfit` the program quits with the following output
(I have replaced the output concerning the arguments passed to
initcomputeMLE by
<<<<< arguments passed to initcomputeMLE>>>>> so that the description of
the
output wouldn't be too long):
*** caught segfault ***
address 0xc, cause 'memory not mapped'
Traceback:
1: .Call("ComputeMLEForR", R, B, max.inner, max.outer, tol)
2: computeMLE(R, B, max.inner = max.inner, max.outer = max.outer, tol
=
tol)
3: initcomputeMLE(<<<<< arguments passed to initcomputeMLE>>>>>)
4: do.call(initfit, args = list(L = L, R = R, Lin = Lin, Rin = Rin, A
=
A))
5: doTryCatch(return(expr), name, parentenv, handler)
6: tryCatchOne(expr, names, parentenv, handlers[[1L]])
7: tryCatchList(expr, classes, parentenv, handlers)
8: tryCatch(expr, error = function(e) { call<- conditionCall(e) if
(!is.null(call)) { if (identical(call[[1L]],
quote(doTryCatch))) call<- sys.call(-4L) dcall<-
deparse(call)[1L] prefix<- paste("Error in", dcall, ": ")
LONG<- 75L msg<- conditionMessage(e) sm<- strsplit(msg,
"\n")[[1L]] w<- 14L + nchar(dcall, type = "w") + nchar(sm[1L], type
= "w") if (is.na(w)) w<- 14L + nchar(dcall, type =
"b")
+ nchar(sm[1L], type = "b") if (w> LONG)
prefix<- paste(prefix, "\n ", sep = "") } else prefix<- "Error :
" msg<- paste(prefix, conditionMessage(e), "\n", sep = "")
.Internal(seterrmessage(msg[1L])) if (!silent&&
identical(getOption("show.error.messages"), TRUE)) {
cat(msg,
file = stderr()) .Internal(printDeferredWarnings()) }
invisible(structure(msg, class = "try-error"))})
9: try(do.call(initfit, args = list(L = L, R = R, Lin = Lin, Rin = Rin,
A = A)))
10: icfit.default(L = left, R = right)
11: icfit(L = left, R = right)
Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
I would greatly appreciate any help provided.
Sincerely yours,
Yuliya Matveyeva.
______________________________________________
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.
