similar to: p-values for the dip test

Hola a tod en s, estoy intentando averiguar el número de modas en una distribución. Para ello utilizo diptest. Mi duda es que no acabo de entender cuando la información suministrada por los test suponen la existencia o no de unimodalidad/multimodalidad. Una parte de la salidad de diptest es la que pego a continuación (el resto esta en el fichero adjunto con las distribuciones kernels y las

>From library(diptest): Shouldn't the following almost always be non-significant for Hartigan's dip test? dip(x = rnorm(1000)) I get dip scores of around 0.0008 which based on p values taken from the table (at N=1000), using the command: qDiptab, are 0.02 < p < 0.05. Anyone familiar with Hartigan's dip test and what I may not be understanding? Thanks, kbrownk

Hi, I just got a value for the dip test out of my data of 0.074 for a sample size of 33. I'm trying to work out what this actually means though? Could someone help me relate this to a p-value? Thanks James

I'm not sure where to begin with this, but I was wondering if someone could refer me to an R package that would test to see if a distribution fits a bimodal distribution better than a unimodal distribution. Thanks, Andrew [[alternative HTML version deleted]]

Dear R users, I am interested in uni- bi- multimodality tests, for analysing reaction times data. I was lead to Hartigan's dip test (Ann. Statistics, 13, 1985, pp. 70-84, Applied Statistics, 34, 1985, 320-325). Not being a programmer I am unable to translate the Fortran code given in ref. 2 into a R function. I'd be glad to learn that someone already did it, or has devised a better

Has a test for bimodality been implemented in R? Thanks, John NIWA is the trading name of the National Institute of Water & Atmospheric Research Ltd.

Hi STATISTICA has a function called "Two-way joining" (see http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the reference material states that this is based on the method as published by Hartigan (found this paper: http://www.jstor.org/pss/2284710 through wikipedia). What is the relationship (if any) between the "heatmap" function in R and this technique? Is there an

Hi I have two questions, the first perhaps dumber than the second. Firstly, I have a data set, and when I plot a histogram it looks like a normal distribution. So I want to overlay a bell-shaped normal distribution on top of it, to demonstrate how similar it is to the normal distribution. I have read the help on dnorm(), rnorm(), pnorm() etc but still can't figure out how to plot a normal

hi r-help, i am doing kmeans clustering in stats. i tried for five clusters clustering using: kcl1 <- kmeans(as1[,c("contlife","somlife","agglife","sexlife",                         "rellife","hordlife","doutlife","symtlife","washlife",                       

Hello. This is not a question about a bug or even best practices; rather I'm trying to understand the philosophy or theory as to why certain portions of the R codebase are written as they are. If this question is better posed elsewhere, please point me in the proper direction. In the thread about the issues with the Tukey line, Martin said [1]: > when this topic came up last (for me) in

Dear Readers: I have two issues in nonparametric statistical analysis that i need help: First, does R have a package that can implement the multimodality test, e.g., the Silverman test, DIP test, MAP test or Runt test. I have seen an earlier thread (sometime in 2003) where someone was trying to write a code for the Silverman test of multimodality. Is there any other tests that can enable me to

Hi, I am new to bioconductor, trying to install KEGGSOAP package, but got warnings() when installing and error message when trying to load the package, can anyone suggest what went wrong? many thanks John > source("http://bioconductor.org/biocLite.R") Bioconductor version 2.11 (BiocInstaller 1.8.3), ?biocLite for help > biocLite("KEGGSOAP") BioC_mirror:

Dear experts, I am encountering a version-dependent issue. My laptop runs Ubuntu 12.04 LTS 64-bit, R 2.14.1; the issue explained below never occurred with this version of R My desktop runs Ubuntu 11.10 64-bit, R 2.13.2; what follows applies to this setup. The data I'm clustering is constituted by the rows of a 320 x 6 matrix containing integers ranging from 1 to 7, no missing data. I applied

The variable local_user was allocated by xstrdup and is not freed or pointed to in this branch. This patch adds the xfree. This entire set of patches passed the regression tests on my system. Bug found by Coverity. Signed-off-by: Kylene Hall <kjhall at us.ibm.com> --- sshconnect.c | 1 + 1 files changed, 1 insertion(+) diff -uprN openssh-4.3p2/sshconnect.c

Hi Kylie, Is it a list with only numerics in it? (I only see REALSXPs there, but obviously inspect isn't showing all of them). If so, you could load it up into one big vector and then also keep partitioning information around. Bioconductor does this (see ?IRanges::CompressedList ). The potential benefit here being that the underlying large vector could then be a big out-of-memory altrep. How

Hello, I was wondering if there were any plans for ALTREP lists (VECSXP)? It seems to me that they could be supported in a similar way to how ALTSTRING works, with Elt() and Set_elt() methods, or would there be some problems with that I?m not seeing due to lists not being atomic vectors? I was taking an approach of converting each list element (of a file-based list data structure) to an ALTREP

I can work on this. Thanks Luke. ~G On Wed, Jul 24, 2019 at 8:25 AM Tierney, Luke <luke-tierney at uiowa.edu> wrote: > If one of you wanted to try to create a patch to support ALTREP > generic vectors here are some notes: > > The main challenge I am aware of (there might be others): Allowing > DATAPTR to return a writable pointer would be too dangerous because > the GC

Hello, here is a working reproducible example which crashes R using kmeans or gives empty clusters using the nstart option with R 15.2. library(cluster) kmeans(ruspini,4) kmeans(ruspini,4,nstart=2) kmeans(ruspini,4,nstart=4) kmeans(ruspini,4,nstart=10) ?kmeans either we got empty always clusters and or, after some further commands an segfault. regards, Detlef Groth ------------ [R] Empty

Hello, I am fitting a coxph model with factors. I am running into problems when using 'survfit'. I am unsure how R is treating the factors when I fit, say: > DATA<-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo) > final<-coxph(Surv(time.sec,done)~f.pom*vo+po,data=DATA) > final.surv<-survfit((final), individual=T,conf.type="log-log")

Hi Jiefei and Kylie, Great to see people engaging with the ALTREP framework and identifying places we may need more tooling. Comments inline. On Thu, Jul 18, 2019 at 12:22 PM King Jiefei <szwjf08 at gmail.com> wrote: > > If that is the case and you are 100% sure the reference number should be 1 > for your variable *y*, my solution is to call *SET_NAMED *in C++ to reset > the