![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
R
默默小艺
乐
展开
-
R与RStudio的安装与第一个hello world程序
呵,今天学点R(默认有梯子)安装R与RStudio博主滴系统是Windows10 64版本(非ubuntu与mac) R下载 这里是清华镜像下载的 R安装 RStudio下载 这里是官网下载 RStudio安装 这里因为安装好了,提示不能再次打开安装包,所以就没有安装图片了,大家就基本默认下一步下一步那样(此处无坑)打开RS一般会出...原创 2018-07-29 11:26:41 · 4288 阅读 · 1 评论 -
R学习之统计算法与R优化包(极大似然估计) --(R语言编程)-----数模
library(MASS);attach(geyser);hist(waiting,freq = F)#mnf<-function(pa,data){ x<-dnorm(data,pa[2],sqrt(pa[4])) y<-dnorm(data,pa[3],sqrt(pa[5])) pdf=pa[1]*x+(1-pa[1])*y l=sum(log(...原创 2018-08-06 15:42:17 · 19733 阅读 · 3 评论 -
R学习之统计算法与R优化包(矩估计法也即Newton) --(R语言编程)-----数模
矩估计说明:R代码结果DONE!!!矩估计说明:就是把求方程解的Newton法过渡到统计学里面,达到参数估计的目的,其实也是对于方程的求解。 直接上例题 R代码#矩估计#下面是之前的Newtons代码Newtons<-function(fun,x,ep=1e-5,it_max=100){#itmax表示最大迭代次数 index<...原创 2018-08-06 10:25:11 · 3852 阅读 · 2 评论 -
R学习之统计算法与R优化包(Newton法) --(R语言编程)-----数模
Newton法说明:Newton法R程序:其中部分说明练习题1代码结果练习题2结果DONE!!!Newton法说明: Newton法R程序:例 {x21+x22−5=0(x1+1)∗x2−(3∗x1+1)=0{x12+x22−5=0(x1+1)∗x2−(3∗x1+1)=0\left\{\begin{array...原创 2018-08-05 17:59:21 · 11836 阅读 · 2 评论 -
R学习之统计算法与R优化包(二分法) --(R语言编程)-----数模
二分法其中函数说明结果作业题:结果DONE!!!二分法 R语言编程:#(二分法)dichotomy<-function(f,a,b,eps=1e-5){ if(f(a)*f(b)>0) list(fail="find root is fail!") else{ repeat{ if(abs(b-a)&...原创 2018-08-05 11:19:50 · 7875 阅读 · 2 评论 -
R学习之蒙特卡罗积分 --(R语言编程)-----数模
蒙特卡罗方法介绍问题1其中函数说明结果问题2结果问题3结果问题4结果问题5结果问题6结果DONE!!!蒙特卡罗方法介绍 问题1∫10sin(x)xdx∫01sin(x)xdx \int_{0}^{1} \frac {sin(x)}xdxR语言编程:#0到1 sin(x)/x积分(投点法...原创 2018-08-04 19:42:49 · 20774 阅读 · 6 评论 -
R学习之统计实验(十)--新乒乓问题(R语言编程)-----数模
问题盒中有12个乒乓球.其中9个是新的,第一次比赛时从中选出3个,比赛后仍放回盒中,第二次比赛时再从盒中选出3个.(1)求第二次取出的球都是新球的概率;(2)又已知第二次取出的球都是新球,求第一次取出的球都是新球的概率.分别作理论和模拟分析. 直接R模拟(大数定律)newpingpang<-function(n){ m<-0;z<-0 for(i in 1...原创 2018-07-30 20:03:10 · 1063 阅读 · 0 评论 -
R学习之统计实验(九)--和的偶数概率问题(R语言编程)-----数模
问题从0-9中随机不放回选4个数,求组成偶数的概率.分别作理论和模拟分析.直接R模拟(大数定律)eventest<-function(n){ m<-0 for(i in 1:n){ x<-sample(x = 0:9,size = 4,replace = F) if(sum(x)%%2==0){m=m+1} } rt<-c...原创 2018-07-30 18:09:31 · 1248 阅读 · 0 评论 -
R学习之统计实验(八)--统计次数概率问题(R语言编程)-----数模
问题从1~10随机有放回选7个数,求10恰好出现2次的概率.分别作理论和模拟分析直接R模拟(大数定律)doubleten<-function(n){ m<-0 for (i in 1:n) { u<-runif(7) x<-floor(u*10)+1 if(sum(x==10)==2){m=m+1} } rt<...原创 2018-07-30 17:36:01 · 6268 阅读 · 0 评论 -
R学习之统计实验(七)--大于1问题(R语言编程)-----数模
问题令U1U1U_1,U2U2U_2,⋅⋅⋅······是一个i.i.d序列,且UiUiU_i~UUU(0,1).令N =min {n:∑ni=1Ui∑i=1nUi\sum_{i=1}^n U_i>1} 用统计模拟方法求E(N)直接R模拟(大数定律)注意括号各种英文输入,R语言编写时候看不出明显区别morethanone <-function(n){ ...原创 2018-07-30 13:50:00 · 3631 阅读 · 5 评论 -
R学习之统计实验(六)--生日问题(R语言编程)-----数模
问题美国数学家伯格米尼曾经做过一个别开生面的实验:在个盛况空前、人山人海的世界杯赛场上,他随机地在某号看台上召唤了22个球迷,请他们分别写下自自己的生日,结果竟发现其中有两人同生日.怎么会这么揍巧呢?通过分析求出理论概率,并运用模拟方法求出概率的近似值.人数分别为40, 50, 64时,又如何? 直接R模拟(大数定律)注意括号各种英文输入,R语言编写时候看不出明显区别birth...原创 2018-07-30 11:45:03 · 4747 阅读 · 1 评论 -
R学习之统计实验(五)--赶火车问题(R语言编程)-----数模
问题问题一列火车从A站开往B站,某人每天赶往B站赶火车.他己了解到火车从A站到B站的运行时间是服从均值为30分钟.标准差为2分钟的 正态分布.火车大约13点离开A站.此人大约13:30到达B站.火车离开A站 的时刻及概率和此人到达B站的时刻及率如表所示.用模拟实验的方法 求他赶上火车的概率. 直接R模拟(大数定律)train<-function(n){ m=0 ...原创 2018-07-30 11:22:26 · 3551 阅读 · 2 评论 -
R学习之统计实验(四)--蒲丰投针(R语言编程)-----数模
问题法国数学家蒲丰于1777年完成的投针实验是历史上第一个统计模拟实验,可以利用实验结果估计圆周率二.该实验可以表述为:向一簇距离为2a的平行线构成的平面中投掷一根长度为2L的针,求针与直线相交的概率.这里假设L直接R模拟(大数定律)注意括号各种英文输入,R语言编写时候看不出明显区别...原创 2018-07-30 10:05:55 · 7535 阅读 · 2 评论 -
R学习之统计实验(三)--报童问题(R语言编程)-----数模
问题设某报每日的潜在卖报数数服从参数为lambda的Poisson分布.每卖出一份报可得报酬a 元,卖不掉退回则侮每份赔偿b元.若某日该报童买 进n份报,试求其期望所得,并对a=1.5, b=0.6, lambda=120,买进报数n1= 100.n2 = 140分别作模拟分析. 直接R模拟(大数定律)newspaper&amp;lt;-function(n,lambda,m,a,b){...原创 2018-07-30 09:02:42 · 2510 阅读 · 2 评论 -
R学习之统计实验(二)--掷色子问题(R语言编程)-----数模
问题问题一颗股子投4次至少得到一个六点与两颗般子投24次至少得到 一个双六点,这两个事件哪一个概率更大? 直接R模拟(大数定律)注意括号各种英文输入,R语言编写时候看不出明显区别sice&lt;-function(n){ p1=0;p2=0; for (i in 1:n) { u1&lt;-runif(4);x1&lt;-floor(u1*6)+1; ...原创 2018-07-29 21:30:00 · 7532 阅读 · 2 评论 -
R学习之统计实验(一)--电梯问题(R语言编程)-----数模
问题有r个人在一楼进入电梯楼上共有n层设每个乘客在任何层楼出电梯的概率相同.试建立一个概率模型,求直到电梯中的人下完为止,电梯需停次数的数学期望并对r=10,n=7进行模拟验证直接R模拟(大数定律)注意括号各种英文输入,R语言编写时候看不出明显区别elevator &lt;- function(r,n,m){ #r个人,n为1楼以上楼层数,m为数据个数 y=0 for...原创 2018-07-29 16:37:40 · 2126 阅读 · 0 评论 -
R学习之统计算法与R优化包(极大似然估计练习) --(R语言编程)-----数模
题目R代码#l=0data=rcauchy(n=1000,location = 0,scale = 10)f<-function(pa){ result= -log(pi*pa*(1+((data)/pa)^2)) -sum(result)}ML.op=optimize(f,interval = c(0,100))#不需要参数初值,interval定义参数...原创 2018-08-06 19:15:13 · 4323 阅读 · 2 评论