Roteiro de Aula Econometria - Testes F e LM, Elasticidades e Operaçoes Matriciais

Melhor visualizado com fonte 16 edit GUI preferences size 16 O objetivo deste roteiro de aula comparar os testes F e LM para o teste de restri o de exclus o apresentar a interpreta o de elasticidades e semielasticidades e apresentar opera es matriciais S o utilizados os mesmos exemplos do laborat rio 1 como segue abaixo Regress o que ser estimada equa o de Mincer Sal rio hor rio real 0 1educa o 2experi ncia 3experi ncia2 4etnia 5sexo 6setor 7ocupa o 8regi o 9urbano u Fontes GIUBERTI Ana Carolina MENEZESFILHO Na rcio Discrimina o de rendimentos por g nero uma compara o entre o Brasil e os Estados Unidos Econ Apl Ribeir o Preto v 9 n 3 p 369384 Set 2005 Meara K Pastore F Webster A 2020 The gender pay gap in the USA a matching study Journal of Population Economics 331 271305 Limpando a rea de trabalho acionar o comando abaixo com ctrlr rm list ls Carregando os dados Para carregar os dados que ja est o dispon veis em seu computador siga os passos A Procedimento manual 1 Clique no topo da janela R Console acionandoa 2 Selecione file load workspace 3 V at o diret rio da desktop em que voc salvou a pasta lista1 descompactada 4 Selecione o arquivo lista1Rdata B Procedimento autom tico Copiar o endere o do diret rio em que se encontra o arquivo dadosRdata Substituir por no endere o como separador das pastas Exemplo setwd DUFABCEnsinoEconometriaII2023Rlistaslab Acionar o comando abaixo com ctrlr load lista1 dadosRdata 1 Regress o Minceriana Relembrando vari veis da regress o Vari vel dependente Sal rio hor rio Lista de explicativas independentes Bloco 1 Educa o e experi ncia 1 Anos de estudo anoseduc 2 Experi ncia experidadetrab idade idade em que come ou a trabalhar 3 Experi ncia ao quadrado experidadetrabsq Bloco 2 Caracter sticas pessoais 4 Etnia preta ou parda etniaprepar 5 Sexo masculino dhom Bloco 3 Setor e ocupa o 6 Setor de atividade agropecu ria ind stria constru o Educa o Sa de Social Outros servi os administra o p blica dativ 7 Ocupa o administra o agricultura com rcio outros servi os dirigente ci ncia tecnologia t cnico de n vel m dio reparo manuten o for as armadas docup Bloco 4 Regional 8 Grande regi o brasileiras Nordeste 2 Sudeste 3 Sul 4 Centro Oeste 5 dreg2dreg5 9 rea urbana durbano Sum rio estat stico das vari veis O nome do objeto com o conjunto de dados dbmin2015 Para observar o resumo estat stico das vari veis do modelo Acione a linha abaixo selecionandoa e apertando ctrlR summarysubsetdbmin2015selectcsalhoranoseducexperidadetrabexperidadetrabsqdetniaprepardhomdativagrodativindusdativconstdativsocialdativservdativapubdocupadmdocupagrodocupcomdocupservdocupdirigentedocupcieartdocuptecmeddocuprepmanutdocuparmadadreg2dreg3dreg4dreg5durbano Quest o 1Quest o 1Quest o 1Quest o 1 Sob as hip teses do modelo cl ssico de regress o linear a matriz de vari nciacovari ncia tem forma simples Tratase de vcov sigma 2 XX1 Sendo sigma 2 a estimativa pontual para a vari ncia do termo de perturba o Esta corresponde soma dos quadrados dos res duos dividida pelo respectivo n mero de graus de liberdade NK1 Vamos calcular a partir de tr s opera es Antes disso preciso gerar matriz X Definindo a matriz X Inicialmentese selecionase estritamente as vari veis explicativas do modelo Para isso basta utilizar o comando subset com op o select matxsubsetdbmin2015selectcanoseducexperidadetrabexperidadetrabsqdetniaprepardhomdativagrodativindusdativconstdativsocialdativservdativapubdocupadmdocupagrodocupcomdocupservdocupdirigentedocupcieartdocuptecmeddocuprepmanutdocuparmadadreg2dreg3dreg4dreg5durbano Agora preciso adicionar o vetor com valores unit rios Gerando um vetor 1 x N de 1s vec1 matrix 1nrow nrow dbmin2015 ncol 1 Acoplando horizontalmente matx e vec1 importante vec1 tem de ser a primeira coluna da esquerda para a direita obrigatoriamente matx cbind vec 1mat x Visualizando View matx De fato o vetor de 1s est na primeira coluna Verificando dimens o da matriz X dim matx Notar que h exatamente 26 colunas exatamente o n mero de par metros do modelo at ent o estimado contendo o intercepto summary matx necess rio transfomar matx em um objeto matricial para poder fazer as opera es matriciais matx asmatrix matx Definindo o vetor Y vecy asmatrix dbmin2015salhornrownrowdbmin2015ncol1 Verificando dimens o do vetor Y dim vecy Notar que o n mero de linhas o mesmo da matriz X Resumo estat stico acionar da primeira terceira linha que seguem abaixo rbind c mean vecy sd vecy tapplysubsetmat xselect canoseducexperidadetrabdetniaprepardativagrodreg2durbano2functionxcmeanxnarmTsdxnarmT Quest o 1bQuest o 1bQuest o 1bQuest o 1b A opera o matricial que gera os estimadores de MQO bmqo XX 1 XY Sendo bmqo o vetor N x 1 de estimadores para todos os par metros da FRP intercepto e coeficientes X matriz N x K1 de vari veis explicativas com uma primeira coluna contendo apenas n meros unit rios para com isso incorporar o intercepto estima o Y vetor N x 1 de valores da vari vel dependente poss vel obter diretamente o vetor de estimativas pontuais realizando a opera o matricial anterior no R o que feito a seguir Opera o 1 XX Notar que a matriz esquerda a transposi o da matriz X A fun o t transp e no R O operador aplica o produto matricial N o se trata do produto escalar em que s o multiplicados os elementos em posi es equivalentes das duas matrizes op1t mat x matx Opera o 2 XX1 A fun o solve inverte matrizes no R Essa sem d vida a opera o que mais exige do computador op2solveop1 Opera o 3 XY Tratase do produto matricial abaixo Aten o nunca alterar a ordem dos fatores de uma multiplica o matricial op3t matx vecy Importante notar a transposi o da matriz X sem a qual a multiplica o seria imposs vel Basta observar as dimens es de matx e vecy abaixo dim matx dim vecy Agora com matx transposta dim t matx dim vecy O n mero de colunas de matx transposta equivalente ao n mero de linhas de vecy como necess rio Opera o 4 opera o 2 x opera o 3 Recordando bmqo XX1XY op 2 op3 Novamente tratase de produto matricial op4op2 op3 Visualizando os resultados op4 View op4 Quest o 2aQuest o 2aQuest o 2aQuest o 2a Calculando a SQR Passo 1 res duos da regress o opres opres vecy matx op4 Passo 2 elevando os res duos ao quadrado e somando o resultado SQRsumopres2 Estimador da vari ncia dos termos de perturba o sigma 2 SQRNK1 Tamanho da amostra e n mero de explicativas N91363 K25 SQRgl SQRNK1 Calculando XX1 J foi calculado anteriormente e resultado se encontra no objeto op2 Multiplicando SQRgl e XX1 Notar que nesse caso como SQRgl um escalar devese utilizar o produto escalar varcovmat SQRgl op2 Dimens o da matriz gerada dim op3 Est correta a dimens o O objeto varcovmat cont m a matriz de vari nciacovari ncia dos estimadores de MQO As vari ncias est o na diagonal principal a qual extra da com o comando diag varmat diag varcovmat Quest o 2bQuest o 2bQuest o 2bQuest o 2b 1F rmula da regress o acionar linha abaixo formulasalhoranoseducexperidadetrabexperidadetrabsqdetniaprepardhomdativagrodativindusdativconstdativsocialdativservdativapubdocupadmdocupagrodocupcomdocupservdocupdirigentedocupcieartdocuptecmeddocuprepmanutdocuparmadadreg2dreg3dreg4dreg5durbano 2Regress o acionar linha abaixo mqo lm formuladata dbmin2015 summary mqo Comparando os vetores de vari ncia matricial e lm O vetor matricial se encontra no objeto varmat J o vetor de vari ncias gerado pelo comando lm corresponde diagonal principal da matriz vcov mqo E isso pois vcov mqo a matriz de vari nciacovari ncia gerada pelo comando lm Passo 1 subtra o elemento a elemento dos dois vetores varmat diag vcov mqo Visualizando apropriadamente View varmat diag vcov mqo notar que as diferen as s o desprez veis com ordem de magnitude de 1011 ou menor De fato a diferen a absoluta m xima max abs varmat diag vcov mqo Um n mero de ordem de magnitude 12 e pois desprez vel Com certeza a opera o matricial realizada equivalente utilizada pelo comando lm A soma das diferen as elemento a elemento atesta isso sum abs varmat diag vcov mqo Ou seja as diferen as somam um n mero desprez vel Passo 2 comparar a diagonal da matriz op3 com os errospadr o dos estimadores segundo comando lm Tabelaresposta ao item 2b cbinditem2avar matautom ticadiagvcovmqodiferen avarmatdiagvcovmqo Fica n tido pois a aus ncia de diferen a entre as vari ncia dos estimadores obtidas com as opera es manuais e via lm Quest o 3aQuest o 3aQuest o 3aQuest o 3a Calculando as previs es As estimativas para os par metros est o no objeto op4 Passo 1 Previs es para todos os indiv duos da amostra prev matx op4 Passo 2 estat sticas das previs es cbind M dia mean prev Mediana median prev DP sd prev Passo 3 m dia para Y observado mean dbmin2015salhor Quest o 4aQuest o 4aQuest o 4aQuest o 4a Procedimento de Wooldrige Passo 1 gerar vetor X base para a previs o vecx0cbindanoseduc10experidadetrab20experidadetrabsq202detniaprepar1dhom1dativagro0dativindus0dativconst0dativsocial0dativserv1dativapub0docupadm0docupagro0docupcom0docupserv1docupdirigente0docupcieart0docuptecmed0docuprepmanut0docuparmada0dreg20dreg31dreg40dreg50durbano1 vecx0 Passo 2 modifica o a matrix X pela subtra o de vecx0 Passo 2a em primeiro lugar preciso criar uma matriz que tem em todas as linhas o vetor x0 todas as linhas s o iguais e equivalentes ao vetor x0 Antes de continuar pense como voc geraria essa matriz a partir de uma opera o matricial A opera o matricial a multiplica o de um vetor vertical de n meros um com N linhas pela vetor x0 este ltimo com 1 linha e K1 colunas Cabe notar que tal opera o dimensionalmente correta representando o vetor de n meros um por V V N x 1 x0 1 x K1 resultado N x K1 matx0 matrix 1nrow nrow matx ncol 1vecx0 Passo 2b em segundo lugar subtrair a matriz X pela matriz x0 por m preciso acrescentar uma primeira coluna de zeros na ltima pois n o se trata de subtrair o intercepto matxmod matx cbind 0mat x0 Passo 2c gerando a base de dados com Y original e X modificada dbmin2015modas dataframe cbindsalhordbmin2015salhormatxmod Passo 3 estimando a regress o modificada mqomod lm formuladata dbmin2015mod Passo 4 calculando o IC Passo 4a extraindo o desviopadr o do intercepto da regress o modificada A vari ncia do intercepto se encontra na primeira linha e primeira coluna da matriz de vari nciacovari ncia do objeto mqomoq A f rmula abaixo calcula a raiz de tal elemento sdintmod vcov mqo mod 1112 Passo 4b coeficiente de confian a na distribui o t Tal coeficiente obtido utilizando a fun o qt pgl a qual informa o n mero esquerda do qual h uma probabilidade p na distribui o t de Student com NK1 graus de liberdade tconf qt 0025NK1 Passo 4c c lculo do IC conforme comando abaixo ccoefficientsmqo mod 1tconfsdintmodcoefficientsmqomod1tconfsdintmod Quest o 4bQuest o 4bQuest o 4bQuest o 4b Procedimento alg brico Passo 1 c lculo da vari ncia da previs o Tratase de realizar a seguinte opera o matricial x0Vb x0 Observar as dimens es dos elementos multiplicados vecx0 1 x K1 vcov mqo K1 x K1 vecx0 K1 x 1 A opera o feita abaixo varprev cbind 1vec x0 vcov mqo t cbind 1vecx0 A opera o acima gerou a vari ncia agora preciso calcular a raiz para com isso obter o desviopadr o sdprev varprev 12 Passo 2 c lculo do IC conforme comando abaixo ccoefficientsmqo mod 1tconfsdprevcoefficientsmqomod1tconfsdprev Notar que no comando acima sdintmod foi substitu do por sdprev Comparando apenas os desviopadr o Come ando pelo obtido no item 4a sdintmod Agora o obtido neste item sdprev Obtendo a diferen a sdintmodsdprev SCRIPT LISTA 2 rm list ls DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS Antes de come ar a rodar este script alterar as instru es que especificam os diret rios Diret rio e arquivo de dados setwd DUFABCEnsinoEconometriaII2023Rlistaslab load lista2 dadosRdata EXPLICATIVAS 1 Escolaridade 30 anoseduc 2 Experiencia 61 exper 3 Experiencia ao quadrado 62 expersq 4 Horas trabalhadas 24 horastrabprisref 5 Etnia Brancoa 43 dbca 6 Regiao brasileira 464850 dnodse dsu e dco 7 Setor de atividade5254 dindudapubdserv 8 Emprego formal 55 dformal 9 Numero de filhos 2629 filhosmdom filhosfdom filhosmndom filhosfndom 10 Renda domiciliar 14 rendadompc 11 Reside em area rural 59 drural Um pouco sobre a rationale destas variaveis a renda domiciliar e uma proxy para a renda gerada pela familia como um todo e portanto nao apenas pela mulher Quanto maior a renda gerada por outros membros da familia na maioria dos casos menor tende a ser a probabilidade da mulher trabalhar pois menor e a necessidade de complementar a rendaE claro que este efeito tende a se verificar apenas para mulheres casadas nao prevalecendo para familias em que a mulher e o unico chefe Para selecionar as variaveis rodar o comando abaixo Quest o 1a Passo 1 rodando regress o de Mincer formularendatrabpri10anoseducexperexpersqhorastrabprisref dbca dindu dapub dserv dno dse dsu dco dformalnfilhos rendadompc drural mqo lm formuladata dmicrohfe summary mqo dim dmicrohfe Incorporando o vetor com o quadrado dos res duos no conjunto de dados dmicrohferessq mqoresiduals2 Passo 2 quest o 1a exame indicativo correla es com algumas vari veis altamente significativas cortest dmicrohferes sqdmicro hfeanoseduc pvalue 5 cortest dmicrohferes sqdmicro hfeexper pvalue 5 cortest dmicrohferes sqdmicro hfehorastrabprisref pvalue 5 Quest o 1b Passo 3 teste de BreuschPagan Passo 3a rodando regress o auxiliar formula ressqanoseducexperexpersqhorastrabprisref dbca dindu dapub dserv dno dse dsu dco dformalnfilhos rendadompc drural auxbp lm formuladata dmicrohfe Passo 3b observando os resultados do teste de signific ncia global teste F Este teste reportando automaticamente na ltima linha do sum rio do objeto com a regress o Basta acionar o comando abaixo summary auxbp A hip tese nula de aus ncia de signific ncia global equivalente hip tese de homocedasticidade Observe a estat stica F e o pvalor dela e conclua v lida a hip tese de homocedasticidade Quest o 1c Passo 4 teste de White Passo 4a o desafio gerar todas os termos quadr ticos e mistos para as 16 vari veis explicativas Tratase de um desafio computacional pois manualmente uma opera o repetitiva e sujeita a erro formularessq2Ianoseduc2Iexper2Iexpersq2Ihorastrabprisref2Idbca2Idindu2Idapub2Idserv2Idno2Idse2Idsu2Idco2Idformal2Infilhos2Irendadompc2Idrural2 dbwhite subset dmicrohfeselectcressqanoseducexperexpersqhorastrabprisref dbca dindu dapub dserv dno dse dsu dco dformalnfilhos rendadompc drural auxwhite lm formuladata dbwhite Passo 4b observando o resultado Basta novamente examinar o resultado do teste F na ltima linha summary auxwhite Quest o 2 Estimador de White Passo 1 instalar o pacote AER installpackages AER Passo 2 habilitar o pacote library AER Passo 3 estimar o modelo sem corrigir os estimadores dos desviospadr o para a heterocedastidade formularendatrabpri10anoseducexperexpersqhorastrabprisref dbca dindu dapub dserv dno dse dsu dco dformalnfilhos rendadompc drural mqo lm formuladata dmicrohfe OBS este passo j foi feito na quest o 1 Passo 4 com base nos resultados de estima o aplicar a corre o de White Para obter resultados de infer ncia robustos para a heterocedasticidade necess rio informar ao R que desejamos utilizar para o c lculo dos valores da estat stica t a matriz de vari nciacovari ncia dos estimadores que incorpora heterocedasticidade O que requer o emprego da instruc o vcov vcovHC s vcov significando matriz de var covar dos estimadores e vcovHC s especificando que desejamos utilizar uma matriz robusta para heterocedasticidade heteroskedasticity robust ou HC O termo s indica o objeto no qual os resultados da regress o foram armazenados no caso mqo Teste de signific ncia individual teste t coeftest mqo vcov vcovHC mqo Passo 3 compara o dos erros padr o com e sem corre o de White cbindsemcorrsqrtdiagvcovmqo comcorr sqrtdiagvcovHCmqodifsqrtdiagvcovmqosqrtdiagvcovHCmqo Quest o 3 multicolinearidade Passo 1 teste de correla o n mero de filhos e experi ncia cortest dmicrohfen filhosdmicro hfeexper Passo 2 teste de correla o n mero de filhos e educa o cortest dmicrohfen filhosdmicro hfeanoseduc Quest o 4 autocorrela o Passo 1 Gerando base de dados em s ries de tempo dmacrots ts dmacro start 1974 end 2010 frequency 1 Sum rio estat stico e visualizac o dos dados macro summary dmacrots plot dmacrots Passo 2 Instalando e habilitando pacote de teste Para rodar a regress o para s ries de tempo preciso instalar e habilitar a pack dynlm o que j foi feito no in cio deste script installpackages dynlm library dynlm O comando que aciona a regress o para s ries de tempo o dynlm Dentro dele podemos especificar diretamente que as vari veis PIB e investimento ser o consideradas em primeiras diferencas e n o em n vel Passo 3 Rodando regress o mqo dynlm d inv d pib taxa data dmacrots summary mqo Passo 4 Realizando o teste gerando o valor da estat stica de Durbin Watson library lmtest dwtest mqo OBS o pvalor deve ser desconsiderado Procurar o pvalor da tabela 4 dim dmacrots N n mero de observa es 201019751 K n mero de vari veis explicativas no modelo SCRIPT LISTA 3 rm list ls ls DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS DIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOSDIRET RIO DE ARQUIVOS Antes de come ar a rodar este script alterar as instru es que especificam os diret rios 1 Obter o endere o do diret rio em que os arquivo dmicrolab2Rdata e dmacropennRdata est o armazenados em seu computador Sugiro que este arquivo seja colocado na pasta metrialab e esta por sua vez seja armazenada na desktop para que assim todos os arquivos que vamos trabalhar durante o curso fiquem armazenados no mesmo diret rio O endere o deve ser algo como o exemplo a seguir diret rio do meu computador CUsersfonsecamorelloDesktopEconometriaIIFluxoUFABCcasa170314ufabcmetrialab 2a O R n o entende o endere o se ele estiver grafado com barras simples ie com preciso substituir as barras simples do diret rio do computador de voc s por barras duplas Al m disso preciso fazer mais dois ajustes 2b o endere o deve terminar com barras duplas 2c o endere o deve ser grafado entre aspas duplas Fazendo os tr s ajustes no endere o do exemplo acima temos CUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa170314ufabcmetrialab 3 Nesta janela em que estamos agora script do R tecle control H e no menu que aparecer digite 4 no campo Localizar CUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa200314ufabcmetrialabusar 5 no campo Substituir por digite o endere o do diret rio corrigido obtido no passo 2 Os dados processados no R ser o salvos neste diret rio Por favor zelar por seus arquivos armazenando o diret rio ao final da aula de laborat rio em pen drives ou sites de hospedagem de dados Google Drive Dropbox etc Eles ser o utilizados nas aulas futuras Antes de continuar preciso fazer um teste para verificar se a opera o de defini o do diret rio foi bem sucedida Ela n o ser a depender do diret rio que est sendo utilizado especificamente se o caminho para ele for muito longo setwdCUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa200314ufabcmetrialabusar setwd GfluxoUFABC2023EconometriaII2023Rlistaslab load lista3 dadosRdata Instalando e carregando packages com as rotinas que ser o utilizadas installpackages AER Importante voc muito provavelmente ainda n o instalou o pacote sampleSelection Por favor fa a isso acionando a linha abaixo installpackages sampleSelection Acionar pacotes library AER library sampleSelection MODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIA MODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIA MODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIAMODELOS VAR DEP BINARIA EXPLICATIVAS 1 Escolaridade 30 anoseduc 2 Experiencia 61 exper 3 Experiencia ao quadrado 62 expersq 4 Etnia Brancoa 43 dbca 5 Regiao brasileira 464850 dnodse dsu e dco 6 Numero de filhos 2629 filhosmdom filhosfdom filhosmndom filhosfndom 7 Renda domiciliar 14 rendadompc 8 Reside em area rural 59 drural MLP probit e logit para o emprego mlphfe lm formula dmicrohfedemp dmicrohfeanoseduc dmicrohfeexper dmicrohfeexpersq dmicrohfedbca dmicrohfedno dmicrohfedse dmicrohfedsu dmicrohfedco dmicrohfenfilhos dmicrohferendadompc dmicrohfedrural prohfe glm formula dmicrohfedemp dmicrohfeanoseduc dmicrohfeexper dmicrohfeexpersq dmicrohfedbca dmicrohfedno dmicrohfedse dmicrohfedsu dmicrohfedco dmicrohfenfilhos dmicrohferendadompc dmicrohfedruralfamily binomiallink probit OBS abaixo acrescente o modelo logit modificando a linha abaixo mudando apenas a link function para logit loghfe mlphferob coeftest mlphfe vcov vcovHC mlphfe prohfe coeftest prohfe loghfe coeftest loghfe setwdCUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa080414 Exportando resultados writetable mlphferob file mlphfetxt sep writetable prohfe file prohfetxt sep writetable loghfe file loghfetxt sep VAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUA VAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUA VAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUAVAR DEP DESCONTINUA Contagem por niveis salariais nrow dmicrohff dmicro hff 10 465 isna dmicrohff 10FALSE nrow dmicrohff dmicro hff 10 465 isna dmicrohff 10FALSE nrow dmicrohff dmicro hff 10 465 isna dmicrohff 10FALSE hist dmicro hff 10 plot density dmicrohff 10 narm T summary dmicro hff 10 nrow dmicrohff Estimando a equacao de Mincer com um modelo tobit descontinuo em y 465 tobit tobit dmicrohferendatrabpri10 dmicrohfeanoseduc dmicrohfeexper dmicrohfeexpersq dmicrohfedbca dmicrohfedno dmicrohfedse dmicrohfedsu dmicrohfedco dmicrohfedindu dmicrohfedapub dmicrohfedserv dmicrohfehorastrabprisref dmicrohfedformal left 465 OBS aqui voc deve acrescentar o modelo MQO tobit coeftest tobit Exportando resultados setwdCUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa080414 writetable tobit file tobithfftxt sep AMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADA AMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADA AMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADAAMOSTRA SELECIONADA Mincer MQO mqomincer Preencher acima o comando atentando para 1 A lista de vari veis incluindo a dependente a mesma do modelo tobit 2 Por m a base de dados outra sendo ela dmicrohfj Sum rio dos resultados de estima o summary mqomincer Corrigindo matriz de var covar mqomincerrob coeftest mqomincer vcov vcovHC mqomincer Visualizando mqomincerrob heckit Modelo de Heckman Rodando Aprendendo sobre o comando selection acionar linha abaixo selection Rodando hecmincer selection dmicrohfedemp dmicrohfeanoseduc dmicrohfeexper dmicrohfeexpersq dmicrohfedbca dmicrohfedno dmicrohfedse dmicrohfedsu dmicrohfedco dmicrohfenfilhos dmicrohferendadompc dmicrohfedrural dmicrohferendatrabpri10 dmicrohfeanoseduc dmicrohfeexper dmicrohfeexpersq dmicrohfedbca dmicrohfedno dmicrohfedse dmicrohfedsu dmicrohfedco dmicrohfedindu dmicrohfedapub dmicrohfedserv dmicrohfehorastrabprisref dmicrohfedformal method 2step Visualizando hecmincer coeftest hecmincer hecmincer Exportando setwdCUsersJesusDocumentsProjetosUFABCEconometriaIIFluxoUFABCcasa080414 writetable mqomincerrob file mqomincerrobtxt sep writetable hecmincer file hecmincertxt sep