similar to: install package mixdist

Hi [R-2.10.0; suse linux] I am having difficulty creating S4 objects in the data directory of a package. I want to create a bunch of simple S4 objects for use in the examples section of the Rd files. It says in R-exts that: " R code should be ?self-sufficient? and not make use of extra functionality provided by the package, so that the data file can also be used without having to load

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

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,

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

Just a report for those who want to try this. (Executive summary: It eventually worked.) The RPM for Fedora Core 1 did not work because it wanted libtcl8.3.so and libtk8.3.so, while this distribution had 8.4 instead of 8.3. I didn't try making soft links, thus pretending that it had 8.3. (I'm planning to install the final version of Core 2 when it comes out soon. This is a computer I

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

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

Hello, I would like to fit a mixture of two Weibull distributions to my data, estimate the model parameters, and compare the fit of the model to that of a single Weibull distribution. I have used the mix() function in the 'mixdist' package to fit the mixed distribution, and have got the parameter estimates, however, I have not been able to get the log-likelihood for the fit of this model

Hello, What I am trying to do is: Generate a density plot of a population of data. This data has a bimodal distribution so I've isolated a couple of possible sub-populations and I want to overlay these two density plots over the first to see whether they are contributing to the bimodal population. I can do this fine with plot(density(...)) and lines(density(...)) . But the resulting plots

I am trying to use Rcmdr 1.5_4 with R-2.11.1 (order to run to run the new version of misdist-0.5.3 which is built under R-2.11.1). however relimp is required for Rcmdr and the version of relimp_1.0.1 downloaded from CRAN will not work with the latest version of Rcmdr (I get error message telling me to reload it). Is there any way round this problem or will there be a new version relimp that is

Hi,all with my data,there are more than 1000 quantitative results of antibody concentrations, there may be 2 components(positive and negative), or 3 components (may be strong positive, positive, and negative), or 4-6 components. Could you tell me how to determine the number of components of the mixture model? the "anova.mix" in mixdist of the R software seems not work. my

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

It's dead code, the core handles all this directly now. The only special case is nouveau and tda988x which used one function for both legacy modeset code and -nv50 atomic world instead of 2 vtables. But amounts to exactly the same. v2: Rebase over the panel/brideg refactorings in stm/ltdc. Signed-off-by: Daniel Vetter <daniel.vetter at intel.com> Cc: Archit Taneja <architt at

It's dead code, the core handles all this directly now. The only special case is nouveau and tda988x which used one function for both legacy modeset code and -nv50 atomic world instead of 2 vtables. But amounts to exactly the same. v2: Rebase over the panel/brideg refactorings in stm/ltdc. Signed-off-by: Daniel Vetter <daniel.vetter at intel.com> Cc: Archit Taneja <architt at

It's dead code, the core handles all this directly now. The only special case is nouveau and tda988x which used one function for both legacy modeset code and -nv50 atomic world instead of 2 vtables. But amounts to exactly the same. v2: Rebase over the panel/brideg refactorings in stm/ltdc. Signed-off-by: Daniel Vetter <daniel.vetter at intel.com> Cc: Archit Taneja <architt at