Как записать статистику хи-квадрат как функцию в R и чем использовать ее в optim для определения параметров распределения вероятностей.Любая помощь будет оценена, спасибо.
n<-28
prob<-function(K, M){
(choose(M+x-2, x-1)*choose(K-M+n-x,n-x+1))/choose(K+n-1,n)}
# prob schould be a vector of functions for x from 1 to 29
# something like that
# prob<-c()
# for(x in 1:(n+1)){prob[x]<-function(K, M){
# (choose(M+x-2, x-1)*choose(K-M+n-x,n-x+1))/choose(K+n-1,n)}}
# but this is not work.
N<-3754
chi<-function(K, M){sum((frequency[,2]-N*prob)^2/(N*prob))}
optim( c(-5.54, 0.37), chi, hessian=TRUE)