求助！！！！求置信区间程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

R或者sas中，怎样求得联合置信区间，请各位高手指点！指点

在此谢过！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：置信区间 高手指点 求助 程序 置信区间

xbar<-c(95.52,164.38,55.69,93.39,17.98,31.13)

  s<-matrix(c(3266.46,  1343.97, 731.54,  1175.50, 162.68, 238.37,  1343.97,  721.91, 324.25, 537.35, 80.17, 117.73,  731.54,  324.25,  179.28, 281.17, 39.15, 56.8,  1175.50,  537.35,  281.17,  474.98,  63.73, 94.85,  162.68,  80.17,  39.15,  63.73,  9.95, 13.88, 238.37,  117.73,  56.80,  94.85,  13.88,  21.26),6,6,byrow=T)

# part(a)

n<-61

p<-6

c<-sqrt(p*(n-1)/n/(n-p)*qf(0.95,p,n-p))

for (i in 1:6) {

  lower<-xbar-c*sqrt(s[i,i])

  upper<-xbar+c*sqrt(s[i,i])

  cat("(",lower,",",upper,")","\n")

}