联立方程组中变量联合显著性检验

根据你的提示，是5个方程，每个方程包括截距项是13个参数，因此，估计完有5*13个参数
因此，你按照前面的例子，构造方程进行估计

约束条件写成下面的形式就可以了

R1 <- matrix( 0, nrow = 5, ncol = 5*13 )
R1[1, 11 ] <- 1
R1[2, 24 ] <- 1
R1[3, 37 ] <- 1
R1[4, 50 ] <- 1
R1[5, 63 ] <- 1
linearHypothesis( fitsur, R1, test = "Chisq" )   #检验

蓝色 发表于 2015-1-14 16:59
根据你的提示，是5个方程，每个方程包括截距项是13个参数，因此，估计完有5*13个参数
因此，你按照前面的例 ...

谢谢老师，这个思路我昨天就试过了，不行

第一个检验没有问题啊


> linearHypothesis(fitsur, R1, test = "Chisq")
Linear hypothesis test (Chi^2 statistic of a Wald test)

Hypothesis:
hrearn_married = 0
hrvac_married = 0
hrsick_married = 0
hrins_married = 0
hrpens_married = 0

Model 1: restricted model
Model 2: fitsur

  Res.Df Df  Chisq Pr(>Chisq)  
1   3020                       
2   3015  5 14.177    0.01452 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>

蓝色 发表于 2015-1-15 07:56
第一个检验没有问题啊

太谢谢老师了。再次感谢。

请问约束条件中有常数项怎么写矩阵呢？比如我想写两个变量系数之和为1 的假设