估計

統計軟體 R

簡介

安裝

操作方式

變數與運算

有序數列

向量

矩陣

多維陣列

複數

因子

串列

資料框

時間數列

流程控制

輸出入

呼叫

函數

2D 繪圖

3D 繪圖

互動介面

套件列表

其他語言呼叫

R 的應用

集合

邏輯推論

模糊邏輯

機率邏輯

檢定

搜尋

優化算法

線性代數

決策樹

人工智慧

分群分類

SVM 向量機

神經網路

遺傳演算法

資料採礦

訊號處理

影像處理

語音處理

自然語言

機器學習

機器人

生物統計

數位訊號處理

方程式求解

數值分析

微積分

微分方程

線性規劃

圖形理論

統計推論

字串處理

正規表示式

視窗程式

網頁程式

文件格式

貝氏網路

訊息

機率統計書

相關網站

參考文獻

最新修改

簡體版

English

均值估計

點估計:

# x 为均值5, 方差1 的总体中抽取的10 个样本
> x=rnorm(10,5)
> x
[1] 4.927264 4.067237 6.136822 5.722123 6.286754 3.266601 4.443779 3.630787
[9] 4.874269 3.748306
# z 值为qnorm(0.025)=-1.959964, qnorm(0.975)=1.959964
> mean(x)+qnorm(0.025)*1/sqrt(10)
[1] 4.090599
> mean(x)+qnorm(0.975)*1/sqrt(10)
[1] 5.330189

區間估計

> 3+qt(p=0.025,df=20)*5/sqrt(20)
[1] 0.667822
> 3+qt(p=0.975,df=20)*5/sqrt(20)
[1] 5.332178

區間估計:總體方差已知

# x 为均值5, 方差1 的总体中抽取的10 个样本
> x=rnorm(10,5)
> x
[1] 4.927264 4.067237 6.136822 5.722123 6.286754 3.266601 4.443779 3.630787
[9] 4.874269 3.748306
# z 值为qnorm(0.025)=-1.959964, qnorm(0.975)=1.959964
> mean(x)+qnorm(0.025)*1/sqrt(10)
[1] 4.090599
> mean(x)+qnorm(0.975)*1/sqrt(10)
[1] 5.330189

區間估計:總體方差未知

> 3+qt(p=0.025,df=20)*5/sqrt(20)
[1] 0.667822
> 3+qt(p=0.975,df=20)*5/sqrt(20)
[1] 5.332178

方差估计 (變異數估計)

點估計:

> x
[1] -5 -4 -3 -2 -1 0 1 2 3 4 5
> var(x)
[1] 11
> sum((x-mean(x))^2)/(length(x)-1)
[1] 11

區間估計

> x
[1] -5 -4 -3 -2 -1 0 1 2 3 4 5
> (10-1)*var(x)/qchisq(0.025,10-1)
[1] 36.66138
> (10-1)*var(x)/qchisq(0.975,10-1)
[1] 5.20429

二項分布的估計

参数p及标准误差的点估计

> x=rbinom(10,1,0.5)
> x
[1] 1 1 0 1 1 1 0 0 1 0
> t=table(x)
> t
x
0 1
4 6
> t['1']/length(x) # 此即为p的点估计, 还可以使用binom.test(table(x))得到.
1
0.6
> sqrt(t['1']*t['0']/length(x)) # 此为标准误差的点估计
1
1.549193

p的区间估计

> binom.test(table(x))
Exact binomial test
data: table(x)
number of successes = 4, number of trials = 10, p-value = 0.7539
alternative hypothesis: true probability of success is not equal to 0.5
95 percent confidence interval:
0.1215523 0.7376219
sample estimates:
probability of success
0.4
> b=binom.test(table(x))
> str(b)
List of 9
$ statistic : Named int 4
..- attr(*, "names")= chr "number of successes"
$ parameter : Named int 10
..- attr(*, "names")= chr "number of trials"
$ p.value : Named num 0.754
..- attr(*, "names")= chr "0"
$ conf.int : atomic [1:2] 0.122 0.738
..- attr(*, "conf.level")= num 0.95
$ estimate : Named num 0.4
..- attr(*, "names")= chr "probability of success"
$ null.value : Named num 0.5
..- attr(*, "names")= chr "probability of success"
$ alternative: chr "two.sided"
$ method : chr "Exact binomial test"
$ data.name : chr "table(x)"
- attr(*, "class")= chr "htest"

Unless otherwise stated, the content of this page is licensed under Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License