Displaying 3 results from an estimated 3 matches for "cube_root".
2020 Oct 15
0
package(moments) issue
...9 * a))
In R, the cube root (or any non-integral power) of a negative real number
is NaN and ((1 - 2/a)/(1 + xx * sqrt(2/(a - 4))))
is c. -928 in this example. I think the intent was for the cube root to be
a negative real, c. -9.76, in this case.
If that is the intent, then a fix is to define a cube_root function for
reals
cube_root <- function (x) sign(x) * abs(x)^(1/3)
and change the z<-... line to use cube_root() instead of ^(1/3).
z <- (1 - 2/(9 * a) - cube_root(((1 - 2/a)/(1 + xx * sqrt(2/(a -
4))))))/sqrt(2/(9 * a))
-Bill
On Thu, Oct 15, 2020 at 3:15 PM Sania Wadud <sania...
2020 Oct 15
2
package(moments) issue
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
2020 Oct 15
0
package(moments) issue
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