Just read your email about mixreg. The problem is the assignment y <- as.matrix(LRINTER) The y argument has to be a ***vector***, not a matrix. Deep in the bowels of mixreg, yhat is formed as a matrix of fitted values for all of the models in the mixture. (The first column of yhat consists of the fitted values from model 1, the second of the fitted values from model 2, und so weiter.) If y is a vector then y - yhat is the matrix whose columns are y - (the respective columns of yhat). This is standard R syntax. If y is a matrix, R won't touch the subtraction with a sterilized barge-pole, and rightly so. Remember that in R there is a big difference between a vector and 1-dimensional array, much as the two may seem to resemble each other. cheers, Rolf Turner rolf at math.unb.ca P. S. The usual protocol for asking questions about ***contributed*** packages is to attempt firstly to contact the author of the package. (In this case me.) Only if the author is a recalcitrant guttersnipe, and fails to respond, should you bug the r-help list about the problem. R. T. ===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===+===
Massimiliano Marinucci wrote:
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)
The data set has about 1200 observations.
I checked out that x and y have the same number of obs.
Also examples work fine.
Any hints on what is going on? Thanks
best regards
M
**********************************************
Massimiliano Marinucci
http://personales.ya.com/max_mar/
Ph.D Candidate in Economics
Fundamentos del Analisis Econ?mico II
(Econom?a Cuantitativa)
Facultad de CC.EE.
Universidad Complutense Madrid
Campus de Somosaguas
Madrid - Spain
**********************************************