similar to: Anscombe-Glynn, Bonett-Seier, D'Agostino

I have a question about the D'Agostino skewness test and the Anscombe-Glynn kurtosis test. agostino.test(x, alternative = c("two.sided", "less", "greater")) anscombe.test(x, alternative = c("two.sided", "less", "greater")) The option "alternative" in those two functions seems to be the null hypothesis. In the output, the

Dear R-users, I am comparing differences in variance, skew, and kurtosis between two groups. For variance the comparison is easy: just var.test(group1, group2) I am using agostino.test() for skew, and anscombe.test() for kurtosis. However, I can't find an equivalent of the F.test or Mood.test for comparing kurtosis or skewness between two samples. Would the test just be a 1 df test on

Another bad case is > moments::anscombe.test(rep(c(1,1.1),length=35)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed I haven't checked the formulas carefully, but I suspect the problem is from taking the cube root of a negative number in z <- (1 - 2/(9 * a) - ((1 - 2/a)/(1 + xx * sqrt(2/(a - 4))))^(1/3))/sqrt(2/(9 * a)) In R, the

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

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

Hi, Does anyone know if the D'agostino test is available with R ? Alex

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

Hi, Is there an implementation of the D'agostino K2 omnibus test in an R package. I've found the D test in the fBasics package, but this doesn't give the same answer as my algorithm. I've already coded a version of the K2 and need to check the results. I have a set of spectra with gaussian distributed noise and I want a statistical test to separate those with no signal (so the

Juraj17 wrote: > > Do I have to write my own, or it exists yet? How name has it, or how can I > use it. > Try the R-function search. It return the function you are looking for as the first match. Dieter -- View this message in context: http://r.789695.n4.nabble.com/D-Agostino-test-tp3424952p3425833.html Sent from the R help mailing list archive at Nabble.com.

On Thu, Aug 20, 2020 at 7:02 AM Christoph Hellwig <hch at lst.de> wrote: > > On Wed, Aug 19, 2020 at 03:07:04PM +0100, Robin Murphy wrote: > >> FWIW, I asked back in time what the plan is for non-coherent > >> allocations and it seemed like DMA_ATTR_NON_CONSISTENT and > >> dma_sync_*() was supposed to be the right thing to go with. [2] The > >> same

On Wed, Aug 19, 2020 at 4:07 PM Robin Murphy <robin.murphy at arm.com> wrote: > > On 2020-08-19 13:49, Tomasz Figa wrote: > > On Wed, Aug 19, 2020 at 1:51 PM Robin Murphy <robin.murphy at arm.com> wrote: > >> > >> Hi Tomasz, > >> > >> On 2020-08-19 12:16, Tomasz Figa wrote: > >>> Hi Christoph, > >>> >

I searched the archives and did not find a solution, so I pose this question to those well-versed in the use of plotmath and expressions. I have a list of strings in an external CSV file which I wish to use sometimes as plot axis labels and sometimes as plot titles. These strings combine plaintext and a few mathematical expressions (Greek letters, subscripts). Moreover, I sometimes need to

On 2020-08-19 13:49, Tomasz Figa wrote: > On Wed, Aug 19, 2020 at 1:51 PM Robin Murphy <robin.murphy at arm.com> wrote: >> >> Hi Tomasz, >> >> On 2020-08-19 12:16, Tomasz Figa wrote: >>> Hi Christoph, >>> >>> On Wed, Aug 19, 2020 at 8:56 AM Christoph Hellwig <hch at lst.de> wrote: >>>> >>>> The

On Wed, Aug 19, 2020 at 03:07:04PM +0100, Robin Murphy wrote: >> FWIW, I asked back in time what the plan is for non-coherent >> allocations and it seemed like DMA_ATTR_NON_CONSISTENT and >> dma_sync_*() was supposed to be the right thing to go with. [2] The >> same thread also explains why dma_alloc_pages() isn't suitable for the >> users of dma_alloc_attrs() and

HI Philippos, Try this: dat1<- read.csv("Validation_data_set3.csv",sep=",",stringsAsFactors=FALSE) #converted to csv str(dat1) #'data.frame':??? 12573 obs. of? 17 variables: # $ Removed.AGC????????????????????????????? : num? 65.67 46.17 41.26 14.09 5.38 ... # $ Removed.SST????????????????????????????? : chr? "" "46.1658" "41.2566"

Hi, I'm a newcomer to R, having previously used SPSS. One problem I have run into is computing kurtosis. A test dataset is here: http://www.whinlatter.ukfsn.org/2401.dat > library(moments) > data <- read.table("2401.dat", header=T) > attach(data) > loglen <- log10(Length) With SPSS, I get Skewness -0.320 Kurtosis -1.138 With R: > skewness(loglen) [1]

Greetings Everybody: I generated a 1.2MB dta file based on the general social survey with Stata8 for linux. The file can be re-opened with Stata, but when I bring it into R, it says all the values are missing for most of the variables. This dataset is called "morgen.dta" and I dropped a copy online in case you are interested http://www.ku.edu/~pauljohn/R/morgen.dta looks like this

Hello I do not know very much about statistics (and English language too :-( ), then I come in search of a clarification (explanation): I found two distinct results on KURTOSIS and I do not know which of them is the correct one. Any aid will be welcome! klebyn ################ CODE rnorm(1000) -> x library(moments) kurtosis(x) skewness(x) detach("package:moments")

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

>From the R FAQ Bug reports on contributed packages should be sent first to the package maintainer, and only submitted to the R-bugs repository by package maintainers, mentioning the package in the subject line. and you are NOT the listed maintainer. Do learn to read FAQs before causing unnecessary work for other people (as the R posting guide asks). Please do as the FAQ asks. On