I wonder if there are some function or package in R which can do the iterative proportional fitting.
In the exponential model f(y1,...yn)=exp(a'yi+b'(yi*yj)+.....+c'(y1*...*yn)+constance),
instead of the canonical parameters, I use maginal probability instead of a and log odds ratio instead of b. and for the order higher than 3, I use the canonical way or even assume them to be 0.
It is good if there is direct function or package in R to do the job, which will save a lot of time of me. Thank you!
Rgs,
Zhen
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help