抽样、模拟、实验设计
qq1323362960
这个作者很懒,什么都没留下…
展开
-
Getting Your Random Sample in Proc SQL
<br />代码摘录自:《Getting Your Random Sample in Proc SQL》<br />libname yugao "C:/BOOKS" ;<br />data yugao.sample;<br /> format rec_no 7.0 age 4.0 sex $6. ;<br /> label rec_no = "Record Number / ID"<br /> age = "Age in Years"<br /> sex = "Gen原创 2011-05-13 20:55:00 · 770 阅读 · 0 评论 -
哑变量矩阵生成dummy matrix with sas
FROM:http://support.sas.com/kb/23/217.htmlFor a specified model, there are several procedures that allow you to save the design matrix to a data set:For models using GLM parameterization (also转载 2012-11-28 00:02:56 · 4423 阅读 · 0 评论 -
Gibbs sampling & R
维基百科的定义:http://en.wikipedia.org/wiki/Gibbs_sampling我们假设一个随机的二元变量(x,y),然后计算其中一个或全部的边缘分布p(x),p(y)。这种抽样思想是考虑条件分布p(xly),p(ylx),比通过联合密度p(xy)来计算求解简单,例如首先,我们要进行初始化y0,然后通过条件概率P(xIy=y0)计算得出x0,然后再基于x0的条原创 2011-12-13 21:28:19 · 4806 阅读 · 1 评论 -
Generating Multivariate Normal Data by Using PROC IML & R
Generating Multivariate Normal Data by Using PROC IML : http://analytics.ncsu.edu/sesug/2006/CC15_06.PDF1.Generate the bivariate normal data/* Generate the bivariate normal data */data one;m原创 2011-12-16 21:51:47 · 1483 阅读 · 0 评论 -
卡方随机数
1.利用RANGAM,这种方法的原理是利用概率的逆运算。我们知道卡方分布是特殊的伽马分布,伽马分布的形状参数alpha=n/2,尺度参数l=0.5时,它就是自由度为n的卡方分布,故可以通过RANGAR生成卡方随机数例如:我们生成自由度为3的卡方分布,100个随机数data a; do i=1 to 100; Ch = RanGam(121212,3/2)/(0原创 2011-11-03 22:23:53 · 4393 阅读 · 0 评论 -
简单随机抽样sas实现
简单随机抽样,就是以相同概率从N个总体单位中抽取n个单位。根据有无放回,它又分为无放回的简单随机抽样和有放回的简单随机抽样。支持以上观点的书籍很多,如《Elements of sampling theory and methods》、《抽样调查的理论及其应用方法》、《抽样技术》等原创 2011-03-09 20:55:00 · 4362 阅读 · 0 评论 -
sas中的随机数发生器
http://blog.renren.com/blog/220869953/729109705?from=fanyeNew Random Number Generator Functions Available in S原创 2011-10-17 20:09:11 · 5528 阅读 · 0 评论 -
Monte Carlo Simulation
A.STEP,(FROM《SAS For Monte Carlo Studies》)1. Designing the system (What are the parameters of the system? What are the relationships among原创 2011-10-16 22:34:47 · 1065 阅读 · 0 评论 -
常用的试验设计方法和SAS程序
摘自:http://www.mysas.net/forum/viewtopic.php?f=4&t=16171.完全随机设计 例1:将10例实验对象随机地分配到A、B两组,先将实验对象按某项特征编号,然后从随机分配(数字)表中任意指定一行,舍掉10―19的数字,规定随转载 2011-03-20 14:52:00 · 4098 阅读 · 0 评论 -
BOOTSTRAP RESAMPLING DISTRIBUTIONS OF THE MEAN WHEN THE CENTRAL LOCATION PARAMETER IS KNOWN
——— Peter M. Saama作者分别利用Mathematica、sas/iml、splus及stat完成该bootstrap实验,这里摘录仅sas/iml代码。值得强调的是,它属于percentile-bootstrap。1.主要算法步骤a) Define p转载 2011-10-07 17:18:47 · 662 阅读 · 0 评论 -
缺失值简单插补方法
1.利用均值/最小值/最大值等进行插补PROC SQL noprint; create table sample as select * ,CASE _value_ when . then MEAN(_val原创 2011-09-08 10:39:55 · 3345 阅读 · 0 评论 -
抽样之proc surveyselect
现摘录SAS文档中的代码:/* Generate Sampling Frame -------------------------------------*/data Customer1; input CustomerID State$ Type$ Usage; format CustomerID SSN11.; datalines;416874322 AL New 839288139763 GA Old 224339008654 GA Old 2451118980542 GA New 3494原创 2011-04-05 20:33:00 · 7899 阅读 · 0 评论 -
bootstraping
之前一位同学问及bootstrap,由此我查阅了几篇文献,初步知晓个皮毛:它是一种非参检验方法,利用重复抽样理论,来减少偏差、控制方差、得到有效置信区间等统计方法。国内bootstrap研究比较少,这里摘录了国外研究成果。1.what is bootstrappingThe method of bootstrapping was first introduced by Efron as a method to derive the estimate of standard error of an arbit原创 2011-04-02 21:37:00 · 3647 阅读 · 0 评论 -
分层抽样Selecting a Stratified Sample with PROC SURVEYSELECT
声明:该文作者是Diana Suhr, University of Northern Colorado,摘自SAS Global Forum 2009. AbstractStratified random sampling is simple and efficient using PROC FREQ and PROC SURVEYSELECT. A routine wasde转载 2011-03-12 21:53:00 · 2248 阅读 · 0 评论