I see it in three different Mac builds, including a quite recent local R-devel build. It boils down to this:> dtruncnorm(numeric(0), mean=6.7, sd=1.38, a=-Inf, b=9)Floating point exception: 8 which looks like a bug in the truncnorm package, where dtruncnorm() is unprepared for a zero-length argument. (The indirect cause is fitdistrplus:::test1fun, which makes calls like the above.) -pd> On 14 Mar 2020, at 19:42 , Bert Gunter <bgunter.4567 at gmail.com> wrote: > > Inline. > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along and > sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > > On Sat, Mar 14, 2020 at 10:36 AM |Juergen Hedderich <j.hedderich at t-online.de> > wrote: > >> Dear R-help list members, >> >> the R Session aborted without any 'comment' for the following 'small >> example': >> >> /library(fitdistrplus) >> library(truncnorm) >> >> filter <- c(4.98, 8.60, 6.37, 4.37, 8.03, 7.43, 6.83, 5.64, 5.43, 6.88, >> 4.57, 7.50, 5.69, 7.88, 8.98, 6.79, 8.61, 6.70, 5.14, 7.29) >> >> fit <- fitdist(filter, "truncnorm", fix.arg=list(a=-Inf, b=9), >> start=list(mean=mean(filter), sd=sd(filter)), >> optim.method="L-BFGS-B", >> lower=c(-0.1, -0.1), upper=c(Inf, Inf))/ >> >> R worked fine in an 'older' R-version (environment). Can anyone help me? >> > > I can't. But without including the specifics of the older and newer > software (including OS version, maybe), maybe no one can. > See ?sessionInfo > > > >> >> Many thanks in advance. >> >> Best regards >> >> J. Hedderich >> >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> 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. >> > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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.-- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com
On Linux it says "Program received signal SIGFPE, Arithmetic
exception". I
think the only way to get a SIGFPE (floating point exception) any more (on
machines with IEEE floating point arithmetic) is taking an integer modulo
zero, which do_druncnorm does when length(x) is 0:
const double cx = x[i % n_x];
Should R catch SIGFPE and turn it into an error condition?
Bill Dunlap
TIBCO Software
wdunlap tibco.com
On Sat, Mar 14, 2020 at 12:19 PM peter dalgaard <pdalgd at gmail.com>
wrote:
> I see it in three different Mac builds, including a quite recent local
> R-devel build.
>
> It boils down to this:
>
> > dtruncnorm(numeric(0), mean=6.7, sd=1.38, a=-Inf, b=9)
> Floating point exception: 8
>
> which looks like a bug in the truncnorm package, where dtruncnorm() is
> unprepared for a zero-length argument.
>
> (The indirect cause is fitdistrplus:::test1fun, which makes calls like the
> above.)
>
> -pd
>
> > On 14 Mar 2020, at 19:42 , Bert Gunter <bgunter.4567 at
gmail.com> wrote:
> >
> > Inline.
> >
> > Bert Gunter
> >
> > "The trouble with having an open mind is that people keep coming
along
> and
> > sticking things into it."
> > -- Opus (aka Berkeley Breathed in his "Bloom County" comic
strip )
> >
> >
> > On Sat, Mar 14, 2020 at 10:36 AM |Juergen Hedderich <
> j.hedderich at t-online.de>
> > wrote:
> >
> >> Dear R-help list members,
> >>
> >> the R Session aborted without any 'comment' for the
following 'small
> >> example':
> >>
> >> /library(fitdistrplus)
> >> library(truncnorm)
> >>
> >> filter <- c(4.98, 8.60, 6.37, 4.37, 8.03, 7.43, 6.83, 5.64,
5.43, 6.88,
> >> 4.57, 7.50, 5.69, 7.88, 8.98, 6.79, 8.61, 6.70, 5.14,
7.29)
> >>
> >> fit <- fitdist(filter, "truncnorm",
fix.arg=list(a=-Inf, b=9),
> >> start=list(mean=mean(filter), sd=sd(filter)),
> >> optim.method="L-BFGS-B",
> >> lower=c(-0.1, -0.1), upper=c(Inf, Inf))/
> >>
> >> R worked fine in an 'older' R-version (environment). Can
anyone help
> me?
> >>
> >
> > I can't. But without including the specifics of the older and
newer
> > software (including OS version, maybe), maybe no one can.
> > See ?sessionInfo
> >
> >
> >
> >>
> >> Many thanks in advance.
> >>
> >> Best regards
> >>
> >> J. Hedderich
> >>
> >>
> >> [[alternative HTML version deleted]]
> >>
> >> ______________________________________________
> >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more,
see
> >> 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.
> >>
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > 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.
>
> --
> Peter Dalgaard, Professor,
> Center for Statistics, Copenhagen Business School
> Solbjerg Plads 3, 2000 Frederiksberg, Denmark
> Phone: (+45)38153501
> Office: A 4.23
> Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.
>
[[alternative HTML version deleted]]
Uh... yes?! On March 14, 2020 2:22:14 PM PDT, William Dunlap via R-help <r-help at r-project.org> wrote:>On Linux it says "Program received signal SIGFPE, Arithmetic >exception". I >think the only way to get a SIGFPE (floating point exception) any more >(on >machines with IEEE floating point arithmetic) is taking an integer >modulo >zero, which do_druncnorm does when length(x) is 0: > const double cx = x[i % n_x]; >Should R catch SIGFPE and turn it into an error condition? > >Bill Dunlap >TIBCO Software >wdunlap tibco.com > > >On Sat, Mar 14, 2020 at 12:19 PM peter dalgaard <pdalgd at gmail.com> >wrote: > >> I see it in three different Mac builds, including a quite recent >local >> R-devel build. >> >> It boils down to this: >> >> > dtruncnorm(numeric(0), mean=6.7, sd=1.38, a=-Inf, b=9) >> Floating point exception: 8 >> >> which looks like a bug in the truncnorm package, where dtruncnorm() >is >> unprepared for a zero-length argument. >> >> (The indirect cause is fitdistrplus:::test1fun, which makes calls >like the >> above.) >> >> -pd >> >> > On 14 Mar 2020, at 19:42 , Bert Gunter <bgunter.4567 at gmail.com> >wrote: >> > >> > Inline. >> > >> > Bert Gunter >> > >> > "The trouble with having an open mind is that people keep coming >along >> and >> > sticking things into it." >> > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) >> > >> > >> > On Sat, Mar 14, 2020 at 10:36 AM |Juergen Hedderich < >> j.hedderich at t-online.de> >> > wrote: >> > >> >> Dear R-help list members, >> >> >> >> the R Session aborted without any 'comment' for the following >'small >> >> example': >> >> >> >> /library(fitdistrplus) >> >> library(truncnorm) >> >> >> >> filter <- c(4.98, 8.60, 6.37, 4.37, 8.03, 7.43, 6.83, 5.64, 5.43, >6.88, >> >> 4.57, 7.50, 5.69, 7.88, 8.98, 6.79, 8.61, 6.70, 5.14, >7.29) >> >> >> >> fit <- fitdist(filter, "truncnorm", fix.arg=list(a=-Inf, b=9), >> >> start=list(mean=mean(filter), sd=sd(filter)), >> >> optim.method="L-BFGS-B", >> >> lower=c(-0.1, -0.1), upper=c(Inf, Inf))/ >> >> >> >> R worked fine in an 'older' R-version (environment). Can anyone >help >> me? >> >> >> > >> > I can't. But without including the specifics of the older and newer >> > software (including OS version, maybe), maybe no one can. >> > See ?sessionInfo >> > >> > >> > >> >> >> >> Many thanks in advance. >> >> >> >> Best regards >> >> >> >> J. Hedderich >> >> >> >> >> >> [[alternative HTML version deleted]] >> >> >> >> ______________________________________________ >> >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> >> 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. >> >> >> > >> > [[alternative HTML version deleted]] >> > >> > ______________________________________________ >> > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> > 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. >> >> -- >> Peter Dalgaard, Professor, >> Center for Statistics, Copenhagen Business School >> Solbjerg Plads 3, 2000 Frederiksberg, Denmark >> Phone: (+45)38153501 >> Office: A 4.23 >> Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com >> >> ______________________________________________ >> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >> 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. >> > > [[alternative HTML version deleted]] > >______________________________________________ >R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see >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.-- Sent from my phone. Please excuse my brevity.