作者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/m.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