Displaying 1 result from an estimated 1 matches for "invltd".
2009 Jun 19
1
a difficulty in boot package
Hi,
I have a problem in programming for bootstrapping.
I don't know why it show the error message.
Please see my code below:
#'st' is my original dataset.
#functions of 'fml.mlogl','pcopula.fam4','ltd','invltd' are already defined
boot.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&...