Displaying 1 result from an estimated 1 matches for "mlebinorm".
2008 Feb 24
0
problem with ML estimation
...ng results
which are very inaccurate.
I am sending my programme. It would be of great help if someone can point
me, where I am going wrong.
############################################################
library(MASS)
Sigma <- matrix(c(4,2.8,2.8,4),2,2)
y<-mvrnorm(n=1000, rep(0, 2), Sigma)
mlebinorm<- function(startval, y)# startval = Initial Values to be passed ,
y= data
{
startval<-as.numeric(startval)
oneside=matrix(c(0,0,0,0,-1,0,0,0,0,1),5,2)
otherside = c(-0.9999999,-0.999999)
n<- ncol(y)
nf<- function(x)
{
mu1<-x[1]
mu2<-x[2]
sig1<-x[3]
sig2<-x[4]...