怎么编写R语言程序

参考文献已经给出


A simple algorithm for generating random variates with a log-concave density.pdf (599.14 KB)

2013-11-4 16:57:00 上传

帮顶，求高手

有程序包...
1) http://cran.r-project.org/web/packages/logcondens/vignettes/logcondens.pdf
2) http://cran.r-project.org/web/packages/LogConcDEAD/vignettes/LogConcDEAD.pdf

ryusukekenji 发表于 2013-11-14 20:58
有程序包...
1) http://cran.r-project.org/web/packages/logcondens/vignettes/logcondens.pdf
2) http: ...

请问你是做这方面的吗?你怎么找到这包?

ym2009 发表于 2013-11-15 01:08
请问你是做这方面的吗?你怎么找到这包?

http://www.rseek.org

ryusukekenji 发表于 2013-11-18 20:30
http://www.rseek.org

你搜索的是什么？

ym2009 发表于 2013-11-18 22:53
你搜索的是什么？

A simple algorithm for generating random variates with a log-concave density.pdf

不错的资源，谢谢分享，收了

1. n <- length(y);m <- 100
   
2. a <- c(10,20,30);b <- c(10,20,30);c <- c(10,20,30);q <- c(10,20,30)
   
3. out_mat <- matrix(numeric(3*(m+1)),ncol=3)
   
4. out_mat[1,] <- c(2,2,2)#initial values
   
5. #step1-4
   
6. times <- 1
   
7. 
   
8. ##par1 =lambda,par2=y,par3=theta,par4=beta
   
9. subfunc <- function(par1,par2,par3,par4){
   
10.     subs1 <- prod((1-exp(par1*par2))^par3)
    
11.     subs2 <- prod((1-subs1)^par4)
    
12.     return(subs1*subs2)
    
13. }
    
14. 
    
15. dfunc <- function(x1){
    
16.     denomionator <- subfunc(x1,y,out_mat[1,2],out_mat[1,3])
    
17.     numerator <- subfunc(out_mat[1,1],y,out_mat[1,2],out_mat[1,3])
    
18.     return(numerator/denomionator)
    
19. }
    
20. 
    
21. 
    
22. while (times<m+1){
    
23.     theta <- rdevroye(param)
    
24.     x <- rglambda(n+a[2],b[2]+sum(y))
    
25.     u <- runif(1,0,1)
    
26.     d_compare<- dfunc(times,x,theta)
    
27.     if (u<d_compare) {
    
28.         times<<-times+1
    
29.         out_mat[times,1] <<-x
    
30.         out_mat[times,2] <<-theta
    
31.         sub1 <- (1-exp(-out_mat[1,1]*y))^out_mat[1,2])
    
32.         out_mat[times,3] <<-rgbeta(n-sum(d)+a[3],b[3]-sum(log(1-sub1))
    
33. 
    
34.     }
    
35. 
    
36. }
    
37. 
    
38. #step5
    
39. comp_mat <- out_mat[-1,]
    
40. theta_head <- (-1/c[1])*log(sum(exp(-c[1]*comp_mat[,2]))/m)
    
41. lambda_head <- (-1/c[2])*log(sum(exp(-c[2]*comp_mat[,1]))/m)
    
42. beta_head <- (-1/c[3])*log(sum(exp(-c[3]*comp_mat[,3]))/m)
    
43. #step6
    
44. theta_ge <- (sum(comp_mat[,2]^(-q1))/m)^(-1/q[1])
    
45. lambda_ge <- (sum(comp_mat[,1]^(-q2))/m)^(-1/q[2])
    
46. beta_ge <- (sum(comp_mat[,3]^(-q3))/m)^(-1/q[3])

复制代码

先看看思路吧，改天再写log-concave 的随机数生成函数

楼上 你的怎么粘贴上来的