日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

sas教程

發布時間:2023/12/20 编程问答 41 豆豆
生活随笔 收集整理的這篇文章主要介紹了 sas教程 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

?

http://web5.pku.edu.cn/pucssr/SASbiancheng.pdf?

?

本教程中的主題將向您介紹 SAS Enterprise Guide。您最好依次瀏覽這些主題。

?

概述
啟動項目并瀏覽主要窗口
向項目添加 SAS 數據
從文本文件導入數據
關于 SAS 任務
創建和修改列表報表
創建條形圖
關于“查詢生成器”
使用查詢連接表
將計算列添加至查詢
從查詢生成“匯總表”
創建餅圖
執行線性模型分析
將多個報表組合成單個文檔
使用過程流
更多內容

?

?

簡明SAS教程

by LongFei?
version:0.11

SAS?統計軟件


前言

簡明致敬郇爺,細節留作習題。(霧)

本教程的內容主要是SAS的簡單基本知識和SAS在多元統計分析中的應用。?
側重在SAS的應用,會介紹如何根據SAS的輸出結果得到結論,但不會深入講解其中的統計原理。?
在每一章的開始,會簡要概括本章內容和學習目標,以幫助同學們有針對性地進行學習。

代碼講解將由一條條的statement講解來構成。

作者的話

這篇文章是為了整理自己最近復習SAS的心得筆記而寫,旨在梳理自己的學習內容。

我以前修讀多元統計分析課程時,SAS上機作業常常是照著課本和課件做的,對代碼的含義一知半解,一旦提交了作業,幾乎就會忘掉代碼了。如今復習,由于已積累了較多的編程實踐經驗,對代碼的理解也加深了。我把自己的理解寫出來,希望能幫助正在或打算學習SAS的同學們。

原本只是打算寫一個小教程,但寫作過程中有點貪多求全,現在越寫越多,只能強行“簡明”了。

注意事項

本教程使用的SAS版本是9.4,SAS版本最低不低于9.3。


快速入門

本章介紹SAS編程的基本概念和邏輯,并寫一個簡單的示例程序,以期讓初學者快速上手SAS編程。

學習目標

  • 了解DataSet(數據集)、Procedure(過程),Statement和Option(選項)的概念
  • 了解SAS程序執行流程和SAS代碼基本語法結構。
  • 了解簡單日志信息的含義
  • 寫一個Hello World程序

基本概念

DataSet(數據集)
一張表,每列代表一個變量,每行代表一個觀測。數據集是SAS進行統計分析的基本單位。
Procedure(過程)
簡稱為proc,SAS集成的各種統計分析功能的總稱。可以使用不同的過程來完成不同的分析。
Statement(語句)

