在正文之前先分享一點生活故事：
上周天傍晚，我閃著腰了邑闲，躺了一晚上不見好算行，發(fā)朋友圈尋醫(yī)問藥了：
因為做的是生信分享，我朋友圈里醫(yī)學(xué)生和醫(yī)生非常多苫耸，我想看他們怎么說州邢。
我收到了幾十條回復(fù)，分享一下給大家鲸阔，就當拿我的教訓(xùn)給大家漲漲經(jīng)驗吧偷霉，有備無患...
1.關(guān)于為什么痛
根據(jù)我的癥狀描述，多位醫(yī)生/醫(yī)學(xué)生指出這是腰肌勞損褐筛，多半是因為久坐姿勢不對引起的类少，如果嚴重一點就會發(fā)展成為小關(guān)節(jié)紊亂，再不注意渔扎，過幾年會發(fā)展成腰椎間盤突出(啊嚇死了硫狞。)
2.關(guān)于去不去醫(yī)院
大多數(shù)醫(yī)生/醫(yī)學(xué)生說得去醫(yī)院，也有幾位說不用去醫(yī)院晃痴，在家躺著休息残吩，感覺有麻木、牽拉感再去醫(yī)院倘核。
3.關(guān)于掛什么科
先要去骨科或者脊柱外科泣侮，要拍個片子確認一下沒有傷到骨頭，也有人說做MR（核磁共振成像）紧唱；然后去針灸/推拿/理療/疼痛/康復(fù)科活尊，針灸和推拿都是有效的，有針灸科的醫(yī)學(xué)生說漏益，這個一針見效蛹锰。不能針灸的話，那就是中醫(yī)正骨推拿绰疤，也能搞定铜犬，不過大家都說推拿不是誰都可以，要有經(jīng)驗的老中醫(yī)才行轻庆。
關(guān)于拍片子癣猾，有朋友指出我這個是彎腰岔氣，不是摔的或者碰的余爆，沒有傷到骨頭煎谍，不拍也行，拍了放心龙屉。
4.病例本例現(xiàn)身說法
知道了推拿可以治這回事呐粘，我查了一下家長，附近有個中醫(yī)推拿館转捕，我就打了電話過去作岖，他們詳細詢問了我的癥狀，問是不是久坐五芝，以前有沒有出現(xiàn)過輕的腰部不適癥狀痘儡，我說是的。她說可以過去枢步，推拿能搞定沉删。最后沒有拍片子渐尿，一個小時的推拿基本緩解了疼痛，本來無法彎腰矾瑰，坐下就痛砖茸，按完可以彎腰了，還有一些不適的感覺殴穴，但是明顯緩解凉夯。
發(fā)現(xiàn)了沒，我的朋友圈可太優(yōu)秀了吧采幌！劇透一波劲够，雙十一我要在公眾號為我朋友圈里的單身男女們發(fā)布一下征友信息，敬請關(guān)注哈哈
以下是正文
1.p值是不是概率休傍？
P value is the probility that rendom chance generated the data, or something else that is equal or rarer.
也就是說征绎，p值由三部分概率組成：目標事件概率+(與目標事件概率相同的事件)的概率+（比目標事件概率更小的事件）的概率。
順便補充兩個R語言的函數(shù)磨取，choose和combn炒瘸。
choose(5,1)
## [1] 5
choose(5,2)
## [1] 10
choose(5,3)
## [1] 10
choose(5,4)
## [1] 5
choose(5,5)
## [1] 1
choose(5,3)表示從五個數(shù)中抽取三個，有多少種可能寝衫，返回具體數(shù)值顷扩。
combn(1:5,1)
## [,1] [,2] [,3] [,4] [,5]
## [1,] 1 2 3 4 5
combn(1:5,2)
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
## [1,] 1 1 1 1 2 2 2 3 3 4
## [2,] 2 3 4 5 3 4 5 4 5 5
combn(1:5,3)
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
## [1,] 1 1 1 1 1 1 2 2 2 3
## [2,] 2 2 2 3 3 4 3 3 4 4
## [3,] 3 4 5 4 5 5 4 5 5 5
combn(1:5,4)
## [,1] [,2] [,3] [,4] [,5]
## [1,] 1 1 1 1 2
## [2,] 2 2 2 3 3
## [3,] 3 3 4 4 4
## [4,] 4 5 5 5 5
combn(1:5,5)
## [,1]
## [1,] 1
## [2,] 2
## [3,] 3
## [4,] 4
## [5,] 5
combn(1:5,3)表示從五個數(shù)中抽取3個，列出所有可能(矩陣形式)
2.計算實例一慰毅，扔兩個硬幣
結(jié)果為HH的概率和p值
（H：head隘截，T：tail）
隨機拋兩次硬幣，HH汹胃、TT的概率都是1/4婶芭，HT的概率是1/2。
根據(jù)三個組成部分着饥，可以計算出HH的p值：1/4 + 1/4 + 0 = 1/2 = 0.5
結(jié)果為HT的概率和p值
順序通常沒什么用犀农，在這里，不考慮順序宰掉，HT和TH是一樣的呵哨。概率為1/2。
p值是 1/2 + 1/4 + 1/4 = 1
3.計算實例二轨奄，扔五個硬幣
五個硬幣孟害，總共有32種可能:
除了五正、五反兩種情況挪拟，其他的30種是：
library(stringr)
lapply(1:4, function(n){
k = combn(1:5,n)
sapply(1:ncol(k), function(x){
y = rep("H",5)
y[k[,x]] = "T"
return(paste(y,collapse = ""))
})
})
## [[1]]
## [1] "THHHH" "HTHHH" "HHTHH" "HHHTH" "HHHHT"
##
## [[2]]
## [1] "TTHHH" "THTHH" "THHTH" "THHHT" "HTTHH" "HTHTH" "HTHHT" "HHTTH" "HHTHT"
## [10] "HHHTT"
##
## [[3]]
## [1] "TTTHH" "TTHTH" "TTHHT" "THTTH" "THTHT" "THHTT" "HTTTH" "HTTHT" "HTHTT"
## [10] "HHTTT"
##
## [[4]]
## [1] "TTTTH" "TTTHT" "TTHTT" "THTTT" "HTTTT"
可以看到上面生成的結(jié)果是一個列表挨务，分了四組，分別代表有1、2谎柄、3丁侄、4個T的所有情況。
五個正面朝上的概率和p值
HHHHH的概率：(1/2)^5 = 1/32朝巫，TTTTT與之概率相等鸿摇，不存在比1/32概率更小的事件。
HHHHH的p值：1/32 + 1/32 + 0 = 0.0625
雖然五個H的幾率很小捍歪，但是p值>0.05，不認為是非常罕見的事件鸵钝。
四正一反的概率和p值
四正一反的概率：choose(5,1)/(2^5) = 5/32
與目標事件相等的事件是四反一正糙臼，比目標事件概率更小的是五正和五反，所以p值是：
5/32 + 5/32 + 1/32 + 1/32 = 0.375
4.連續(xù)型的數(shù)據(jù)如何計算p值
例如身高恩商，無法窮舉所有的可能性变逃。
密度圖曲線下的面積代表橫坐標對應(yīng)范圍事件發(fā)生的概率
生成示例圖，身高的正態(tài)分布曲線怠堪，均值約為155.5
library(ggplot2)
x = seq(130,180,0.1)
df = data.frame(x = x,
y = dnorm(x,mean = 155.5,sd = 7),
gro = ifelse(x< 142,"A",ifelse(x >= 169,"C","B")))
ggplot(df, aes(x = x,y = y)) +
geom_line() + theme_bw()
R語言如何計算曲線下的面積揽乱？我找到了一個積分函數(shù)，奇怪的知識增加了
連續(xù)型數(shù)據(jù)的概率
小于142cm:
integrate(function(x){dnorm(x=x,mean = 155.5,sd = 7)}, -Inf, 142)
## 0.02689204 with absolute error < 8.3e-06
大于169cm:
integrate(function(x){dnorm(x=x,mean = 155.5,sd = 7)}, 169, Inf)
## 0.02689204 with absolute error < 8.3e-06
介于142~169cm之間
integrate(function(x){dnorm(x=x,mean = 155.5,sd = 7)} , 142 , 169)
## 0.9462159 with absolute error < 6.7e-12
R語言是個好玩具粟矿，還可以畫出相應(yīng)的圖（面積圖）
ggplot(df, aes(x = x)) +
geom_line(aes(y = y)) +
geom_area(data = df[df$gro == "A",],aes(y = y), fill = 'blue',color = NA,alpha = 0.5)+
geom_area(data = df[df$gro == "B",],aes(y = y), fill = 'red',color = NA,alpha = 0.5)+
geom_area(data = df[df$gro == "C",],aes(y = y), fill = 'blue',color = NA,alpha = 0.5)+
theme_bw()
重申一下凰棉，面積就是概率啊。
連續(xù)型數(shù)據(jù)的p值
再來看p值陌粹，身高小于142的p值如何計算撒犀？
身高>169和<142概率相等，沒有比<142更加小概率的事件掏秩。所以p值為0.027*2 = 0.054
身高為142~169之間的p值：
0.946 + 0.027*2 = 1
0.946 是目標事件的概率或舞，兩個0.027是(比目標事件概率更小的事件)的概率。
區(qū)分p值和概率的極端例子
身高介于155~156之間的概率是：
integrate(function(x){dnorm(x=x,mean = 155.5,sd = 7)} , 155 , 156)
## 0.05694333 with absolute error < 6.3e-16
面積圖：
df2 = data.frame(x = x,
y = dnorm(x,mean = 155.5,sd = 7),
gro = ifelse(x< 155,"A",ifelse(x >= 156,"C","B")))
ggplot(df2, aes(x = x)) +
geom_line(aes(y = y)) +
geom_area(data = df2[df2$gro == "B",],aes(y = y), fill = 'red',color = NA,alpha = 0.5)+
theme_bw()
概率只為0.057蒙幻，但p值呢?
目標事件概率是0.057映凳，比目標事件概率更小的事件概率？因為155~156 處在正態(tài)分布的中位數(shù)附近邮破，所以他的概率是最大的诈豌，身高落在所有的其他范圍內(nèi)的概率都比155~156要低，所以p值等于0.057+（1-0.057）=1
