[版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規]
1.不回覆私底下email的問題。有問題的請在此留言。
2.不保證能夠回答所有問題。
3.不協助debug。
4.SAS軟體本身的問題請自行聯繫SAS客服部門。
[版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規][版規]
2015年9月12日 星期六
A SAS® Macro for Computing Point Estimates and Confidence Intervals of Effect Sizes Associated with Mediation Analysis
這篇技術文件主要是在介紹在進行中介分析(mediator analysis)時，如何用拔靴法(bootstrap)計算effect size及其信賴區間。作者並提供了一個巨集程式供其他使用者方便使用。首先簡單描述一下中介效應。
2015年9月7日 星期一
A Generalized Approach to Estimating Sample Sizes
原文載點：http://support.sas.com/resources/papers/proceedings12/336-2012.pdf
幾年前曾經寫過幾篇關於介紹PROC POWER和PROC GLMPOWER功能的技術文件。本篇特地再引述一篇計算樣本的技術文件的用意是這一篇特別介紹了如何用模擬的方式來做檢定力分析(Power Analysis)，尤其當模型中涉及到交互作用項時。想看原文的可以直接跳到原文的第四頁來看。
幾年前曾經寫過幾篇關於介紹PROC POWER和PROC GLMPOWER功能的技術文件。本篇特地再引述一篇計算樣本的技術文件的用意是這一篇特別介紹了如何用模擬的方式來做檢定力分析(Power Analysis)，尤其當模型中涉及到交互作用項時。想看原文的可以直接跳到原文的第四頁來看。
2015年8月30日 星期日
Colon(:)izing My Programs
原文載點：http://support.sas.com/resources/papers/proceedings15/3367-2015.pdf
冒號( : )在SAS的資料處理中扮演了一個鮮為人知的角色。這一篇2015年在SAS Global Forum發表的技術文件介紹了如何適當的使用冒號在SAS程式裡面。我列了幾個一般使用者比較會用到的功能。
冒號( : )在SAS的資料處理中扮演了一個鮮為人知的角色。這一篇2015年在SAS Global Forum發表的技術文件介紹了如何適當的使用冒號在SAS程式裡面。我列了幾個一般使用者比較會用到的功能。
2015年8月24日 星期一
A SAS Macro to Compute Effect Size (Cohen’s d) and its Confidence Interval from Raw Survey Data
原文載點：http://analytics.ncsu.edu/sesug/2012/SD-06.pdf
Effect size是一種量化兩個或多個變數間關係的測度。在許多臨床試驗裡面都需要這個數據。不同的研究設計需要不同種類的Effect size，其中又以Cohen's d統計量最為流行，因為他是專門用在下面這三種的研究：(1)兩個獨立群組的試後比較; (2)一個群組內的試前試後比較; (3)兩個獨立群組的試前試後差異比較。由於SAS沒有內建任何的程序來計算Effect size，所以通常需要去找線上的計算器或者自己用Excel試算表去計算。但其實Cohen's d的公式相當簡單，裡面只有用到三種元素：各組樣本，各組平均數，各組標準差。所以我在算Cohen's d時都是自己先用PROC MEANS去算出那三個元素，然後再用PROC IML去算。2012年的SESUG發表了一篇技術文件，由兩個南佛大的人公布了一個計算Cohen's d的SAS巨集程式。
Effect size是一種量化兩個或多個變數間關係的測度。在許多臨床試驗裡面都需要這個數據。不同的研究設計需要不同種類的Effect size，其中又以Cohen's d統計量最為流行，因為他是專門用在下面這三種的研究：(1)兩個獨立群組的試後比較; (2)一個群組內的試前試後比較; (3)兩個獨立群組的試前試後差異比較。由於SAS沒有內建任何的程序來計算Effect size，所以通常需要去找線上的計算器或者自己用Excel試算表去計算。但其實Cohen's d的公式相當簡單，裡面只有用到三種元素：各組樣本，各組平均數，各組標準差。所以我在算Cohen's d時都是自己先用PROC MEANS去算出那三個元素，然後再用PROC IML去算。2012年的SESUG發表了一篇技術文件，由兩個南佛大的人公布了一個計算Cohen's d的SAS巨集程式。
2014年12月18日 星期四
Generalized McNemar’s Test for Homogeneity of the Marginal Distributions
原文載點：http://www2.sas.com/proceedings/forum2008/382-2008.pdf
在配對分析裡面，McNemar's test一直是主要被拿來使用的統計分析，但這個方法有個先天上的限制是只支援二項資料。一旦資料成為多項資料時，就必須使用Generalized McNemar's test (或稱Stuart-Maxwell test)或者是Bhapkar's test。其中，Bhapkar's test是有被SAS/STAT內的程序所支援，但Generalized McNemar's test則無。Xuezheng Sun和Zhao Yang在2008年發表了一篇技術文件裡面，敘述了如何使用既有的SAS程序來執行Bhapkar's test，並且公佈了一個自製的巨集程式來執行Generalized McNemar's test。
在配對分析裡面，McNemar's test一直是主要被拿來使用的統計分析，但這個方法有個先天上的限制是只支援二項資料。一旦資料成為多項資料時，就必須使用Generalized McNemar's test (或稱Stuart-Maxwell test)或者是Bhapkar's test。其中，Bhapkar's test是有被SAS/STAT內的程序所支援，但Generalized McNemar's test則無。Xuezheng Sun和Zhao Yang在2008年發表了一篇技術文件裡面，敘述了如何使用既有的SAS程序來執行Bhapkar's test，並且公佈了一個自製的巨集程式來執行Generalized McNemar's test。
2013年7月22日 星期一
Creating ZIP Files with ODS
原文載點：http://support.sas.com/resources/papers/proceedings13/131-2013.pdf
SAS雖然可以處理龐大的資料，可是有時候過於龐大的資料沒有辦法用電子郵件寄送，直接放在網路硬碟又很佔空間，上傳下載也花時間，所以能夠壓縮成ZIP檔相信是許多人會做的動作。一般人都是額外安裝winzip, winrar或7zip來額外處理壓縮檔案的事情，但其實在SAS內部整理完資料後就可以直接進行壓縮。以下這篇發表於SAS Global Forum 2013的技術文件就教導了如何用SAS ODS程序來製作ZIP壓縮檔案。
SAS雖然可以處理龐大的資料，可是有時候過於龐大的資料沒有辦法用電子郵件寄送，直接放在網路硬碟又很佔空間，上傳下載也花時間，所以能夠壓縮成ZIP檔相信是許多人會做的動作。一般人都是額外安裝winzip, winrar或7zip來額外處理壓縮檔案的事情，但其實在SAS內部整理完資料後就可以直接進行壓縮。以下這篇發表於SAS Global Forum 2013的技術文件就教導了如何用SAS ODS程序來製作ZIP壓縮檔案。
2013年7月19日 星期五
CREATING MACRO VARIABLES VIA PROC SQL
文章來源：http://www.ats.ucla.edu/stat/sas/library/nesug99/cc107.pdf
在 SAS 裡面製作巨集變數，常用的方法不外乎是用 %let 或者是 %symput 來宣告，但功能就僅僅於此了。不過，若使用 PROC SQL 程序來宣告巨集變數的話，可以做些變化來提升程式寫作效率。
2012年5月5日 星期六
Everything in Its Place: Efficient Geostatistical Analysis with SAS/STAT ® Spatial Procedures
原文載點：http://support.sas.com/resources/papers/proceedings10/337-2010.pdf
SAS 在處理空間統計和空間模型上一向相當薄弱，根據這幾年全球最大的統計會議 Joint Statistical Meeting 所顯示，越來越多人在做空間統計分析（我自己也是）。所以特別利用這一篇文章來介紹一份 2010 年發表的 SAS 技術文件，來介紹在 SAS 裡面到底開怎麼樣進行空間統計分析。
SAS 在處理空間統計和空間模型上一向相當薄弱，根據這幾年全球最大的統計會議 Joint Statistical Meeting 所顯示，越來越多人在做空間統計分析（我自己也是）。所以特別利用這一篇文章來介紹一份 2010 年發表的 SAS 技術文件，來介紹在 SAS 裡面到底開怎麼樣進行空間統計分析。
2012年4月27日 星期五
ODS Graphics Editor
原文載點：http://www2.sas.com/proceedings/forum2008/235-2008.pdf
SAS 雖然在作圖上比較複雜，但有一些利用ODS直接生成的圖形是還真的不錯看。不過有些細微的調整可能就還要再用Photoshop或其他圖形編輯軟體去改，不過SAS/GRAPH提供的ODS Graphics Editor已經可以讓使用者直接在圖上修改，裡面提供了許多統計圖形可能會用到的功能，的確比直接用程式去寫快了許多。這一篇發表在 SAS Global Forum 2008 的技術文件就對這個SAS的圖形編輯器做了詳盡的介紹。
SAS 雖然在作圖上比較複雜，但有一些利用ODS直接生成的圖形是還真的不錯看。不過有些細微的調整可能就還要再用Photoshop或其他圖形編輯軟體去改，不過SAS/GRAPH提供的ODS Graphics Editor已經可以讓使用者直接在圖上修改，裡面提供了許多統計圖形可能會用到的功能，的確比直接用程式去寫快了許多。這一篇發表在 SAS Global Forum 2008 的技術文件就對這個SAS的圖形編輯器做了詳盡的介紹。
2012年3月31日 星期六
Computation of Correlation Coefficient and Its Confidence Interval in SAS®
原文載點：http://www2.sas.com/proceedings/sugi31/170-31.pdf
在一些特殊的情況下，相關係數需要搭配信賴區間。PROC CORR程序不像其他程序語法一樣在後面的選項加上一個 CL 就可以自動生出信賴區間，而是需要用別的選項語法才能產生。而且，這個功能是只有 SAS v9.1.3 之後的版本才有，舊版本則需要靠使用者自己寫程式去算。David Shen 和 Zaizai Lu 在 SUGI31 發表了一篇技術文件來解說如何在新舊版本下算出相關係數的信賴區間。
在一些特殊的情況下，相關係數需要搭配信賴區間。PROC CORR程序不像其他程序語法一樣在後面的選項加上一個 CL 就可以自動生出信賴區間，而是需要用別的選項語法才能產生。而且，這個功能是只有 SAS v9.1.3 之後的版本才有，舊版本則需要靠使用者自己寫程式去算。David Shen 和 Zaizai Lu 在 SUGI31 發表了一篇技術文件來解說如何在新舊版本下算出相關係數的信賴區間。
2012年3月24日 星期六
Mixed Model Influence Diagnostics
原文載點：http://www2.sas.com/proceedings/sugi29/189-29.pdf
在進行混合模型的估計中，模型檢定也是不可忽略的一個重要步驟，雖然 SAS 有現成的語法可以執行，但因為所產生的報表和圖型落落長，不熟悉的人根本不會知道哪些圖或報表是代表什麼意義。SAS 內部工作人員利用官方的資料和程式發表了一篇技術文件，詳盡地解說每一張圖表在混合模型檢定中的意義。
在進行混合模型的估計中，模型檢定也是不可忽略的一個重要步驟，雖然 SAS 有現成的語法可以執行，但因為所產生的報表和圖型落落長，不熟悉的人根本不會知道哪些圖或報表是代表什麼意義。SAS 內部工作人員利用官方的資料和程式發表了一篇技術文件，詳盡地解說每一張圖表在混合模型檢定中的意義。
2012年3月13日 星期二
Count Data Models in SAS®
原文載點：[LINK]
在統計應用領域，我們經常會遇到所謂的"計數資料"(count data)，比方說就診人數，死亡人數，顧客人數等等。在分析這些計數資料時，最常使用的模型莫過於卜瓦松模型（Poisson model）。但在實際的應用上，使用者經常需要面對兩大問題：一是過度離散（overdispersion）問題，二是觀測值有太多零（excess zeros）的問題。這些問題通常都會導致卜瓦松模型估計的結果有偏誤。因此 WenSui Liu 和 Jimmy Cela 在 2008 年 SAS Global Forum 發表了一篇技術文件，介紹一些其他比較適合處理計數資料的模型。
在統計應用領域，我們經常會遇到所謂的"計數資料"(count data)，比方說就診人數，死亡人數，顧客人數等等。在分析這些計數資料時，最常使用的模型莫過於卜瓦松模型（Poisson model）。但在實際的應用上，使用者經常需要面對兩大問題：一是過度離散（overdispersion）問題，二是觀測值有太多零（excess zeros）的問題。這些問題通常都會導致卜瓦松模型估計的結果有偏誤。因此 WenSui Liu 和 Jimmy Cela 在 2008 年 SAS Global Forum 發表了一篇技術文件，介紹一些其他比較適合處理計數資料的模型。
2011年10月14日 星期五
Fitting Cox Model Using PROC PHREG and Beyond in SAS
http://support.sas.com/resources/papers/proceedings09/236-2009.pdf
Cox PH 模型在倖存分析（survival analysis）被廣泛的使用。一篇發表在SAS GLOBAL FORUM 2009的技術文件解說了在 SAS 底下如何用 Cox PH 模型來計算一些重要估計量的方法，並且提供一個方便的巨集程式來簡化繁複的程式寫作和運行。
Cox PH 模型在倖存分析（survival analysis）被廣泛的使用。一篇發表在SAS GLOBAL FORUM 2009的技術文件解說了在 SAS 底下如何用 Cox PH 模型來計算一些重要估計量的方法，並且提供一個方便的巨集程式來簡化繁複的程式寫作和運行。
2011年8月16日 星期二
Wow! You Did That Map with SAS/GRAPH®?
原文載點：http://www.nesug.info/Proceedings/nesug07/np/np01.pdf
身為一個專門從事空間時間序列分析的學者我來講，每天跟不同的地圖奮鬥是很稀鬆平常的是。目前市面上最熱門的地圖分析軟體莫過於ArcGIS，但這套軟體僅提供有限的分析方法，而其主要的功能僅僅就是繪製地圖而已。如果你用了別的軟體進行分析，就必須要將結果另存新檔，然後把結果跟地圖檔一起輸入到ArcGIS裡面才能開始繪圖。想當然爾，這裡面也沒有太多強大的資料整理功能，所以一切的一切都還是要經過另外的軟體將資料結果整理好才行。如果分析和繪製地圖的動作可以一次在SAS裡面完成的話將會省下許多時間和精力。NESUG 2007年有一篇教學文件，提供一些在SAS繪製地圖的資訊。
身為一個專門從事空間時間序列分析的學者我來講，每天跟不同的地圖奮鬥是很稀鬆平常的是。目前市面上最熱門的地圖分析軟體莫過於ArcGIS，但這套軟體僅提供有限的分析方法，而其主要的功能僅僅就是繪製地圖而已。如果你用了別的軟體進行分析，就必須要將結果另存新檔，然後把結果跟地圖檔一起輸入到ArcGIS裡面才能開始繪圖。想當然爾，這裡面也沒有太多強大的資料整理功能，所以一切的一切都還是要經過另外的軟體將資料結果整理好才行。如果分析和繪製地圖的動作可以一次在SAS裡面完成的話將會省下許多時間和精力。NESUG 2007年有一篇教學文件，提供一些在SAS繪製地圖的資訊。
2011年7月9日 星期六
MCPOWER: a Flexible Macro Suite for Generating Monte Carlo Power Estimates for Linear, Logistic, and Poisson Regression Models
原文載點：http://support.sas.com/resources/papers/proceedings11/430-2011.pdf
哈佛的一位統計分析師分享了一個名為 %mcpower 的巨集程式，能夠輕鬆的完成在線性迴歸、羅吉斯迴歸以及卜瓦松迴歸的蒙地卡羅檢定力模擬與估計。
哈佛的一位統計分析師分享了一個名為 %mcpower 的巨集程式，能夠輕鬆的完成在線性迴歸、羅吉斯迴歸以及卜瓦松迴歸的蒙地卡羅檢定力模擬與估計。
2011年7月2日 星期六
Fundamental Diagnostics for Two-Level Mixed Models: The SAS ® Macro MIXED_DX
原文載點：http://support.sas.com/resources/papers/proceedings10/201-2010.pdf
在進行 mixed model 分析時，模型的檢定是相當重要但也是相當煩人的一環。幾位南卡大和南佛大的學者寫了一個巨集程式來簡化相關程式的寫作，雖然目前這個程式只能使用在 2-level model，但也是具有相當貢獻。
在進行 mixed model 分析時，模型的檢定是相當重要但也是相當煩人的一環。幾位南卡大和南佛大的學者寫了一個巨集程式來簡化相關程式的寫作，雖然目前這個程式只能使用在 2-level model，但也是具有相當貢獻。
2011年6月25日 星期六
Expanding PROC EXPAND
原文載點：http://support.sas.com/resources/papers/proceedings11/417-2011.pdf
PROC EXPAND 程序是 SAS/ETS 裡面其中一個語法，提供許多跟時間相關資料的簡便處理。比方說，將日資料累加成月資料，月資料累加成季資料或年資料，也可以倒過來把季/年資料變成月資料，月資料變成日資料。若遇到缺漏值，該程序可用內插法將資料補齊。也可進行諸如moving average的資料轉換。Bruce Gilsen 在 2011 年的 SAS Global Forum 發表了一篇技術文件，介紹了兩個用 PROC EXPAND 運作的例子。
PROC EXPAND 程序是 SAS/ETS 裡面其中一個語法，提供許多跟時間相關資料的簡便處理。比方說，將日資料累加成月資料，月資料累加成季資料或年資料，也可以倒過來把季/年資料變成月資料，月資料變成日資料。若遇到缺漏值，該程序可用內插法將資料補齊。也可進行諸如moving average的資料轉換。Bruce Gilsen 在 2011 年的 SAS Global Forum 發表了一篇技術文件，介紹了兩個用 PROC EXPAND 運作的例子。
2011年6月18日 星期六
MORE CUSTOMIZATION?: CREATING SYMBOLS IN RTF FILES USING ODS
原文載點：http://www.lexjansen.com/pharmasug/2004/technicaltechniques/tt12.pdf
有時候想要在 SAS 輸出報表裡面寫一些標題或註解，需要加上一些科學符號、希臘字母、甚至是邏輯運算子等特殊符號時，得使用特別的方法顯示出來。這一篇 2004 年的 PHARMASUG 技術文件可供參考。
有時候想要在 SAS 輸出報表裡面寫一些標題或註解，需要加上一些科學符號、希臘字母、甚至是邏輯運算子等特殊符號時，得使用特別的方法顯示出來。這一篇 2004 年的 PHARMASUG 技術文件可供參考。
2011年6月4日 星期六
Creating Forest Plots from Pre-computed Data using PROC SGPLOT and Graph Template Language
原文載點：http://support.sas.com/resources/papers/proceedings10/195-2010.pdf
Forest plot 是研究環境流行病學的人經常會用到的一種圖表，主要是可以將許多變數或樣本的參數估計值以及信賴區間畫成一條一條的橫線，以供其他人易於判讀相對大小以及有無顯著。我在自己的博士論文裡面也放了相當多這種圖型，但之前我都是自己用 R 寫程式畫的，過程挺複雜的，沒想到 SAS 裡面有現成的語法可以畫這種圖。
Forest plot 是研究環境流行病學的人經常會用到的一種圖表，主要是可以將許多變數或樣本的參數估計值以及信賴區間畫成一條一條的橫線，以供其他人易於判讀相對大小以及有無顯著。我在自己的博士論文裡面也放了相當多這種圖型，但之前我都是自己用 R 寫程式畫的，過程挺複雜的，沒想到 SAS 裡面有現成的語法可以畫這種圖。
訂閱：
文章 (Atom)