【正文】 。 STDERR給出 LSM的標準誤 。 cards。 run。 回歸分析依自變數(shù)個數(shù)的多少分為一元回歸和多元回歸；依依變數(shù)和自變數(shù)之間關系的性質分為線性回歸和非線性回歸 。 默認情況下 ， CORR過程所進行的相關分析將給出分析變量的描述性統(tǒng)計結果 、 Pearson相關系數(shù)以及每個分析變量所對應的概率值 。 Corr過程－一個實例 許多害蟲的發(fā)生都和氣象條件有一定的關系 。 run。 Plot y軸變量 *x軸變量 。 model y=x。 三元線性回歸模型； Model y=x1 x2 x1*x2 x1*x1 x2*x2。 run。 cards。 output。 10 非線性回歸分析 在農學和生物學研究中 ， 大多數(shù)變數(shù)之間的關系不是簡單的線性關系 ， 而是各種各樣的非線性關系 ， 非線性回歸分析比較復雜 ， 主要有三個方面的內容：一是選擇合適的模型；二是非線性方程的參數(shù)估計；三是非線性方程和參數(shù)的假設測驗 。 ＝偏導表達式 。 DER語句 給出非線性回歸方程對參數(shù)的一階或二階偏導 。 Model 依變量＝模型表達式 。 model y=x1x4 /selection=forward。 do i=1 to 15。 應用實例 1 data new2。 proc glm。 回歸分析模型一般有以下形式 Model y=x。 180 175 28 40 25 160 117 120 165 80 。 Weight 變量 。 proc corr。 設定 partial變數(shù)時進行偏相關分析 。 PARTIAL 變量 。 從而能依據自變數(shù)x的數(shù)據對目標變數(shù) y進行預測或插值 。 means tr/t。 output。 LSMEANS 語句計算效應的最小二乘估計的平均數(shù)（ LSM） 。 若 x和 y之間存在線性回歸關系 ， 說明基礎生產力對來年產量有影響 ， 不能用原有的 y值進行方差分析 ， 必須消去 x的不同對 y帶來的影響 ， 即通過求y依 x的線性回歸方程 ， 將各處理的 yi都矯正到 x在同一水平時的值 。 要真正反映目標變數(shù) y的處理效應 ， 應先將不易控制的自變數(shù) x對目標變數(shù) y的影響剔除 ， 再進行方差分析 ， 這種分析即是協(xié)方差分析 。 數(shù)據塊結束 proc gplot data= wt5。 nn越大越平滑 。 VAXIS=數(shù)值 定義縱坐標的刻度。 與 P LOT 配對。 proc print。 將 x全距分為 10組 ,繪出直方圖 RUN。 LEVELS＝ n：將區(qū)間變量的數(shù)據分為 n組。 腳注內容 NOTE 選項 11 ’ 注解 ’ 。 垂直條形圖 BLOCK 變量表 / 選項 3。 而在 WINDOWS時代 ， 美觀漂亮也成為了軟件最重要的性能指標之一 ， 因此 SAS公司推出了高分辨率圖形 。dist39。 proc print。 do blk=1 to 2。 class i j。 output。 run。 應用舉例 2— 單因素隨機區(qū)組試驗 data new。 cards。 do a=1 to 3。 proc print。 Do a=1 to 5。 proc print。 do j=1 to 11。 Class語句 指定分類變數(shù)，此語句必須設定，并且應出現(xiàn)在 model語句之前。 與 ANOVA過程相比 ， GLM過程運行時一般占用較多的計算機內存和時間 。 24 30 28 26 27 24 21 26 31 28 25 30 32 33 33 28 21 22 16 21 。 do i=1 to 5。 模型定義中可用 |和 n簡化表達方法 水稻施肥試驗 ， 5個處理 ABCDE分別代表氨水 、廢氨水 、 碳銨 、 尿素和對照 （ 不施肥 ） ， 每處理 4盆（ 完全隨機試驗 ） ， 產量 （ 克 /盆 ） 結果如下： A: 24 30 28 26 B: 27 24 21 26 C: 31 28 25 30 D: 32 33 33 28 E: 21 22 16 21 請作方差分析。因而應依據試驗設計，給出正確的線性數(shù)學模型。 Means 效應表 /選項 。 當處理效應為隨機模型時 ， 估計方差分量 。 [eg13] NEANS過程可用于推斷單個樣本是否抽自平均數(shù)為 0的總體 。 data new。 run。 run。 結果輸出包括測驗平均數(shù)與 0差異顯著性的 t 值和概率 p。 clear。 Univariate 過程的重點在于描述變量的分布 。 未指定時 ， 輸出的統(tǒng)計數(shù)只有 N, MEAN, STD, MAX和 MIN。 指定要分析的數(shù)據集名及 一些選項 VAR 變量表 。 By語句用于指定分組變數(shù) ， 如果你需要分組處理數(shù)據 ， 例如要按 （ x） 類別分組輸出統(tǒng)計結果 ， 采用 by語句即可 。 Run。 幾點注意： 1） 該命令只能讀入 txt文本或 dbf文件； 2） 文本文件必須是純文本格式 （ 不帶逗號 、 不帶制表符或 TAB符 ） ； 3） 也不宜有表頭等格式 。 Run。 *數(shù)據集為 ； input x y 。 data 語句 若 libname mylib ‘e:\work’。 對于復雜的數(shù)據準備 ， 擬在稍后再作交代 ，本節(jié)只討論相對簡單的數(shù)據準備 。 其基本格式為： 注釋語句：使用星號 “ *” 開始 ， 可占多行 ， 以分號 “ 。 Do … end。 WHERE 條件表達式 。 用 SAS軟件 ， 你不必清楚怎么做 ， 但需清楚做什么 ！ 只需將數(shù)據按一定格式準備 ， 交由過程步即可實現(xiàn)相應的統(tǒng)計分析 。 Run(運行 )－用于執(zhí)行或調用程序。 在該窗口中將 SAS系統(tǒng)的所有輸出結果依次按照目錄樹的結構加以排列 ， 每一個過程步的結果被表示為一個節(jié)點 ， 展開該節(jié)點就可進一步看到表示不同輸出內容的子節(jié)點 ， 使用鼠標右鍵單擊每個節(jié)點 ， 就可對輸出結果的各部分進行察看 、 存儲 、 打印 、 刪除等操作 。 警告 ， 綠色 ， 以 Warning開始 ， 一般在程序中含有系統(tǒng)可以自動更正的小錯誤時出現(xiàn) ， 此時會提供錯誤序列號 。 SAS是第四代計算機語言的代表性軟件 ， 用戶只需清楚干什么 ， 不必清楚怎么做 ！ 使用相對方便 ， 用戶可以完成所有統(tǒng)計分析 、 預測 、 建模和模擬抽樣等工作 。 公司大事記（續(xù)） 2022： SAS CEO Jim Goodnight 在 2022美國商業(yè)大獎中榮獲StevieTM最佳企業(yè)管理人獎；新推出的 SAS Enterprise ETL Server在性能方面無人能及； 2022： SAS實現(xiàn)年銷售收入 19億美元； 2022： Ann Goodnight 進入北卡羅萊那大學董事會； 2022：銷售收入為 ； SAS在全球約有 45,000家客戶； 《 財富 》 全球 500強企業(yè)前 100家企業(yè)中有 91家是 SAS 客戶； 2022年銷售收入的 22%用于研發(fā)投入； SAS在全球設有 400多個辦事處 。 隨著統(tǒng)計軟件需求市場增大 ， 于1976年成立了 SAS軟件研究所 ， 主要從事統(tǒng)計計算軟件開發(fā)和銷售 ， 逐步發(fā)展成為全球一流的商業(yè)智能軟件和服務提供商 。 SAS（發(fā)音 sass） 是“統(tǒng)計分析軟件” （ statistical analysis software） 的首字母縮寫 ， 由 Jim Goodnight 和北卡羅萊那州立大學的同事 （ 包括 John Sall） 于上世紀 70年代初開始從事農業(yè)及生物試驗研究數(shù)據分析 。 公司大事記（續(xù)） 2022： SAS 在 Intel Developer Forum 中將電子智能集成到 64位計算中；SAS 將軟件移植到 Linux； 2022： SAS 慶賀公司成立 25周年；用戶評選 SAS 為分析軟件冠軍； 2022： SAS 連續(xù)第 6年榮登 Fortune的 “ 美國 100強公司 ” 名單； SAS 成立 Rapid Risk Deployment； Sun、 FEA初期解決方案合作伙伴； 2022： SAS 收購 Marketmax； Bank of America Foundation、 DPI、 SAS in School 向北卡羅萊納學校捐贈禮物；收購 OpRisk Analytics；IBM和 SAS成立戰(zhàn)略聯(lián)盟；推出 ； 2022： SAS CEO Jim Goodnight榮獲 StevieTM Award for Best Executive獎； SAS歷史上最重要的發(fā)布新的 SAS9軟件對商業(yè)智能行業(yè)帶來革命性的影響；這是一個功能強大的分析軟件 ， 能夠進行數(shù)百萬個高質量的 、 自動化預測 。 在 SAS/BASE的基礎上 ， 用戶還可以增加各種模塊而增加不同的功能 ， 如 SAS/STAT(統(tǒng)計分析模塊 )、SAS/GRAPH(繪圖模塊 )、 SAS/OR(運籌學模塊 )、SAS/IML(交互式矩陣程序設計語言模塊 )等 。 提示 ， 藍色 ， 以 Note開始 ， 提供系統(tǒng)或程序運行的一些常規(guī)信息 ， 大多數(shù)時候我們可以視而不見 。 3) 結果輸出（ Output）窗口 結果窗口幫助用戶瀏覽和管理所提交 SAS程序的輸出結果 。 Tools(工具 )－提供對各種輸出結構進行編輯的工具 ， 如：表格 、 圖形 、 報告等 ， 并支持進行系統(tǒng)環(huán)境和狀態(tài)的設置 。 很多統(tǒng)計分析難度較大 ， 常令非統(tǒng)計專業(yè)人士頭疼 。 VAR 變數(shù) 。 else if … then statement 2； else statement 3。 一行多句與多行一句 為使長程序清晰易讀 ， 在相應程序段間可以使用空行分隔 ， 也可以使用注釋加以說明 。 但統(tǒng)計分析大多比較復雜 ， 不同的過程有不同的數(shù)據格式要求 ， 數(shù)據準備亦有一定的難度 。 若原來就有該數(shù)據集 ， 則把其中的內容清除 。 *確定 x1x10 十個變數(shù)； Input x$ y ； *$指明變量 x為字符變數(shù)， 表明數(shù)據是連續(xù)讀入的； Input 語句 data temp。 Proc print。 即可將 e:\work\ x1x4讀入到aaa數(shù)據集中 [eg6]。 By 變數(shù)序列 。 如果條件變?yōu)?x、 y兩個變數(shù)的值都要大于50， 則 where語句改為： Where x50 and y50。 描述性統(tǒng)計過程 MEANS過程格式 PROC means data=數(shù)據 集名 。 指定統(tǒng)計量對應的新變量名 MEANS語句的選項 設定結果輸出的統(tǒng)計數(shù) ， 如 STD, MEAN, VAR等 （ 可選項有 CSS, CV, KURTOSIS, MAX, MEAN, MEDIAN, MIN, N, P1, P99, PROBT, Q1, Q3, QRANGE, RANGE, SKEWNESS, STDDEV, STDERR, SUM, SUMWGT, T, VAR等共 38項 ） 。 SUMMARY過程 univariate過程除可計算基本統(tǒng)計數(shù)之外 ， 主要用于計算分位數(shù) （ Quartile， 中位數(shù)為一特例 ， 即 50％ 分位數(shù) ） 、 繪制分布圖 、 次數(shù)表及正態(tài)分布測驗等 。 log。 RUN。 proc means n mean stderr t prt。 var y。 測驗兩品種千粒重差異的顯著性 。 run。 當處理效應為固定模型時 ， 對各個平均數(shù)進行多重比較 。 *用于指明分類變數(shù)； Model 依變數(shù)＝效應表 /選項 。同一試驗資料，分析結果依模型不同而異。 Model y=a b(a)； 嵌套（ nested）模型，用于系統(tǒng)分組資料。 Data new。 cards。 當用于方差分析時 ， 主要處理非平衡資料 ， 可進行多個變量的對比檢驗 。 語句說明 Proc glm語句 設定分析數(shù)據集和輸出數(shù)據集等。 do i=1 to 3。 . . 。 應用舉例 1－兩向分組資料方差分析 Data new。 13 14 14 12 12 13 3 3 3 10 9 10 2 5 4 。 二因素有重復資料方差分析表 ANOVA的例子 data new。 end。 [eganova22] F A E D C G B 區(qū)組 1 20 24 22 18 21 20 20 A D C F B E G 區(qū)組 2 26 16 19 21 19 20 19 F B G E C A D 區(qū)組 3 21 21 21 23 20 23 19