similar to: testing for bimodal distribution

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

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

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.

>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

Hi, I was looking for a package that would help with outlier detection for bimodal distributions. I have tried 'outliers' and 'extremevalues' packages, but am not sure if they are ok for bimodal distribution. Any help would be highly appreciated! thanks, [[alternative HTML version deleted]]

Hello everyone, I'm trying to use mle from package stats4 to fit a bi/multi-modal distribution to some data, but I have some problems with it. Here's what I'm doing (for a bimodal distribution): # Build some fake binormally distributed data, the procedure fails also with real data, so the problem isn't here data = c(rnorm(1000, 3, 0.5), rnorm(500, 5, 0.3)) # Just to check

We currently ship a PAE 32-bit domU that we can trivially make bimodal, except that if we set it to "bimodal", then older Xens will default to thinking the domU is not PAE: 353 dsi->pae_kernel = PAEKERN_no; 354 if ( dsi->__elfnote_section ) 355 { 356 p = xen_elfnote_string(dsi, XEN_ELFNOTE_PAE_MODE); 357 if ( p != NULL && strncmp(p,

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

For several years, I have been using Splus to analyze an ongoing series of datasets that have a bimodal distribution. I have used the following functions, in particular the ms() function, to estimate the parameters: two means, two standard deviations, and one proportion. Here is the code I've been using in S: btmp.bi <- function(vec, p, m1, m2, sd1, sd2) {

Hi, Does anybody have a good tip of how to treat bimodal data to perform statistical analyses? My data set ranges from -1 to 1 (any values are posssible in between) and most data are either close to -1 or close to 1. They are the results of a two choice experiment where individuals could choose more than once in either direction and scores were calculated. Simone Simone Immler

Hi, Suppose I have the following vector (data points): > x [1] 36.0 57.3 73.3 92.0 300.4 80.9 19.8 31.4 85.8 44.9 24.6 48.0 [13] 28.0 38.3 85.2 103.6 154.4 128.5 38.3 72.4 122.7 123.1 41.8 21.7 [25] 143.6 120.2 46.6 29.2 44.8 25.0 57.3 96.4 29.4 62.9 66.4 30.0 [37] 24.1 14.8 56.6 102.4 117.5 90.4 37.2 79.6 27.8 17.1 26.6 16.3 [49] 41.4 48.9 24.1

Hi, is there any package/function that can tell if a numeric vector (continuous data) has a bimodal or trimodal distribution and caluclate the location of the corresponding modes? Thanks

Hello, Is there any test  for bimodality in R that x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) hist(x,nclass=100) Thank you in advance for any help. Regards, Samor [[alternative HTML version deleted]]

Hello R Users, I am doing a Latin Hypercube type simulation. I have found the improvedLHS function and have used it to generate a bunch of properly distributed uniform probabilities. Now I am using functions like qlnorm to transform that into the appropriately lognormal or triangularly distributed parameters for my modes. However I have a parameter which I believe is bimodally distributed,

Hi Everyone- After searching through posts and my favorite R-help websites I'm still confused about a problem. I have data which is bimodal in nature, but there is no clearly obvious separation between the two peaks. In programs such as Origin, I can deconvolute the two distributions and have it generate a "best guess" as to what the two subpopulations are which make up my

Hi All, I have a parameter that is bimodal, and I want to get some sort of linear model done with it results = some.linear.function(bimodal.param ~ factor1 + some other stuff, mydata) I want to see if factor 1 matters (it has 3 levels, of of which can be taken as baseline), i.e: summary(results) Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -2.108522

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

Hi all, I need to transform classic 32bit Fortran code to 64bit Fortran code, see the discussion [R-SIG-Mac] rnorm.halton. But I'm clearly a beginner in Fortran... Does someone already do this for his package? From here, http://techpubs.sgi.com/library/tpl/cgi-bin/getdoc.cgi?coll=linux&db=bks&fname=/SGI_Developer/Porting_Guide/ch03.html , I identify the following changes

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to