similar to: Help with flexmix

Hi, Marion, I believe the package FlexMix provides a more generalized version of finite mixture modeling than is found in mclust/mclust02. Please see: http://cran.r-project.org/doc/vignettes/flexmix/flexmix-intro.pdf Karen --- Karen M. Green, Ph.D. Karen.Green@sanofi-aventis.com Research Investigator Drug Design Group Sanofi Aventis Pharmaceuticals -----Original Message-----

Hello, I see that "mclust" is a pacakge that handles fitting mixtures of normals. Are there any other packages out that that can handle mixtures of gammas or other exponentials? Additionally, are there any packages out there that can fit bimodal distributions without mixtures? i.e., Cobb et al. 1983 using moment recursion relations? Thank you! mw Marion Wittmann, Ph.D.

I've been trying with no success to model mixtures of Gamma distributions using the package flexmix (see examples below). Can anyone help me get it to model better? Thanks very much. -Ben ## ## Please help me get flexmix to correctly model mixtures of ## Gamma distributions. See examples below. ## library('flexmix') ## ## Plot a histogram of dat and the Gamma mixture model given

Hello, I'm trying out flexmix and mclust for the first time on some univariate data which is typically best described as lognormal, but can sometimes be gamma distributed as well. I first tried using EM on mclust assuming the data was lognormally distributed and could only get it to work in "E" mode, i.e. the equal variance mode. I could never get it to work on "V" mode [

I am trying to estimate a latent class multinomial logit regression with flexmix. I am not sure if I should do it as follows: m4<-flexmix(cbind(y,1-y)~x1+x2|id,model=FLXMRglm(family="binomial"),data=NPreg,k=2) , where id links each row with the corresponding respondent. Each respondent has 4 alternatives to chose from. y takes the value 1 only for the alternative chosen; x1 and

Hi, I have a data matrix in which there are 1000 rows by 30 columns. The first column of each row is a numeric indicating the number of trips taken to a particular location with location attributes in the following column entries for that row. I want to repeat each row based on the number of trips taken (as indicated by the number in the first column)...i.e., if 1,1 indicates 4 trips, I want

I'm trying to perform ANCOVAs in R 1.14, on a Mac OS X, but I can't figure out what I am doing wrong. Essentially, I'm testing whether a number of quantitative dental measurements (the response variables in each ANCOVA) show sexual dimorphism (the sexes are the groups) independently of the animal's size (the concomitant variable). I have attached a 13-column matrix as a data

Dear all, As I resume my dissertation work next month, I'd like to actually start an R package this time around. I haven't done so because I update my code very often (still in development phase), so running the skeleton function, running checks, building, and re-installing the package onto the system seemed like a long and tedious process. I would like to hear your experience on how

Greetings, I am using the nlm function in a minimization exercise but have consistently received the code: "last global step failed to locate a point lower than 'estimate'...." So it seems that this function is only finding a local minimum, not necessarily a global minimum. Does anyone have a suggestion or recommendation for an alternative package or function within R

Hello, Sorry if this is a repeat, having trouble finding what I need on archive and am a novice with graphics in R! I have plotted the results of a logit regression by factor, where "city" is the factor set name. The problem is city names are too long and I need to rotate them 90 degrees along the xaxis. This is what I've tried... plot(city, predicted.probability, xlab =

Greetings I'm a relatively new R user and I'm trying to build a latent class model. I've used the 'R Site Search' and it appears there's not much dialogue on these packages On mmlcr, I've gotten it working, but not sure if I'm using it correctly. On flexmix, I can only seem to get results for one class. I'm attaching my code below - if anyone

We want to exclude virtio-gpu since it's irrelevant, large, and pulls in other unwanted dependencies (modeswitching, drm). Add a second list of kmods which is a blacklist, applied after the first. This reduces the libguestfs initrd size from 39M down to 17M, with concomitant small reductions in boot time. --- src/ext2_initrd.ml | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-)

Eric, What do you think about renaming ‘src’ to ‘server’? I think it would help when newcomers navigate the source code. (And I hope that my FOSDEM talk will bring lots of new users and developers, hence why I've been working on this over the holiday.) If you agree then I'll push a commit which changes the directory name and makes the concomitant changes to makefiles, docs etc. Rich.

Dear list, I have been comparing the outputs of two packages for latent class regression, namely 'flexmix', and 'mmlcr'. What I have noticed is that the flexmix package appears to come up with a much better fit than the mmlcr package (based on logLik, AIC, BIC, and visual inspection). Has anyone else observed such behaviour? Has anyone else been successful in using the mmlcr

Dear R community, I'm hoping someone out there has perhaps done this and can share their code and/or expertise with me. I need to pull recession periods out of a hydrograph - can anyone help me with this? I want to create a subset from streamflow data that consists of just the recession curves - the decreasing runoff after the passage of a peak flow. would really appreciate any help on

Dear R users, i am wondering if anyone has some hints for this problem (i have not found a clear answer after searching the R-mailing list archive, 'help.search' in R, and R-Wiki, and the like...): let's assume that i have 4 periods compositional time series data: t=1, A=0.1; B=0.5; C=0.4 t=2, A=0.2; B=0.4; C=0.4 t=3, A=0.5; B=0.3; C=0.2 t=4, A=0.4; B=0.3;

> On Mon, 1 Apr 2019 21:49:06 +0000 (UTC) > Billy Bob via samba <samba at lists.samba.org> wrote: > > > > While 'make test' is intended to> pass generally, the reference > > > environment is Ubuntu 14.04 (soon to be> upgraded to 18.04 for the > > > next release).  > > > > By "next release," do you mean 4.10.1 or 4.11?

Dear R-help, A colleague of mine was running some code on two of our boxes, and noticed a rather large difference in running time. We've so far isolated the problem to the difference between R 1.7.1 and 1.8.0, but not more than that. The exact same code took 933.5 seconds in 1.7.1, and 3594.4 seconds in 1.8.1, on the same box. Basically, the code calls boot() to bootstrap fitting mixture

On Tue, 2 Apr 2019 15:48:41 +0000 (UTC) Billy Bob <billysbobs at yahoo.com> wrote: > > On Mon, 1 Apr 2019 21:49:06 +0000 (UTC) > > Billy Bob via samba <samba at lists.samba.org> wrote: > > > > > > While 'make test' is intended to> pass generally, the reference > > > > environment is Ubuntu 14.04 (soon to be> upgraded to 18.04

Hello, I have a data frame with 1425 observations, 539 of which are zeros. I am trying to fit the following ZINB: f3<-formula(Nbr_Abs~ Zone * Year + Source) ZINB2<-zeroinfl(f3, dist="negbin", link= "logit", data=TheData, offset=log(trans.area), trace=TRUE) Zone is a factor with 4 levels, Year a factor with 27 levels, and Source a factor with 3 levels. Nbr_Abs is counts