R语言实现正太总体下的区间估计_天财研院1410班团支部

http://blog.sina.com.cn/u/2033416165

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

R语言实现正太总体下的区间估计

(2014-12-13 21:32:38)

标签：

it

一、均值的区间估计的R程序

当总体方差已知的时候公式

http://s7/bmiddle/002dC0qVty6OmcqHCPs56&690

当总体方差未知的时候公式

http://s4/mw690/002dC0qVty6OmcrfD1N43&690

代码如下：

interval_estimate1<-funct ion (x, sigma=-1, alpha=0.05){

n<-length(x); xb<-mean(x)

if (sigma>=0){

tmp<-sigma/sqrt(n)*qnorm(1 -alp ha/2 ); df<-n

}

else{

tmp<-sd(x)/sqrt(n)*qt(1-al pha/ 2,n- 1); df<-n-1

}

data.frame(mean=xb, df=df, a=xb-tmp, b=xb+tmp)

}

另外,R软件的t.test(x)检验函数可以完成相应的区间估计任务。

二、方差的区间估计

当均值已知时公式为

http://s15/mw690/002dC0qVty6OmcQVv7Ude&690

当均值未知时公式为
 http://s7/mw690/002dC0qVty6OmcQX94Ob6&690

根据这两种情况，方差的区间估计R程序

i2<-function(x,mu=inf,alpha=0.05){

n<-length(x)

if(mu

s2<-sum((x-mu)^2)/n;df<-n

}

else{

s2<-var(x);df<-n-1

}

a<-df*s2/qchisq(1-alpha/2,df)

b<-df*s2/qchisq(alpha/2,df)

data.frame(var=s2,df=df,a=a,b=b)

}

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：关于dnorm(),pnorm(),qnorm(),rnorm()的解释

后一篇：两正态总体的均值比较

新浪BLOG意见反馈留言板　欢迎批评指正