求助一个简单的程序编写

除了p,b其他的都设定为一个常数。

1. x <- rf(10,2,3)
   
2. m <- 10;r <- 5;n <- length(x);k <- 0.2
   
3. alpha <- c(1,2)
   
4. beta <- c(3,4)
   
5. shape <- r+alpha
   
6. ##step2
   
7. scale_b <- beta[2]-sum(log(x[1:r]))
   
8. b <- rgamma(m,shape[2],scale=scale_b)
   
9. ##step1
   
10. scale_p <-beta[1]+sum(log(1+x^b)+(n-r)*log(1+x[r]^b))
    
11. p <- rgamma(m,shape[1],scale=scale_p)
    
12. 
    
13. ##step 3
    
14. pb <- matrix(cbind(p,b),ncol=2,byrow=F)
    
15. 
    
16. ##step4
    
17. h3_func<-function(b){
    
18.   return(scale_p^shape[1]*exp(-1*sum(log(1+x[1:r]^b))))
    
19. }
    
20. 
    
21. g_func <- function(z,y){
    
22.   return(sqrt(z^2+y^2))
    
23. }
    
24. 
    
25. 
    
26. gsel <- sum(g_func(pb[1:m,1],pb[1:m,2])*h3_func(pb[1:m,2]))/sum(h3_func(pb[1:m,2]))
    
27. 
    
28. glinex <--1*(1/k)*log(
    
29.    sum(h3_func(pb[1:m,2])/exp(k*g_func(pb[1:m,1],pb[1:m,2]))
    
30.        )/sum(h3_func(pb[1:m,2]))
    
31.     )

复制代码

What is your function g(p,b) and h_3(p,b) ?

Bayesian.var 发表于 2013-11-20 01:28
What is your function g(p,b) and h_3(p,b) ?

已经更新了g和h_3函数了，

nuomin 发表于 2013-11-19 10:36

我想问一下，在这个程序中, r表示能够观察的数量，n-r表示删失的，这个怎么表示