Displaying 1 result from an estimated 1 matches for "pcdif".
Did you mean:
pcdf
2008 Sep 19
0
panel data analysis possible with mle2 (bbmle)?
...oss] <- atech[1,1:cross]/cfrontr1[1,1:cross]
for (i in 2:years) {
agwth[i,1:cross]=(rr1[i-1,1:cross]-alpha1[i-1,1:cross])*beta1*(1-rr1[i-1,1:cross])+cfgwth[i,1:cross]
atech[i,1:cross]=atech[i-1,1:cross]*(1+agwth[i,1:cross])
rr1[i,1:cross]=atech[i,1:cross]/cfrontr1[i,1:cross]}
pcdiffun <- function(a,b) (a-b)/b
pcdif <- pcdiffun(atech[2:years,],aobs[2:years,]);sd1 = sd(pcdif)
rt1 = -sum(dnorm(pcdif,mean=0,sd=sd1,log=TRUE))}
##############here is where the function ends
library(bbmle)
maxlk <- mle2(exfunc,
start=list(c1 = -3.74466306894815,c2 = -0.2401359...