search for: nor1mix

...also provides the 16 "Marron-Wand densities" (known in the (1d) density estimation business). Erik J?rgensen has provided "p" and "q" functions -- such that there's now the full range of "dpqr" functionality. I have recently *renamed* the package to 'nor1mix' (now on CRAN where 'normix' has gone). Note that *estimation* (i.e. fitting normal mixtures to data) is not part of this package, but rather in (the more general framework of) the package 'mclust'. Martin Maechler <maechler at stat.math.ethz.ch> http://stat.ethz.ch/~m...

...also provides the 16 "Marron-Wand densities" (known in the (1d) density estimation business). Erik J?rgensen has provided "p" and "q" functions -- such that there's now the full range of "dpqr" functionality. I have recently *renamed* the package to 'nor1mix' (now on CRAN where 'normix' has gone). Note that *estimation* (i.e. fitting normal mixtures to data) is not part of this package, but rather in (the more general framework of) the package 'mclust'. Martin Maechler <maechler at stat.math.ethz.ch> http://stat.ethz.ch/~m...

Dear R Users?? I woould like to generate mixture model data. I try two type method as followings two method, can anyone tell me which is right? or better? 1. generate two vectors data y1 and y2 from f1 and f2 seperately, and calculate y=alpha1*y1+alpha2*y2, 2. generate prob=unif(1), if prob <alpha1 y=y1, else y=y2. Can anyone give me some idea? Thanks so much! Goeland goeland at

Hi! I need to fit a data that shows up as two gaussians partially superimposed to the corresponding gaussian distributions, i.e. data=c(rnorm(100,5,2),rnorm(100,-6,1)) I figured it out how to do it with mle or fitdistr when only one gaussian is necessary, but not with two or more. Is there a function in R to do this? Thank you very much in advance, Luis

I've had a need for multinomial "random number generation" occasionally. And other people too. The following code is currently in the (very small ``not very high importance'') CRAN package normix --- which I will rename to "nor1mix" very seen because of a ``name registration'' problem I want to add "this" (well the functionality) to a standard package -- "mva" probably. The reason for my post is to ask about importance for C code (and a an official API to that). 1) I think we have no...

I have not followed the history of this thread, but I am quite flummoxed as to why the OP is rewriting code to estimate parameters from an univariate Gaussian mixture model when alternatives such as EMCluster (which generally appears to handle initialization better than MClust) exist. Or perhaps there is more to it in which case I apologize. But I thought that I would make the OP aware of the

Recently, I have been working with some data that look like two overlapping gaussian distributions. I would like to either 1) determine the mean and SD for each of the two distributions OR 2) get some (bayesian ?) statistic that estimates how likely an observation is to belong to the left-hand or right-hand distribution In case I'm using the wrong language, my data looks something like

...different number of modes, with and without infinite peaks and also some standard densities like the normal and the uniform. There is also a small intersection (e.g. the claw density) with the collection of normal mixtures introduced by Marron and Wand and implemented in R in the package "nor1mix". Similar to the test bed functions by Donoho and Johnstone (Blocks, Bumps etc.) commonly used in regression or the "Peppers" and "Lena" images popular in image analysis, the densities in this collection should be useful for testing and comparing new and existing density e...

Hi R Users, Sorry if its out of topic. I would like to ask you about twin peaks - normal distribution. How R can handle it, any example to explain it in R. Thanks, regards, Unung

Hello, I am trying to find out if R can do the following: I have a mixture of normals say f = 0.2*Normal(2, 5) + 0.8*Normal(3,2) How do I find the difference in the densities at any particular point of f and at Normal(2,5)? -- Thanks, Jim. [[alternative HTML version deleted]]

...different number of modes, with and without infinite peaks and also some standard densities like the normal and the uniform. There is also a small intersection (e.g. the claw density) with the collection of normal mixtures introduced by Marron and Wand and implemented in R in the package "nor1mix". Similar to the test bed functions by Donoho and Johnstone (Blocks, Bumps etc.) commonly used in regression or the "Peppers" and "Lena" images popular in image analysis, the densities in this collection should be useful for testing and comparing new and existing density e...

...orrect ?interval = c(lower,upper)? if ?sign(f(x))? does not satisfy the requirements at the interval end points; see the ?Details? section. We've had this function, called safeUroot() in our package copula for a while now, where an earlier not-exported version has been in my package nor1mix even longer. When I was tempted to introduce it into yet another CRAN package I maintain, I decided that this may be a sign that such a simple [ utility for / generalization of ] uniroot() should probably rather be added to R itself. The function definition, also visible in R's devel.sources,...

I am trying to fit a normal linear model with response y and predictor x and two factors sex and group. I would like each combination of sex and group to have individual slopes and then subsequently have parallel slopes. I tried the model y ~ x*sex*group and it seemed to work for the first case.. Is this how it is supposed to be done? For the second the model y ~ sex + group seems to work. In

...t","mvoutlier","mvpart","mvtnorm","NADA","ncomplete","ncvar","negenes","network","neural","nFDR","NISTnls","nlme","nlmeODE","nlreg","nnet","nor1mix","norm","normalp","NORMT3","nortest","noverlap","npmc","numDeriv","nws","Oarray","odesolve","onion","optmatch","orientlib","ouch","outliers",&qu...

....5), MAc (1.0.6), MAd (0.2), mboost (2.0-3), mc2d (0.1-7), mclust (3.4.3), MCMCglmm (2.03), mhsmm (0.3.2), mixAK (0.8), MKmisc (0.6.1), monomvn (1.8-1), MplusAutomation (0.2-2), MplusAutomation (0.2-3), multcomp (1.1-6), mvpart (1.3-1), ncvreg (2.0), neuralnet (1.3), nlstools (0.0-10), noia (0.94), nor1mix (1.1-2), party (0.9-9993), pcaPP (1.8), penalized (0.9-30), pgirmess (1.4.4), phyclust (0.1-5), picante (1.0-1), plotrix (2.8-3), plsdof (0.2-0), pmg (0.9-42), portfolio (0.4-5), portfolioSim (0.2-6), prefmod (0.8-19), PropCIs (0.1-5), PwrGSD (1.16), pyramid (1.2), qpcR (1.2-6), QRMlib (1.4.5), qvc...