Displaying 1 result from an estimated 1 matches for "xpar".
Did you mean:
par
2005 Mar 02
1
Warning: number of items to replace is not a multiple of replacement length
...t)
{
library(mvtnorm)
#============================
Makeham<-function(tt)
{
a2=0.030386513
a3=0.006688287
b3=0.039047537
t<-tt-20
h.t<-a2+a3*exp(b3*t)
S.t<-exp(-a2*t+a3/b3*(1-exp(b3*t)))
return(S.t*h.t)
}
#===========================
trans<-function (age)
{
indic=c(2,2)
lage=log(age)
xpars=c(2.1862908,7.5528077,8.5806697,2.3319461,8.8959507,9.1380187,0.3720293)
beta<-c(0,0)
alpha1<-c(0,0)
alpha2<-c(0,0)
beta[1]<-xpars[1]
alpha1[1]<-xpars[2]
alpha2[1]<-xpars[3]
beta[2]<-xpars[4]
alpha1[2]<-xpars[5]
alpha2[2]<-xpars[6]
corr<-matrix(rep(0,4),nc=2)
corr[1,1...