similar to: Mixture Analysis

Dear All I am trying to estimate a mixture of regression and get the following error using the mixreg package: Error in y - yhat : non-conformable arrays The instruction I used were: x <- as.matrix(LRHUN) y <- as.matrix(LRINTER) TS <- list(list(beta=c(3.0,1.0),sigsq=1,lambda=0.4), list(beta=c(0.0,1.0),sigsq=1,lambda=0.6)) prova <- mixreg(x,y, ncomp=2, theta.start=TS)

Hi, I woul like to know if it is possible to have a "R code" to generate EM Algorithm for a normal bivariate mixture. Best regard, S.F.

Dear R users it's my first time using R and I am specially interested in doing analyses with the "mixreg" package. I use W2000 but it seems the code for R is available only for Linux Os. Anyway I I downloaded the *tar.gz file of the package for R and the converted it in a *zip file. Apparently the installation process was ok, but when loading it I got the following error: Error in

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to

Hello! I am working on a manuscript on sexual dimorphism in an aquatic invertebrate, where we have estimated sexual dimorphism (SD) for 7 different traits in four populations (a total of 28 SD-estimates). We have used the following formula for estimating SD: 100 * (mean male trait value - mean female trait value)/overall trait mean). Then, we have used these SD-estimates to perform a GLM against

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

Hello all, I am trying to apply a conjoint analysis in order to determine the best profile that captures the most preferred combination of levels of given categorical factors. For this a set of factors is given and initially a fractional factorial design has to be produced as a subset of all possible factor levels combinations, sufficient to estimate the main effects utilities. Then the

Hi, Does anyone know if the D'agostino test is available with R ? Alex

Some students of mine are trying to use gnls, the generalized non-linear least squares function within the nlme library, to study evolutionary questions where correlations between traits at the species level are non-independent because of the evolutionary relatedness of the species. Specifically, they're using a non-linear function (log(sexual dimorphism) ~ log(a + b*variation in mating

Hi, I would like to minimize the value of x1-x2, x2 is a fixed value of 0.01, x1 is the quantile of normal distribution (0.0032,x) with probability of 0.7, and the changing value should be x. Initial value for x is 0.0207. I am using the following codes, but it does not work. fr <- function(x) { x1<-qnorm(0.7,0.0032,x) x2=0.01 x1-x2 } xsd <- optim(0.0207, fr,

Hi all, I have problems with years of dates using "chron" package. I don't understand why R by this istruction: > dates("01/02/29",out.format="d/m/year") [1] 02/Jan/2029 > dates("01/02/30",out.format="d/m/year") [1] 02/Jan/1930 reads "29" as 2029 and "30" as 1930. How could I change to read "00" to

Dear all, It's possible to obtain the residuals with AR time series function. Thanks in advance, Massimiliano Talarico Pianificazione Commerciale - Area Crm Unicredit Xelion Banca S.p.A. Via Pirelli 32 - 20124 Milano Tel.: 02 67360 525 Fax: 02 67738 525 www.xelion.it [[alternative HTML version deleted]]

Hi, I have a quite "simple" setup for a particular customer that loves redundancy and failover. PDC + BDC with LDAP Passwords on two 389-ds in multimaster node + several samba member servers Actually pointing singularly on both the systems everything works great. As soon as I modify my passdb backend line from the single form to the form containing both backends that is from passdb

Hello all, I would like to fit a mixture model whose components are normal distributions confined in a closed interval. Since there are already several packages for EM, I would like to extend one of these instead of writing a new script from scratch. What would be the best way to customize such an existing package for doing that? Which one is the best in terms of extensibility? (e.g mixtools,

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 sirs, I'm using the Splancs package to compute standard errors of the estimates of a spatio-temporal k function. When I use as spatial and temporal distances too long vectors (respectively 60 and 80 entries) for a dataset of 1000 observations, R gives me the message Error: cannot allocate vector of size 180000 Kb Reached total allocation of 512 Mb. I ran the function memory.size()

Dear all, I have the following R dataframe: set.seed(11) (df <- data.frame(ID=rep(1:10,1:10),a=factor(sample(1:4,55,rep=T)) )) where ID is an identification code. I need to create a new variable "b" in which I would paste the full group of "a" variable according to ID variable. For example for ID=2 "b" should be 13, for ID=3 "b" should

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

Dear R-help, I am doing clustering via finite mixture model. Please suggest some packages in R to find clusters via finite mixture model with continuous variables. And also I wish to verify the distributional properties of the mixture distributions by fitting the model with lognormal, gamma, exponentials etc,. Thanks in advance,  warm regards,Ms.Karunambigai M PhD Scholar Dept. of Biostatistics