环境与生态统计||置信区间的估计
经典统计学的核心思想就是用样本（我们所能采集到的数据，如总磷浓度，TP）去估计总体（关于总体情况的描述）；总体的参数是未知的，不可测度或难以测度，注意它是固定了的数值；我们采集了有限的TP，而且在前人研究的基础上已知一个生态系统中TP可以被近似为对数正态分布，那么剩下的工作就是要估计分布的对数均值和对数标准差。估计总体的均值和标准差，最自然的方法就是用样本的均值和标准差去估计。
为了估计这个总体的参数，我们就要通过样本来构造统计量，注意它是一个随机变量。因为我们不能仅用一次采样的数据来估计总体，我们每一次采样都会得到一个均值和标准差。这样这个用于估计的值就是一个随机变量。
随机变量的意思是随着你样本选取的不同，具体到每一个样本的统计量的统计值也不尽相同。这个随机变量的统计值就是对总体参数的点估计，由于样本估计总体总是会存在一定的偏差，所以我们为了更好的估计总体参数，于是用到了置信区间。
95％的置信度的意思是如果你从总体中抽取100个不同样本，每个样本都用相同的统计量构造的置信区间(注意:由于样本不相同，这些置信区间的范围也不尽相同)，那么有95个置信区间包含了总体参数的真值。换句话说，95%的置信区间意味着真值 落在区间内的概率是0.95。
最关键的是要理解统计量是随机变量而总体的参数是一个实实在在的数值。
置信区间的计算在知道方差和不知道方差的情况下，计算公式是不一样的。
如果不知道方差，则confint（x,alpha） 知道方差，则confint（x,sigma,alpha)
按照我的简要理解，置信椭圆基本上是对置信区域的描述方式，其长轴和短轴分别为置信区域的参数，置信椭圆的长短半轴，分别表示二维位置坐标分量的标准差（如经度的σλ和纬度的σφ）。一倍标准差（1σ）的概率值是68.3%，二倍标准差（2σ）的概率值为95.5%；三倍标准差（3σ）的概率值是99.7%。
μ的置信度为1−α的置信区域为
当p=1时，它是一个区间；当p=2时，它是一个椭圆，这时可将其在坐标平面上画出；当p=3时，它是一个椭球；当p＞3时，它是一个超椭球；它们均以 为中心。
同置信区间与假设检验的关系一样，置信区域与假设检验之间也有着同样的密切关系。一般来说，μ0包含在上述置信区域内，当且仅当原假设 H0：μ=μ0在显著性水平α下被接受。因此，可以通过构造置信区域的方法来进行假设检验。
我们用 rnorm(20, 0.2, 1)来构建TP的浓度，计算95% 和50% 的置信区间：
在这个案例中，我们假设TP浓度对数值的真实分布是N(2.05,0.34)，并且利用计算机从这个分布中采取30个随机数来模仿采样过程，然后计算置信区间。当多次重复，我就会期望区间会包含2.05.
以上均是对均值的推断，统计推断还有第二部分，即是对标准差的推断。同样是用样本的标准差来估计总体的标准差。
通过对总体均值和标准差的估计，可以确定模型的参数。但是，Everglades湿地研究背后的问题是要设定TP的环境标准：使用背景浓度的75百分点作为标准。因此接下来的问题是如何估计0.75分位数。如果我们知道总体分布是正态分布的，且均值和标准差的真值是已知的，那么可以直接估计0.75分位数。假设均值=2.05，标准差=0.34.
TP浓度分布的0.75分位数就是 =9.77pp。但是，如何描述0.75分位数的不确定性呢？
参考：
如何理解 95% 置信区间？
用R语言求置信区间
基础知识：置信椭圆confidence ellipse
绛旓細闅忔満鍙橀噺鐨勬剰鎬濇槸闅忕潃浣犳牱鏈夊彇鐨勪笉鍚岋紝鍏蜂綋鍒版瘡涓涓牱鏈殑缁熻閲忕殑缁熻鍊间篃涓嶅敖鐩稿悓銆傝繖涓殢鏈哄彉閲忕殑缁熻鍊煎氨鏄鎬讳綋鍙傛暟鐨勭偣浼拌锛岀敱浜庢牱鏈及璁℃讳綋鎬绘槸浼氬瓨鍦ㄤ竴瀹氱殑鍋忓樊锛屾墍浠ユ垜浠负浜嗘洿濂鐨勪及璁℃讳綋鍙傛暟锛屼簬鏄敤鍒颁簡缃俊鍖洪棿銆95锛呯殑缃俊搴︾殑鎰忔濇槸濡傛灉浣犱粠鎬讳綋涓娊鍙100涓笉鍚屾牱鏈紝姣忎釜鏍锋湰閮界敤...
绛旓細鎺ㄦ柇缁熻鍒嗘瀽锛 鎺ㄦ柇缁熻鍒嗘瀽娑夊強浠庢牱鏈暟鎹腑鎺ㄦ柇鎬讳綋鐗瑰緛銆傝繖鍖呮嫭鍋囪妫楠屻缃俊鍖洪棿浼拌銆佹柟宸垎鏋愮瓑鏂规硶锛岀敤浜庡垽鏂瀵熷埌鐨勫樊寮傛槸鍚﹀叿鏈夌粺璁″鎰忎箟銆傚洖褰掑垎鏋愶細 鍥炲綊鍒嗘瀽鐢ㄤ簬寤虹珛鍙橀噺涔嬮棿鐨勫叧绯绘ā鍨嬨傚湪鐢熺墿瀛︿腑锛屽彲浠ョ敤鍥炲綊鍒嗘瀽鏉ユ帰璁ㄧ敓鐗╁鍙橀噺涔嬮棿鐨勫叧绯伙紝渚嬪鐢熼暱鐜囦笌鐜鍥犵礌涔嬮棿鐨勫叧绯汇傜敓瀛樺垎鏋愶細 鐢熷瓨...
绛旓細绗竴,浠庡悇绫诲尯鐨勫湡鍦伴潰绉強鍒嗗竷鏉ョ湅,鍦熷湴瀹夊叏鍖虹浉褰撲簬鍦熷￥鐜璐ㄩ噺鐨勨厾绫诲拰鈪＄被鍖,鍩烘湰瀹夊叏鍖虹浉褰撲簬鈪㈢被鍖,璀︽垝鍖哄拰鍗遍櫓鍖虹浉褰撲簬瓒呪參绫诲尯銆備粎浠庤繖涓鐐规潵璇,鐩墠璇勪环宸ヤ綔涓,灏嗙敤鍦熷￥鐜璐ㄩ噺鏍囧噯璇勪环寰楀嚭鐨勨厾绫汇佲叀绫汇佲參绫汇佽秴鈪㈢被鐨勫湡澹ょ幆澧冭川閲忓垎绾,鍒嗗埆瀹氫箟涓烘竻娲併佽交搴︽薄鏌撱佷腑搴︽薄鏌撱侀噸搴︽薄鏌,涓嶇鍚堝疄闄...
绛旓細搴旂敤缁熻纭曞＋涓撲笟瀛︿綅灞炰簬缁忔祹瀛︾澹傚湪搴旂敤缁熻瀛︾殑瀛︿範涓紝瀛︾敓灏嗗涔犵粺璁″鐨勫熀鏈悊璁恒佹蹇靛拰鏂规硶锛屽寘鎷暟鎹噰闆嗐佺粺璁℃帹鏂佸洖褰掑垎鏋愩佺粺璁″缓妯＄瓑鍐呭銆備粬浠篃浼氬涔犲浣曚娇鐢ㄨ绠楁満鍜岀粺璁杞欢鏉ュ鐞嗗拰鍒嗘瀽澶ч噺鐨勬暟鎹傚簲鐢ㄧ粺璁″閫氬父娑夊強鍒板疄璇佺爺绌躲佹暟鎹垎鏋愩侀娴嬪拰鍐崇瓥绛夋柟闈紝鍥犳鍏峰鑹ソ鐨勬暟瀛﹀熀纭鍜...
绛旓細1 鏈変竴澶ф壒绯栨灉,璁炬瘡琚嬮噸閲忔湇浠庢爣鍑嗗樊涓6鍏嬬殑姝ｆ佸垎甯,鐜颁粠涓殢鏈烘娊鍙栦簡16琚,缁忕О閲嶅緱鐭,骞冲潎姣忚閲嶉噺涓500鍏,璇曚互95%鐨缃俊姘村钩浼璁¤鎵圭硸鏋滃钩鍧囬噸閲忕殑缃俊鍖洪棿銆2 鏌愪紒涓氬墠涓ゅ勾鍚勫搴﹂攢鍞噺璧勬枡濡備笅(鍗曚綅:涓囧厓): 骞翠唤 涓瀛ｅ害 浜屽搴 涓夊搴 鍥涘搴 0 5 4 2 2 5 0 6 6 4 6 13 瑕佹眰: (1)閲囩敤...
绛旓細鎵浠,涓瀹氳鍗冩柟鐧捐鍦板缓璁句汉姘村拰璋愮殑缁胯壊姘村伐绋,鍧氬喅浠庢簮澶存潨缁濅笉鑹」鐩殑涓婇┈寤鸿,杩欏氨闇瑕佸鎷熷缓鐨勬按宸ョ▼椤圭洰杩涜鍏呭垎璁鸿瘉銆佽鐪熸瘮閫夈佺簿蹇冭璁;鍧氭寔缁忔祹鏁堢泭銆佺ぞ浼氭晥鐩涓庣敓鎬鏁堢泭骞堕噸;鍧氭寔浠ョ瀛﹀彂灞曡鐨勮姹,搴﹂噺鍜屾妸鎻″伐绋嬬洰鏍囧疄鐜扮殑缃俊鍖洪棿鍜屽亸宸 閮芥睙鍫版槸涓崕姘戞棌绉戝鏅烘収涓庢枃鏄庣殑缁撴櫠,閮芥睙鍫扮殑寤鸿鍜岃繍鐢ㄦ墍...
绛旓細绯荤粺浠跨湡鐨勬蹇垫槸鏍规嵁绯荤粺鍒嗘瀽鐨勭洰鐨勶紝鍦ㄥ垎鏋愮郴缁熷悇瑕佺礌鎬ц川鍙婂叾鐩镐簰鍏崇郴鐨勫熀纭涓婏紝寤虹珛鑳芥弿杩扮郴缁熺粨鏋勬垨琛屼负杩囩▼鐨勩佷笖鍏锋湁涓瀹氶昏緫鍏崇郴鎴栨暟閲忓叧绯荤殑浠跨湡妯″瀷锛屾嵁姝よ繘琛岃瘯楠屾垨瀹氶噺鍒嗘瀽锛屼互鑾峰緱姝ｇ‘鍐崇瓥鎵闇鐨勫悇绉嶄俊鎭傝绠楁満璇曢獙甯歌鐢ㄦ潵鐮旂┒浠跨湡妯″瀷锛坰imulation model锛夈備豢鐪熶篃琚敤浜庡鑷劧绯荤粺鎴栦汉閫犵郴缁熺殑...
绛旓細搴旂敤缁熻瀛︾殑搴旂敤棰嗗煙闈炲父骞挎硾锛屽寘鎷絾涓嶉檺浜庡競鍦鸿皟鐮斻侀噾铻嶉闄╄瘎浼般佸尰瀛︾爺绌躲鐢熸佺幆澧鍒嗘瀽銆佺ぞ浼氳皟鏌ャ佽川閲忔帶鍒躲佹斂搴滄斂绛栬瘎浼扮瓑銆備笓涓氭瘯涓氱敓閫氬父鑳藉鍦ㄤ紒涓氥佹斂搴滈儴闂ㄣ佺爺绌舵満鏋勩佸尰鐤楁満鏋勭瓑棰嗗煙浠庝簨鏁版嵁鍒嗘瀽銆佺粺璁″缓妯°侀闄╄瘎浼般佸喅绛栨敮鎸佺瓑宸ヤ綔銆傛讳箣锛屽簲鐢ㄧ粺璁″鏄皢缁熻瀛﹀簲鐢ㄤ簬瀹為檯闂鐨勪竴闂ㄤ笓涓氥傚畠...
绛旓細4.鍔姏寤鸿浜烘按鍜岃皭銆佺敓鎬佸弸濂界殑缁胯壊姘村伐绋嬨 姘村伐绋嬬殑寤鸿,鐗瑰埆鏄ぇ鍨嬫按宸ョ▼鐨勫缓璁,鍏充箮缁忔祹绀句細鍙戝睍銆鐢熸佺幆澧鏀瑰杽浠ュ強浜轰滑鐢熸椿鍜鐢熷瓨鏂瑰紡鐨勫彉鍖栥備互寰鏈夌殑姘村伐绋嬪缓璁,铏界粡杩囦簡鍓嶆湡宸ヤ綔鈥滃崈閿ょ櫨鐐尖濈殑鎵撻,浣嗗缓鎴愬悗浠嶆棫鏄晥鐩婁笉楂樸佺幆澧冨奖鍝嶇獊鍑恒佽繍琛岃竟闄呮垚鏈繃澶с傜幇鍦ㄧ湅,杩欐牱涓浜涙按宸ョ▼鐨勪慨寤,鍙槸浣挎棦寰楀埄鐩婅...
