gretl version 1.9.12 Copyright Ramu Ramanathan, Allin Cottrell and Riccardo "Jack" Lucchetti This is free software with ABSOLUTELY NO WARRANTY Current session: 2014-06-16 11:14 ? 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 Successive criterion values within tolerance (1e-007) Convergence achieved after 9 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.41645 3.29456 0.4299 0.6672 mod_AGE_log 1.26453 0.466693 2.710 0.0067 *** mod_BirthCount~_2 -1.57710 0.335183 -4.705 2.54e-06 *** mod_BirthCount~_3 -0.186075 0.131108 -1.419 0.1558 mod_CITIZEN_2 0.694781 0.756394 0.9185 0.3583 mod_CITIZEN_3 -0.297287 0.131505 -2.261 0.0238 ** mod_coord_x 0.0799813 0.302393 0.2645 0.7914 mod_coord_y -0.417008 0.343800 -1.213 0.2252 mod_DEG_URB_1 -0.320124 0.149825 -2.137 0.0326 ** mod_DEG_URB_2 -0.422517 0.121690 -3.472 0.0005 *** mod_FatCountry_1 0.440811 0.304247 1.449 0.1474 mod_FatCountry_2 0.435019 0.378258 1.150 0.2501 mod_FatCountry_3 0.339932 0.310094 1.096 0.2730 mod_FatCountry_4 0.853548 0.513382 1.663 0.0964 * mod_HATLEVEL_0 3.45769 2.43213 1.422 0.1551 mod_HATLEVEL_1 2.96436 2.40620 1.232 0.2180 mod_HATLEVEL_2 2.26377 2.39630 0.9447 0.3448 mod_HATLEVEL_3 2.81118 2.39551 1.174 0.2406 mod_HATLEVEL_4 2.56400 2.39923 1.069 0.2852 mod_HATLEVEL_5 4.28110 2.39631 1.787 0.0740 * mod_HATLEVEL_6 6.66804 2.43308 2.741 0.0061 *** mod_HHLABOUR_EMP~ 1.77494 0.136298 13.02 9.11e-039 *** mod_HHLABOUR_NEM~ -0.423910 0.0774400 -5.474 4.40e-08 *** mod_HHNBPERS_0_4~ -0.271204 0.159276 -1.703 0.0886 * mod_HHNBPERS_14_~ 0.0387695 0.190035 0.2040 0.8383 mod_HHNBPERS_16_~ -0.214398 0.148256 -1.446 0.1481 mod_HHNBPERS_25_~ 0.150883 0.211607 0.7130 0.4758 mod_HHNBPERS_5_1~ -0.127727 0.158582 -0.8054 0.4206 mod_HHNBPERS_65p~ 1.00878 0.163382 6.174 6.64e-010 *** mod_HHNBPERS_sq 1.14638 0.382680 2.996 0.0027 *** mod_HHTYPE_1 -0.654556 0.180461 -3.627 0.0003 *** mod_HHTYPE_2 -0.656514 0.225397 -2.913 0.0036 *** mod_HHTYPE_3 0.780312 0.127393 6.125 9.05e-010 *** mod_HHTYPE_4 0.379911 0.141605 2.683 0.0073 *** mod_HHTYPE_5 0.128512 0.182886 0.7027 0.4822 mod_ien_adm_log 0.0899976 0.0136540 6.591 4.36e-011 *** mod_JOBSTAT_10 -0.475185 0.469934 -1.011 0.3119 mod_JOBSTAT_21 -1.16884 0.456024 -2.563 0.0104 ** mod_JOBSTAT_22 -2.70888 0.553073 -4.898 9.69e-07 *** mod_MAINSTAT_11 1.39176 0.171396 8.120 4.66e-016 *** mod_MAINSTAT_20 -2.67110 0.443714 -6.020 1.75e-09 *** mod_MAINSTAT_31 -0.298274 0.537071 -0.5554 0.5786 mod_MAINSTAT_32 0.529301 0.454533 1.164 0.2442 mod_MAINSTAT_33 0.638843 0.512178 1.247 0.2123 mod_MAINSTAT_35 -1.42882 0.457348 -3.124 0.0018 *** mod_MARSTADEFA~_2 0.344734 0.139032 2.480 0.0132 ** mod_MARSTALEGAL_1 0.510691 0.166273 3.071 0.0021 *** mod_MARSTALEGAL_2 0.695475 0.154483 4.502 6.73e-06 *** mod_MARSTALEGAL_3 0.450898 0.149363 3.019 0.0025 *** mod_MEHMHS1_1 0.329761 0.290499 1.135 0.2563 mod_MEHMHS1_2 0.498487 0.262063 1.902 0.0571 * mod_MEHMHS1_3 -0.205904 0.244046 -0.8437 0.3988 mod_MEHMHS1_4 -0.101133 0.238846 -0.4234 0.6720 mod_MEHMHS2_1 0.341278 0.0992990 3.437 0.0006 *** mod_MEHMHS3_1 -0.614606 0.163974 -3.748 0.0002 *** mod_MEHMHS3_2 -0.801920 0.102015 -7.861 3.82e-015 *** mod_mode_1_CAPI -0.149414 0.0804638 -1.857 0.0633 * mod_MotCountry_1 -0.542419 0.471768 -1.150 0.2502 mod_MotCountry_2 0.00719308 0.515575 0.01395 0.9889 mod_MotCountry_3 -0.632689 0.479735 -1.319 0.1872 mod_MotCountry_4 0.138608 0.641596 0.2160 0.8290 mod_PSUstrata_1 1.13982 0.195226 5.838 5.27e-09 *** mod_PSUstrata_2 0.633040 0.146076 4.334 1.47e-05 *** mod_PSUstrata_3 0.113762 0.123229 0.9232 0.3559 mod_reg_code_LV0~ 0.195513 0.196994 0.9925 0.3210 mod_reg_code_LV0~ -1.07957 0.180572 -5.979 2.25e-09 *** mod_reg_code_LV0~ 0.687220 0.153410 4.480 7.48e-06 *** mod_reg_code_LV0~ 0.358410 0.171419 2.091 0.0365 ** mod_SEX_1 0.534638 0.0754601 7.085 1.39e-012 *** mod_SSUstrata_11 1.77168 4.11376 0.4307 0.6667 mod_SSUstrata_12 -0.213577 0.520293 -0.4105 0.6814 mod_SSUstrata_21 2.29726 0.619723 3.707 0.0002 *** mod_vcgr3_2 1.58065 0.447488 3.532 0.0004 *** mod_vcgr3_3 1.42191 0.339212 4.192 2.77e-05 *** mod_vcgr3_4 0.304754 0.269307 1.132 0.2578 mod_vcgr3_5 0.343228 0.213083 1.611 0.1072 Log-likelihood 13047.23 Akaike criterion -25940.46 Schwarz criterion -25422.45 Hannan-Quinn -25760.82 sigma = 0.000622456 (5.87071e-006) Left-unbounded observations: 0 Right-unbounded observations: 12 Bounded observations: 2724 Point observations: 3434 Test for normality of residual - Null hypothesis: error is normally distributed Test statistic: Chi-square(2) = 168.761 with p-value = 2.25897e-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 store: using filename T:/!Vadiba/MatNodrD/!Info/Metodologija/Int_Reg/EHSIS2012_example/data/ymod.txt Data written OK. Done