similar to: smartpred depends on fitted() in flexmix?

I am fitting a regression model with a bs term and then making predictions based on the model. According to some info on the internet at http://www.stat.auckland.ac.nz/~yee/smartpred/DummiesGuide.txt there are some problems with using predict.lm when you have a model with terms such as bs,ns,or poly. However when I used one of the examples they said would illustrate the problems I get virtually

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

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

Hi, I am trying to understand how to use the flexmix package, I have read the Leisch paper but am very unclear what is needed for the M-step driver. I am just fitting a simple linear regression model. The documentation is far from clear what the FLXmclust function does, but, it in principle could do all I need, however, I do not get sentible results as if I try the following the result is poor:

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'

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

Dear R-helpers & Friedrich Leisch, I'm trying to fit a mixture of regression models on attached data set by doing the following: agl <- read.table("agl.txt") mod1 <- flexmix(resp~gng+csc|Subject,data=agl, model=FLXglm(family="binomial"),k=2) The result is a (varying) number of the following warnings: Warning messages: 1: non-integer #successes in a

I am receiving the following error while executing KLdiv method: kl<-KLdiv(y) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "KLdiv", for signature "data.frame" the code is as follows: require(flexmix) KLdiv(y) the object y (snippet, actual object contains 188 entries) is as follows: s1p5 s2p5 1 0.6 0.8 2 0.8

I've been trying to install some R packages such as mclust and flexmix on linux but have had the following error messages. > I've been trying to install mclust on my notebook which has linpus linux lite os and I have installed R as well as some packages all right. However, when I tried to install mclust, it gave me the following messages. Any suggestions? Tim ++++++++++++++++

Hi, Does anyone know how to do a coord_trans() in which the y-axis is tranformed into (for example) -1000/y? Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

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

Dear r-helpers, Does anyone have a straightforward example of putting together three unrelated (expect for a common y-axis) xyplot() figures in what would be in base graphics a par(mfrow = c(1, 3)) arrangement? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102

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

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

Hi, I would like to use the effect() function (actually a slightly modified version of it) on the output of the lmer() function in the lme4 package. But the effects package requires the nlme pacvkage, which is incompatible with lme4. Workaround? ______________________________________________ Professor Michael Kubovy University of Virginia Department of Psychology for mail add: for FedEx or

Dear R-helpers, How do I transform v <- c('insd-otsd', 'sppr-unsp') into c('insd--otsd', 'sppr--unsp') ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903

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)