gretl version 1.9.13 Copyright Ramu Ramanathan, Allin Cottrell and Riccardo "Jack" Lucchetti This is free software with ABSOLUTELY NO WARRANTY Current session: 2014-06-16 11:15 ? run rungretl.inp T:\!Vadiba\MatNodrD\!Info\Metodologija\Int_Reg\EHSIS2012_example\data\rungretl.inp # EHSIS 2012 # Martins Liberts # Generated file from R ? open dat.mod.csv parsing T:\!Vadiba\MatNodrD\!Info\Metodologija\Int_Reg\EHSIS2012_example\data\dat.mod.csv... using delimiter ',' longest line: 1385 characters first field: 'ID' number of columns = 81 number of variables: 81 number of non-blank lines: 6171 scanning for variable names... line: ID,WEIGHT,IR,lb,ub,mod_AGE_log,mod_BirthCountry_2,mod_BirthCountry_3,mod_CITIZEN_2,mod_CITIZEN_3,mod_coord_x,mod_co... scanning for row labels and data... warning: missing value for variable 5, obs 9 warning: missing value for variable 5, obs 288 warning: missing value for variable 5, obs 728 warning: missing value for variable 5, obs 1129 warning: missing value for variable 5, obs 1444 treating these as undated data Listing 82 variables: 0) const 1) ID 2) WEIGHT 3) IR 4) lb 5) ub 6) mod_AGE_log 7) mod_BirthCountry_2 8) mod_BirthCountry_3 9) mod_CITIZEN_2 10) mod_CITIZEN_3 11) mod_coord_x 12) mod_coord_y 13) mod_DEG_URB_1 14) mod_DEG_URB_2 15) mod_FatCountry_1 16) mod_FatCountry_2 17) mod_FatCountry_3 18) mod_FatCountry_4 19) mod_HATLEVEL_0 20) mod_HATLEVEL_1 21) mod_HATLEVEL_2 22) mod_HATLEVEL_3 23) mod_HATLEVEL_4 24) mod_HATLEVEL_5 25) mod_HATLEVEL_6 26) mod_HHLABOUR_EMP_sq 27) mod_HHLABOUR_NEMP_sq 28) mod_HHNBPERS_0_4_sq 29) mod_HHNBPERS_14_15_sq 30) mod_HHNBPERS_16_24_sq 31) mod_HHNBPERS_25_64_sq 32) mod_HHNBPERS_5_13_sq 33) mod_HHNBPERS_65plus_sq 34) mod_HHNBPERS_sq 35) mod_HHTYPE_1 36) mod_HHTYPE_2 37) mod_HHTYPE_3 38) mod_HHTYPE_4 39) mod_HHTYPE_5 40) mod_ien_adm_log 41) mod_JOBSTAT_10 42) mod_JOBSTAT_21 43) mod_JOBSTAT_22 44) mod_MAINSTAT_11 45) mod_MAINSTAT_20 46) mod_MAINSTAT_31 47) mod_MAINSTAT_32 48) mod_MAINSTAT_33 49) mod_MAINSTAT_35 50) mod_MARSTADEFACTO_2 51) mod_MARSTALEGAL_1 52) mod_MARSTALEGAL_2 53) mod_MARSTALEGAL_3 54) mod_MEHMHS1_1 55) mod_MEHMHS1_2 56) mod_MEHMHS1_3 57) mod_MEHMHS1_4 58) mod_MEHMHS2_1 59) mod_MEHMHS3_1 60) mod_MEHMHS3_2 61) mod_mode_1_CAPI 62) mod_MotCountry_1 63) mod_MotCountry_2 64) mod_MotCountry_3 65) mod_MotCountry_4 66) mod_PSUstrata_1 67) mod_PSUstrata_2 68) mod_PSUstrata_3 69) mod_reg_code_LV003 70) mod_reg_code_LV005 71) mod_reg_code_LV007 72) mod_reg_code_LV008 73) mod_SEX_1 74) mod_SSUstrata_11 75) mod_SSUstrata_12 76) mod_SSUstrata_21 77) mod_vcgr3_2 78) mod_vcgr3_3 79) mod_vcgr3_4 80) mod_vcgr3_5 81) w ? intreg lb ub w mod_AGE_log mod_BirthCountry_2 mod_BirthCountry_3 \ mod_CITIZEN_2 mod_CITIZEN_3 mod_coord_x mod_coord_y mod_DEG_URB_1 \ mod_DEG_URB_2 mod_FatCountry_1 mod_FatCountry_2 mod_FatCountry_3 \ mod_FatCountry_4 mod_HATLEVEL_0 mod_HATLEVEL_1 mod_HATLEVEL_2 \ mod_HATLEVEL_3 mod_HATLEVEL_4 mod_HATLEVEL_5 mod_HATLEVEL_6 mod_HHLABOUR_EMP_sq \ mod_HHLABOUR_NEMP_sq mod_HHNBPERS_0_4_sq mod_HHNBPERS_14_15_sq \ mod_HHNBPERS_16_24_sq mod_HHNBPERS_25_64_sq mod_HHNBPERS_5_13_sq \ mod_HHNBPERS_65plus_sq mod_HHNBPERS_sq mod_HHTYPE_1 mod_HHTYPE_2 \ mod_HHTYPE_3 mod_HHTYPE_4 mod_HHTYPE_5 mod_ien_adm_log mod_JOBSTAT_10 \ mod_JOBSTAT_21 mod_JOBSTAT_22 mod_MAINSTAT_11 mod_MAINSTAT_20 \ mod_MAINSTAT_31 mod_MAINSTAT_32 mod_MAINSTAT_33 mod_MAINSTAT_35 \ mod_MARSTADEFACTO_2 mod_MARSTALEGAL_1 mod_MARSTALEGAL_2 mod_MARSTALEGAL_3 \ mod_MEHMHS1_1 mod_MEHMHS1_2 mod_MEHMHS1_3 mod_MEHMHS1_4 mod_MEHMHS2_1 \ mod_MEHMHS3_1 mod_MEHMHS3_2 mod_mode_1_CAPI mod_MotCountry_1 mod_MotCountry_2 \ mod_MotCountry_3 mod_MotCountry_4 mod_PSUstrata_1 mod_PSUstrata_2 \ mod_PSUstrata_3 mod_reg_code_LV003 mod_reg_code_LV005 mod_reg_code_LV007 \ mod_reg_code_LV008 mod_SEX_1 mod_SSUstrata_11 mod_SSUstrata_12 \ mod_SSUstrata_21 mod_vcgr3_2 mod_vcgr3_3 mod_vcgr3_4 mod_vcgr3_5 Warning: couldn't improve criterion (gradient = 14.3017) Convergence achieved after 14 iterations Model 1: Interval estimates, using observations 1-6170 Lower limit: lb, Upper limit: ub Standard errors based on Hessian coefficient std. error z p-value ---------------------------------------------------------------- w 1.44937 3.29022 0.4405 0.6596 mod_AGE_log 1.25609 0.466043 2.695 0.0070 *** mod_BirthCount~_2 -1.56769 0.334717 -4.684 2.82e-06 *** mod_BirthCount~_3 -0.186957 0.130927 -1.428 0.1533 mod_CITIZEN_2 0.688109 0.755141 0.9112 0.3622 mod_CITIZEN_3 -0.296486 0.131323 -2.258 0.0240 ** mod_coord_x 0.0814391 0.301973 0.2697 0.7874 mod_coord_y -0.413898 0.343321 -1.206 0.2280 mod_DEG_URB_1 -0.317988 0.149615 -2.125 0.0336 ** mod_DEG_URB_2 -0.422131 0.121522 -3.474 0.0005 *** mod_FatCountry_1 0.438014 0.303823 1.442 0.1494 mod_FatCountry_2 0.431428 0.377731 1.142 0.2534 mod_FatCountry_3 0.337332 0.309661 1.089 0.2760 mod_FatCountry_4 0.851586 0.512679 1.661 0.0967 * mod_HATLEVEL_0 3.45562 2.42910 1.423 0.1549 mod_HATLEVEL_1 2.96324 2.40321 1.233 0.2176 mod_HATLEVEL_2 2.26163 2.39333 0.9450 0.3447 mod_HATLEVEL_3 2.80804 2.39254 1.174 0.2405 mod_HATLEVEL_4 2.55969 2.39626 1.068 0.2854 mod_HATLEVEL_5 4.27765 2.39334 1.787 0.0739 * mod_HATLEVEL_6 6.67154 2.43004 2.745 0.0060 *** mod_HHLABOUR_EMP~ 1.77375 0.136110 13.03 8.08e-039 *** mod_HHLABOUR_NEM~ -0.423035 0.0773339 -5.470 4.49e-08 *** mod_HHNBPERS_0_4~ -0.272164 0.159058 -1.711 0.0871 * mod_HHNBPERS_14_~ 0.0378624 0.189772 0.1995 0.8419 mod_HHNBPERS_16_~ -0.217054 0.148052 -1.466 0.1426 mod_HHNBPERS_25_~ 0.150488 0.211317 0.7121 0.4764 mod_HHNBPERS_5_1~ -0.129119 0.158365 -0.8153 0.4149 mod_HHNBPERS_65p~ 1.00793 0.163159 6.178 6.51e-010 *** mod_HHNBPERS_sq 1.14892 0.382157 3.006 0.0026 *** mod_HHTYPE_1 -0.654304 0.180208 -3.631 0.0003 *** mod_HHTYPE_2 -0.656100 0.225087 -2.915 0.0036 *** mod_HHTYPE_3 0.780548 0.127216 6.136 8.48e-010 *** mod_HHTYPE_4 0.381043 0.141411 2.695 0.0070 *** mod_HHTYPE_5 0.128856 0.182634 0.7055 0.4805 mod_ien_adm_log 0.0900776 0.0136355 6.606 3.95e-011 *** mod_JOBSTAT_10 -0.473545 0.469232 -1.009 0.3129 mod_JOBSTAT_21 -1.16563 0.455336 -2.560 0.0105 ** mod_JOBSTAT_22 -2.70560 0.552263 -4.899 9.63e-07 *** mod_MAINSTAT_11 1.39014 0.171161 8.122 4.59e-016 *** mod_MAINSTAT_20 -2.66690 0.443041 -6.020 1.75e-09 *** mod_MAINSTAT_31 -0.296127 0.536293 -0.5522 0.5808 mod_MAINSTAT_32 0.530521 0.453849 1.169 0.2424 mod_MAINSTAT_33 0.638550 0.511423 1.249 0.2118 mod_MAINSTAT_35 -1.42941 0.456656 -3.130 0.0017 *** mod_MARSTADEFA~_2 0.344014 0.138838 2.478 0.0132 ** mod_MARSTALEGAL_1 0.510722 0.166040 3.076 0.0021 *** mod_MARSTALEGAL_2 0.694639 0.154271 4.503 6.71e-06 *** mod_MARSTALEGAL_3 0.450344 0.149157 3.019 0.0025 *** mod_MEHMHS1_1 0.328696 0.290101 1.133 0.2572 mod_MEHMHS1_2 0.499013 0.261706 1.907 0.0566 * mod_MEHMHS1_3 -0.205434 0.243714 -0.8429 0.3993 mod_MEHMHS1_4 -0.100799 0.238520 -0.4226 0.6726 mod_MEHMHS2_1 0.342953 0.0991617 3.459 0.0005 *** mod_MEHMHS3_1 -0.615486 0.163750 -3.759 0.0002 *** mod_MEHMHS3_2 -0.801992 0.101874 -7.872 3.48e-015 *** mod_mode_1_CAPI -0.151349 0.0803523 -1.884 0.0596 * mod_MotCountry_1 -0.538746 0.470991 -1.144 0.2527 mod_MotCountry_2 0.00954809 0.514755 0.01855 0.9852 mod_MotCountry_3 -0.628205 0.478951 -1.312 0.1896 mod_MotCountry_4 0.141909 0.640618 0.2215 0.8247 mod_PSUstrata_1 1.13533 0.194952 5.824 5.76e-09 *** mod_PSUstrata_2 0.630457 0.145871 4.322 1.55e-05 *** mod_PSUstrata_3 0.112854 0.123059 0.9171 0.3591 mod_reg_code_LV0~ 0.195213 0.196723 0.9923 0.3210 mod_reg_code_LV0~ -1.08058 0.180324 -5.992 2.07e-09 *** mod_reg_code_LV0~ 0.685181 0.153197 4.473 7.73e-06 *** mod_reg_code_LV0~ 0.357613 0.171183 2.089 0.0367 ** mod_SEX_1 0.534240 0.0753557 7.090 1.35e-012 *** mod_SSUstrata_11 1.76300 4.10792 0.4292 0.6678 mod_SSUstrata_12 -0.214340 0.519556 -0.4125 0.6799 mod_SSUstrata_21 2.28746 0.618856 3.696 0.0002 *** mod_vcgr3_2 1.57136 0.446858 3.516 0.0004 *** mod_vcgr3_3 1.41904 0.338744 4.189 2.80e-05 *** mod_vcgr3_4 0.301911 0.268935 1.123 0.2616 mod_vcgr3_5 0.341756 0.212789 1.606 0.1083 Log-likelihood 13056.57 Akaike criterion -25959.13 Schwarz criterion -25441.12 Hannan-Quinn -25779.50 sigma = 0.000621572 (5.86566e-006) Left-unbounded observations: 0 Right-unbounded observations: 12 Bounded observations: 2723 Point observations: 3434 Pseudo-point observations: 1 Test for normality of residual - Null hypothesis: error is normally distributed Test statistic: Chi-square(2) = 168.303 with p-value = 2.84097e-037 ? yhat1 = $yhat Generated series yhat1 (ID 82) ? yres1 = $uhat Generated series yres1 (ID 83) ? matrix b = $coeff Generated matrix b ? err = mwrite(b, \ "T:/!Vadiba/MatNodrD/!Info/Metodologija/Int_Reg/EHSIS2012_example/data/coeff \ .txt") Generated scalar err = 0 ? store T:/!Vadiba/MatNodrD/!Info/Metodologija/Int_Reg/EHSIS2012_example/data/ymod.txt ID \ w yhat1 yres1 wrote T:/!Vadiba/MatNodrD/!Info/Metodologija/Int_Reg/EHSIS2012_example/data/ymod.txt Done