加载中…
个人资料
独闲居士
独闲居士
  • 博客等级:
  • 博客积分:0
  • 博客访问:2,623
  • 关注人气:126
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
博文
(2020-07-20 20:26)
标签:

教育

搞定省级、城市、制造业细分行业的 出口技术复杂度计算,2000-2019年。

需要联系qq:773721498
(2013-02-23 10:01)
Stata

2011-10-02 08:58:28 分类: 默认分类  标签: |字号 订阅

Using STATA 11:

Define the x and y matrices.

: b=invsym(x’x)

(2012-12-17 22:01)
标签:

it

描述统计
describe
describe命令可以描述数据文件的整体,包括观测总数,变量总数,生成日期,每个变量的存储类型(storage type),标签(label)等。

list [varlist] [if exp] [in range]

summarize [varlist] [weight] [if exp] [in range] [,detail]
summarize可以提供varlist指定变量(可以不止一个)的如下统计量:Percentiles(分位数),四大最大的数和四个最小的数,Variance(方差),Std. Dev.(标准差),Skewness(偏度),Kurtosis(斜度)

tabstat
tabstat varlist [weight] [if exp] [in range] [, stats(statname [...]) ]
tabstat提供[, stats(statname [...]) ]指定的统计量,可供选择的有mean(均值),count(非缺失观测值个数),sum(总和),max(最大值),min(最小值),range(最大值-最小值),sd(标准差),var(方差),cv(变易系数=标准差/均值),skewness(偏度),kurtosis(斜度),median(中位数),p1(1%分位数,类似地有p5, p10, p25,

标签:

杂谈

 一、带标注的散点图 sysuse lifeexp, clear keep if region==2 | region==3 replace gnppc = gnppc / 1000 label var gnppc 'GNP per capita (thousands of dollars)' gen lgnp = log(gnp) qui reg lexp lgnp predict hat  label var hat 'Linear prediction' replace country = 'Trinidad' if country=='Trinidad and Tobago' replace country = 'Para' if country == 'Paraguay' gen pos = 3  replace pos = 9 if lexp > hat  replace pos = 3 if country == 'Colombia' replace pos = 3 if country == 'Para' replace pos = 3 if country == 'Trinidad' replace pos = 9 if country == 'United States'  #delimit ;  twoway (scatter lexp gnppc, mlabel(country) mlabv(pos)) (line hat gnppc, sort) , xsca(log) xlabel(.5 5 10 15 20 25 30, grid) legend(off)    title('Life expectancy vs. GNP per capita')    subtitle('
标签:

教育

表中第③栏是我国1978-2002年全社会客运量的资料,据期绘制散点图,见下图,可以看出,各年的客运量资料基本呈线性趋势,但在几个不同的时期直线有不同的斜率,因此考虑用变参数线性趋势模型进行预测。具体步骤如下:

  表 我国1978-2002年全社会客运量及预测值 单位:万人

年份 时间t 全社会客运量y 各期的一次指数平滑值http://wiki.mbalib.com/w/images/math/5/c/2/5c2f495ad0869994a3e7dcaf10908650.png 各期的二次指数平滑值http://wiki.mbalib.com/w/images/math/3/1/a/31a3aac99dd5668600b6a1f9435a661b.png at
标签:

it

MATLAB画时间序列的自相关函数和偏自相关函数图

 

autocorr(Series)       %画出自相关图,图中上下两条横线分别表示自相关系数的上下界,超出边界的部分表示存在相关关系。
[a,b] = autocorr(Series)   %a 为各阶的相关系数,b 为滞后阶数
parcorr(Series)  %画出偏自相关图
[c,d] = parcorr(Series)    %c 为各阶的偏自相关系数,d 为滞后阶数

子相关系数     偏相关系数     模型定阶
拖尾                 p阶截尾           AR(p)模型
q阶截尾           拖尾                 MA(q)模型
拖尾                 拖尾                 ARMA(p,q)模型
标签:

教育

控制变量在工具书中的解释:除自变量之外,一切能使因变量发生变化的变量。这类变量是应该加以控制的,如果不加控制,它也会造成因变量的变化,即自变量和一些未加控制的因素共同造成了因变量的变化,这叫自变量的混淆。因此,只有将自变量以外一切能引起因变量变化的变量控制好,才能弄清实验中的因果关系。

  控制变量——在实验中,不仅仅只有自变量才是和因变量有关的,于自变量之外往往存在额外相关变量,此类变量简称额外变量,因其必须被想办法控制,在实验中保持恒定不变,又称其为控制变量。

 

 

某一个变量模型中随机解释变量高度相关,但却不与

标签:

教育

注:JB统计量对应的p大于0.05,则表明非正态,这点跟sktest和swilk检验刚好相反;
dta为数据文件;
gph为图文件;
do为程序文件;
注意stata要区别大小写;
不得用作用户变量名:
_all _n _N _skip _b _coef _cons _pi _pred _rc _weight double
float long int in if using with
命令:
读入数据一种方式
input  y
        
       5.5 
       6.2 
       7.7 
       8.5 
end

su/summarise/sum x 或 su/summarise/sum x,d
对分组的描述:
sort group
by group:su x
%%%%%
tabstat economy,stats(max)  %返回变量economy的最大值
%%stats括号里可以是:mean,count(非缺失观测值个数),sum(总和),max,min,range,
%%      s

标签:

it

原文地址:http://yinlianqian.blog.163.com/blog/static/6488756720082268510559/

作者:思想的成长轨迹

 

写在这里保存下来,我自己以后要看看:)

 Bugs 网站:www.mrc-bsu.cam.ac.uk/bugs

1. 点选 图示,进入WinBugs,开启 (1)新的程式编辑视窗or (2)已存在的程式功能表列:File-New 功能表列:File-Open

2. 撰写程式:程式包含三部份.

(1) Model:建构贝氏统计模式,设定各参数的prior distribution及各参数间的关系等.

(2) Data:以list指令起始,列出各参数的样本观察值及样本个数(N= ).

(3) Initial value:同样运用list指令,列出各参数的起始值.

ps:上述三大

标签:

it

Matlab工具箱命令汇总
学习   2008-07-01 11:37   阅读62   评论1   字号: 大大  中中  小小 Ⅰ.1  统计工具箱函数
表Ⅰ-1  概率密度函数
函数名
对应分布的概率密度函数
betapdf
贝塔分布的概率密度函数
binopdf
二项分布的概率密度函数
chi2pdf
卡方分布的概率密度函数
exppdf
指数分布的概率密度函数
fpdf
f分布的概率密度函数
gampdf
伽玛分布的概率密度函数
geopdf
几何分布的概率密度函数
hygepdf
超几何分布的概率密度函数
normpdf
正态(高斯)分布的概率密度函数
lognpdf
对数正态分布的概率密度函数
nbinpdf
负二项分布的概率密度函数
ncfpdf
非中心f分布的概率密度函数
nctpdf
非中心t分布的概率密度函数
ncx2pdf
非中心卡方分布的概率密度函数
poisspdf
泊松分布的概率密度函数
raylpdf
雷利分布的概率密度函数
tpdf
学生氏
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有