Displaying 2 results from an estimated 2 matches for "svec2".
Did you mean:
vec2
2010 Mar 24
0
optimize a joint lieklihood with mle2
...the correct argument for the
mle2 function. Here is what I've tried:
##################################
the lieklihood function to be optimized is called mfun(logN, s, h=h,
cohort=cohort). I want to maximize the vectors logN and s. I define the
starting values for logN as svec1 and for s as svec2
then I write for the optimization:
m1 = mle2(mfun,start=list(svec1=logN,svec2=s),
method="L-BFGS-B",lower=lower,upper=10,
vecpar=TRUE,
data=list(h=h,cohort=cohort))
but I get the error:
Error in mle2(mfun, start = list(svec1 = logN, svec2 = s), method =
"L-BFGS-B"...
2011 Sep 01
3
betareg question - keeping the mean fixed?
...link.phi="log",
control=betareg.control(start=svec))
I understood that y~-1 could be used to give a fixed mean of 0.5 however I
get the following error:
Error in linkinv(x %*% beta + offset) :
Argument eta must be a nonempty numeric vector
I think I can work round this by using:
svec2<-c(qlogis(mean(data1$scaled)),0,0,0,0,0)
f2<-betareg(scaled ~ expt_label + grouped_hpi | expt_label + grouped_hpi,
data=data1,
+ link.phi="log",control=betareg.control(start=svec2))
This appears to work (ie doesn't return errors), but given i know the mean
to be f...