R语言.如何写code
问题描述：
R语言.如何写code
Q1：
>xmean ( exp(-x+1)^2) / dexp(x)#能解释这个dexp(x)什么意思,然后整条式的意思.
Q2：如何在R里打出pi, 不能的话就用3.1415926代替?
Q3： 用monte carlo integration 有什么要注意的地方.
Q4： ∫ e^(-x) dx 上限： 无限, 下限：0 .如何用mote carlo integration 算这个? 那∫ 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)