问题标题：
R语言.如何写codeQ1：>xmean(exp(-x+1)^2)/dexp(x)#能解释这个dexp(x)什么意思,然后整条式的意思.Q2：如何在R里打出pi,不能的话就用3.1415926代替?Q3：用montecarlointegration有什么要注意的地方.Q4：
问题描述：
R语言.如何写code
Q1：
>xmean(exp(-x+1)^2)/dexp(x)#能解释这个dexp(x)什么意思,然后整条式的意思.
Q2：如何在R里打出pi,不能的话就用3.1415926代替?
Q3：用montecarlointegration有什么要注意的地方.
Q4：∫e^(-x)dx上限：无限,下限：0.如何用motecarlointegration算这个?那∫e^(-x^3)dx呢?
要是前面有常数又怎样处理,例如∫(1/根号2pi)*e^(-x)dx上限1,下限0
问题比较多,非常感谢先
黄山回答：
q1,dexp是指数分布的概率密度； r*是生成*分布的随机数,p是*分布的概率分布,d是概率密度,还有q是d的反函数,貌似 R可向量操作,数乘向量是用数乘向量内...
田帆回答：
要是前面有常数又怎样处理，例如∫(1/根号2pi)*e^(-x)dx上限1，下限0
黄山回答：
常数可以提到积分符号前面；也可以直接放到f里面，一次积分出来；
比如 ∫ (1/根号2pi)*e^(-x)dx
f <- function(x) { exp(-x)/sqrt(2*pi) } #常数就是那个sqrt（2*pi）
integrate(f, 0, 1)
查看更多
