2011年12月_qq1323362960

12月 11月 10月 09月 08月 07月 06月 05月 04月 03月 02月

转载 R学习网址

转自：清风驿栈http://www.stat.auckland.ac.nz/~paul/RGraphics/rgraphics.htmlhttp://had.co.nz/ggplot2/http://addictedtor.free.fr/graphiques/http://zoonek2.free.fr/UNIX/48_R/all.html 《Statistics with R》

2011-12-30 19:45:34 3459

转载 R 基础

#======================R-demonstration script=============================#--------------------------start and quit---------------------------------Rquit()save.image() #save data, although I w

2011-12-30 19:25:02 637

原创判别分析、logistic在信用风险评估中的应用

利用判别分析、logistic模型来分析信用风险的文章很多，这类文章应该传入国内最早的方法之一了。但是看期刊文章可以发现，作者都喜欢直接给出结果，至于如何得到结果的过程通常并不是非常的明晰。《Application of Proc Discrim and Proc Logistic in Credit Risk Modeling》一文给出SAS软件的具体操作：1.判别分析之前转载过一篇判别

2011-12-24 21:07:19 2099

原创多元统计学习资源

传统的多元统计主要涉及多元分布与多元回归、均值向量、协方差矩阵、聚类判别、主成分与因子分析、对应分析、典型相关等等。理论学习一直都会是长期的过程，至于如何实现上述多元统计方法，推荐如下sas:http://psych.colorado.edu/~carey/Courses/PSYC7291/ExampleCode.htm,这个链接主要是用procs完成http://www.psyc

2011-12-22 22:22:49 882

转载信用风险理论、模型及应用研究（综述）

作者：李兴法国内外研究历史和现状一、国外研究（一）均值一方差理论综述。Harry Markowitz(1952)引入了均值一方差框架用以科学计量风险与收益问题，从而为风险的定量研究建立了数学基础。Fiseher Black和Myron Seholes(1973)推导出股票的欧式期权的价值。Robert Merton(1974)采用Black-Seholes模型解出期权的价值。Ge

2011-12-21 21:44:45 15149

原创 PROC SQL & DATA Step do the same job

sas中，很多情况下 data setp能够搞定的，proc sql也可以搞定。论坛里面经常碰到一题多解，例如汇总、重复数据、数据合并等问题。很多网友也曾作过一些汇总，且巧浏览网页时发现NESUG 2011的一篇文章总结的相当齐全，特此摘录文中例子。感兴趣的网友可以自己google，原文题目《PROC SQL for DATA Step Die-hards》，作者：Christianna S.

2011-12-19 17:22:50 2459

原创 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

转载牛顿迭代例子Newton-Raphson Method

Aim: Find oˆ such thatProblem: Analytic solution of likelihood equations not always available. Example: Censored exponentially distributed observations Suppose that and that the censored time

2011-12-15 20:29:50 3742

原创 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 4807 1

原创初识EM algorithm

很多书阐述了EM algorithm的原理，直接点可以在维基百科上找到它的基本情况：http://en.wikipedia.org/wiki/Expectation%E2%80%93maximization_algorithm。EM算法分两步走即E步和M步，和很多书本和文章中一样，维基中只给出了核心结论，至于其中的数理逻辑关系并没有给出，需要探究EM算法的推导过程可以参考JerryLead

2011-12-10 21:15:36 766

原创 iml数据的处理和输入输出

之前在“IML 编程的基本函数概况”博文中有过简单的概述，这里的例子更详尽。http://www.stat.ucl.ac.be/SMCS/serveur/SAS/samples/iml/samp_iml.htmlods html body='datademo.html' ; run ; /*--create filename for class data--*/%macro fil

2011-12-09 21:05:14 1197

原创一个Statistical Computing主页

#Plotting in preparation for the ImpSamp function. Suggests ideal g* would be an unusual distributionx=seq(-3,3,by=.1)gx=sqrt(abs(x/(1-x)))*dt(x,5)plot(x,gx,type="l")lines(x,dt(x,5),lty=2,col

2011-12-06 22:07:01 739

原创 SAS读取复杂的csv文件

最初接触读取excel文件时候，学习过几种简单的读取方式，见: http://blog.csdn.net/yugao1986/article/details/6359406但最近接触了一些稍微复杂的数据，再简单无变化的调用前面讲到的方法，发现sas总是报错。自己摸索了下，总结两条：首先，数据形式大致形如下面两种情况： 1.方法来自《Reading Delimited Text

2011-12-05 19:22:56 9528

原创利用filename and pipe读取文件

最近在网上看到一些关于批量读取文件的帖子，想到今后定会遇到，特作记录。一.Filename1.filename function语法a.FILENAME(fileref )帮助文档里面有着经典的例子：**Example 1: Assigning a Fileref to an External File ;%let filrf=myfile;%let rc=%sy

2011-12-03 14:16:26 12479 1

转载 CALL SYMPUT Routine

CALL SYMPUT(macro-variable,value);macro-variable can be one of the following items:a character string that is a SAS name, enclosed in quotation marks. For example, to assign the

2011-12-01 20:06:19 2340