similar to: Avoiding choosing parameters with mix[mixdist]

Dear all, I am trying to fit a double gaussian to some data using the mixdist package: --- begin code --- ? library(mixdist) ? time <- seq(673,723) ??counts <- c(3,12,8,12,18,24,39,48,64,88,101,132,198,253,331,419,563,781,1134,1423,1842,2505,374,6099,9343,13009, ? ? ? ? ? ? ? ? ? ?15097,13712,9969,6785,4742,3626,3794,4737,5494,5656,4806,3474,2165,1290,799,431,213,137,66,57,41,35,27,27,27) ?

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

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,

Dear R users, I tried to find a solution in the search list, but I cannot find it. I would like to read a .txt file with, let say, three variables, with two of which have repeated values in a number a columns. An example: The variables: Treat, x1, x2. The values: A 2.5 3.4 2.7 5.6 5.7 5.4 10.1 9.4 ... B 5.3 5.4 6.5 7.5 1.3 4.5 10.5 4.1 ... ... In the first column, the letters represent the

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, I try to use a very large file (~3 Gib) with the filehash package. The length of the dataset is around 4,000,000 obs. I get this message from R while trying to "load" the dataset (named "cc084.csv"): > dumpDF(read.csv("cc084.csv", header=T), dbName="db01") Erreur : impossible d'allouer un vecteur de taille 15.6 Mo (French) Error:

Thanks for telling me that you could not get my message, I hope this work better... so my question was: I built a population matrix to which I applied the fonction eigen in order to find the main parameters about my population. I know that the first eigen value correspond to lambda or exponential growth rate of my population. My problem is that I want to have the 95% confidence interval of the

Hi, I would like to extract the fitted values from a model using LMER but only for the fix portion of the model and not for the fix and random portion (e.g it is the procedure outpm or outp in SAS). I am aware of the procedure fitted() but I not sure it give the fitted values both for the fixed and random or only the fixed. I looked in the r help and the r list and I haven?t not found much

Hi, Quick question. With an old phone system a receptionist receiving a call has 1 button to push to transfer calls to a specific extension, with Asterisk, a receptionist would actually put the caller on hold, pick up another line, call the extension, ask if the person is available, hang up pick up the caller again and transfer. To me it's seems a long way to simply do a receptionist

I've been beating my head against this one and just can't figure it out. I hope someone here may have an answer. The employees using the workstations on this network are getting increasingly upset with this problem. The problem is wildly varying logon and logoff times over the network. This is definitely not a matter of long profile transfers. An individual can log onto a workstation

Dear R community, I'm trying to remove a loop from my code but I'm stock and I can't find a good way to do it. Hopefully one of you will have something clever to propose. Here is a simplified example: I have a squared matrix: > nom.plac2 <- c("102", "103", "301", "303","304", "403") > poids2 <- matrix(NA, 6,6,

Full_Name: Richard Simard Version: OS: Linux Submission from: (NULL) (132.204.25.139) In the random number generator MARSAGLIA_MULTICARRY, your algorithm is different than the one in the original reference that you give in the R manual: Marsaglia in his post to the mailing list {\it sci.stat.math} on September 29, 1997. The last line in the R program has a ^ while Marsaglia's algorithm has

I hope someone can enlighten me on this. Situation: NT network, Samba PDC, about 20 NT 4.0 workstations. log(s).smbd are created per machine for easier analysis (as log.smbd.<NetBios Name>). Not often, but often enough to be of concern, are errors in request_oplock_break that seem to indicate that another smbd process that should be listening for break requests on a UDP socket

Hi, Does anybody know where I can get an untouched smb.conf file? When I used samba-swat, it rewrote smb.conf and deleted out the comments that was in the original file that came with samba. Minh --- samba-request@lists.samba.org wrote: > Send samba mailing list submissions to > samba@lists.samba.org > > To subscribe or unsubscribe via the World Wide Web, > visit >

Hi, All, I want to fit a normal mixture model. Which package in R is best for this? I was using the package 'mixdist', but I need to group the data into groups before fitting model, and different groupings seem to lead to different results. What other package can I use which is stable? And are there packages that can automatically determine the number of components? Thank you, Cindy

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

Dear all, I am attempting to model some one-dimensional data using a mixture model of non-central Student's t distributions. However, I haven't been able to find any R package that provides this functionality. Could there be a way to "manipulate" the EM algorithms from the mixdist or mixtools package to fit the model, or do you have any other suggestions? If anyone could help