子程序
中相应的哑元S确设置成了real*8,这可能就是不一致的原因。(不知道你把uz_int设置成什么?)
再有
kxy_integral(k_end,k_start,uz_int)
应该是
call
kxy_integral(k_end,k_start,uz_int)
Fisher判别分析,即 LDA相应的R实现为:
MASS包中的 lad() 函数,qda() 函数
lad(x, grouping, prior = proportions ,tol = 1.0e-4, method , CV = FALSE, nu, .....)
lad(formula, data, .... ,subset , na.action )