[R] fit for truncated normal distribution

[peter dalgaard]

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 using 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 using 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 using 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 using 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 using cbs.dk  Priv: PDalgd using gmail.com