similar to: Bimodal distribution

Hi all, Is there any grid search function to compute the maximum likelihood for a random variable? thanks, Samor --------------------------------- [[alternative HTML version deleted]]

Hello everybody, I am trying to install the library stepplr under windows (http://www.maths.bris.ac.uk/R/web/packages/stepPlr/index.html), in order to use the function plr, but I still have problem to find the right link for this purpose! I am very thankful for your help! Samor

Dear all, I am very thankful, if you could tell wheather it is possible to write paste("Mus., 10 ", expression(mu)," g", sep="") Thank you in advance, Samor --------------------------------- [[alternative HTML version deleted]]

Dear R users, I'm looking for a test of bimodality in order to make some decisions about how to procede with an analysis algorithm. I have not come across any such tests in my readings and discussions apart from the Rao which appears to be applicable to cyclic data. The data I'm interested in characterizing as uni- or bimodal are frequency x amplitude spectra of consonant speech sounds,

Hello, Is it possible to have two meta-plots in one graph (not par(mfrow=c(2,1))? But somthing like  library(metafor)  library("igraph")    if (interactive()) {     forest(dat.Treat$RR, ci.lb=dat.Treat$lower, ci.ub=dat.Treat$upper, xlab="Relative Risk",slab=dat.Treat$ID,refline=1)     forest(dat.Control$RR, ci.lb=dat.Control$lower, ci.ub=dat.Control$upper, xlab="Relative

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.

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,

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

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

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

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

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

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 am calling the following code with the loop! It makes 3 out graphs but empty! COuld you help me plase on that? Thank you in advance z1 <- dbConnect(MyData, "something", "A1", "A2") for (tt in c("xyz", "abc", "m1")) { message(paste("Here", tt, "!!!")) mydata <- dbReadTable(z1, tt)

I am creating a CDF plot function more user-friendly than any default r function. Depending upon the bimodality of the data (it is often bimodal), or any other strange data trends, the points can end up gathering in just about any corner of the plot. So, when I add a legend, whether I choose to add it in the bottom right, top left, or wherever, it will sometimes end up putting the legend right

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

Hi, I am running the following loop, but it takes hours to run as n is big. Is there any way "apply" can be used? Thanks. ### Start nclass <- dim(data)[[2]] - 1 z <- matrix(0, ncol = nclass, nrow = nclass) n <- dim(data)[[1]] x <- c(1:nclass) # loop starts for(loop in 1:n) { r <- data[loop, 1:nclass] classified <- x[r == max(r)]