Displaying 1 result from an estimated 1 matches for "int_3".
Did you mean:
int3
2011 Oct 20
1
R code Error : Hybrid Censored Weibull Distribution
...n,j) *
(pweibull(x2,shape=p,scale=lamda))^(j-1) * (1 -
pweibull(x2,shape=p,scale=lamda))^(n-j) *
dweibull(x2,shape=p,scale=lamda)
int_2
}
v2[j]<-integrate(int2,lower=0,upper=t)$value
}
sum(v2)
}
#calculation for t3
fb3<-function(t,r){
v3<-numeric(0)
for(j in 1:r){
int3<-function(x3){
int_3<- (p/lamda)^2 * j * choose(n,j) *
(pweibull(x3,shape=p,scale=lamda))^(j-1) * (1 -
pweibull(x3,shape=p,scale=lamda))^(n-j) *
dweibull(x3,shape=p,scale=lamda)
int_3
}
v3[j]<-integrate(int3,lower=0,upper=t)$value
}
sum(v3)
}
a<-c(fb1(t,r),fb2(t,r),fb2(t,r),fb3(t,r))
A<-matrix(a,nrow=2,nco...