Displaying 1 result from an estimated 1 matches for "fml1".
Did you mean:
3ml1
2009 Jun 19
1
a difficulty in boot package
...ot.OR<-function(data,i)
{
E=data[i,]
ml1<-glm(c_VAsex90_bf ~ trt,family=binomial,data=E)
ml2<-glm(c_VAsex90_bm ~ trt,family=binomial,data=E)
marg.covariates<-cbind(rep(1,length(E$trt)),E$trt)
dep.covariates<-cbind(rep(1,length(E$age_avr)),E$age_avr)
start<-c(ml1$coef,ml2$coef,0,0)
fml1<-optim(start,fml.mlogl,control=c(maxit=10000),hessian=F)
x<-(1+exp(fml1$par[1]))^(-1)
y<-(1+exp(fml1$par[3]))^(-1)
b<-exp(fml1$par[5]+fml1$par[6]*43)+1
p00<-ltd(-log(pcopula.fam4(exp(-invltd(x,b)),exp(-invltd(y,b)),b)),b)
p1<-exp(fml1$par[1])/(1+exp(fml1$par[1]))
p2<-exp(fml1$p...