similar to: unimodality test

Dear R users, I've written the following functions to implement Silverman's modality test ("Using kernel density estimates to investigate multimodality", J.R. Stat. Soc. B, 43, (1981), 97-99.), and I tested them on the Chondrite data set (Good & Gaskins, J. Amer. Stat. Ass., 75, (1980), 42-56). Values for the critical window width seem OK, which is not the case for the

Hi all, I am using Hartigan & Hartigan's [1] "dip test" of unimodality via the diptest package in R. The function dip() returns the value of the test statistic but I am having problems calculating the p-value associated with that value. I'm hoping someone here is familiar with this process and can explain it. In the original article there is an example using n=63 and a

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

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

>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

Eigen seems to be giving different results in 0.64 : eigen(matrix( c(0.00000000, 0.0000000, 0.00000000, 1.000000000, 0.01262442, -0.1141633, 0.04989433, 0.006112607, 0.04781087, 0.4177869, -0.15569656, 0.074242237, 0.20167331, 0.2531109, -0.04973573, 0.034046467) ,4,4))$values R: [1] 0.35519256 -0.29294606 -0.26637369 -0.03168619 S: [1] 0.4689288205 -0.4246063863 -0.2810596564

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

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]]

Hello, I downloaded R today because I was told it has very good bootstrapping abilities. What I need to do is to program it (or use an existing program) to bootstrapping test for multimodality using nonparametric kernel density estimates as proposed by Efron and Tibshirani (1993). If anyone can get me started, I will be immensely grateful. Thanks. Sangeeta

Hi, Running Asterisk CVS-Head latest on a Dual P3 800 1Gb Ram. For some odd reason now that I have the asterisk box almost to the stage I want it, I hit a problem. I have a te405p in the system, Zap/g1 is connected to the telco as an ISDN 30, Zap/g4 is connected by ISDN Primary Rate to an Ericcson BP250 phone system. When calls come in on g1 they go straight through instantaneously to the

Dear list, I am trying to choose the scale parameter for the cgt transform but I don't know how to do it. In time I would like to be able to separate points 30 samples apart, and in frequency I would like to separate bands 0.04 Hz apart. I tried the two approaches described below and they gave me different results. I would appreciate advise on how to do this. The Rwave Gabor transform uses

>>>>> Serguei Sokol <sokol at insa-toulouse.fr> >>>>> on Tue, 30 May 2017 16:01:17 +0200 writes: > Le 30/05/2017 ? 09:33, Martin Maechler a ?crit : ... >> However, even after the patch, The example from the SO >> post differs from the result of Richie Cotton's >> function... > The explanation is quite simple.

>>>>> Serguei Sokol <sokol at insa-toulouse.fr> >>>>> on Mon, 29 May 2017 15:28:12 +0200 writes: > Sorry, I have seen it too late that we had different tab > width in the original file and my editor. Here is the > patch with all white spaces instead of mixing tabs and > white spaces. thank you - it still gives quite a few

Hi, Urgently response would be wonderful, system is a Fedora Core 2. I have a Ericsson BP250 connected to 1 port on the TE405P and another connected to a local telco ISDN30. I have been running CVS-HEAD from about a 2 months ago and upgraded it again just in cause it was a version issue (didn't fix it) but this is what I am getting. When a person calls out from an extension on the BP250 to

Dear All, A method for the estimation is univariate unimodal densities (with unknown mode) is described in "Statistical Inference under Order Restrictions" by Barlow et al.. Would anyone know whether there is an R-implementation (preferably with reference) for the estimation of univariate discrete unimodal densities (with unknown mode)? Thanks in advance for your help. Kind

Hi there, I'm doing some clustering analysis and try to find all the algorithms related to clustering in R. Here is the list of the algorithms I found. But I'm not sure if It's the complete list. Could you please check it and see if there're other ones? Thank you very much! P.S.: List of the algorithms related to clustering: (1) Hierarchical methods hclust

Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan _________________________________________________________________ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009 [[alternative HTML version deleted]]

Dear All, Does anybody know if there is a functionality in R to break histograms that show a clear bi-modal (or multi-modal) distribution into a series of unimodal histograms that added up result in the original histogram? I was thinking of using QQ-plots (for which tools are available in R), and then observing the number of times the observed quantiles cross the 1:1 line, but this only gives an

Here are three reasons for converting Fortran code, especially older Fortran code, to C: 1. The C-Fortran interface is not standardized. Various Fortran compilers pass logical and character arguments in various ways. Various Fortran compilers mangle function and common block names in variousl ways. You can avoid that problem by restricting R to using a certain Fortran compiler, but that can

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