Displaying 1 result from an estimated 1 matches for "ms22".
Did you mean:
ms02
2012 Aug 28
1
Optim Problem
...he code command for optim i want it to contain both the t and q data so i can obtain the correct estimate. Is there any suggestion as to how this can be done. I have tried h<-c(t,q) but it is not working because q lies within an interval.
rate<-15;n<-100;a<-40;b<-60;rr<-1000
ms11=ms22=0
bia11=bia22=0
t<-rexp(a,rate)
for(i in 1:rr){
C1<-runif(b,0,rate)
C2<-rexp(b,rate)
f2 <- function(C1, C2) {
? r <- pmax(C1 , C2 + C1)
? cbind(C1, r)
}?
m<-f2(C1,C2)
x[1:b]<-(m[,1])
u<-x[1:b]
x[1:b]<-(m[,2])
v<-x[1:b]
q<-cbind(u,v)
h<-c(t,q)
z<-function(dat...