similar to: finite mixture model, using flexmix

Dear R-users, Is there a way to get p-values for a one-sided hypothesis test about a poisson mean? Thanks, Jan Wijffels University Center for Statistics W. de Croylaan 54 3001 Heverlee Belgium tel: +32 (0)16 322784 fax: +32 (0)16 322831 <http://www.kuleuven.be/ucs> http://www.kuleuven.be/ucs Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm [[alternative HTML version

Dear useRs, I very much like the effect display of the proportional odds model on page 29 (Figure 8) of the following paper by John Fox: http://socserv.mcmaster.ca/jfox/Papers/logit-effect-displays.pdf It really gives a very concise overview of the model. I would like to use it to illustrate the proportional odds mixed models we fit here for a project on Diabetes but I can't seem to reproduce

Dear R-users, Does anyone of you have binary files for the packages rpart.permutation, snow and rsprng. I would like to use them in my classification tree. I know they are still at the 0.x development stage, though. Where can I get information on how to compile tar.gz files? I'm using windows XP with R 2.1.0. Thanks, Jan --------------------------------------------------------------------

Dear R-users, I tried to install the R(D)COM server version 1.35, to connect MS Excel and R. I am using R-version 2.1.0 on a Windows environment. And I am using MS Excel 2000 in a dutch version. I get the following error when I open Excel: "Compile error, can not find the project or the library" (translated from dutch). And it indicates LEFT in blue bold. Does anyone of you know what

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

Dear Listers, My post might be somewhat OT. Currently, I am trying to use flexmix to build a finite mixture model. For instance, I am getting the prior probability and coefficients for each latent class from training data. Is there a way to get the posterior probablity and prediction of a new dataset? What I am thinking is to apply the prior prob and coefficient from training set to testing data

Dear useRs, FlexMix version 1.0-0 has been released on CRAN. FlexMix implements a general framework for finite mixtures of regression models using the EM algorithm. FlexMix provides the E-step and all data handling, while the M-step can be supplied by the user to easily define new models. Existing drivers implement mixtures of standard linear models, generalized linear models and model-based

Dear useRs, FlexMix version 1.0-0 has been released on CRAN. FlexMix implements a general framework for finite mixtures of regression models using the EM algorithm. FlexMix provides the E-step and all data handling, while the M-step can be supplied by the user to easily define new models. Existing drivers implement mixtures of standard linear models, generalized linear models and model-based

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 [

Hello list, Has anyone had any luck creating an M-step driver for negative binomial regression for use with package flexmix? I've had a look here: http://cran.r-project.org/web/packages/flexmix/vignettes/flexmix-intro.pdf as well as poking around in the flexmix source, but I haven't had much luck getting anything to work. I can't figure out how to a) come up with an initial estimate

Hi, I am using the package "flexmix" and would like get some assistance. I am trying to run two equations jointly Y1=X1B+E1 Y2=X2G+E2 So that I have X and Y in a matrix format and would like to run the latent class model using flexmix. Though, my problem here is that Flexmix automatically generates intercepts. I have intercept for both of the equations that my X matrix looks like

Dear R Users, R Core Team, I currently wonder how to predict the probability of an event with new data resulting from a finite mixture. I read the documentation of the flexmix package and the examples of applications provided on CRAN but I could not find how to predict (except "manually" but I am looking for a simpler solution) the final probability of the mixture (for each individual)

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

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

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

Problem on flexmix when trying to apply signature developed in one model to a new sample. Dear R Users, R Core Team, I have a problem when trying to know the classification of the tested cases using two variables with the function of flexmix: After importing the database and creating a matrix: BM<-cbind(Data$var1,Data$var2) I see that the best model has 2 groups and use: ex2

Dear list: I would like to recreate how the artificial data set BregFix was generated in package flexmix (thanks Bettina and Friedrich). The data set is thoroughly described in Grun and Leisch's Computational Statistics & Data Analysis 51(11) :5247-5252 but references to the appropriate seed number(s) are missing (providing these details was certainly beyond the scope of the authors'

Hi, I was going through the examples in smartpred. It seems there's an unstated dependency on the fitted() function in package flexmix. n = 20 set.seed(86) x = sort(runif(n)) y = sort(runif(n)) library(splines) fit = lm(y ~ ns(x, df=5)) plot(x, y) lines(x, fitted(fit)) # won't work w/o prior loading of the flexmix package. newx = seq(0, 1, len=n) points(newx, predict(fit,

Hello all, I am trying to fit a truncated mixture model and I wrote a driver for flexmix following the example in the vignette, but it doesn't work for me: it assigns all data points to one component only, e.g.: > > source('bugged.R') > > Call: > flexmix(formula = x ~ 1, k = 2, model = truncatedmodel(lower = -4, > upper = 4)) > > prior size