[R] parametric bootstrap

2008-07-28 Thread laila khalfan
Hi I am using the boot function in boot package. I am facing a problem in getting the values of bias and st.error in the output. r<-36 n<-40 shape<-2 theta11<-exp(1) # (=2.718282) theta21<-exp(.5) #( =1.648721) m0<- function(Ti,Tj) #a function that generates the MLestimates { loglik<-function(ti

[R] parametric bootstrap

2008-07-26 Thread laila khalfan
HiI am using the boot function in boot package. I am facing a problem in getting the values of bias and st.error in the output. r<-36n<-40shape<-2theta11<-exp(1) # (=2.718282)theta21<-exp(.5) #( =1.648721)m0<- function(Ti,Tj) #a function that generates the MLestimates{ loglik<-function(ti,t

Re: [R] parametric bootstrap

2008-07-26 Thread Prof Brian Ripley
On Sat, 26 Jul 2008, laila khalfan wrote: Hi I am trying to find a parametric bootstrap confidence interval and when I used the boot function I get zero bias and zero st.error? What could be my mistake? The way you simulated the 'parametric bootstrap'. Thank you and take care. Laila

[R] parametric bootstrap

2008-07-25 Thread laila khalfan
Hi I am using the boot function in boot package. I am facing a problem in getting the values of bias and st.error in the output. r<-36n<-40shape<-2theta11<-exp(1) # (=2.718282)theta21<-exp(.5) #( =1.648721) m0<- function(Ti,Tj) #a function that generates the estimates{ loglik<-function(ti,tj,

[R] parametric bootstrap

2008-07-25 Thread laila khalfan
Hi I am trying to find a parametric bootstrap confidence interval and when I used the boot function I get zero bias and zero st.error? What could be my mistake? Thank you and take care. Laila [[alternative HTML version deleted]] __ R-help@r-proj