作者mbaplayer ( #@#,^^)
看板Statistics
标题[程式] 请问一题R的问题
时间Wed Dec 27 23:35:18 2017
[软体程式类别]:
R
[程式问题]:
请问有关区间估计、信赖区间的问题
[软体熟悉度]:
新手
[问题叙述]:
题目如下:
假设有一个超级大的班级，已知他们期末考的分数服从常态分配平均值为65，
标准差为9（上帝视角），用以上资讯模拟出资料并回答下列问题：
（注意set.seed=1,2,3,…,100，不需要画图）
(1)若某学者从班级里随机取25 人，且标准差未知，重复次数为100 次，
请算出所有信赖区间有涵盖真实平均数65的比例。
小弟是R新手，想请问
(1)第8行那个x<-rt(n,n)，自由度是n吗?还是n-1?
(2)这一题的信赖区间的alpha要怎麽用程式求出来啊??
也就是第9行的qt(,n-1)要怎麽填啊?
(3)我的程式码如下，请问还有哪里有错误吗??
跪谢各位大神!!
mu<-65;sigma<-9;n<-25;no.rep<-100
l<-rep(NA,no.rep)
r<-rep(NA,no.rep)
for(i in 1:no.rep){
print(i)
set.seed(i)
x<-rt(n,n)
S^2<-1/(n-1)*sum((x-mean(x))^2)
l<-mean(x)-qt(,n-1)*sqrt(S^2/n)
r<-mean(x)+qt(,n-1)*sqrt(S^2/n)
}
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.164.14.155
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1514388920.A.9EE.html
※ 编辑: mbaplayer (1.164.14.155), 12/27/2017 23:38:41
1F：→ andrew43: 1. df = n-1 2. alpha是人定的，配合1-alpha信赖区间 12/28 00:33
2F：→ andrew43: 3. 直接模拟样本，而不是模拟t-value 12/28 00:36
谢谢!
※ 编辑: mbaplayer (1.164.14.155), 12/28/2017 18:19:40