對於統計製圖有興趣的人，絕對不能錯過。
吃早餐時，滑手機的當下看到這篇文章。
回到家後「立刻」開電腦，把文章下載下來(我一定要備份，這是我日後個人統計著作講到error bar如何下語法畫出時，可以拿出來講的負面案例)，按我閱讀文章(圖9)。
刊登此文章的期刊
文章於官網被警示(顯示目前正調查及驗證中，它隨時會被下架)
期刊文章官網(隨時會被下架致連結失效)
來看看大作，也就是圖9
聚焦重點部位，Ｔ竟然站在高崗上！(有些T甚至跳起來了)
然後每個T的身高都一樣，這也太巧合了。
然後灰色高崗上的T似乎不在置中位置上。
然後然後，這篇文章其它圖也很可疑，但我就不再點出了。
這種事後加工(小畫家程式？)行為本來就不被允許(用人工拉的無法精準到位)，這就是統計圖造假，作者八成得罪了統計圖繪圖專家，所以吃了閉門羹。
輕鬆點 ^___^ 來唱首歌！站在高崗上！會唱的一起唱。
"連綿的青山百里長呀，蘶蘶聳起像屏障呀喂"
如果資料量不大(通常是千筆以下)，建議不要用這種長條圖+這種符號(姑且叫它T)，應使用散布(佈？)圖和盒形圖的疊圖，見這篇文章。
這種T符號通常可用來描述信賴區間(較常以森林圖方式呈現)、標準差或標準誤。它可以正負很大(即大區間)，也可以正負很小(即小區間)。白話講，就是給你的估計值留點餘地，話不要講得太滿，說多少就是多少，而是有一個合理的估計區間(因為人畢竟不是神)。
SAS在這裡有對應的語法可讓你使用（見官網），例如要畫垂直長條圖，就須用到 Vbar Statement．
可以指定一個變項，這個變項有100個數值(僅舉例)，而Vbar要長多高，端看這100個數值的何描述性統計值，你下SUM，SAS就用總和呈現其高度，你下MEDIAN，SAS就用中位數呈現高度，不須額外再計算總和值或中位數值各別寄放在不同變數中待提取。最後，再分別用Limitstat和Limits控制上述的T符號(笑)。如果想呈現該論文應有圖形，下Limitstat=Stderr Limits=Upper語法就能達成。同理，若用使STDERR，SAS一樣會自動計算標準誤，再調整T身高(好一個淺顯易懂的名詞)至對應的刻度，不須額外計算。
下面兩張圖出自SAS官網文件，已清楚說明 Limitstat option和Limits option用法。
下面這張SAS統計圖，是我多年前畫的。若要同時呈現站T和倒T，下Limits=both語法即可達成。除了p值為事後加工打上去，其餘部份全由SAS語法產出。
至於有個財團法人中X發展中心常辦的品管圈大賽，說什麼繪圖時，軸項上的刻度線(具體的中文名稱我不知叫什麼，但英文名稱肯定叫tick，有major和minor之分)要朝內不能朝外不然會扣分，這樣的規定不覺得很無理取鬧嗎？SAS繪圖語法 Xaxis statement中，有朝內或朝外的語法，也可以不要有tick，但我就是不朝內，故意朝外跟中X作對！
沒有留言:
張貼留言