SAS常用的随机数函数简介文档
SAS函数汇总范文
SAS函数汇总范文SAS（Statistical Analysis System）是一种统计分析系统，它提供了一套完整的分析工具和编程语言，用于数据的处理、管理、分析和展示。
SAS提供了多种函数，用于实现各种统计分析方法和数据处理操作。
以下是一些常用的SAS函数汇总：1.数据管理函数-FORMAT函数：用于为数据指定格式，如日期格式、货币格式等。
-LENGTH函数：用于定义变量的长度。
-SORT函数：用于对数据集进行排序。
-DROP函数：用于删除数据集中的一些变量。
-MERGE函数：用于将两个或多个数据集按照相同的变量进行合并。
-DELETE函数：用于删除数据集中的一些观测。
-RENAME函数：用于重命名数据集或变量。
2.数据转换函数-SUBSTR函数：用于提取字符串的子串。
-TRANSLATE函数：用于将字符串中的一些字符替换为其他字符。
-UPCASE函数：将字符串转换为大写。
-LOWCASE函数：将字符串转换为小写。
-PROPCASE函数：将字符串转换为首字母大写，其余字母小写。
-COMPRESS函数：用于去除字符串中的空格或指定字符。
-CAT函数：用于将多个字符串连接成一个字符串。
-INPUT函数：将字符型变量转换为数值型变量。
-PUT函数：将数值型变量转换为字符型变量。
3.数据统计函数-MEAN函数：计算变量的平均值。
-N函数：计算变量的非缺失观测数。
-NMISS函数：计算变量的缺失观测数。
-STD函数：计算变量的标准差。
-VAR函数：计算变量的方差。
-MIN函数：计算变量的最小值。
-MAX函数：计算变量的最大值。
-SUM函数：计算变量的总和。
-FREQ函数：计算变量的频数。
-PROB函数：计算变量的概率密度或累积分布。
4.数据分析函数-REG函数：进行回归分析。
-TTEST函数：进行独立样本t检验。
-ANOVA函数：进行方差分析。
-CORR函数：计算变量之间的相关系数。
-FACTOR函数：进行因子分析。
SAS函数介绍
Functions and CALL Routines by CategoryCategories and Descriptions of Functions Category Function DescriptionArray DIM Returns the number of elements in an arrayHBOUND Returns the upper bound of an arrayLBOUND Returns the lower bound of an arrayBitwise Logical Operations BAND Returns the bitwise logical AND of two argumentsBLSHIFT Returns the bitwise logical left shift of two argumentsBNOT Returns the bitwise logical NOT of an argumentBOR Returns the bitwise logical OR of two argumentsBRSHIFT Returns the bitwise logical right shift of two argumentsBXOR Returns the bitwise logical EXCLUSIVE OR of two argumentsCharacter String Matching CALL RXCHANGE Changes one or more substrings that match a patternCALL RXFREE Frees memory allocated by other regular expression_r(RX) functions and CALL routinesCALL RXSUBSTR Finds the position, length, and score of a substring that matches a patternRXMA TCH Finds the beginning of a substring that matches a pattern and returns a valueRXPARSE Parses a pattern and returns a valueCharacter BYTE Returns one character in the ASCII or the EBCDIC collating sequenceCOLLATE Returns an ASCII or EBCDIC collating sequence character stringCOMPBL Removes multiple blanks from a character stringCOMPRESS Removes specific characters from a character stringDEQUOTE Removes quotation marks from a character valueINDEX Searches a character expression for a string of charactersINDEXC Searches a character expression for specific charactersINDEXW Searches a character expression for a specified string as a wordLEFT Left aligns a SAS character expressionLENGTH Returns the length of an argumentLOWCASE Converts all letters in an argument to lowercaseMISSING Returns a numeric result that indicates whether the argument contains a missing valueQUOTE Adds double quotation marks to a character valueRANK Returns the position of a character in the ASCII or EBCDIC collating sequenceREPEAT Repeats a character expressionREVERSE Reverses a character expressionRIGHT Right aligns a character expressionSCAN Selects a given word from a character expressionSOUNDEX Encodes a string to facilitate searchingSPEDIS Determines the likelihood of two words matching, expressed as the asymmetric spelling distance between the two wordsSUBSTR (left of =) Replaces character value contentsSUBSTR (right of =) Extracts a substring from an argumentTRANSLATE Replaces specific characters in a character expressionTRANWRD Replaces or removes all occurrences of a word in a character stringTRIM Removes trailing blanks from character expressions and returns one blank if the expression is missingTRIMN Removes trailing blanks from character expressions and returns a null string (zero blanks) if the expression is missingUPCASE Converts all letters in an argument to uppercaseVERIFY Returns the position of the first character that is unique to an expressionDBCS KCOMPARE Returns the result of a comparison of character stringsKCOMPRESS Removes specific characters from a character stringKCOUNT Returns the number of double-byte characters in a stringKINDEX Searches a character expression for a string of charactersKINDEXC Searches a character expression for specific charactersKLEFT Left aligns a SAS character expression by removing unnecessary leading DBCS blanks and SO/SIKLENGTH Returns the length of an argumentKLOWCASE Converts all letters in an argument to lowercaseKREVERSE Reverses a character expressionKRIGHT Right aligns a character expression by trimming trailing DBCS blanks and SO/SIKSCAN Selects a given word from a character expressionKSTRCA T Concatenates two or more character stringsKSUBSTR Extracts a substring from an argumentKSUBSTRB Extracts a substring from an argument based on byte positionKTRANSLATE Replaces specific characters in a character expressionKTRIM Removes trailing DBCS blanks and SO/SI from character expressionsKTRUNCATE Truncates a numeric value to a specified lengthKUPCASE Converts all single-byte letters in an argument to uppercaseKUPDATE Inserts, deletes, and replaces character value contentsKUPDATEB Inserts, deletes, and replaces character value contents based on byte unitKVERIFY Returns the position of the first character that is unique to an expressionDate and Time DATDIF Returns the number of days between two datesDA TE Returns the current date as a SAS date valueDA TEJUL Converts a Julian date to a SAS date valueDA TEPART Extracts the date from a SAS datetime valueDA TETIME Returns the current date and time of day as a SAS datetime valueDAY Returns the day of the month from a SAS date valueDHMS Returns a SAS datetime value from date, hour, minute, and secondHMS Returns a SAS time value from hour, minute, and second valuesHOUR Returns the hour from a SAS time or datetime valueINTCK Returns the integer number of time intervals in a given time spanINTNX Advances a date, time, or datetime value by a given interval, and returns a date, time, or datetime valueJULDATE Returns the Julian date from a SAS date valueJULDATE7 Returns a seven-digit Julian date from a SAS date valueMDY Returns a SAS date value from month, day, and year valuesMINUTE Returns the minute from a SAS time or datetime valueMONTH Returns the month from a SAS date valueQTR Returns the quarter of the year from a SAS date valueSECOND Returns the second from a SAS time or datetime valueTIME Returns the current time of dayTIMEPART Extracts a time value from a SAS datetime valueTODAY Returns the current date as a SAS date valueWEEKDAY Returns the day of the week from a SAS date valueYEAR Returns the year from a SAS date valueYRDIF Returns the difference in years between two datesYYQ Returns a SAS date value from the year and quarterDescriptive Statistics CSS Returns the corrected sum of squaresCV Returns the coefficient of variationKURTOSIS Returns the kurtosisMAX Returns the largest valueMEAN Returns the arithmetic mean (average)MIN Returns the smallest valueMISSING Returns a numeric result that indicates whether the argument contains a missing valueN Returns the number of nonmissing valuesNMISS Returns the number of missing valuesORDINAL Returns any specified order statisticRANGE Returns the range of valuesSKEWNESS Returns the skewnessSTD Returns the standard deviationSTDERR Returns the standard error of the meanSUM Returns the sum of the nonmissing argumentsUSS Returns the uncorrected sum of squaresV AR Returns the varianceExternal Files DCLOSE Closes a directory that was opened by the DOPEN function and returns a valueDINFO Returns information about a directoryDNUM Returns the number of members in a directoryDOPEN Opens a directory and returns a directory identifier valueDOPTNAME Returns directory attribute informationDOPTNUM Returns the number of information items that are available for a directoryDREAD Returns the name of a directory memberDROPNOTE Deletes a note marker from a SAS data set or an external file and returns a value FAPPEND Appends the current record to the end of an external file and returns a value FCLOSE Closes an external file, directory, or directory member, and returns a valueFCOL Returns the current column position in the File Data Buffer (FDB)FDELETE Deletes an external file or an empty directoryFEXIST Verifies the existence of an external file associated with a fileref and returns a value FGET Copies data from the File Data Buffer (FDB) into a variable and returns a value FILEEXIST Verifies the existence of an external file by its physical name and returns a valueFILENAME Assigns or deassigns a fileref for an external file, directory, or output device and returns a valueFILEREF Verifies that a fileref has been assigned for the current SAS session and returns a value FINFO Returns the value of a file information itemFNOTE Identifies the last record that was read and returns a value that FPOINT can useFOPEN Opens an external file and returns a file identifier valueFOPTNAME Returns the name of an item of information about a fileFOPTNUM Returns the number of information items that are available for an external file FPOINT Positions the read pointer on the next record to be read and returns a valueFPOS Sets the position of the column pointer in the File Data Buffer (FDB) and returns a valueFPUT Moves data to the File Data Buffer (FDB) of an external file, starting at the FDB's current column position, and returns a valueFREAD Reads a record from an external file into the File Data Buffer (FDB) and returns a valueFREWIND Positions the file pointer to the start of the file and returns a valueFRLEN Returns the size of the last record read, or, if the file is opened for output, returns the current record sizeFSEP Sets the token delimiters for the FGET function and returns a valueFWRITE Writes a record to an external file and returns a valueMOPEN Opens a file by directory id and member name, and returns the file identifier or a 0PATHNAME Returns the physical name of a SAS data library or of an external file, or returns a blankSYSMSG Returns the text of error messages or warning messages from the last data set or external file function executionSYSRC Returns a system error numberExternal Routines CALL MODULE Calls the external routine without any return codeCALL MODULEI Calls the external routine without any return code (in IML environment only)MODULEC Calls an external routine and returns a character valueMODULEIC Calls an external routine and returns a character value (in IML environment only)MODULEIN Calls an external routine and returns a numeric value (in IML environment only)MODULEN Calls an external routine and returns a numeric valueFinancial COMPOUND Returns compound interest parametersCONVX Returns the convexity for an enumerated cashflowCONVXP Returns the convexity for a periodic cashflow stream, such as a bondDACCDB Returns the accumulated declining balance depreciationDACCDBSL Returns the accumulated declining balance with conversion to a straight-line depreciationDACCSL Returns the accumulated straight-line depreciationDACCSYD Returns the accumulated sum-of-years-digits depreciationDACCTAB Returns the accumulated depreciation from specified tablesDEPDB Returns the declining balance depreciationDEPDBSL Returns the declining balance with conversion to a straight-line depreciationDEPSL Returns the straight-line depreciationDEPSYD Returns the sum-of-years-digits depreciationDEPTAB Returns the depreciation from specified tablesDUR Returns the modified duration for an enumerated cashflowDURP Returns the modified duration for a periodic cashflow stream, such as a bondINTRR Returns the internal rate of return as a fractionIRR Returns the internal rate of return as a percentageMORT Returns amortization parametersNETPV Returns the net present value as a fractionNPV Returns the net present value with the rate expressed as a percentagePVP Returns the present value for a periodic cashflow stream, such as a bondSA VING Returns the future value of a periodic savingYIELDP Returns the yield-to-maturity for a periodic cashflow stream, such as a bond Hyperbolic COSH Returns the hyperbolic cosineSINH Returns the hyperbolic sineTANH Returns the hyperbolic tangentMacro CALL EXECUTE Resolves an argument and issues the resolved value for executionCALL SYMPUT Assigns DATA step information to a macro variableRESOLVE Returns the resolved value of an argument after it has been processed by the macrofacilitySYMGET Returns the value of a macro variable during DATA step executionMathematical ABS Returns the absolute valueAIRY Returns the value of the airy functionCNONCT Returns the noncentrality parameter from a chi-squared distributionCOMB Computes the number of combinations of n elements taken r at a time and returns a value CONSTANT Computes some machine and mathematical constants and returns a valueDAIRY Returns the derivative of the airy functionDEVIANCE Computes the deviance and returns a valueDIGAMMA Returns the value of the DIGAMMA functionERF Returns the value of the (normal) error functionERFC Returns the value of the complementary (normal) error functionEXP Returns the value of the exponential functionFACT Computes a factorial and returns a valueFNONCT Returns the value of the noncentrality parameter of an F distributionGAMMA Returns the value of the Gamma functionIBESSEL Returns the value of the modified bessel functionJBESSEL Returns the value of the bessel functionLGAMMA Returns the natural logarithm of the Gamma functionLOG Returns the natural (base e) logarithmLOG10 Returns the logarithm to the base 10LOG2 Returns the logarithm to the base 2MOD Returns the remainder valuePERM Computes the number of permutations of n items taken r at a time and returns a valueSIGN Returns the sign of a valueSQRT Returns the square root of a valueTNONCT Returns the value of the noncentrality parameter from the student's t distributionTRIGAMMA Returns the value of the TRIGAMMA functionProbability CDF Computes cumulative distribution functionsLOGPDF Computes the logarithm of a probability (mass) functionLOGSDF Computes the logarithm of a survival functionPDF Computes probability density (mass) functionsPOISSON Returns the probability from a Poisson distributionPROBBETA Returns the probability from a beta distributionPROBBNML Returns the probability from a binomial distributionPROBBNRM Computes a probability from the bivariate normal distribution and returns a value PROBCHI Returns the probability from a chi-squared distributionPROBF Returns the probability from an F distributionPROBGAM Returns the probability from a gamma distributionPROBHYPR Returns the probability from a hypergeometric distributionPROBMC Computes a probability or a quantile from various distributions for multiple comparisons of means, and returns a valuePROBNEGB Returns the probability from a negative binomial distributionPROBNORM Returns the probability from the standard normal distributionPROBT Returns the probability from a t distributionSDF Computes a survival functionQuantile BETAINV Returns a quantile from the beta distributionCINV Returns a quantile from the chi-squared distributionFINV Returns a quantile from the F distributionGAMINV Returns a quantile from the gamma distributionPROBIT Returns a quantile from the standard normal distributionTINV Returns a quantile from the t distributionRandom Number CALL RANBIN Returns a random variate from a binomial distribution CALL RANCAU Returns a random variate from a Cauchy distributionCALL RANEXP Returns a random variate from an exponential distributionCALL RANGAM Returns a random variate from a gamma distributionCALL RANNOR Returns a random variate from a normal distributionCALL RANPOI Returns a random variate from a Poisson distributionCALL RANTBL Returns a random variate from a tabled probability distributionCALL RANTRI Returns a random variate from a triangular distributionCALL RANUNI Returns a random variate from a uniform distributionNORMAL Returns a random variate from a normal distributionRANBIN Returns a random variate from a binomial distributionRANCAU Returns a random variate from a Cauchy distributionRANEXP Returns a random variate from an exponential distributionRANGAM Returns a random variate from a gamma distributionRANNOR Returns a random variate from a normal distributionRANPOI Returns a random variate from a Poisson distributionRANTBL Returns a random variate from a tabled probabilityRANTRI Random variate from a triangular distributionRANUNI Returns a random variate from a uniform distributionUNIFORM Random variate from a uniform distributionSAS File I/O ATTRC Returns the value of a character attribute for a SAS data setATTRN Returns the value of a numeric attribute for the specified SAS data setCEXIST Verifies the existence of a SAS catalog or SAS catalog entry and returns a valueCLOSE Closes a SAS data set and returns a valueCUROBS Returns the observation number of the current observationDROPNOTE Deletes a note marker from a SAS data set or an external file and returns a valueDSNAME Returns the SAS data set name that is associated with a data set identifierEXIST Verifies the existence of a SAS data library memberFETCH Reads the next nondeleted observation from a SAS data set into the Data Set Data Vector (DDV) and returns a valueFETCHOBS Reads a specified observation from a SAS data set into the Data Set Data Vector (DDV) and returns a valueGETV ARC Returns the value of a SAS data set character variableGETV ARN Returns the value of a SAS data set numeric variableIORCMSG Returns a formatted error message for _IORC_LIBNAME Assigns or deassigns a libref for a SAS data library and returns a valueLIBREF Verifies that a libref has been assigned and returns a valueNOTE Returns an observation ID for the current observation of a SAS data setOPEN Opens a SAS data set and returns a valuePATHNAME Returns the physical name of a SAS data library or of an external file, or returns a blankPOINT Locates an observation identified by the NOTE function and returns a valueREWIND Positions the data set pointer at the beginning of a SAS data set and returns a valueSYSMSG Returns the text of error messages or warning messages from the last data set or external file function executionSYSRC Returns a system error numberV ARFMT Returns the format assigned to a SAS data set variableV ARINFMT Returns the informat assigned to a SAS data set variableV ARLABEL Returns the label assigned to a SAS data set variableV ARLEN Returns the length of a SAS data set variableV ARNAME Returns the name of a SAS data set variableV ARNUM Returns the number of a variable's position in a SAS data setV ARTYPE Returns the data type of a SAS data set variableSpecial ADDR Returns the memory address of a variableCALL POKE Writes a value directly into memoryCALL SYSTEM Submits an operating environment command for executionDIF Returns differences between the argument and its nth lagGETOPTION Returns the value of a SAS system or graphics optionINPUT Returns the value produced when a SAS expression that uses a specified informat expression is readINPUTC Enables you to specify a character informat at run timeINPUTN Enables you to specify a numeric informat at run timeLAG Returns values from a queuePEEK Stores the contents of a memory address into a numeric variablePEEKC Stores the contents of a memory address into a character variablePOKE Writes a value directly into memoryPUT Returns a value using a specified formatPUTC Enables you to specify a character format at run timePUTN Enables you to specify a numeric format at run timeSYSGET Returns the value of the specified operating environment variableSYSPARM Returns the system parameter stringSYSPROD Determines if a product is licensedSYSTEM Issues an operating environment command during a SAS session State and ZIP Code FIPNAME Converts FIPS codes to uppercase state namesFIPNAMEL Converts FIPS codes to mixed case state namesFIPSTATE Converts FIPS codes to two-character postal codesSTFIPS Converts state postal codes to FIPS state codesSTNAME Converts state postal codes to uppercase state namesSTNAMEL Converts state postal codes to mixed case state namesZIPFIPS Converts ZIP codes to FIPS state codesZIPNAME Converts ZIP codes to uppercase state namesZIPNAMEL Converts ZIP codes to mixed case state namesZIPSTA TE Converts ZIP codes to state postal codesTrigonometric ARCOS Returns the arccosineARSIN Returns the arcsineATAN Returns the arctangentCOS Returns the cosineSIN Returns the sineTAN Returns the tangentTruncation CEIL Returns the smallest integer that is greater than or equal to the argumentFLOOR Returns the largest integer that is less than or equal to the argumentFUZZ Returns the nearest integer if the argument is within 1E-12INT Returns the integer valueROUND Rounds to the nearest round-off unitTRUNC Truncates a numeric value to a specified lengthVariable Control CALL LABEL Assigns a variable label to a specified character variableCALL SET Links SAS data set variables to DATA step or macro variables that have the same name and data typeCALL VNAME Assigns a variable name as the value of a specified variableVariable Information V ARRAY Returns a value that indicates whether the specified name is an arrayV ARRAYX Returns a value that indicates whether the value of the specified argument is an arrayVFORMAT Returns the format that is associated with the specified variableVFORMA TD Returns the format decimal value that is associated with the specified variableVFORMA TDX Returns the format decimal value that is associated with the value of the specified argumentVFORMA TN Returns the format name that is associated with the specified variableVFORMA TNX Returns the format name that is associated with the value of the specified argumentVFORMA TW Returns the format width that is associated with the specified variableVFORMA TWX Returns the format width that is associated with the value of the specified argumentVFORMA TX Returns the format that is associated with the value of the specified argumentVINARRAY Returns a value that indicates whether the specified variable is a member of an arrayVINARRAYX Returns a value that indicates whether the value of the specified argument is a member of an arrayVINFORMAT Returns the informat that is associated with the specified variable VINFORMATD Returns the informat decimal value that is associated with the specified variableVINFORMATDX Returns the informat decimal value that is associated with the value of the specified argumentVINFORMATN Returns the informat name that is associated with the specified variableVINFORMATNX Returns the informat name that is associated with the value of the specified argumentVINFORMATW Returns the informat width that is associated with the specified variableVINFORMATWX Returns the informat width that is associated with the value of the specified argumentVINFORMATX Returns the informat that is associated with the value of the specified argument VLABEL Returns the label that is associated with the specified variableVLABELX Returns the variable label for the value of a specified argumentVLENGTH Returns the compile-time (allocated) size of the specified variableVLENGTHX Returns the compile-time (allocated) size for the value of the specified argument VNAME Returns the name of the specified variableVNAMEX Validates the value of the specified argument as a variable nameVTYPE Returns the type (character or numeric) of the specified variableVTYPEX Returns the type (character or numeric) for the value of the specified argumentWeb Tools HTMLDECODE Decodes a string containing HTML numeric character references or HTML character entity references and returns the decoded stringHTMLENCODE Encodes characters using HTML character entity references and returns the encoded stringURLDECODE Returns a string that was decoded using the URL escape syntax URLENCODE Returns a string that was encoded using the URL escape syntax。
SAS产生随机数的方法随机数函数和CALL子程序
SAS产生随机数的方法随机数函数和CALL 子程序原文地址:SAS产生随机数的方法:随机数函数和CALL子程序作者:supersasmacro运用SAS进行MonteCarlo蒙特卡罗模拟(第四弹):SAS产生随机数的方法:随机数函数和CALL子程序本文未经作者同意严禁转载1随机数函数产生随机数序列随机数函数产生随机数序列的语法:var=name(seed,arg)，我们在前面的文章里使用的都是此类方法，变量var记录了由随机数种子为seed的随机函数name产生的一个随机数。
我们举两个例子来深入说明随机数函数产生随机数序列的原理。
程序1:DATATEMP1(DROP=I);DOI=1TO10;RUNI=RANUNI(123);SEED=RUNI*(2*31-1);OUTPUT;END;RUN;PROCPRINTDATA=TEMP1;RUN;程序2:DATATEMP2(DROP=I);DOI=1TO5;RUNI1=RANUNI(123);RUNI2=RANUNI(456);*这里我们虽然指定了另一个随机数种子的值，但其实并不起作用;OUTPUT;END;RUN;PROCPRINTDATA=TEMP2;RUN;程序1的结果:程序2的结果:将程序1的结果和程序2的结果进行比较，程序2中的RUNI1是程序1中的RUNI的第1,3,5,7,9条数据，而程序2中的RUNI2是程序1中的RUNI的第2,4,6,8,10条数据。
我们可以看到，虽然RUNI2=RANUNI(456);，但是这条语句中的随机数种子456并没有起作用。
RANUNI等随机函数的值产生后，SAS系统会隐含地将其转化为(0,1)的区间内，如果要看到这些随机数序列原始的值，可以通过程序1中的SEED=RUNI*(2*31-1);语句，结果如程序1输出结果中的seed所示。
2CALL子程序产生随机数序列CALL子程序产生随机数序列的语法如下:callname(seed,arg,var)，与随机数函数产生随机数序列的语法类似。
SAS常用函数介绍
SAS常用函数介绍
SAS常用函数介绍
SAS的功能非常强大和全面，这一点在函数上得到了极为充分的体现，SAS 8.2中共有25大类函数，功能涉及到文件管理、金触、数学计算、统计计算、字符串处理、货币转换等各个方面。
因篇幅限制，这里无法将数百种函数一一解释，只列举出饺为常用的十大类函数，关于这些函数的详细信息，以及术列出的货币转换函数、金融函数、地区以及邮政编码函数等相关信息请大家参见SAS帮助。
1.算术函数
2.字符函数
3.数学函数
4．概率与密度函数
5．分位数函数
6．随机函数
7．样本统计函数
8．三角函数
9．截断函数
10．日期时间函数。
第50章 SAS函数简介
第50章 SAS函数简介SAS函数是一个子程序，它由0个或几个自变量返回一个结果值。
每个SAS函数都有一个关键词名字。
为了引用函数，要写出它的名字，然后写出一个自变量或几个自变量，将其用放入括号中，而后这个函数对这些自变量进行某种运算。
SAS函数形式为 functionname (argument ,…)。
本章列举了常用SAS函数10类，以及关于这些函数的信息。
50.1 日期时间函数应用举例【例50-1】计算两个日期之间的天数。
具体SAS程序（程序名为SASTJFX50_1）：程序说明data _null_;sdate=’01jan2008’d;edate=’01jan2009’d;days=datdif(sdate,edate,’act/act’); days360=datdif(sdate,edate,’30/360’); put days= days360=;run; 调用datdif函数按每个月实际天数计算按每个月30天计算结果显示为：days=366 days360=360【例50-2】计算两个日期之间的年数。
具体SAS程序（程序名为SASTJFX50_2）： 程序说明data _null_;sdate=’01jan1999’d;edate=’01jan2009’d;year1= yrdif(sdate,edate,’30/360’); year2= yrdif(sdate,edate,’act/act’); year3=yrdif(sdate,edate,'act/360'); year4=yrdif(sdate,edate,'act/365'); put year1= year2= year3= year4= ; run; 调用yrdif函数按一年360天一个月30天计算按实际天数计算按一年360天每个月实际天数计算按一年365天每个月实际天数计算结果显示为：year1=10.083333333 year2=10.084931507 year3=10.233333333 year4=10.093150685【例50-3】计算当前时间的小时数和分钟数。
常用SAS函数简介
152.3 146.7
150.6 148.1
5. 截取函数简介
SAS截取函数应用举例
例5：请把STD=1.92354和CV=0.601106用 3个不同的截取函数(ceil、int、round)进行运 算，并分析这些函数的作用。
q1=2.9199855804 q2=13.894376071
利用HELP查找SAS系统提供的各种函数
Help下拉菜单
SAS 帮助 与文件
索引中键入： SAS Function
SAS函数种类
在SAS系统9.1版本,SAS函数分为二十五种类型, 分别是算术函数、截取函数、数学函数、概率 和密度函数、分位数函数、非中心函数、样本 统计函数、随机数函数、财政金融函数、日期 和时间函数、换算函数和特殊函数等。 在SAS系统6.12版本，SAS函数分为十七种类型 共有178个,SAS系统8.1版本，SAS函数分为二十 二种类型共有319个函数。
例1.2: 以日为单位计算当前日期的天数。
data; x=date (); y=today(); format x y yymmdd10.; put x= y=; run;
例1.3: 将时间值换算为以秒计数的 计算方法。
data; hms=hms (21,50,51); put hms=; x=21*60*60+50*60+51; put x=; mdy=mdy (08,18,2001); put mdy=; format mdy yymmdd10.; run;
练习1
计算以下数据表中日期值所在的年季月以及 处于某月的第几天、某周的第几天。
Date 1991-01-02 1991-01-03 1991-01-04 1991-01-05 1991-01-07 1991-01-08 1991-01-09 1991-01-10 1991-01-11 1991-01-12
2.SAS基础_表达式及函数简介
变量类型转换
字符自动转换为数值（无法转换为数值时，取缺失值）
data; X=1; /*X为数值变量 */ Y=’10’; /*Y为字符变量 */ X=Y; /*将字符变量赋值给一个数值变量时， 自动将字符变量Y转换为数值变量 */ run; data; X=1; /*X为数值变量 */ Y=’10’; /*Y为字符变量 */ Z=X+Y; /*算术表达式中有字符变量时，自动将字符变量Y 转换为数值变量 */ run;
Slide 10
创建变量
使用赋值语句
Data a; x=1; run;/*数据集a中有个变量x，值为1，为数值型*/
使用INPUT语句
Data b; Input x $ y z; Datalines; Aaa 23 45 ; Run; /*创建三个变量，x为字符型，y、z为数值型*/
Slide 11
Slide 25
二 SAS函数的参数
参数类型 • 变量名 • 常数 • 函数
• 表达式 函数及其参数类型举例。
函数举例 Max(X,Y); X=Repeat(‘---‘, 20); Least＝Min(sum(of x1-x10),y) ； 参数类型 变量名X，Y 常数 变量名和函数
SAS表达式及函数简介
Slide 1
SAS表达式
SAS表达式是由一系列算符和运算对象形成的一个指令集，它 被执行后产生一个目标值。 运算对象是SAS变量和SAS常数；算符是特殊的运算符、函数 和括号。
表达式分为简单表达式和复合表达式。 下列式子都是表达式。 X+1 3 LOG(X) LOG10(X) P/A*100 1-EXP(N/(N-1)) A=B=C STATE=‘CA’
2.SAS基础_表达式及函数简介
3.逻辑算符：
逻辑算符也称为布尔(Boolean)算符，在表达式里通常用来连接 一系列比较式。逻辑算符包括： &(AND) 与； | (OR) 或； ^ (NOT) 非 例如，age=25 AND sex=‘男’ 或 age=25 & sex=‘男’ age=25 OR sex=‘男’
4.运算次序：
SAS表达式及函数简介
Slide 1
SAS表达式
SAS表达式是由一系列算符和运算对象形成的一个指令集，它 被执行后产生一个目标值。 运算对象是SAS变量和SAS常数；算符是特殊的运算符、函数 和括号。
表达式分为简单表达式和复合表达式。 下列式子都是表达式。 X+1 3 LOG(X) LOG10(X) P/A*100 1-EXP(N/(N-1)) A=B=C STATE=‘CA’
AGE<100
在DATA步中，为了对变量进行变换，建立新变量、条件处理、计 算新值及指定新值，都需要使用SAS表达式。
Slide 2
一、SAS常数 SAS常数用来表示固定的值，它或者是一个数字， 或者是用引号引起来的字符串、或者是其它特殊记号。 SAS使用的常数有五种类型：数值常量、字符常量、 日期、时间和日期时间常量、十六进制数值常量、十六 进制字符常量。我们主要介绍数值常量、字符常量和日 期常量。 1.数值常量： 一个数值常量就是出现在SAS语句里的数字。数值 常量可以包括小数点，负号和E记号。 例如：1 -5 1.23 1.2E23 0.1 0.5E-10。 用E表示法时，如2E4，它表示 2 10 4 ，或者20000。
前缀算符用于变量、常数、函数或者用括号括起来的表达式 。例如： +Y -25 -COS(A) -(x1+x2-100)
SAS Function(函数)
SAS的函数对重要的函数加以介绍SAS提供了比一般程序设计语言多几倍的标准函数可以直接用在数据步的计算中，其中包括所有语言都有的数学函数、字符串函数，还包括特有的统计分布函数、分位数函数、随机数函数、日期时间函数、财政金融函数，等等。
这些函数的调用方法类似其它语言，比如求x1，x2，x3三个自变量的和可以用函数SUM(x1,x2,x3) 。
另外，SAS还提供了函数调用的另一种语法以便于把多个数据集变量作为函数自变量，其格式为“函数名(OF 变量名列表)”，其中变量名列表可以是任何合法的变量名列表，比如x1 ，x2，x3的和等价地可以用SUM(OF x1 x2 x3)或SUM(OF x1-x3)表示。
注意两种写法不能混在一起，比如SUM(OF x1,x2,x3)和SUM(x1-x3)都是错的。
本小节对重要的函数加以介绍，其它详见《SAS软件：Base SAS软件使用手册》（高惠璇等编译，中国统计出版社出版）。
一、数学函数ABS(x) 求x的绝对值。
MAX(x1,x2,…,xn) 求所有自变量中的最大一个。
MIN(x1,x2,…,xn) 求所有自变量中的最小一个。
MOD(x,y) 求x除以y的余数。
SQRT(x) 求x的平方根。
ROUND(x,eps) 求x按照eps指定的精度四舍五入后的结果，比如ROUND(5654.5654,0.01) 结果为5654.57，CEIL(x) 求大于等于x的最小整数。
当x为整数时就是x本身，否则为x右边最近的整数。
FLOOR(x) 求小于等于x的最大整数。
当x为整数时就是x本身，否则为x左边最近的整数。
INT(x) 求x扔掉小数部分后的结果。
FUZZ(x) 当x与其四舍五入整数值相差小于1E-12时取四舍五入。
LOG(x) 求x的自然对数。
LOG10(x) 求x的常用对数。
EXP(x) 指数函数。
SIN(x), COS(x), TAN(x) 求x的正弦、余弦、正切函数。
SAS常用函数介绍
一、数学函数ABS(x) 求x的绝对值。
MAX(x1,x2,…,xn) 求所有自变量中的最大一个。
MIN(x1,x2,…,xn) 求所有自变量中的最小一个。
MOD(x,y) 求x除以y的余数。
SQRT(x) 求x的平方根。
ROUND(x,eps) 求x按照eps指定的精度四舍五入后的结果，比如ROUND(5654.5654,0.01) 结果为5654.57，ROUND(5654.5654,10)结果为5650。
CEIL(x) 求大于等于x的最小整数。
当x为整数时就是x本身，否则为x右边最近的整数。
FLOOR(x) 求小于等于x的最大整数。
当x为整数时就是x本身，否则为x左边最近的整数。
INT(x) 求x扔掉小数部分后的结果。
FUZZ(x) 当x与其四舍五入整数值相差小于1E-12时取四舍五入。
LOG(x) 求x的自然对数。
LOG10(x) 求x的常用对数。
EXP(x) 指数函数。
SIN(x), COS(x), TAN(x) 求x的正弦、余弦、正切函数。
ARSIN(y) 计算函数y=sin(x)在区间的反函数，y取[-1,1]间值。
ARCOS(y) 计算函数y=cos(x)在的反函数，y取[-1,1]间值。
ATAN(y) 计算函数y=tan(x)在的反函数，y取间值。
SINH(x), COSH(x), TANH(x) 双曲正弦、余弦、正切ERF(x) 误差函数GAMMA(x) 完全函数此外还有符号函数SIGN，函数一阶导数函数DIGAMMA，二阶导数函数TRIGAMMA ，误差函数余函数ERFC，函数自然对数LGAMMA，ORDINAL函数，AIRY 函数，DAIRY函数，Bessel函数JBESSEL，修正的Bessel函数IBESSEL，等等。
二、数组函数数组函数计算数组的维数、上下界，有利于写出可移植的程序。
数组函数包括：DIM(x) 求数组x第一维的元素的个数（注意当下界为1时元素个数与上界相同，否则元素个数不一定与上界相同）。
SAS 函数简介
SAS 函数简介SAS 软件中共有13 类147 个S AS 函数, 现将其中最常用者的用法作一扼要介绍。
１．算术函数(1)求一个数的绝对值的函数A BS 和求一个正数的算术平根的函数S QRT如求-250 的绝对值的语句为x=ABS(-250)、求625 的算术平根的语句为y=SQRT(625)。
(2)求若干个数的最大值的函数MAX 和最小值的函数MIN DATA;a1=2; a2=-9; a3=15; a4=-23; x=MAX(of a1-a4);y=MIN(of a1-a4); FILE PRINT; PUT x y; RUN;程序运行结果为∶x=15、y=-23。
２．截取函数(1)取一个数的整数部分的函数I NT、按规定舍入精度取一个数的近似值的函数R OUND; (2)取≥自变量的最小整数的函数C EIL、取≤自变量的最大整数的函数F LOORDATA; 运行结果a=INT(223.456); b=ROUND(223.456, 100); a=223 b=200 c=ROUND(223.456,10); d=ROUND(223.456,1); c=220 d=223 e=ROUND(223.456,0.1); f=ROUND(223.456,0.01); e=223.5f=223.46g=CEIL(2.4); h=CEIL(-2.4); i=FLOOR(2.4); j=FLOOR(-2.4); g=3 h=-2FILE PRINT; PUT a b c d e f g h i j; RUN; i=2 j=-3３．数学函数(1)以ｅ为底的指数E XP、自然对数L OG(以e为底)、以2为底的对数L OG2 和常用对数函数LOG10(以10 为底)DATA; a=2; b=2.718281828; c=10; w=EXP(a); x=LOG(b);y=LOG2(a); z=LOG10(c); FILE PRINT; PUT w x y z; RUN;程序运行结果为∶w≈7.389056、x≈1、y=1、z=1。
sas中的常用函数（数学函数、字符…
一、数学函数1. ABS（X）：求X的绝对值；2. MAX(A1，A2，...，An)：求最大；3. MIN(A1，A2，...，An)：求最小；4. MOD(x,y)：x除以y的余数；5. sort(x)：x的平方根；6. round(x,eps)：x按照eps的精度四舍五入的结果，如round(5.777,2)=5.78；7. ceil(x)：大于等于x的最小整数；8. floor(x)：小于等于x的最大整数；9. int(x)：扔掉小数部分后的结果；10.fuzz(x)：当x与四舍五入整数值相差1E-12时取四舍五入；11.log(x)：取x的自然对数；12.log10(x)：取x的常用对数；13.exp(x)：指数函数；14.sin(x)，cos(x)，tan(x)：三角函数；15.arsin(x)：计算y=sin(x)的反函数，在[-1,1]区间；16.arcos(x)：计算y=cos(x)的反函数；17.atan(x)：计算y=tan(x)的反函数；18.sinh(x)/cosh(x)/tanh(x)：双曲正弦、余弦、正切；19.erf(x)：误差函数；20.gamma(x)：完全函数；二、字符函数1.trim(s)返回去掉字符串s的尾随空格的结果；2.upcase(s)把字符串s中的所有小写字母转换成大写字母的结果；3.lowcase(s)把字符串s中的所有大写字母转换成小写字母的结果；4.index(s,s1)查找s1在s中出现的位置。
找不到时返回0；5.rank(s)字符s的ASCII码；6.byte(n)第n个ASCII码值的对应字符；7.repeat(s,n)字符串s重复n次；8.substr(s,p,n)从字符串s中的第p个字符开始抽取n个字符长的子串；9.tranwrd(s,s1,s2)从字符串s中把所有字符串s1替换成字符串s2后的结果；其他字符函数还有collate、compress、indexc、left、length、reverse、right、scan、translate、verify、compbl、dequote、indexw、quote、soundex、trimn、indexw；三、日期和时间函数1.mdy(m,d,yr)生成yr年m月d日的sas日期值；2.year(date)由sas日期值data得到年；3.mouth(date)由sas日期值data得到月；4.day(date)由sas日期值data得到日；5.weekday(date)由sas日期值data得到星期几；6.qtr(date)由sas日期值data得到季度值；7.hms(h,m,s)由小时h，分钟m，秒s生成sas日期时间值；8.dhms(d,h,m,s)由sas日期值d，小时h，分钟m，秒s生成sas日期时间值；9.datepart(dt)求sas日期值dt的日期部分；10.index(interval,from,n)计算从from开始经过n个间隔后的sas日期。
spss随机数字函数说明
随机函数：RV.BERNOULLI(prob)。
数值。
返回具有指定概率参数prob 的Bernoulli 分布中的随机值。
RV.BETA(shape1, shape2)。
数值。
返回具有指定形状参数的Beta 分布中的随机值。
RV.BINOM(n, prob)。
数值。
返回具有指定试验次数和概率参数的二项式分布中的随机值。
RV.CAUCHY(loc, scale)。
数值。
返回具有指定位置和刻度参数的Cauchy 分布中的随机值。
RV.CHISQ(df)。
数值。
返回具有指定自由度df 的卡方分布中的随机值。
RV.EXP(scale)。
数值。
返回具有指定刻度参数的指数分布中的随机值。
RV.F(df1, df2)。
数值。
返回具有指定自由度df1 和df2 的F 分布中的随机值。
RV.GAMMA(shape, scale)。
数值。
返回具有指定形状和刻度参数的Gamma 分布中的随机值。
RV.GEOM(prob)。
数值。
以指定的概率参数从几何分布中返回一个随机值。
RV.HALFNRM(mean, stddev)。
数值。
以指定的均值和标准差从半正态分布中返回一个随机值。
RV.HYPER(total, sample, hits)。
数值。
以指定的参数从超几何分布中返回一个随机值。
RV.IGAUSS(loc, scale)。
数值。
以指定的位置和刻度参数从逆高斯分布中返回一个随机值。
PLACE(mean, scale)。
数值。
返回具有指定均值和刻度参数的Laplace 分布中的随机值。
RV.LNORMAL(a, b)。
数值。
返回具有指定参数的对数正态分布中的随机值。
RV.LOGISTIC(mean, scale)。
数值。
返回具有指定均值和刻度参数的logistic 分布中的随机值。
RV.NEGBIN(threshold, prob)。
数值。
返回具有指定阈值和概率参数的负二项式分布中的随机值。
RV.NORMAL(mean, stddev)。
SAS常用函数参考
SAS常⽤函数参考SAS常⽤函数参考这⼀节分类列出常⽤的函数，需要时可以参看帮助。
基本⼀、数据管理vector：向量 numeric：数值型向量 logical：逻辑型向量character；字符型向量 list：列表 data.frame：数据框c：连接为向量或列表 length：求长度 subset：求⼦集seq，from:to，sequence：等差序列rep：重复 NA：缺失值 NULL：空对象sort，order，unique，rev：排序unlist：展平列表attr，attributes：对象属性mode，typeof：对象存储模式与类型names：对象的名字属性⼆、字符串处理character：字符型向量 nchar：字符数 substr：取⼦串format，formatC：把对象⽤格式转换为字符串paste，strsplit：连接或拆分charmatch，pmatch：字符串匹配grep，sub，gsub：模式匹配与替换三、复数complex，Re，Im，Mod，Arg，Conj：复数函数四、因⼦factor：因⼦ codes：因⼦的编码 levels：因⼦的各⽔平的名字nlevels：因⼦的⽔平个数 cut：把数值型对象分区间转换为因⼦table：交叉频数表 split：按因⼦分组aggregate：计算各数据⼦集的概括统计量tapply：对“不规则”数组应⽤函数数学⼀、计算+, -, *, /, ^, %%, %/%：四则运算ceiling，floor，round，signif，trunc，zapsmall：舍⼊max，min，pmax，pmin：最⼤最⼩值range：最⼤值和最⼩值sum，prod：向量元素和，积cumsum，cumprod，cummax，cummin：累加、累乘sort：排序approx和approx fun：插值diff：差分sign：符号函数⼆、数学函数abs，sqrt：绝对值，平⽅根log, exp, log10, log2：对数与指数函数sin，cos，tan，asin，acos，atan，atan2：三⾓函数sinh，cosh，tanh，asinh，acosh，atanh：双曲函数beta，lbeta，gamma，lgamma，digamma，trigamma，tetragamma，pentagamma，choose ，lchoose：与贝塔函数、伽玛函数、组合数有关的特殊函数fft，mvfft，convolve：富利叶变换及卷积polyroot：多项式求根poly：正交多项式spline，splinefun：样条差值besselI，besselK，besselJ，besselY，gammaCody：Bessel函数deriv：简单表达式的符号微分或算法微分三、数组array：建⽴数组 matrix：⽣成矩阵data.matrix：把数据框转换为数值型矩阵lower.tri：矩阵的下三⾓部分 mat.or.vec：⽣成矩阵或向量t：矩阵转置 cbind：把列合并为矩阵 rbind：把⾏合并为矩阵diag：矩阵对⾓元素向量或⽣成对⾓矩阵aperm：数组转置 nrow, ncol：计算数组的⾏数和列数dim：对象的维向量 dimnames：对象的维名row/colnames：⾏名或列名 %*%：矩阵乘法crossprod：矩阵交叉乘积（内积） outer：数组外积kronecker：数组的Kronecker积 apply：对数组的某些维应⽤函数tapply：对“不规则”数组应⽤函数 sweep：计算数组的概括统计量aggregate：计算数据⼦集的概括统计量 scale：矩阵标准化matplot：对矩阵各列绘图 cor：相关阵或协差阵Contrast：对照矩阵 row：矩阵的⾏下标集col：求列下标集四、线性代数solve：解线性⽅程组或求逆 eigen：矩阵的特征值分解svd：矩阵的奇异值分解 backsolve：解上三⾓或下三⾓⽅程组chol：Choleski分解 qr：矩阵的QR分解chol2inv：由Choleski分解求逆五、逻辑运算<，>，<=，>=，==，!=：⽐较运算符!，&，&&，|，||，xor()：逻辑运算符logical：⽣成逻辑向量 all，any：逻辑向量都为真或存在真ifelse()：⼆者择⼀ match，%in%：查找unique：找出互不相同的元素 which：找到真值下标集合duplicated：找到重复元素六、优化及求根optimize，uniroot，polyroot：⼀维优化与求根程序设计⼀、控制结构if，else，ifelse，switch：分⽀for，while，repeat，break，next：循环apply，lapply，sapply，tapply，sweep：替代循环的函数。
SAS常用的随机数函数简介文档
运用SAS进行Monte Carlo蒙特卡罗模拟（第五弹）：SAS常用的随机数函数简介前一篇文章我们介绍了两种产生随机数序列的方法，即随机数函数产生随机数序列，其语法为：var = name(seed,<arg>)和CALL子程序产生随机数序列，其语法为：callname(seed,<arg>,var)。
本节我们将介绍SAS常用的随机数函数（其概率函数我们这里就不作详细介绍，感兴趣的话请查阅相关文献；SAS随机数函数中的seed均为随机数种子）：SAS随机数函数分布情况参数说明RANBIN(seed,n,p) 二项分布n：独立实验的次数，p：成功的概率RANCAU(seed) 柯西分布RANEXP(seed) 指数分布RANGAM(seed,a) 伽玛分布 a：a>0，形状参数RANNOR(seed) 正态分布NORMAL(seed) 正态分布RANPOI(seed,m) 泊松分布m：m>0，均值RANTBL(seed,p1,p2,...p n) 离散分布p(i)：p(i) >0，且Σp(i)=1，概率RANTRI(seed,h) 三角分布h：0<=h<=1，斜边RANUNI(seed) 均匀分布UNIFORM(seed) 均匀分布这里要注意：Functions RANUNI and UNIFORM are identical. Function UNIFORM cannot be utilized as a CALL routine.文章中还举例说明了用上述基础的SAS随机数函数通过变换，可以产生很多有趣的分布，本人对此没有研究，请大家查看相关文献。
所有的SAS随机数函数都是通过RANUNI随机数函数变换得到的，例如我们通过就可以得到一个正态分布，通过e=-ln(u3)就可以得到指数分布。
通过下面的例子我们可以证明刚才的结论：程序一：DATA TEMP5(DROP=I);DO I=1 TO 12;RUNI=RANUNI(123);OUTPUT;END;RUN;PROC PRINT DATA=TEMP5;RUN;程序二：DATA TEMP6(DROP=I);DO I=1 TO 3;RUNI=RANUNI(123);RNOR=RANNOR(456);REXP=RANEXP(789);OUTPUT;END;RUN;PROC PRINT DATA=TEMP6;RUN;程序一结果：OBS RUNI1 0.750402 0.320913 0.178394 0.906035 0.357126 0.221117 0.786448 0.398089 0.1246710 0.1876911 0.7761812 0.43607程序二结果：OBS RUNI RNOR REXP1 0.75040 0.65572 0.098682 0.35712 0.39428 0.921103 0.12467 0.29958 0.82994从上面的结果我们可以看到，程序二的结果中，变量RUNI的三个值分别对应程序一结果里的程序二结果中的变量RUNI的第一、五、九个值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
运用SAS进行Monte Carlo蒙特卡罗模拟（第五弹）：
SAS常用的随机数函数简介
前一篇文章我们介绍了两种产生随机数序列的方法，即随机数函数产生随机数序列，其语法为：var = name(seed,<arg>)和CALL子程序产生随机数序列，其语法为：call
name(seed,<arg>,var)。
本节我们将介绍SAS常用的随机数函数（其概率函数我们这里就不作详细介绍，感兴趣的话请查阅相关文献；SAS随机数函数中的seed均为随机数种子）：
SAS随机数函数分布情况参数说明
RANBIN(seed,n,p) 二项分布n：独立实验的次数，p：成功的概率
RANCAU(seed) 柯西分布
RANEXP(seed) 指数分布
RANGAM(seed,a) 伽玛分布 a：a>0，形状参数
RANNOR(seed) 正态分布
NORMAL(seed) 正态分布
RANPOI(seed,m) 泊松分布m：m>0，均值
RANTBL(seed,p1,p2,...p n) 离散分布p(i)：p(i) >0，且Σp(i)=1，概率
RANTRI(seed,h) 三角分布h：0<=h<=1，斜边
RANUNI(seed) 均匀分布
UNIFORM(seed) 均匀分布
这里要注意：Functions RANUNI and UNIFORM are identical. Function UNIFORM cannot be utilized as a CALL routine.
文章中还举例说明了用上述基础的SAS随机数函数通过变换，可以产生很多有趣的分布，本人对此没有研究，请大家查看相关文献。
所有的SAS随机数函数都是通过RANUNI随机
数函数变换得到的，例如我们通过就可以得到一个正态分布，通过e=-ln(u3)就可以得到指数分布。
通过下面的例子我们可以证明刚才的结论：程序一：
DATA TEMP5(DROP=I);
DO I=1 TO 12;
RUNI=RANUNI(123);
OUTPUT;
END;
RUN;
PROC PRINT DATA=TEMP5;
RUN;
程序二：
DATA TEMP6(DROP=I);
DO I=1 TO 3;
RUNI=RANUNI(123);
RNOR=RANNOR(456);
REXP=RANEXP(789);
OUTPUT;
END;
RUN;
PROC PRINT DATA=TEMP6;
RUN;
程序一结果：
OBS RUNI
1 0.75040
2 0.32091
3 0.17839
4 0.90603
5 0.35712
6 0.22111
7 0.78644
8 0.39808
9 0.12467
10 0.18769
11 0.77618
12 0.43607
程序二结果：
OBS RUNI RNOR REXP
1 0.75040 0.6557
2 0.09868
2 0.35712 0.39428 0.92110
3 0.12467 0.29958 0.82994
从上面的结果我们可以看到，程序二的结果中，变量RUNI的三个值分别对应程序一结果里的程序二结果中的变量RUNI的第一、五、九个值。
为什么中间只隔了两个函数，结果却隔了三个数值呢，原因如上面公式里提到的，产生一个正态分布需要两个均匀分布函数，而产生一个指数分布只需要一个均匀分布函数。
为了进一步说明刚才提到的正态分布和指数分布是由均匀分布函数转换而来的，我们可以看下面这个例子：
程序三：
DATA TEMP7(DROP=I);
DO I=1 TO 3;
RUNI=RANUNI(123);
U1=RANUNI(123);***三个均匀分布随机数;
U2=RANUNI(123);
U3=RANUNI(123);
RNOR=SQRT(-2*LOG(U1))*COS(4*ARSIN(1)*U2); ***正态分布公式;
REXP=-LOG(U3); ***指数分布公式;
OUTPUT;
END;
RUN;
PROC PRINT DATA=TEMP7;
RUN;
程序三结果：
OBS RUNI U1 U2 U3 RNOR REXP
1 0.75040 0.32091 0.17839 0.90603 0.6557
2 0.09868
2 0.35712 0.22111 0.78644 0.39808 0.39428 0.92110
3 0.12467 0.18769 0.77618 0.43607 0.29958 0.82994
这个结果要结合程序二的结果看，我们就会发现，两个程序的RNOR和REXP值是一样的，这就说明了这两个分布确实是通过上述公式由均匀分布转换而来的。
参考资料
Xitao Fan, etc..Monte Carlo Studies: A Guide for Quantitative Researchers. SAS Institute
Inc.,2002
本文章只用于学习，请不要用于商业目的，否则后果自负。