[R] package(moments) issue

[Sania Wadud]

Hi Bill,

Thanks for prompt reply and letting me know a way around it.

I have more than 1200 observations and not all the values are the same.
However, my data points are quite similar, for example,
0.079275, 0.078867, 0.070716  in millions and etc. I have run the data
without converting it to millions and I still get the same error
message. As I have kurtosis value, it should be fine for the time being.
However, in case if you have any other explanation behind the
error, I will highly appreciate it.

Best wishes,
Sania

On Thu, 15 Oct 2020 at 20:43, Bill Dunlap <williamwdunlap using gmail.com> wrote:

> moments::anscombe.test(x) does give errors when x has too few values or if
> all the values in x are the same
>
> > moments::anscombe.test(c(1,2,6))
> Error in if (pval > 1) pval <- 2 - pval :
>   missing value where TRUE/FALSE needed
> > moments::anscombe.test(c(2,2,2,2,2,2,2,2))
> Error in if (pval > 1) pval <- 2 - pval :
>   missing value where TRUE/FALSE needed
>
> You can use tryCatch() to return some special value where these errors
> occur.  E.g., use
>
> tryCatch(anscombe.test(x), error=function(e) list(p.value=NA))
>
> instead of ascombe.test (this assumes you will only look at p.value later
> - add other entries as needed.).
>
> -Bill
>
> On Thu, Oct 15, 2020 at 12:09 PM Sania Wadud <sania.wadud using gmail.com>
> wrote:
>
>> Hi all,
>>
>> While running the anscombe.test in R, I'm getting an error of *Error in if
>> (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed* for a
>> few time series columns whereas for most of the series the function is
>> working fine. I have checked for those specific columns for missing
>> values.
>> However, there is no NA/NAN value in the dataset.
>>
>> I have also run kurtosis for those time series and I can get results. I
>> have noticed the same issue was reported in the following post for python
>> (scipy.stat) and the report mentions that there is a bug in moments
>> package
>> in R. I was wondering whether the issue has been resolved.
>>
>> https://github.com/scipy/scipy/issues/1950
>>
>> Could anyone please help me with the error?
>>
>> I look forward to hearing from you.
>>
>> Regards,
>> Sania
>>
>>         [[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]]