Displaying 1 result from an estimated 1 matches for "vpsi".
Did you mean:
vps
2010 May 24
2
How to set parameters constraints in a function?
...pa *<- exp(param[5]);
for (i in 2:n) {psi[i] <- a0 + a1*x[i-1] + b1*psi[i-1]}
lam <- gamma(kappa)/gamma(kappa+(1/alpha))
ll <-
n*log(alpha/gamma(kappa))+kappa*alpha*sum(log(x))-n*kappa*alpha*log(lam)-kappa*alpha*sum(log(psi))-lam^(-alpha)*sum((x/psi)^alpha)
return(list(maxl=-ll,vpsi=psi))
}
Thanks in advance!!
Cheers,
Carol
[[alternative HTML version deleted]]