similar to: 'nor1mix' for 1-dimensional normal mixture distributions

The new package "benchden" 1.0.0 implements 28 benchmark densities for nonparametric density estimation that were introduced by A. Berlinet and L. Devroye ("A Comparison of Kernel Density Estimates", Pub. Inst. Stat. Univ. Paris, XXXVIII, fasc. 3, 1994, 3-59, http://cg.scs.carleton.ca/~luc/devs.html ). This collection includes a variety of densities with different degrees of

The new package "benchden" 1.0.0 implements 28 benchmark densities for nonparametric density estimation that were introduced by A. Berlinet and L. Devroye ("A Comparison of Kernel Density Estimates", Pub. Inst. Stat. Univ. Paris, XXXVIII, fasc. 3, 1994, 3-59, http://cg.scs.carleton.ca/~luc/devs.html ). This collection includes a variety of densities with different degrees of

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

>>>>> "Jacob" == Jacob van Wyk <jlvw at rau.ac.za> writes: Jacob> Could anybody please help: I simply want to generate Jacob> random samples from a multinomial distribution with Jacob> fixed n and given probability vector p. Jacob> How can I do this? I've recently written the following {for the R package "normix"}. The

Hello, I have a concrete statistical question: I have a sample of an univariate mixture of an unknown number (k) of normal distributions, each time with an unknown mean `m_i' and a standard deviation `k * m_i', where k is known factor constant for all the normal distributions. (The `i' is a subscript.) Is there a function in R that can estimate the number of normal distributions k

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

Hi, I have some very noisy, relatively sparse data; a biological response of roughly ~8 subjects at ~8 times points). I've been following the data trend using a lowess line, over-plotted with several values of bandwidth, 'f <- seq(0.3, 0.9, by=0.1)'. At this point, we have no models for these data. I wonder if there is any way under R to assign some sort of confidence interval to

Hi all, once more about the old subj :-) My data has too much various distribution families and for every particular experiment I need just to decide whether the data is "quite homogeneous" or it has two or more clusters. I've revisited the following libraries: amap, clust, cclust, mclust, multiv, normix, survey. And I didn't find any ready-to-use general

It is a pleasure for me to announce the new built for Rmetrics Version 191.0057. The source files and Windows binary packages can be downloaded from www.rmetrics.org . The new built has also been submitted to the CRAN server. Some new functions and example files have been added. Unfortunately the user guides and reference guides are not yet updated, they have still the status of Version

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

I am attempting to install a anchoring vingettes package for R and I seem to have a problem with unzipping the files in R. I use the command that can be found here <http://wand.stanford.edu/anchors/>. Here is the dialog I get after entering the command: > install.packages("anchors", dependencies = TRUE, +

Hi, Does anyone know how to compute the quantile of a mixture of four bivariate normal distriutions? Many thanks! Hannah [[alternative HTML version deleted]]

For those R user's who don't subscribe to R-sig-db and are having troubles with the latest RMySQL binary on CRAN, please read the email thread at the end of this message. RMySQL 0.7-2 does work with MySQL 5.1, however the CRAN binary is linked against the 5.0 version. Best, Jeff -------- Original Message -------- Subject: Re: [R-sig-DB] Problems with RMySQL and MySQL server version

Can someone explain to me in what way the new (dpqr)gamma parameter can be interpreted as a rate (when shape != 1)? The only gamma rate that I am aware of is the hazard rate given by dgamma/(1-pgamma), the log of which is returned by my hgamma function (event library). Jim -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello I am trying to find an automated way of fitting a mixture of normal and log-normal distributions to data which is clearly bimodal. Here's a simulated example: x.1<-rnorm(6000, 2.4, 0.6)x.2<-rlnorm(10000, 1.3,0.1)X<-c(x.1, x.2) hist(X,100,freq=FALSE, ylim=c(0,1.5))lines(density(x.1), lty=2, lwd=2)lines(density(x.2), lty=2, lwd=2)lines(density(X), lty=4) Currently i am using

Hello, I''m trying to install/update rmagick on a Debian machine using: gem install --local rmagick-2.10.0.gem and get following error: /usr/local/bin/ruby extconf.rb install --local rmagick-2.10.0.gem checking for Ruby version >= 1.8.5... yes checking for gcc... yes checking for Magick-config... yes Warning: Found more than one ImageMagick installation. This could cause problems at

Hi all, How can we use ks.test() to evaluate goodness of fit of mixtures of distributions? For example I have the following dataset: > x [1] 176.1 176.8 259.6 171.6 90.0 234.3 145.7 113.7 105.9 176.2 168.9 136.1 [13] 109.2 110.3 164.3 117.7 131.3 163.7 200.4 196.4 196.2 168.6 190.4 127.5 [25] 136.0 114.2 112.0 91.9 333.4 295.5 172.0 293.3 91.7 289.7 118.8 55.1 [37] 161.9 233.9 197.7

Hi, I am using the KernSmooth package to estimate nonparametrically bivariate density functions. However, it seems that the bandwidths (one for each co-ordinate direction) have to be selected manually. This does not apply for the univariate case, for which dpik (included in KernSmooth) uses up-to-date plug-in rules. Does anyone know about a package, or function, which estimates bandwidths

Dear Ingmar & Dave, Thanks a lot for your help and sorry for the late reply. Finally, I've found a way to separate the mixture of distributions (empirically). But the gamlss package looks great, I'm sure it will help me during my further studies. Kind regards, Susanne On 15 Jun 2009, at 20:09, Ingmar Visser wrote: > Dear Susanne & Dave, > > The gamlss package family

David Ahern <dsahern at gmail.com> writes: > On 4/18/19 8:24 AM, Toke H?iland-J?rgensen wrote: >>>> >>> >>> Understood. Hopefully in March I will get some time to come back to this >>> and propose an idea on what I would like to see - namely, the admin has >>> a config option at load time to enable driver counters versus custom map