Displaying 1 result from an estimated 1 matches for "xmui".
Did you mean:
mui
2012 May 31
0
Problem in a programming
...not running correctly, Can anyone please correct this Please!
dopt<-function(beta)
{
set.covar<-seq(-1,1,0.001)
n<-length(set.covar)
xi<-sample(set.covar,5)
## start of information loop
infor<-function(xi)
{
mu<-exp(t(xi)*beta)/(1+exp(t(xi)*beta))
mui<-mu%*%t(1-mu)
xmui<-mui%*%t(xi)
info<-xmui%*%xi
}
## end of information loop
info<-infor(xi)
##start of loop to replace every element of x (i loop)
for(i in 1:5)
{
## start of a loop to replace each element of xi with every element of set.covar (j loop)
for(j in 1:n)
{
x<-xi
xi<-replace(x,i...