Hi Are there any functions that de-convolute data into a given number of clusters, rather like the NPMLE GLIM macros from Murray Aitkin and Brian Francis? Basically I would like to code into R the same approach but include the possiblility of some data being censored. In principle the formulae are the same (just replace the likelihood function) but I haven't managed to get my head round the model building problems. I am looking for an elegant way of doing the following: Consider a model where the variables are partitioned into two groups - call them fixed and random. I have these as data frames - each can be a mixture in principle of vectors and factors. Then I want to fit: Y~fixed+F:random where Y is the response (possibly including the censoring indicator), F is a factor. The response, the factor F and the random and fixed data frames are stacked N long for N mass points. F indexes the mass point. I have been trying to build up the model and use gnlr() except the complexities are beating me at the moment. Should I use a list or a model frame? Can anyone advise? John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
I'm running a discriminant analysis in R (0.64.1) to compare it with SPlus 4.5R2. The following command line works fine in SPlus but gives an error in R. I've only used R for a little while so I'm not certain here what R (or lda) is complaining about. The dependent variable (sarich.na[,3]) is an alpha categorical variable, if that makes a difference. I'm using version VR5.3 (file name VR5.3pl037.zip). lda.out<-lda(sarich.na[,3]~., data=sarich.na[,4:32]) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type Is this an lda issue or an R issue? Dr. Marc R. Feldesman email: feldesmanm at pdx.edu email: feldesman at ibm.net fax: 503-725-3905 "Math is hard. Let's go to the mall" Barbie Powered by: Monstrochoerus - the 300 MHz Pentium II -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._