
Walmes, segue o CMR de um exemplo que consta no livro Analise de Modelos de Regressao Linear do Charnet. Pag 152. Depois, ja na pagina 243, é refeito com variaveis dummy, usando "0" para mulher e "1"para homem. #####Exemplo Variaveis Dummy ####### #Os dados referem-se a a salarios e tempo de experiencia no cargo #de gerente de agencias bancarias de um grande banco. #Alem disso, gostaria de saber se há diferenças quando #o sexo é levado em conta. y<-c(1.93,3.17,2.27,3.13,2.77,3.09,2.65,2.22,2.85,3.23,2.82,1.90,2.53,2.56,4.22,4.09,3.60,4.70,3.14,2.99,4.74,4.11,2.36,4.09,4.50,2.90,4.48) x<-c(0,17,5,15,9,15,8,5,13,20,11,1,6,7,23,20,18,27,11,10,29,23,4,22,25,9,25) reg1<-lm(y~x) ## Depois fizemos incorporando as dummies, do seguinte modo: z<-c(rep(0,12),rep(1,15)) #até 12 feminino - 15 restante masculino xz<-x*z #incluindo a interação xz reg<-lm(y~x+z+xz) summary(reg) ### Um outro modo: regm<-lm(y[z=="0"]~x[z=="0"]) #modelo para mulheres summary(regm) regh<-lm(y[z=="1"]~x[z=="1"]) #modelo para homens plot(x[z=="1"],regh$residuals) summary(regh) E se eu quiser entao, associar um segundo fator de variaveis dummy, alem de sexo, poderia colocar "pós graduação", 0 se nao possui, 1 se possui.

É só você criar os zeros e uns a um objeto, por exemplo posgrad <- c(0,1,1,0,1,0,1, ...) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Julia Maria
-
Walmes Zeviani