* t å¼çè®¡ç® t = ç³»æ°/æ å误
use B1_consume, clear
regress consume income
dis %4.2f 0.6848014/0.2487525 /*income ç t å¼*/
dis %4.2f 51.89511/80.84397 /*
常æ°é¡¹ ç t å¼*/
* ç©éµè§£æ
mat b0 = diag(b)
mat list b0
mat inv_se_b = inv(se_b)
mat list inv_se_b
mat t = hadamard(b0, inv_se_b)
mat list t
mat t = vecdiag(t)
reg consume income
* p å¼
* H0 : bj = 0 å³ï¼ç³»æ°ä¼°è®¡å¼æ¯å¦æ¾èä¸çäºé¶
* ç±äº tå¼ æä» t åå¸ï¼æ以æ们å¾å®¹æ计ç®å
¶ p å¼
help density functions
mat list t
local p_income = ttail(11-2, 2.75)*2 /*åå°¾*/
local p_cons = ttail(11-2, 0.64)*2
mat pvalue = (`p_income' \ `p_cons')
mat list pvalue