mata: st_matrix("bst",st_matrix("e(b)") :/ sqrt(diagonal(st_matrix("e(V)"))')) /*该行代码具体解释如下: st_matrix()的使用可利用 help mata st_matrix 进行查看 bst为逗号后矩阵命名。st_matrix(name, X)将Stata中的名为name的矩阵新建/替换为X中的内容 st_matrix("e(b)")调用系数矩阵...
比如现在Stata中有一个矩阵A和一个数值型标量a,那么B=st_matrix("A")和b=st_scalar("a")就会在mata中生成对应的矩阵和数值型标量。如果这个标量是字符型的,则需要命令st_strscalar(),如果仍然使用st_numscalar()的话,该标量进入mata之后,就变成一个0行0阶的矩阵。 最后再看看宏如何进入mata: st_local()...
例如,假设一个二次型趋势: mata st_matrix("deltaquad",0.024 *((-4::3):-(-1)):^2)pretrends,numpre(3)b(beta)v(sigma)deltatrue(deltaquad)coefplot image.png 详细的统计量展示可以用: matlistr(results)returnlist
mata:st_matrix("sum", colsum(st_matrix("freq"))) matrix rowjoin freq=freq sum mat per=(freq/sum[1,1],freq/r(N)) mat coljoin multable=freq per matrix rownames multable=`r(list) ' "总计" matrix colnames multable="频率" "响应百分比" "个案百分比" 第四部分:结果导出excel putexcel a...
x=st_data(.,"x*") c=correlation(x) n=rows(c) b=strofreal(lowertriangle(c)+uppertria gle(st_matrix("r(Rho)")),"%9.3f") p=st_matrix("r(P)") for (i=2; i<=n; i++) {for (j=1; j<=i-1; j++) {p[i,j]=2*ttail ...
st_matrix("H",(re, im)) end local j=0 gen yhat=0 forvalues i=-63(1)64{ local j=`j'+1 quietly replace yhat=yhat+H[`j',1]*cos(`i'*x) /* */ +H[`j',2]*sin(`i'*x) } replace yhat=(1/128)*yhat #delimit;
index=st_addvar("long",("X1","X2","X3")) st_addobs(rows(Y)) st_store((1,rows(Y)),index,Y) 其中st_addobs(rows(Y))的作用是将Stata数据中的obs设定为这个矩阵的行数。缺少这一步,Stata可能会报告出错。 输出矩阵就容易多了,st_matrix("r(A)",A)将矩阵输出到r(A)中,st_...
x=st_data(.,“x*”) c=correlation(x) n=rows© b=strofreal(lowertriangle©+uppertriangle(st_matrix(“r(Rho)”)),“%9.3f”) p=st_matrix(“r§”) for (i=2; i<=n; i++) { for (j=1; j<=i-1; j++) { p[i,j]=2ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,j]...
spearmanx*n matax=st_data(.,"x")c=correlation(x)n=rows(c)b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")),"%9.3f")p=st_matrix("r(P)")for (i=2; i<=n; i++) {for (j=1; j<=i-1; j++) {p[i,j]=2ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,...
x=st_data(.,"x*") c=correlation(x) n=rows(c) b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")),"%9.3f") p=st_matrix("r(P)")for(i=2; i<=n; i++) {for(j=1; j<=i-1; j++) { p[i,j]=2*ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,j]^2)/...