SAS代碼,無論分成多少行,都要以一個分號結束,這一條到分號結束的代碼,就被稱為一個statement或語句,這是SAS執行代碼的基本單位。一條條的語句,構成了SAS程序。?
statement之間是有層次關系的。頂層statement是一個獨立執行單元的statement,比如data xxx;,proc print;,run;等等。次級statement則需要依賴于頂層statement——先執行頂層statement,才能執行次級statement,典型的次級statement就是by xx;,input a b;等等。?
次級statement是有使用范圍的,需要在正確的頂層statement下才能被SAS執行,他們的主要作用是在細節上控制主statement的行為,在某些過程(proc)中,某些次級statement是必須設置的。比如說在排序過程(sort)里,就不能直接執行proc sort;,而必須在它后面添加by 變量名;來指定用于排序的變量后,才能執行。?
以下示例代碼的每一行都是一個statement

  • data xxx;
  • set a;
  • proc print data=xxx;
  • run;
  • Option(選項)
    選項是指寫在statement語句后,分號前,以空格間隔的一系列單詞或等式,比如noprint,data=xxx等。它們通常不是必要的,但能控制statement的宏觀行為如:指定用于分析的數據集,不輸出結果,指定繪制直方的數量等等。?
    頂層statement和次級statement的option寫法稍有不同。頂層statement的option是接空格如?proc procname data=xxx noprint;,而次級statement是接/如histogram RM/normal nmidpoints=12;

    程序編寫流程

    dataprocrun data
    代碼以data作為標識符起頭,稱為data步(數據步)。數據步的作用是建立數據集,進行數據的讀取和修改,后續可以提供給過程步進行統計分析。SAS提供了強大的數據輸入功能,能夠把各種格式的數據輸入到SAS中。此外,還能在數據輸入的過程中進行計算,生成新的數據列,甚至能處理缺失數據,功能十分強大。
    proc
    代碼以proc起頭,稱為proc步(過程步)。過程對指定的數據集進行統計分析。SAS提供了豐富多樣的統計分析功能,還有各種選項可供細化。
    run
    寫下一個run;,就可以讓SAS在讀取到這一步時,運行前面的所有語句。如果沒有這一條,則代碼只是被SAS讀取了,但不會運行。

    我們編寫SAS程序,正是按照上述的運行過程來做的:?
    1. data步——獲取和整理數據?
    2. 寫proc步——分析數據?
    3. 提交代碼給SAS,得到output——輸出分析結果

    程序執行流程

    提交代碼逐條讀取語句讀取run后運行并輸出繼續讀讀完代碼

    在提交代碼后,SAS會逐條讀入statement。當讀到run;這條statement時,SAS就會執行前面的所有statement,輸出分析結果。如果我們的SAS程序中沒有run;語句,則程序會處于持續運行狀態,但不會執行前面讀取到的statement。?
    有些proc步的statement是持續性的,會阻止下一個proc的執行。這時就需要加一句quit;來終止該proc步。這是后話,待到用時再說。

    在SAS 9.3以后,在默認設置下,proc步的分析結果會被輸出到一個html文件上,這個html文件會被SAS用結果瀏覽器打開。如果不喜歡這種方式,可以使用ODS(Output Delivery System 輸出交付系統)來改變輸出方式和樣式。這在后續的章節中會講到。

    基本語法

    當初我學習的時候沒有特別考究過SAS的語法,老師也沒有強調,導致我常常需要記憶整段的代碼來解決某個具體問題。這種學習方式效率十分低下,也缺乏拓展性——問題要求稍微變化一點就沒轍了。?
    事實上,SAS的代碼雖然語法簡單而選項繁多,但也有基本的編程語法結構。了解語法結構,能讓我們在記憶時更有針對性,更容易理解代碼含義,從而提高我們解決相似問題的能力。

    從上一節的SAS運行過程可以看出,要編寫一個SAS程序,就需要編寫data步和proc步的代碼,它們有著相似的語法結構。編寫完成之后,再寫一行run;讓程序執行statement即可。

    Data步

    Data步以data作為開頭標識,空格后接自行設定的數據集名,生成一個數據集。下面接上各種statement描述如何處理這個數據集。Data步有自己的次級statement語句,只可在data步中使用。?
    一個完整data步的語法結構如下,以<>括起來的部分都是可選的,不必嚴格遵守下面的順序:

  • DATA datasetname;/* 設定數據集名 */
  • <code that defines the variables in the data set>;/* 定義數據集的變量 */
  • <code to enter data>;/* 輸入數據 */
  • <code to create new variables>;/* 創建新變量 */
  • <code to assign missing values>;/* 設定缺失值 */
  • <code to output data>;/* 設定輸出的數據集 */
  • <code to assign labels to variables>;/* 設定變量的標簽 */
  • <and other data tasks>;/* 其它數據任務 */
  • Proc步

    Proc步以proc?作為開頭標識,空格后接SAS過程名,再一個空格后接該過程的選項,最后仍以分號結束。分號后面接上多句次級statement,。同data步一樣,proc有自己的次級statement語句。不過每一個過程的次級statement語句都不一樣,有些能通用,有些則不能,有些還是必寫的,需要具體問題具體分析。一般的語法結構如下,:

  • proc procname <options>;/* 啟動過程步 對數據進行統計分析 */
  • <statements>/<options>;/* 次級statement,如果要添加選項,必須加/隔開 */
  • run和quit
    run;和quit都是一條statement語句,而且是不依賴于data步和proc步的statement。?
    當SAS執行了run;時,它就會把前面讀取到的所有data步或proc步的statement都執行完。如果一個程序里沒有run;,那么這個程序在提交后也不會被執行。?
    當SAS執行了quit;時,它就會把某些會持續運行,占用SAS資源的的proc停掉,讓SAS能夠執行后面的proc。
    大小寫問題
    在SAS中,各種statement的關鍵字,如data,proc等,都是不區分大小寫的。不過,在字符串中,則會區分大小寫。

    學習心得:?
    “做中學”是以實際問題驅動學習的方式,這在編程語言的學習中非常高效。我認為這一方式同樣適用于SAS的編程學習。在學習中應當多多實踐,盡可能完成課本的所有練習,加深對SAS代碼的理解。?

    對data步,重要的是掌握它的次級statement操作,包括創建、讀取、修改、刪除、輸出等等。處理數據是統計分析的開始,也是重要的一步。在多元統計分析的學習中,我們更關注統計分析的實踐。練習中使用的數據一般都是處理好了的標準數據,可以直接使用。不過,在工作中可能會遇到各種各樣的數據,有很多可能會有缺失值,如果想要在SAS的使用上更進一步,就一定要對data步有更多的了解。?

    對proc步,則需要按不同的過程,逐個學習。對有些功能相似的過程,可以先學習其中一個過程,然后對比另一個過程,進行遷移學習。?
    在學習一個過程時,我喜歡首先了解它的應用背景,也就是了解它是用來做什么的。然后學習它常用的statement,知道如何從細節上控制分析過程。做到這一步就足夠完成一些基本的統計分析了。但有時候有一些特殊需求,比如不輸出分析結果,或者畫直方圖時順帶畫出正態分布曲線,這就需要設置可選項了,這時候再去了解這個過程的可選項,以及相應的某些核心的statement的可選項即可。

    日志(Log)

    在向SAS提交代碼之后,務必檢查日志里的內容,確認程序的執行情況。特別是需要關注錯誤和警告。

    錯誤
    紅色的字表示錯誤(error),說明代碼出錯了,會導致程序在停止在SAS所提示的出錯行,不再繼續往下執行代碼。
    警告
    綠色的字表示警告(warning),說明當前提示的這句代碼可能有問題,有可能會讓你得不到預期的結果,但程序還是能繼續運行下去。

    學會閱讀日志能幫助快速地排除小BUG,避免不必要的時間浪費。?
    日志默認是打開的,如果沒有打開,或者不小心關掉了。可以手動打開:菜單欄-->視圖-->日志

    Hello World

    將以下代碼復制到SAS的代碼編輯器里并提交

  • data hellodata;
  • input ISay $ 1-20;
  • datalines;
  • Hello World One
  • Hello World Two
  • Hello World Three
  • ;
  • proc print;
  • run;
  • 會得到以下輸出:

    ObsISay
    1Hello World One
    2Hello World Two
    3Hello World Thre

    數據輸入和儲存

    本章介紹SAS的數據集的基本知識和data步中的數據輸入方法。

    學習目標

    • 了解SAS數據集類型、結構和變量數據類型
    • 掌握input的使用方法,能夠精細地設定數據輸入格式
    • 掌握infile的使用方法,能夠導入數據文件
    • 掌握datalines的使用方法,能夠手動輸入數據

    SAS數據集基本知識

    數據集類型

    數據集分為臨時數據集和持久數據集兩種類型。

    臨時數據集
    數據集名是一個single name,也就是簡單的名字,如:data example;,data dec2016等。臨時數據集的數據,會暫時儲存在SAS程序里,但一旦關掉SAS,就會被銷毀,想要重新使用數據的話,就需要再次提交輸入數據的代碼。
    持久數據集
    數據集名則是一個two-part name,也就是用一個.連接的復合名字,如:data research.2016,data MyData.0001等。這個復合名字的前半部分是SAS的 library(邏輯庫)名,后半部分是數據集名。這樣的數據集會以文件的形式被儲存到硬盤上,具體的位置是library指定的路徑。即使關掉了SAS,下次再打開,只需要分配邏輯庫名,就能直接使用,無需再次提交數據輸入的代碼。

    實際上,臨時數據集在SAS中也是存在library里的。這個library的名字叫work。每次關閉SAS時,系統會自動銷毀這個library里的數據。?
    在數據量比較小,數據輸入耗時可忍受(幾秒鐘)的情況下,兩種數據集沒有什么使用上的差別。但一旦數據規模擴大到輸入耗時幾十分鐘甚至幾十個小時的時候,持久數據集就要優于臨時數據集了。?
    一般來說,實驗可以使用臨時數據集,做實際工作則應該使用持久數據集。

  • /* Ex 2.1 臨時數據集和持久數據集的區別 */
  • /* 臨時數據集實例
  • 每一次重新打開SAS,要運行代碼中的proc步,就必須先提交data步。
  • */
  • data tempdata;/* 輸入數據集 */
  • input t @@;
  • datalines;
  • 1 2 3 4
  • ;
  • proc print data=tempdata;/* 打印數據集 */
  • run;
  • /* 持久數據集實例
  • 只有第一次需要提交全部的代碼,之后數據集就會被以文件形式儲存在MyData庫中。之后要使用,只需要執行第一句分配邏輯庫名,就可以直接提交proc步,而無需再次提交data步輸入數據。
  • */
  • libname MyData "D:\MyData";/* 分配邏輯庫名 */
  • data MyData.tempdata;/* 輸入數據集 */
  • input t @@;
  • datalines;
  • 1 2 3 4
  • ;
  • proc print data=MyData.tempdata;/* 打印數據集 */
  • run;
  • 如果知道數據集文件的儲存路徑,還可以直接用文件路徑名調用。比如,假設文件`D:\MyData目錄下有一個名為tempdata的SAS數據集文件,我們可以直接使用:

  • proc print data="D:\MyData\tempdata";
  • run;
  • 但是,不建議使用。

    數據集的結構

    數據集的結構就是一張表,每一列表示一個變量,每一行表示一個觀測。以下是一個典型的SAS數據集?
    在第一行中,第一列的Obs表示觀測的序號,后面的a b c分別是變量名

    Obsabc
    1123
    2456
    3789
    變量命名規則

    SAS允許為變量自由命名,但存在一些規則:

    • 可以設定長度為1-32的任意字母、數字或下劃線,不能包含空格。
    • 必須以字母或下劃線開頭
    • 不區分大小寫
    • 最好具有描述性(就是說能從變量名看出這個變量的含義,少用諸如a,b,c這樣無語義的變量名),這項規則是可選的,但建議盡量遵守。

    一些正確的變量名:?
    GENDER?AGE_IN_1999?AGEin1999?_OUTCOME_HEIGHT_IN_CM?WT_IN_LBS?
    一些錯誤的變量名:?
    AGE IN 2000?2000MeaslesCount?S S Number?Question 5?WEIGHT IN KG?AGE-In-2000

    變量數據類型

    學習過編程語言的同學都知道,在編程語言中,變量的數據類型是很重要的。在SAS中,雖然沒有明確指出,但實際上也是存在不同的數據類型的,分別如下:Numberic
    數字型變量,也是SAS的默認變量類型。決定一個變量是不是數字型,只需要看它在邏輯上是否可以進行運算,如身高、體重之類,就可以用于計算,ID號、電話號碼之類,不可用于計算,所以前兩者是數字型變量,后兩者不是。這種變量類型會自動忽略掉高位的0,如果你把ID號設置為數字類型,那么00012就會變成12。
    Character
    字符型變量,也包括Text(文本)或String(字符串)。這種類型的變量不能用于運算,并且會區分大小寫。如果你對性別變量分組,那么M和m會被分成不同的組。SAS默認會把數據看成數字型變量,所以必須要手動把諸如ID號之類的數字形式的變量設置成字符型變量。
    Date
    日期型變量。要輸入SAS的數據集中的日期變量可能會有多種多樣的格式,比如10/15/09、01/05/2010或JAN052010等等。在后面會看到如何處理這多種多樣的日期格式數據,并把他們轉化為SAS能夠用于分析的數據。實際上,在SAS數據集中儲存的日期型變量數據是用整數表示的,這能很方便地進行數據操作。

    變量數據類型實例程序

  • /* Ex 2.2 變量數據類型實例程序 */
  • data vartypes;
  • input id_num id_char $ weight_num @18 birth_date MMDDYY10.;/* 輸入4個變量,分別是id_num(默認-數字型) id_char(字符型) weight_num(默認-數字型) birth_date(日期型,格式為:月/日/年) */
  • datalines;
  • 0001 0001 60 03/07/2016
  • 0002 0002 50 04/02/2016
  • 0003 0003 30 12/31/2016
  • 0004 0004 45 04/12/2016
  • 0005 0005 49 04/07/2016
  • ;
  • proc print data = vartypes;
  • run;
  • 打印的結果如下:

    Obsid_numid_charweight_numbirth_date
    1100016020520
    2200025020546
    3300033020819
    4400044520556
    5500054920551

    可以看到id_num的0被去掉了,而birth_date變成了整數形式

    數據的輸入

    數據的輸入分為兩步:首先,從外部文件或流式數據(in-stream data,可以簡單地理解為直接嵌在代碼里的數據)中讀取原始數據,然后是將原始數據轉化為SAS使用的標準數據集,需要設定一定的轉化規則。?
    SAS提供的數據輸入功能十分強大,但是實際用到的不多,只講簡單實用的部分即可。

    讀取原始數據

    外部文件
    使用infile語句進行讀取。格式為?infile filepath;。可以后接選項詞做進一步的輸入控制。常用的有DLM,用于指定分隔符;FIRSTOBS,用于指定從哪一行開始讀數據;OBS指定到哪一行結束讀取。?
    infile "D:\MyData\MyCSV.csv" DLM=',' FIRSTOBS=2 OBS=10表明從D:\MyData\MyCSV.csv中讀取數據,分割符為逗號,從第2行開始讀,第10行結束。
    流式數據
    以datalines;開頭,換行之后就可以寫數據了,到最后換行,以一個分號結束即可。?
    如果要寫入的數據中有分號,則可以以dalines4;開頭,最后結尾換成三個分號。 不要使用cards,這是很老舊的用法了,現在的datalines;完全可以替代cards;的作用

    input語句用于設定如何將各種格式的數據轉化成SAS使用的標準數據集。

    設定轉化規則

    讀入了原始數據之后,需要設定轉化規則,用數學說原始數據集到SAS數據集的映射,把原始數據轉化成標準的SAS數據集。SAS提供了多種多樣的方法來確定轉化規則,但一般來說,我們只需要使用input語句來設定規則就足夠了。?
    input語句的通用的語法為:input <specification(s)> <@ | @@>;?
    specification(s)表示一個或多個特殊的表述,內容可以包括變量,以圓括號括起來的變量列表,$符號,指針控制,列控制,格式化標識等等,非常復雜,具體可以通過SAS的官方幫助手冊里查input Statement進行了解。下面列舉最常用的幾種輸入格式,能夠解決絕大多數的輸入問題。它們之間可以互相混合使用以完成最佳輸入。

    Freeform List Input

    Freeform List Input(自由輸入)是最簡單也是最常用的。只需要設置變量列表,以及指定某些變量為字符類型即可,SAS會自動讀取數據,數據之間以空格作為間隔符。

    通用語法:input variable <$> ;?
    variable為變量名,$是可選項,表明指定它前面的變量為字符類型。這最簡單的語法,在大多數時候已經能滿足學習中的數據需要了。這主要用于以作為空格間隔的原始數據。

  • /* Ex 2.3 Freeform List Input */
  • data test;
  • input x $ y z; /* x后接一個$,表明x是字符型變量,不省略高位的0 */
  • datalines;
  • 001 2 3
  • 004 5 6
  • 007 8 9
  • ;
  • proc print;
  • run;
  • 輸出結果:

    Obsxyz
    100123
    200456
    300789
    Compact Data Format

    Compact Data Format(緊湊輸入)也是非常常用的,比起自由輸入,只是增加了末尾用于控制數據輸入指針的符號。?
    通用語法:input variable <@|@@>;?
    <@ | @@>表示可以寫@或@@或不寫,這決定著讀取數據的指針的行為。?
    -?不寫:指針立刻換行?
    -?@:指針暫不換行,直到下一個input或一輪循環結束。?
    -?@@:指針不換行。循環結束也不換行。?
    緊湊輸入指的是末尾添加@@符號的寫法,這可以讓SAS連續地讀取所有的數據。

  • /* Ex 2.4 Compact Data Format-input中@數量的區別 */
  • data test;
  • input x @; /*單個@,暫時不換行,直到下一個input或一輪循環結束 */
  • input y; /*沒有@,下一個input語句會去讀新的一行*/
  • input z @@; /*兩個@,下一個input不換行,直接從下一個數據開始讀*/
  • datalines;
  • 1 2 3
  • 4 5 6
  • 7
  • ;
  • proc print;
  • run;
  • 輸出結果:

    Obsxyz
    1124
    2567
    Column Input

    Column Input(列輸入)是以列控制指針來控制輸入的,可以以精確到具體列的形式進行數據的讀取。?
    基本語法:input variable <$> startcol<-endcol> ...;?
    <$>正是前面說的字符類型控制符。startcol指明從哪一行開始讀,<-endcol>是可選項,指明到哪一行讀取結束,不寫的話,就默認取starcol的值。

  • /* Ex 2.5 Column Input 列輸入 */
  • data MyData;
  • input ID $ 1-3 SBP 4-6 DBP 7-9 GENDER $ 10 AGE 11-12 WT 13-15;
  • datalines;
  • 001120 80M15115
  • 002130 70F25180
  • 003140100M89170
  • 004120 80F30150
  • 005125 80F20110
  • ;
  • run;
  • 輸出結果

    ObsIDSBPDBPGENDERAGEWT
    100112080M151
    200213070F251
    3003140100M891
    400412080F301
    500512580F201
    Formatted Input

    Formatted Input(格式化輸入)是用于將特殊的數據格式轉化為SAS可識別的數據類型的,比如日期或者帶小數點的數字。?
    通用語法:input formated-vars ;?
    formated-vars指的是格式化變量,標準格式是:col-pointer variable informat?
    其中,col-pointer指明從第幾列開始讀取數據,用@n來表示,n為指定開始讀取的列,variable指出變量名,informat指出輸入的格式,所有的informat都必須以點號.或點號.后接1個數字的形式作結尾,具體格式表如下:

    informat輸入格式表

    InformatMeaning
    n.n 列數據,數字類型
    $n.n列數據,字符類型,移除頭部的空格
    $CHARn.n列數據,字符類型,保留頭部的空格
    COMMAn.n列數據,數字類型,去除小數點或美元符號(比如$40,000會被讀取為40000)
    COMMAn.mn列數據,數字類型,其中帶有2位小數,去除小數點或美元符號(比如$19,020.22會被讀取為19020.22)
    MMDDYY8.輸入日期,形如01/12/16(月-日-年)(注意Y2K問題)
    MMDDYY10.輸入日期,形如04/07/2016(月-日-年)
    DATE7.輸入日期,形如20JUL16(日月年)
    DATE9.輸入日期,形如12JAN2016(日月年) (沒有Y2K問題)

    需要特別說明地是,日期在SAS中是以整數的形式儲存的,這樣便于進行統計分析。如果需要輸出日期的數據,則需要在輸出時,特別指出輸出格式。

  • /* Ex 2.6 Formatted Input */
  • data MyData;
  • input @1 SBP 3. @4 DBP 3. @7
  • GENDER $1. @8 WT 3. @12 OWE
  • COMMA9. @22 InputDate MMDDYY8.;
  • format InputDate MMDDYY8.; /* 指定InputDate的輸出格式,沒有這句則會輸出整數 */
  • /* 列標尺
  • ----+----1----+----2----+
  • */
  • datalines;
  • 120 80M115 $5,431.00 01/11/16
  • 130 70F180 $12,122 01/12/16
  • 140100M170 7550 01/12/16
  • 120 80F150 4,523.2 01/12/16
  • 125 80F110 $1000.99 01/12/16
  • ;
  • proc print data=MyData;
  • run;
  • 輸出結果

    ObsSBPDBPGENDERWTOWEInputDate
    112080M1155431.0001/11/16
    213070F18012122.0001/12/16
    3140100M1707550.0001/12/16
    412080F1504523.2001/12/16
    512580F1101000.9901/12/16

    判別分析

    判別分析的應用場景:對于總體,各自有自己的觀測集,現在給出一個樣本,需要判斷它屬于哪一個總體。

    根據不同的判別指標和條件,我們可以用不同的衡量方法,來進行不同的判別分析,它們有各自不同的優劣。

    距離判別

    距離無疑是最直觀的衡量方法,樣本距離哪個總體近,就判它屬于哪個總體。?
    具體地說,就是利用觀測集來估計總體(或總體的某些關鍵統計量),然后以某種方式來衡量觀測與總體的遠近程度,然后選擇最近的總體。

    馬氏距離

    要判斷距離,首先得明確計算距離的方法。現在我們希望了解的是給定樣本到某個總體的距離。聯系歐氏空間,很自然地就想到,這個距離可以定義成到總體的中心的距離,而總體的中心可以用它的均值向量來表示,為了方便計算,我們一般取距離的平方。此外,還要考慮總體的集中程度,越集中的,自然就越高,這個在一元下可以用方差表示,在多元下則是用協差陣來表示。由此想法,我們可以定義樣品與總體的馬氏距離為?

    ?

    對于多個總體,我們只需要比較的大小即可

    特別地,在兩總體比較中,我們可以定義一個線性函數,從而得到判別準則:?

    ?

    SAS程序

    在SAS中,我們可以使用discrim過程來完成上述的各種判別分析。

    通用語法是

  • proc discrim <options>;
  • class variable;
  • <other statements>;
  • class variable;是必填的,它決定著要用什么變量來劃分類別。然后我們可以根據不同的需求來添加options和statement完成統計分析。

    馬氏距離判別

    ?

    ?

    ?

    ?

    ?

    ?

    ?

    ?

    ?


    統計軟件SAS使用教程

    2009-07-17??fire216
    百度空間?|?百度首頁? |?登錄 長江的空間 技術的匯集地:人杰地靈. 主頁博客相冊|個人檔案?|好友
    ? 查看文章 ??
    統計軟件SAS使用教程 2007-06-16 23:31
    第1章SAS系統簡介
    §1.1 SAS發展概況
    SAS(Statistical Analysis System)是一個大型的數據管理與數據統計分析處理的軟件包。1966年由美國North Carolina州立大學開始研制,1976年在美國成立了SAS研究所。
    SAS主要用于數據處理和統計分析領域,是一個功能齊全、使用方便靈活。只須要有少量的、簡單的語句、寫出SAS程序,進行運行,就可以滿足擁護要求,一些特殊的計算或處理方式可以通過選項來指定。從而達到且讓用戶將SAS程序在SAS環境下提高運行,及時了解到程序運行情況及出現的錯誤,程序可以方便的編輯修改和運行,直到用戶得出滿意的計算結果。
    目前,SAS已經發展成為一個功能齊全、應用范圍廣泛、使用靈活方便的數據庫管理和數據分析的標準軟件系統。其統計分析部分,在數據處理和統計分析領域,被業界和國際上公認為標準軟件和最權威的統計軟件包。
    SAS應用廣泛,其應用范圍涉及到理、工、農、林、醫、管理、商業、行政事物等各個領域。國際上成立了專門的SAS協會SUGI(SAS User Group International),每年有學術會議討論研究有關SAS的問題。
    SAS在我國的應用。
    SAS的主要版本。
    §1.2 SAS的結構、功能、特點
    一、SAS結構與功能:
    SAS軟件包由多個大的功能模塊組成,用戶可以根據需要,選擇安裝部分或全部SAS功能模塊來組成一個運行系統。
    SAS系統的核心(基本)部分是SAS/BASE模塊,其功能是承擔數據管理,管理用戶使用環境,進行用戶語言的處理,調用其他模塊。在SAS/BASE模塊的基礎上,還可以增加如下不同的模塊、從而實現不同的功能。
    ⑴SAS/BASE?
    ⑵SAS/STAT
    ⑶SAS/AF
    ⑷SAS/FSP
    ⑸SAS/GRAPH
    ⑹SAS/ETS
    ⑺SAS/IML
    ⑻SAS/OR
    ⑼SAS/QC
    二、SAS的特點:
    l 實用性強、功能完善、使用方便、編程簡單、容易學習。
    l 信息存儲簡單,數據交換接口豐富,可以調用多種格式的數據文件。
    l 語言編程能力強大。
    l 對數據連續處理。
    l 統計分析方法豐富,使用簡單。
    l 具有強大的報表輸出和繪圖功能。
    三、SAS的運行環境:
    l 硬件系統:586計算機
    l 操作系統:Windows95/98/2000
    §1.3 SAS for Windows 的操作入門
    SAS for windows系統的操作界面是標準的圖形用戶界面,是標準的多文檔界面,如同Word、Excel等軟件一樣,操作使用簡單、類似。軟件操作可以通過菜單或命令方式完成。
    一、 SAS的安裝與啟動
    ⒈ SAS的安裝:
    ⒉ SAS的啟動:
    ⒊ 系統的設置:
    二、 SAS集成環境介紹:
    ⒈ SAS的集成環境:
    ⒉ 程序編輯窗口(Program Edit):
    ⒊ 日志窗口(Log):
    ⒋ 輸出窗口(Output):
    ⒌ 其它窗口:
    三、 SAS操作入門:
    ⒈ 程序的編輯:
    ⒉ 程序的運行:
    ⒊ 結果查看:
    ⒋ 有關文件的操作。
    §1.4 SAS/ASSIST簡介
    SAS/ASSIST是為初學者學習使用SAS系統提供的一個交互式輔助系統,目的是初學者更容易的掌握SAS的強大功能。
    一、 SAS/ASSIST的啟動:
    二、 SAS/ASSIST的操作入門:
    三、 SAS/ASSIST的的學習資料:
    請訪問:cie.nwsuaf.edu.cn/apptech/sjm/sas?? οnlοad="this.style.overflowX='auto';">第2章SAS語言的基本概念
    §2.1 SAS語言概述
    SAS提供了一種完善的編程語言。如同大多數計算機高級語言一樣,SAS用戶只須要熟悉其命令、語句及簡單的語法規則就可進行數據管理和分析處理工作。因此,掌握SAS編程技術是學習SAS的關鍵環節。
    在SAS中,系統具有眾多的計算過程,把大部分常用的復雜的數據計算的算法作為標準過程調用,用戶僅需要指出過程名及必要的參數即可,這樣使得SAS的變成變的十分簡單。
    SAS程序是SAS語句的有序集合。以程序的功能將SAS程序劃分成多個程序模塊(程序段),這些模塊分為兩類:數據準備模塊——數據步(DATA Step)和過程部(PROC Step)。在一個SAS程序中可以有一個或多個數據步與過程步。
    SAS程序的結構和數據流向示意圖
    數據步的作用是把數據源中的數據作為輸入,經過加工后輸出到一個或多個“SAS數據集”。
    過程步是把由數據步建立和產生的數據集中的數據作為輸入,調用一個或多個標準過程,經標準過程計算處理后,將結果以清晰的表格或圖形方式輸出到OUTPUT窗口中或寫入磁盤文件。
    簡單說:數據步是為過程步準備數據,并將數據放入數據集中;
    過程步是把指定的數據集中的數據計算處理后并輸出結果。
    一個簡單的SAS程序的如下:
    data student;
    ??? input xh$ xm$ yy gs;
    ??? p=yy gs;
    ??? cards;
    ?????? 001 aaa 89 91
    ?????? 002 bbb 70 88
    ?????? 003 ccc 90 92
    ?????? 004 ddd 70 87
    ?????? 005 eee 85 83
    proc print;
    proc means;
    run;
    運行結果如下:
    數據集:
    OBS???? XH????? XM????? YY???? GS????? P
    1????? 001???? aaa???? 89???? 91???? 180
    2????? 002???? bbb???? 70???? 88???? 158
    3????? 003???? ccc???? 90???? 92???? 182
    4????? 004???? ddd???? 70???? 87???? 157
    5????? 005???? eee???? 85???? 83???? 168
    平均數:
    Variable?? N?????????? Mean??????? Std Dev??????? Minimum??????? Maximum
    -------------------------------------------------------------------
    YY???????? 5???? 80.8000000???? 10.0349390???? 70.0000000???? 90.0000000
    GS???????? 5???? 88.2000000????? 3.5637059???? 83.0000000???? 92.0000000
    P????????? 5??? 169.0000000???? 11.7898261??? 157.0000000??? 182.0000000
    §2.2 SAS 常量
    SAS中有多種常量:數值常量,字符常量,日期、時間、日期時間常量。
    ①數值常量:
    ?? 數值常量是一個十進制數值,在書寫上與數學書寫習慣一致。可分為:正負整數、小數和科學記數法。例如:
    ????? 數???????????? SAS書寫
    -3???? -3
    3.12???? 3.12
    18000000?? 1.8E6??? (大于1032-1必須用科學記數法表示)
    -0.0000125?? -1.23E-5
    ②字符常量:
    字符常量是用一對單引號或雙引號括起來的一串字符,其長度不超過200個字符。
    ③日期、時間、日期時間常量:
    這三類常量是將其對應的格式化值用單引號括起來,并且后跟類型符來表示的。
    l 日期值(DATE):表示為’ddmmmyy’D或’ddmmmyyyy’D。
    其中dd表示日,mmm表示月份,yy或yyyy表示年份。
    如:’24MAR03’D或24MAR2003’D
    l 時間值(TIME):表示為’hh:mm:ss’T或’hh:mm:ss. ss’T。
    其中hh表示小時,mm表示分鐘,ss表示分鐘。
    如:’14:30:25’T或’14:30:25. 87’T。
    l 日期時間值(DATATIME):表示為:’ddmmmyy:hh:mm:ss’DT。
    §2.3 SAS 變量
    一、變量及變量的命名
    變量:指具有相同性質的數據值的集合。變量是用變量名表示的。
    變量的命名:以字母或下畫線開頭,后跟字母或數字或下畫線的字符串。長度不多于8個字符。空格和特殊的字符($、@、#)不能出現在變量名中。
    二、變量的屬性
    SAS的變量分為兩類:數值型與字符型。每個變量都有長度、輸入格式、輸出格式和標號四種屬性。變量的屬性可以由首次出現的上下文給出定義,也可以用專用的語句定義明確給出。
    1、變量的取值范圍:
    數值型變量的值是數值,數值的范圍是:±10-307到±10 308。
    字符型變量的值是字符串,其取值范圍最多為200個字符。
    2、變量長度:
    變量長度是指SAS數據集中存儲它的每一個值的字節數。
    數值型變量的長度為:3"8個字節,默認為8個字節。
    字符型變量的的長度為:1"200個字節,默認為8個字節。
    3、變量的格式標號
    由于變量名長度不能超過8個字符,有時變量名受此限制而不能完全表達變量的含義。變量的格式標號是用來給變量定義一個標簽,這樣可以提高輸出變量值時的可讀性。用LABEL語句可以定義變量的格式標號。
    4、數據缺值(Missing Values)
    在數據分析中,經常會遇到數據不全或數據缺少的問題。SAS在數據分析時允許數據缺值,一般用“.”表示。
    5、變量的簡化表示
    SAS程序定義了完整的變量清單,可以在許多語句中使用縮寫變量清單,達到簡化書寫的目的。(參見P9表2.1)
    §2.5 SAS 函數
    如同其他高級語言一樣,SAS提供了豐富的標準函數,為計算帶來極大方便。SAS函數的書寫格式為:
    函數名(參數表)
    例如:SIN(X) ABS(X)
    一、數的參數
    函數的參數可以是常量、變量、表達式。函數的參數可以是一個或多個。如:
    sin(x)?? 一個參數
    max(x1,x2,x3,x4,x5)?? 多個參數
    參數書寫形式有如下三種:
    ① 將所有參數列入括號內,參數間用“,”隔開。如:max(x1,x2,x3,x4,x5)
    ② 將所有參數列入括號內,參數間用“空格”隔開,前面冠以OF,如:
    max(of?? x1?? x2?? x3?? x4?? x5)
    ③ 參數變量名的尾部是有序整數時,可寫成:(OF 變量名1-變量名n),如:
    max(of?? x1-x5)
    二、函數值
    除個別函數外,多數函數值的類型與其參數類型是一致的。
    三、關于函數的幾點說明
    ① 函數不能直接用在PUT語句中。
    ② 函數參數都有特定的取值范圍。
    ③ 除樣本統計函數外,多數函數不允許有缺值。
    ④ 對某些概率函數,若參數選擇不當,可能引起不收斂的問題。
    四、函數分類
    SAS函數的類型有:算術、數學、截尾、三角、樣本統計、概率函數。(略)
    §2.6 SAS表達式
    SAS表達式是由操作符和圓括號將各種操作數連接起來表示運算關系的的式子。表達式的值可能是:數值、字符串或缺失值。
    一、操作符
    ① 算術運算符
    **(乘方)?? *(乘)?? /(除)??? (加)?? -(減)
    ② 比較(關系)運算符
    = 或EQ??? 等于?? ^= 或 NE?? 不等于
    > 或GT??? 大于??? < 或LT??? 小于
    >=或GE??? 大于等于?? <=或LE?? 小于等于
    比較兩個量大小關系,構成比較表達式,運算結果為邏輯值。如:
    5>3 1??? X<2*Y??? ”ab c”>=”XYZ”
    ③ 邏輯運算符
    &;amp; 或 AND?? | 或 OR??? ^ 或NOT
    對邏輯值再進行運算,構成邏輯表達式,運算結果為邏輯值。如:
    5>3 1 AND X<2*Y
    邏輯運算真值表:
    邏輯值A 邏輯值B A AND B A OR B NOT A
    真 真 真 真 假
    真 假 假 真 假
    假 真 假 真 真
    假 假 假 假 真
    ④ 其它運算符
    <>??? 取兩個量中最大值?? 53<>34
    ><??? 取兩個量中最小值?? 53><34
    || 或!聯接兩個字符串?? ”abc” || ”de”
    二、操作符運算優先級
    在一個SAS表達式中,同時可能包含多個操作符,各種操作符的執行優先順序規定如下:
    規則1:括號內優先
    規則2:具有較高優先級的操作先執行,優先級次序如下:
    組1:**,正負號,NOT,<>,><
    組2:*,/
    組3: ,-
    組4:||
    組5:<,<=,=,^=,> ,> 組6:AND
    組7:OR
    規則3:組號小者優先組號高者,優先級先同時,依次從左向右執行。
    §2.7 SAS語句與程序
    SAS語句是以SAS關鍵字(保留字)開頭、后跟SAS名、特殊字符、操作符等組成,以分號結束。一個語句規定了一種操作或為系統提供默寫信息。
    一、SAS關鍵字:
    關鍵字是SAS系統已經賦于特定含義的一個單詞或單詞縮寫,是一個特定的字符串。在SAS語句中,多數語句是以關鍵字作為開頭的。
    二、SAS 名:
    SAS名是以字母或下畫線開頭,后跟字母或數字或下畫線的字符串。長度不多于8個字符。空格和特殊的字符($、@、#)不能出現在變量名中。
    另外,SAS保留了一些特殊的變量名并賦予特定的含義。這些變量都是以下畫線開頭和結尾,如:_N_表示當前數據步已經執行過的次數。
    在SAS語句中,可能出現的SAS名有:變量名、數據集名、過程名、輸出格式名、選項名、數組名、語句標號名、文件邏輯名等。
    三、SAS程序
    SAS程序是由若干SAS語句組成的有序特定的集合,典型的SAS程序都包括兩個部分:數據步和過程步,整個程序以RUN語句結束。
    同其它高級語言一樣,SAS程序由多行語句組成,其書寫格式如下:
    ① 語句可以在一行任意列開始書寫
    ② 一個語句可以寫一行或多行,但要以分號結束
    ③ 多條語句可以寫在一行,語句見必須以分號隔開(英文分號“;”)
    ④ SAS語句對英文字符不區分大小寫。
    ⑤ 建議將程序書寫成鋦齒形,便于程序閱讀與調試。

    總結

    以上是生活随笔為你收集整理的sas教程的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

    av中文在线观看 | 一区二区中文字幕在线观看 | 在线观看一区 | 婷婷中文字幕 | 最近日本韩国中文字幕 | 亚洲精品美女 | 国产精品日韩欧美一区二区 | 香蕉视频免费看 | 国产 一区二区三区 在线 | 91大神免费在线观看 | 日韩精品中文字幕在线 | 天天se天天cao天天干 | 在线观看色网站 | 免费h精品视频在线播放 | 久久99热这里只有精品国产 | 日韩精品免费一区二区在线观看 | 成人午夜免费剧场 | 天堂黄色片| 99爱精品视频 | 五月香视频在线观看 | 24小时日本在线www免费的 | 美女视频又黄又免费 | 国产精品视频内 | 激情综合色播五月 | 特级毛片网 | 精品视频免费久久久看 | 丁香5月婷婷 | 2021国产在线 | 久久蜜臀av | 在线观看中文字幕视频 | 久久久久亚洲精品中文字幕 | 欧美一二在线 | 日本在线成人 | 91麻豆精品久久久久久 | 免费看的av片 | 啪啪免费观看网站 | 日韩三级久久 | 久久免费观看视频 | 精品国产免费观看 | 一区二区三区电影 | 国产资源在线免费观看 | 黄色精品一区二区 | 丁香花在线观看视频在线 | av中文字幕在线播放 | 免费看av在线 | 在线视频 亚洲 | 日韩免费看的电影 | 成人久久久久久久久久 | 欧美最猛性xxxxx亚洲精品 | 不卡视频一区二区三区 | 精品一二三区 | 国产精品一区免费观看 | 999成人网 | 中文字幕色网站 | 欧美日韩不卡在线 | 色播五月激情五月 | 狠狠操欧美 | 成人免费xxx在线观看 | 亚洲激情小视频 | 日本精品久久久久中文字幕 | 日韩在线免费小视频 | 亚洲一级黄色 | 免费看片网址 | 天天舔天天射天天操 | a视频在线观看免费 | 日韩中文字幕免费在线观看 | 成人a免费看 | 成人av免费看 | 91精选在线| 麻豆极品 | 五月婷婷丁香色 | 天天操·夜夜操 | 久久久久国产精品视频 | 欧美日韩二区三区 | av黄色av | 91视频麻豆 | 日韩中文字幕电影 | 国产亚洲成av人片在线观看桃 | 国产成人黄色网址 | 免费av在 | 99九九免费视频 | 亚洲人视频在线 | 久久久久久久久久久久久影院 | 欧美日韩视频免费看 | 99久久久国产精品免费观看 | 国产精品一区二区在线观看免费 | 四虎影视成人永久免费观看亚洲欧美 | 成人av在线资源 | 亚洲欧美日韩在线一区二区 | 日韩av影视在线观看 | 国产亚洲成人精品 | 天天干天天看 | 国产资源在线播放 | 久久影视网 | 丁香六月婷婷开心 | 亚洲成人动漫在线观看 | 午夜视频导航 | 黄色一级大片在线免费看产 | 热热热热热色 | 88av色| 91片在线观看| 日韩av偷拍 | 中文字幕一区二区三区四区久久 | 成人毛片100免费观看 | 国产成人久久av免费高清密臂 | 久久97久久97精品免视看 | 狠狠干婷婷色 | 久久免费久久 | 99精品国产99久久久久久97 | 国产精品a久久 | 日韩三区在线观看 | 96精品视频 | 欧美日韩精品在线视频 | 久久国产精品一二三区 | 超碰日韩在线 | 亚洲精品国产精品国 | 色综合久久66 | 亚洲片在线观看 | 亚洲aⅴ免费在线观看 | 久久久人人爽 | 国产精品久久一区二区无卡 | 最新av观看| 日韩精品在线观看av | 欧美污在线观看 | 欧洲一区二区三区精品 | 国内亚洲精品 | 欧美精品成人在线 | 成全在线视频免费观看 | av超碰免费在线 | 天天爽夜夜爽精品视频婷婷 | 亚洲欧洲一区二区在线观看 | 99久久精品国产网站 | 青青草国产成人99久久 | 一区二区久久久久 | 2021av在线 | 久久成人毛片 | 成人avav| 午夜成人免费电影 | 丝袜少妇在线 | 久久撸在线视频 | 国产精品视频免费在线观看 | 久久久久国产精品视频 | 亚洲一区二区高潮无套美女 | 91九色蝌蚪视频网站 | 色偷偷av男人天堂 | 亚洲另类视频在线观看 | 五月宗合网 | 91中文字幕在线 | 中文字幕在线观看一区二区 | 国内视频在线 | 中文字幕乱码电影 | av一级一片 | 国产免费久久av | 亚洲天堂毛片 | 久操久| 91av视频| 精品久久久久久久久久久院品网 | 亚洲视频高清 | 欧美专区日韩专区 | 精品电影一区 | a天堂在线看 | 国产精品毛片一区二区在线 | 国产精品自产拍在线观看 | 免费的黄色av | 亚洲aⅴ在线观看 | 在线观看成人国产 | 最近免费中文字幕大全高清10 | 91亚洲精品乱码久久久久久蜜桃 | 免费在线电影网址大全 | 日韩二区在线播放 | 国产精品69av | 日本中文字幕在线 | 久久免费播放视频 | 综合久久网站 | 国产精品一区二区在线看 | 美女在线免费视频 | 99免费精品 | 久av电影| 美女视频久久久 | 成人av片免费看 | 国产精品99久久久久人中文网介绍 | 五月天婷婷狠狠 | 88av网站| 天天操天天爽天天干 | 免费高清av在线看 | 欧美极品少妇xxxxⅹ欧美极品少妇xxxx亚洲精品 | 国产精品尤物 | 久久亚洲福利视频 | 99爱视频| 日韩电影在线观看一区二区三区 | 国产视频一区在线播放 | 日韩美视频 | 亚洲天堂网在线观看视频 | 97精品欧美91久久久久久 | 久草视频在| 一区二区成人国产精品 | 91九色国产在线 | 精品在线观看一区二区三区 | 久久久久久久国产精品 | 亚洲国产高清在线观看视频 | 91成人国产 | 亚洲国产日本 | 免费视频色 | 久久久精品国产一区二区 | 九九免费精品视频 | 久久免费在线 | 亚洲国产日韩在线 | 国产精品福利午夜在线观看 | 欧洲一区精品 | avav片| 91成人黄色 | 国产精品久久久久久久久久久免费 | 天堂在线一区二区 | 激情综合交 | 美女免费网视频 | 91麻豆精品一区二区三区 | 免费观看v片在线观看 | 国内精品国产三级国产aⅴ久 | 人人舔人人舔 | 国内精品久久久久久中文字幕 | 国产伦精品一区二区三区免费 | 欧美成人一二区 | 国产精品免费高清 | 四虎在线观看精品视频 | 亚洲aⅴ乱码精品成人区 | 色美女在线 | 午夜精品一区二区三区在线观看 | 欧美一级电影 | 国内少妇自拍视频一区 | 国产香蕉97碰碰碰视频在线观看 | 久久高清国产视频 | 亚洲国产精品电影在线观看 | 国产 字幕 制服 中文 在线 | 在线精品亚洲 | 欧美另类交在线观看 | 天天天天色综合 | 在线观看成人一级片 | 国产精品每日更新 | 日韩在线高清免费视频 | 国产v在线播放 | 九九免费在线观看视频 | 夜色.com | 蜜桃视频精品 | 二区在线播放 | 99久久免费看 | 最近最新中文字幕 | 日日夜夜人人天天 | 国产日产精品一区二区三区四区 | 五月婷婷狠狠 | 国产视频欧美视频 | 我要看黄色一级片 | 91国内在线 | 日韩国产高清在线 | 国产一区观看 | 99视| 69av视频在线 | 天天舔天天搞 | 一区二区视频在线免费观看 | 亚洲一区二区三区毛片 | 久久久久久久久久久久电影 | 果冻av在线 | 亚洲aⅴ在线 | 成人免费看电影 | 91xav| 欧美了一区在线观看 | 欧美日韩高清在线 | 亚洲精品乱码久久久久久高潮 | 美女久久久久久久久久久 | 午夜久久成人 | 色网站在线观看 | 欧美色精品天天在线观看视频 | 久久午夜电影 | 在线黄色免费av | 久久久精品国产一区二区三区 | 日韩免费一级a毛片在线播放一级 | 97色在线观看免费视频 | 久久国产一区二区三区 | 日韩三区在线 | 久久久久久久免费看 | 婷婷丁香国产 | 亚洲国产精品va在线看黑人动漫 | 色婷五月天 | 精品免费视频123区 午夜久久成人 | 国产 成人 久久 | 色婷婷综合在线 | 99久久99| 在线看免费| 久久这里只有精品首页 | 特级西西444www大胆高清无视频 | 国产一区欧美日韩 | 青青草国产在线 | 综合久久久久 | 天天操天天玩 | 日夜夜精品视频 | 国语精品视频 | 最新av在线播放 | 国产精品久久久久久一区二区三区 | 亚洲精品美女免费 | 91精品啪在线观看国产81旧版 | 日韩欧美国产视频 | 欧美精品久久99 | 成人在线你懂得 | 日韩r级电影在线观看 | 天天久久夜夜 | 天天天综合| 狠狠色伊人亚洲综合网站野外 | 亚洲黄色激情小说 | 99在线观看免费视频精品观看 | 91爱爱电影 | 国精产品一二三线999 | 亚洲最快最全在线视频 | 久草视频首页 | 免费av的网站 | 天天爱天天射 | 欧美久草网 | 精品亚洲二区 | 成人一级视频在线观看 | 黄色精品一区 | 国产香蕉av| 久久久久久久久影视 | а中文在线天堂 | 久青草视频在线观看 | 亚洲a成人v | 黄色录像av| 手机av在线免费观看 | 天天干天天干天天色 | 国产无套精品久久久久久 | 最近中文字幕 | 成人久久国产 | 激情综合网婷婷 | 色噜噜狠狠色综合中国 | 不卡的一区二区三区 | 综合久久2023 | 久久免费的精品国产v∧ | 日本精品视频一区 | 久草在线中文视频 | 中文字幕国产一区二区 | 久久伊人91 | 日本精品视频在线 | 五月婷影院 | 五月婷婷视频在线观看 | 色婷婷在线播放 | 亚洲人成在线电影 | 深夜福利视频在线观看 | 日日夜夜精品网站 | 香蕉视频国产在线观看 | 日本大尺码专区mv | 久久久久久蜜av免费网站 | 91成人精品一区在线播放 | 欧美精品免费一区二区 | a视频在线 | 丁香av| 午夜视频福利 | www免费黄色 | 精品亚洲男同gayvideo网站 | 欧美日韩国产xxx | 日韩中文免费视频 | 久久久国产精品视频 | 中文在线字幕免费观看 | 日本最新高清不卡中文字幕 | 日韩精品在线免费播放 | 国产精品一区二区三区在线播放 | 日日夜夜免费精品 | 夜夜躁日日躁狠狠久久av | 亚洲免费黄色 | 在线观看免费黄色 | 中文av一区二区 | 久久精品国产免费 | 黄色三级在线观看 | 中文字幕一区在线 | 不卡日韩av | 亚洲欧洲国产视频 | 国产精品免费一区二区三区在线观看 | 国产精品免费观看在线 | 天天爱综合 | 亚洲婷婷伊人 | 97网站| 在线中文字幕一区二区 | 中文字幕日本特黄aa毛片 | 免费观看av网站 | 国产成人精品一区二区三区网站观看 | 天天干夜夜想 | 91少妇精拍在线播放 | 超碰伊人网 | 久久久国际精品 | 一级α片| 在线色网站 | 国产精品久久久久久久久久免费 | 中文视频在线看 | 国产在线色站 | 免费日韩三级 | 欧美夫妻性生活电影 | 高潮久久久久久久久 | 国产区网址 | 久久精品视频在线 | 欧美精品免费一区二区 | 黄色一级大片在线免费看产 | www久久99 | 国内外激情视频 | 操操操日日日干干干 | 精品一区二区在线免费观看 | 九色视频自拍 | 欧美日韩国产一区二区三区在线观看 | 欧美污网站| 蜜臀久久99静品久久久久久 | av高清影院| 国产高h视频| av一级二级| 亚洲成人精品在线 | av黄色免费在线观看 | 中文字幕av在线播放 | 九草视频在线 | 丁香六月婷婷激情 | 久久久久久久久久久综合 | 国产精品久久片 | 婷婷久久亚洲 | 国产一级高清 | 黄色成人av | 国产一区二区在线免费 | 久久久亚洲精品 | 日本精品视频网站 | 亚州精品国产 | 狠狠躁天天躁 | 国产a级片免费观看 | 天天翘av| 久久综合狠狠综合久久狠狠色综合 | 久久精品国产久精国产 | 欧美日韩一区二区三区在线免费观看 | 亚洲va欧美va国产va黑人 | 亚洲成aⅴ人在线观看 | 国产成人久久77777精品 | 99免费在线视频 | 91精品久久久久久粉嫩 | 激情电影影院 | 中文字幕免费高清在线 | 成人在线播放视频 | 欧美在线aa | 99热官网 | 8090yy亚洲精品久久 | 国产精品久久久久久久久久久久午夜 | 色婷婷激情网 | 中文字幕 在线 一 二 | 国产资源在线观看 | 国产黄在线 | 97精产国品一二三产区在线 | 精品中文字幕在线播放 | 香蕉视频导航 | 午夜电影一区 | 在线免费视频你懂的 | 亚洲欧美日韩精品一区二区 | 日本夜夜草视频网站 | 亚洲成av人片一区二区梦乃 | 国内精品久久久久久中文字幕 | 色老板在线视频 | 香蕉影视| 国产精品久久久久一区 | 国产一区二区在线免费 | a黄色片| 99视频在线 | 欧美综合国产 | www.五月婷婷 | 国产黄色大全 | 黄色一级大片免费看 | 免费中午字幕无吗 | 精品国产一区二区三区日日嗨 | 国产精品手机在线播放 | 天天色天天搞 | 日韩中文字幕免费 | 九九99靖品 | 日批在线观看 | 午夜av在线电影 | 韩国一区二区三区在线观看 | 精品在线你懂的 | 亚洲在线视频免费 | 91视频专区| 中文字幕在线精品 | 久久久精品国产一区二区电影四季 | 免费a v在线| 欧美一级性生活视频 | 中国美女一级看片 | 黄色小网站免费看 | 国产在线精 | 久久久久久伊人 | 国产精品一区二区麻豆 | 奇米先锋 | 91精品在线观看入口 | 亚洲2019精品 | 免费福利视频网站 | 国产精品久久久久久久久软件 | 国产精品亚洲人在线观看 | 国产在线a视频 | freejavvideo日本免费 | 91黄色视屏| 2022中文字幕在线观看 | 欧美精品九九99久久 | 久久久久黄 | 色中文字幕在线观看 | 久久精品视频在线免费观看 | 日韩av三区 | 日韩精品久久一区二区三区 | 成人在线观看免费 | 91av大全| 99草视频 | 青春草免费在线视频 | 亚洲成人999 | av.com在线| 成人综合免费 | 久色免费视频 | 国产裸体永久免费视频网站 | 九九精品久久久 | 在线视频 精品 | 国产免费嫩草影院 | 91福利区一区二区三区 | 久久精品电影院 | 丁香婷五月 | 国产 色 | 亚洲最大激情中文字幕 | 九九免费在线观看 | 日韩69视频 | 成年人精品 | 国产精品激情 | 欧洲一区二区在线观看 | 99精品一级欧美片免费播放 | 日韩精品视频在线观看网址 | 黄在线免费看 | 久久久精品 一区二区三区 国产99视频在线观看 | 免费看搞黄视频网站 | 亚洲综合欧美日韩狠狠色 | 日韩欧美网址 | 91重口视频 | 国产色视频一区二区三区qq号 | av黄色在线| 久久久久亚洲精品成人网小说 | 亚洲午夜av | 日韩av免费观看网站 | 又粗又长又大又爽又黄少妇毛片 | 中文字幕中文字幕在线一区 | 91精品视频在线观看免费 | 狠狠成人 | 亚洲国内精品在线 | 在线视频1卡二卡三卡 | 一本—道久久a久久精品蜜桃 | 久久国产精品色婷婷 | 亚洲 综合 精品 | 中文字幕成人在线观看 | 色综合久 | 中文字幕av在线电影 | 亚洲欧美日韩国产一区二区 | 亚州五月 | 91看片网址| 成年人黄色大全 | 最新一区二区三区 | 国产视频中文字幕 | 国产精品丝袜久久久久久久不卡 | 久精品在线 | 国产91在线播放 | 欧美日韩在线视频一区二区 | 国产精品久久久毛片 | 久久综合色影院 | 在线亚洲激情 | 91精品国产成人www | 日韩av网页 | 色综合久久久久综合体 | 国产精品v欧美精品 | 国产精品成| 久久久久亚洲国产 | 久久综合婷婷综合 | 中文字幕在线字幕中文 | 成年人免费在线看 | 狠狠干天天色 | 99精品热视频只有精品10 | 国产精品久久久久亚洲影视 | 欧美精品一区二区在线播放 | 精品一区二区久久久久久久网站 | 在线午夜电影神马影院 | 中文字幕精品三级久久久 | 精品亚洲欧美一区 | 国产成人精品福利 | 成人免费看视频 | 草免费视频| 亚洲精品国产综合99久久夜夜嗨 | 成人片在线播放 | 中文在线免费观看 | 亚洲欧美日本A∨在线观看 青青河边草观看完整版高清 | 欧美夫妻生活视频 | 日韩三级视频在线看 | 亚洲综合最新在线 | 天天综合久久综合 | 午夜av色 | 亚洲日本国产精品 | 亚洲一级久久 | 日韩免费观看一区二区 | 亚洲国产成人精品在线观看 | 中文字幕在线资源 | 中文字幕国产一区 | 久久久国产一区二区 | 日韩免费一级a毛片在线播放一级 | 在线播放 日韩专区 | 日韩一级网站 | 91在线网站| 欧美 国产 视频 | 国产一区二区综合 | 九九精品在线观看 | 欧美日韩三区二区 | 精品夜夜嗨av一区二区三区 | 欧洲精品久久久久毛片完整版 | 日日夜夜精品 | 日韩av区 | 国产精品久久久毛片 | 国产尤物在线观看 | 精品高清视频 | 欧美一区二区在线免费看 | 爱爱一区 | 日韩久久精品一区二区三区下载 | 国产精品 美女 | 国产无套精品久久久久久 | 亚洲区另类春色综合小说校园片 | 欧美日韩国产一区二 | 五月天婷亚洲天综合网精品偷 | 久久久久欠精品国产毛片国产毛生 | 婷婷六月久久 | 永久免费毛片在线观看 | 香蕉视频日本 | 91在线操| 国产精品美女视频网站 | 欧美福利久久 | 狠狠色狠狠色综合日日小说 | 色偷偷88888欧美精品久久久 | 欧美 日韩 性 | 鲁一鲁影院 | 91九色国产视频 | 在线电影日韩 | 99c视频在线| 久久视频在线视频 | 91精品对白一区国产伦 | 久久er99热精品一区二区三区 | 欧美91在线 | 久久久国产精华液 | 久久免费福利 | 国产一区二区在线免费播放 | 成人a级大片 | 在线视频日韩欧美 | 在线播放 日韩专区 | 中文字幕日韩无 | 97精品久久| 欧美日韩一区二区三区不卡 | www.天天干| 在线观看黄色大片 | 天天草夜夜 | 亚洲欧洲视频 | 精品视频在线免费观看 | 欧美一级日韩免费不卡 | 久久综合桃花 | 日韩高清在线一区二区 | 日韩在线高清免费视频 | 一本一本久久a久久精品牛牛影视 | 天天操夜夜拍 | 国产日韩精品欧美 | 手机看片国产日韩 | 97久久久免费福利网址 | 三级av免费观看 | 国产专区在线播放 | 亚洲欧美怡红院 | 最近中文字幕mv免费高清在线 | 激情久久五月天 | 日韩网站在线 | 91福利视频免费观看 | 人人插人人 | 亚洲aⅴ一区二区三区 | 91久久精品一区二区二区 | 人人搞人人干 | 一区二区在线影院 | 国内揄拍国内精品 | 岛国一区在线 | 精品国产诱惑 | 欧美日韩精品久久久 | 欧美日韩视频一区二区 | 亚洲第一区精品 | 久久一区二区三区超碰国产精品 | 久精品一区 | 丁香九月婷婷 | 国产精品一区二区美女视频免费看 | 丁香五月亚洲综合在线 | 91精品免费在线观看 | 天天操狠狠操 | 色资源二区在线视频 | 综合在线色 | 不卡av电影在线观看 | 狠狠狠色狠狠色综合 | 日韩欧美xxx | 久久免费国产电影 | 欧美日韩国产精品久久 | 国产亚洲精品久久19p | 国产精品一区二区三区久久 | 国产亚洲在 | 久久精品视频2 | 久久久久久久久毛片 | 久草在线观看 | 亚洲精品一区二区三区高潮 | 国产糖心vlog在线观看 | 色狠狠一区二区 | 黄色软件大全网站 | 日韩精品中文字幕久久臀 | 在线中文字幕电影 | 亚洲精品免费观看 | 欧美亚洲另类在线视频 | 99国产精品| 婷婷看片| 91在线影视 | av成人黄色| 香蕉视频一级 | av色图天堂网 | 久久久久五月 | 亚洲精品乱码白浆高清久久久久久 | 久久久久久久久久伊人 | 日韩日韩日韩日韩 | 中文字幕在线专区 | 精品黄色视| 欧美资源 | 久久久久久高清 | 五月激情婷婷丁香 | 亚洲免费视频在线观看 | 亚洲欧美国产精品va在线观看 | 日韩欧美网站 | 91重口视频 | 伊人天天色 | 久久久高清一区二区三区 | 日韩欧美国产视频 | 超碰精品在线观看 | 射久久久| 欧美男男激情videos | 国产精品久久久久久久久久免费 | 成人精品国产免费网站 | 国产精品99久久久久久久久 | 久久国产日韩 | 在线观看免费视频你懂的 | 亚洲午夜久久久久久久久久久 | 久久视频在线视频 | av软件在线观看 | 91爱爱视频 | 国产精品成人a免费观看 | 亚洲欧洲精品视频 | 91网站观看 | 91视频大全 | 久久伊人婷婷 | 亚洲精品999 | 黄色91在线 | 欧美日韩国产一区 | 日韩免费久久 | 精品国产一区二区三区四区vr | 久久不射电影网 | 成人黄色小说视频 | 五月开心综合 | 午夜精品一区二区三区在线视频 | 五月天色中色 | 激情五月婷婷网 | 久久久免费看视频 | 免费观看的黄色片 | 91看片在线播放 | 国产亚洲亚洲 | 国产亚洲精品日韩在线tv黄 | 一区二区三区久久 | 中文乱码视频在线观看 | 国产午夜在线 | 国产成年人av | 91在线视频免费观看 | 色天堂在线视频 | 日韩精品一区二区三区免费视频观看 | 高清av中文在线字幕观看1 | 久久视频免费在线观看 | 久久综合九色综合欧美就去吻 | 久久久国产成人 | av网站在线观看播放 | 国产一级在线 | 国产精品成久久久久三级 | 精品高清美女精品国产区 | 国产精品成人一区二区三区吃奶 | 欧美日韩国语 | 狠狠色噜噜狠狠狠 | av免费高清观看 | 97精品国产91久久久久久 | 日韩精品一区二区三区第95 | 色国产视频 | 日韩视频免费观看高清完整版在线 | 国产精品99久久久久的智能播放 | 免费久久精品视频 | 欧美精品一区二区三区四区在线 | 欧美另类视频 | 99视频网站 | 国产精品毛片一区视频播不卡 | 丁香六月婷 | 成年人在线观看免费视频 | 国产精品高潮呻吟久久av无 | 国产视频一区二区三区在线 | 亚洲国产手机在线 | 人人擦 | 国产成人在线免费观看 | 九九九九精品 | 国产精品99久久久久 | 干天天 | 成人资源在线 | 国产成人亚洲在线观看 | 日本精品中文字幕在线观看 | www.香蕉视频| 日本成人a| 国产69精品久久久久久 | 天天干天天天 | 精品免费久久 | 亚洲综合色激情五月 | 久艹在线播放 | 夜夜骑日日操 | 波多野结衣一区二区三区中文字幕 | 中文字幕免费一区 | 欧美日韩裸体免费视频 | 四虎国产免费 | av电影免费在线播放 | 91精品国产乱码在线观看 | 国产精品一区二区三区在线看 | 日韩免费观看高清 | 婷婷中文字幕在线观看 | 亚洲国产天堂av | 在线国产视频观看 | 国产精品免费在线播放 | 韩国在线一区二区 | 国产综合婷婷 | 奇米影视8888在线观看大全免费 | 美女福利视频在线 | 综合激情婷婷 | 9在线观看免费高清完整 | 亚洲精品美女视频 | 99精品偷拍视频一区二区三区 | 亚洲日本在线一区 | 午夜精品电影一区二区在线 | 国产精品午夜免费福利视频 | 91久久国产自产拍夜夜嗨 | 超碰在线免费福利 | 91视频黄色 | 免费又黄又爽 | 欧美激情精品久久 | 日韩高清一区在线 | 一区二区高清在线 | 视频 天天草| 欧美视频99 | 亚洲国产无 | 日韩高清dvd| 成人久久久电影 | 国产精品永久在线观看 | 久久成人麻豆午夜电影 | 91在线视频网址 | 亚洲va天堂va欧美ⅴa在线 | 免费观看www小视频的软件 | 久久精品三级 | 天天操天天草 | 日韩在线免费视频 | 香蕉色综合 | 69国产精品视频 | 日韩大片免费在线观看 | av免费黄色 | 成人午夜免费福利 | 国产黑丝袜在线 | 特级西西444www大胆高清无视频 | 丁香视频在线观看 | 手机在线日韩视频 | 国产高清小视频 | 国产精品免费观看国产网曝瓜 | 亚洲免费资源 | 午夜精品久久久久久久99婷婷 | 国产高清永久免费 | 不卡av电影在线 | 久久伊人综合 | 中文字幕在线播放一区二区 | 久久99久久99精品免观看软件 | 蜜桃av人人夜夜澡人人爽 | 中文字幕有码在线 | 国产精品网站一区二区三区 | 成年人在线观看 | 黄色三级免费看 | 色吊丝在线永久观看最新版本 | 日日干夜夜草 | 91在线看黄| www色网站| 中文字幕在线字幕中文 | 五月婷婷综合在线视频 | 日韩欧美视频一区 | 成人全视频免费观看在线看 | 婷婷精品国产欧美精品亚洲人人爽 | 欧美专区日韩专区 | 久久精品一二三区白丝高潮 | 激情婷婷久久 | 97碰在线视频 | 欧美午夜性 | 国语精品免费视频 | 国产精品免费在线 | 一级性生活片 | 国产在线a免费观看 | 激情五月婷婷综合 | 丁香六月在线观看 | 欧美做受高潮 | 一区免费视频 | 国产又粗又猛又爽又黄的视频先 | 国产黄色av网站 | 亚州精品天堂中文字幕 | 亚洲女欲精品久久久久久久18 | 97色在线观看免费视频 | 久久久久久久久久影视 | 国产成人香蕉 | 婷婷网站天天婷婷网站 | 精品毛片一区二区免费看 | 国产免费久久精品 | 中文字幕永久免费 | 日本电影久久 | 国产精品久久久久亚洲影视 | 有码中文在线 | 丁香五月亚洲综合在线 | 成人国产精品av | 91九色蝌蚪视频 | 丁香电影小说免费视频观看 | 三级av免费看 | 亚洲国产精品成人va在线观看 | 婷婷丁香在线观看 | 亚洲天堂社区 | 亚洲一级在线观看 | 欧美日产在线观看 | 中文在线www | 在线视频你懂得 | 日韩欧美综合在线视频 | 久久成人在线视频 | 中文字幕高清有码 | 五月婷婷中文字幕 | 成人天堂网 | 香蕉免费 | 日本特黄一级 | 97在线成人 | 日韩av一区二区三区四区 | 毛片在线网 | 国产高清不卡av | 国内精品亚洲 | 国产精品视频专区 | 欧美 国产 视频 | 国产精品高潮呻吟久久久久 | 免费看日韩片 | 高潮毛片无遮挡高清免费 | 美女激情影院 | 天天干天天射天天操 | 99视频在线观看一区三区 | 在线天堂中文在线资源网 | 中文乱幕日产无线码1区 | 免费黄色在线网址 | 午夜成人免费影院 | 成人 国产 在线 | 午夜精品福利一区二区三区蜜桃 | 最新中文字幕在线播放 | 成人a免费看 | av中文在线影视 | 日本久久精品 | 久久精品在线 | 日日夜夜免费精品 | 天天操天| 毛片永久免费 | 亚洲精选视频免费看 | 亚洲三级影院 | 久久乱码卡一卡2卡三卡四 五月婷婷久 | 91精品国产92久久久久 | 一级免费黄视频 | 国产视频不卡一区 | 免费在线观看av网站 | 欧美天堂久久 | 亚洲国产三级在线观看 | 午夜精品视频一区二区三区在线看 | 超碰人人在线观看 | 黄色免费高清视频 | 五月在线 | 国产精品乱码久久久 | 超碰人人99| 亚洲视频综合在线 | 一级黄色片在线观看 | 91视视频在线直接观看在线看网页在线看 | 久久久久久片 | 亚洲国产婷婷 | 日韩精品中文字幕有码 | 亚洲天堂va | 又黄又刺激又爽的视频 | 亚州精品视频 | 最近2019年日本中文免费字幕 | 日韩字幕在线观看 | 色诱亚洲精品久久久久久 | 久久avav| 在线视频手机国产 | 麻豆va一区二区三区久久浪 | 欧美精品小视频 | 永久免费av在线播放 | av片在线观看 | 91精选在线观看 | 亚洲精品视频大全 | av超碰在线 | 在线a视频免费观看 | 国产精品va | 国产精品美女在线观看 |