Dear all,
I have a question on how to set up the starting model in POLYCLASS and
make sure the terms in the starting model retained in the final
POLYCLASS model.
In the function POLYMARS, this can be done using the STARTMODEL option.
See below for example, I started with model
y= b0 + b1*X1 + b2*X2 + b3*X4 + b4*X5 + b5*X2*X5 + e
> m00 <- matrix(c(
1, NA, 0, NA, 1,
2, NA, 0, NA, 1,
4, NA, 0, NA, 1,
5, NA, 0, NA, 1,
2, NA, 5, NA, 1),nrow = 5, ncol=5, byrow=TRUE);
> m2 <- polymars(response=PID2$y, predictors=PID2[,1:7],
startmodel=m00)
> summary(m2)
But I could not figure out how this works for POLYCLASS. There is an
option FIT in POLYCLASS, which needs to be a POLYCLASS object though.
Any suggestion or information is greatly appreciated.
Sincerely,
Xiaogang Su
================================
Xiaogang Su, Assistant Professor
Department of Statistics and Actuarial Science
University of Central Florida
Orlando, FL 32816
(407) 823-2940 [O]
[EMAIL PROTECTED]
http://pegasus.cc.ucf.edu/~xsu/
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
