统计建模与R软件 薛毅 陈立萍 清华大学出版社第四章课后答案

 

用R软件来求解：n=6
（1）矩估计法
x=c(0.1,0.2,0.9,0.8,0.7,0.7)
(2*mean(x)-1)/(1-mean(x))
[1] 0.3076923
#注意分子分母都要加括号，不然结果会出错
(2)极大似然估计

f=function(a) 6/(a+1)+sum(log(x))
uniroot(f,c(0,1))#求根函数
运算结果为：
$root
[1] 0.211182

$f.root
[1] -3.844668e-05

$iter
[1] 5

$init.it
[1] NA

$estim.prec
[1] 6.103516e-05

在这里更正一下，上图中求矩估计法和极大似然函数那里两个错误，一个是求均值那里，一个是对数似然函数求导那里，以R软件中的内容为准，仔细看哦，不要误导了大家。
 

2.

 

用R软件来实现：
x=c(rep(5.365),rep(15,245),rep(25,150),rep(35,100),rep(45,70),rep(55,45),rep(65,25))
sum(x)/1000

3.