R语言-第八章(数学运算与模拟).ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,第八章,数学运算与模拟,1、数学函数,exp(x)；log(x)；log10(x)；sqrt(x),abs(x),sin(x)；cos(x),min(x);max(x),which.min(x);which.max(x),pmin(x);pmax(x)#对多个向量进行逐元素对比，并返回各个最大值,sum(x);prod(x)#求和/求积,cumsum(x);cumprod(x)#前K个元素累计求和/求积,round(x);floor(x);ceiling(x)#四舍五入/向下取整/向上取整,factorial(x)#阶乘,x-c(2,3,4),y-c(3,1,5),z pmin(x,y,z),1 2 0 2,cumprod(y),1 3 3 15,factorial(y),1 6 1 120,案例：求函数的最大最小值,nlm(function(x)return(x2-sin(x),8),$minimum,1-0.2324656,$estimate,1 0.4501831#最小值在改点取得,$gradient,1 4.024558e-09#取得最小估计值时的梯度,$code#最小梯度的调用次数,1 1,$iterations#迭代次数,1 5,optim(8,function(x)return(x2-sin(x),$par,1 0.45#取得最小值的X,$value,1-0.2324655#最小值,$counts,function gradient,32 NA,$convergence#收敛,1 0,$message,NULL,微积分,微分,D(expression(exp(x2),x),exp(x2)*(2*x),积分,integrate(function(x)x2,0,1),0.3333333 with absolute error 3.7e-15,2、统计分布函数,d对应概率密度函数,dnorm()dchisq()dbinom(),p对应累计分布函数,pnorm()pchisq()pbinom(),q对应分布的分位数,qnorm qchisq qbinom,r对应随机数生成函数,rnorm rchisq rbinom runif,order()还可以对数据框排序,y y,v1 v2,1 def 2,2 ab 5,3 zz 1,r r,1 3 1 2,r r,1 2 1 3,rank()排序略有不同,x rank(x),1 4.0 2.5 1.0 2.5,解线性方程组用solve(),其本质是求解 形式的方程组,a b solve(a,b),1 0 1,其它线代函数,t():转置,qr()：QR分解,det():求矩阵行列式,集合运算,x-c(1,2,5),y-c(5,1,8,9),union(x,y)#并集,intersect(x,y)#交集,setdiff(x,y)#差集,setequal(x,y)#检验是否相等,2%in%x#检验是否是集合中元素,choose(n,k)#检验从n元素抽取k元素的方法数,x y union(x,y),1 1 2 5 8 9,intersect(x,y),1 1 5,setdiff(x,y),1 2,setequal(x,y),1 FALSE,2%in%x,1 TRUE,choose(5,2),1 10,5、用R做模拟,例1,x=4)#大于等于四返回1.否则返回0,1 0.1975,例2：求解 ，其中X、Y都服从标准正太分布。,sum-0,nreps-100000,for(i in 1:nreps),xy-rnorm(2),sum sum,1 56267.34,另一种更易理解的写法,emax-function(nreps),x-rnorm(2*nreps),maxxy-pmax(x1:nreps,x(nreps+1):(2*nreps),return(mean(maxxy),emax(100000),1 0.5635426,