r求含参积分

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求问大神怎么用R求含参积分！！

我想求μ = 0, σ = 1, α = 1, β = 1 和μ = 0, σ = 1, α = 1, β = −1的f(x)值，能用r实现嘛

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

1. fun = function(x,mu,sigma,alpha, beta){
   
2.   re = function(y){
   
3.     return(exp(-(y-mu)^2/(2*sigma^2))*exp(alpha+beta*y)/(1+exp(alpha+beta*y)))
   
4.   }
   
5.   return(re(x)/(integrate(re, -100, 100))$value)
   
6. }
   
7. 
   
8. 
   
9. funv = Vectorize(fun)
   
10. 
    
11. x = seq(-10, 10, len=200)
    
12. 
    
13. y1=funv(x,mu=0,sigma=1,alpha=1,beta=1)
    
14. y2=funv(x,mu=0,sigma=1,alpha=1,beta=-1)
    
15. 
    
16. plot(x, y1, type = 'n', xlab = 'x', ylab = 'y')
    
17. lines(x, y1, lty=2, lwd=2)
    
18. lines(x, y2, lty=3, lwd=2)
    
19. legend('topright', c('y1','y2'), lty=2:3, lwd=2, bty='n')

复制代码

zhou1_20 发表于 2019-10-4 08:23

跪谢大神！！！！！

thanks for sharing