similar to: double gaussian with mixdist: what's wrong?

Dear R User, I am an electrical engineering student and have just come across a curve fitting problem. I need to find the constituent Gaussian distribution curves fitting the data attached in Workbook1.txt here. I tried to use Mixdist on R but ran into following problem. Can you suggest me where I am going wrong? > super <-read.table("Workbook1.txt",,sep ="\t") >

Hello there, I am using nls the first time for a non-linear regression with a logistic growth function: startparam <- c(alpha=3e+07,beta=4000,gamma=2) fit <- nls(dataset$V2~(( alpha / ( 1 + exp( beta - gamma * dataset$V1 ) ) ) ),data=dataset,start=startparam) Everything works fine and i get good results. Now I would like to improve the results using my DUMMY Variable (dataset$V6) the

Hi, I am working on a population of an invasive clam. The data are the size of each clam per station (2mm on average). Each station is found at a different distance from a power nuclear station, so at different water temperatures. The fist step I want to do is to identify cohort size at each station or (zone of water temperature). The second step will be to see whether the size or number of

Hi, I am a new ext3 user and I am having problems similar to what many people are facing on this list. I have done everything as the emails suggest. I am using SuSE 7.3 with the the 2.4.10 kernel and ext3 compiled as a module. Because it is 2.4.10 kernel, I created all my journals from a rescue disk and everything seems to be working right for the non-root filesystems. My fstab has all my

Hello, I'm trying to run constrOptim. It returns to me an error about the fact that constraints arguments (ui and ci) are non compatibles: > optout= constrOptim(startparams, f=ImpulseSS, grad=grImpulse, ui=UI, ci=CI, data=gexp[k,], t=t) Error in ui %*% theta : non-conformable arguments I would like to point out that I can calculate that product in the command line: > UI %*%

Dear List, We are using the mix() function in the mixdist package to fit mixture models to some of our data. The package provides a function to compare the fits of nested models using an ANOVA function. However, we were wondering whether there are methods that can be used to compare models that differ in the number of distributions fitted and/or the shape of the distributions. Is there a way,

Hi, All, I fitted a 3-component normal mixture model with the mixdist package in R. How can I get the density of a new data after I fit the model? Is there any function to do it? Thanks, Cindy [[alternative HTML version deleted]]

Dear R-help, my aim is to fit a constrained 2 component gamma mixture using mix within package mixdist. The constraints being fixing both of the means and both of the sigmas. Thus the only unconstrained parameter I wish to fit is pi the proportions of the mixtures. This is the code I used: PP2Pars<-data.frame(c(0.23,0.77), c(9.174,29.54), c(9.983,12.369));

Hi all, I'm trying to install the mixdist package with: install.packages("mixdist") But I'm getting the following errors: * installing *source* package ‘mixdist’ ... ** R ** data ** preparing package for lazy loading ** help *** installing help indices ** building package indices ... Warning in utils::data(list = f, package = package, lib.loc = lib.loc, envir = dataEnv) :

Dear R users, Please help me on using "mix" function under package "mixdist". My data distribution shows there are two components for the mixture distribution: left part is an exponential and right part is a normal. So I plan to use "gamma" mixture distribution to estimate the parameters. Here is what I am using for the "mix" function. Test<-mix(x,

The Wine development release 1.7.10 is now available. What's new in this release (see below for details): - Window maximization in the Mac driver. - AVI compressor implementation. - Thread local storage support in dynamically loaded libraries. - Beginnings of a Task Scheduler implementation. - Extended IPX protocol support. - Various bug fixes. The source is available from the

Dear all, I have a plot with two gaussian distributions: ggplot( tofdf, aes( x = tof, y = counts ) ) + geom_histogram( stat = "identity", position = "stack", fill = 2 ) + geom_smooth( method = "nls", formula = y ~ (a/b)*exp(-(x-c)^2/(2*b^2)) + (d/e)*exp(-(x-f)^2/(2*e^2)), se=F, start=list(a=100, b=1, c=695, d=100, e=1, f=710),