Don't like ads? PRO users don't see any ads ;-)
Guest

Untitled

By: a guest on Apr 3rd, 2012  |  syntax: R  |  size: 1.37 KB  |  hits: 38  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
This paste has a previous version, view the difference. Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. varssociov<-c("c1", "c2", "n.c3", "c4","n.c5","n.c6");varssociov
  2. summary(cdn[ ,varssocio ])
  3.  
  4. # Revisar bloc socioeconomic
  5. attach(cdn)
  6. summary(c3)
  7. table(c3)
  8.  
  9. sel<-which(c3<1);sel
  10. is.na(cdn$c3)<-sel
  11. cdn$c3[ sel ]<-NA
  12. detach(cdn)
  13. cdn$n.c3<-factor(cdn$c3,levels=1:2,labels=c("Si","No"))
  14. summary(cdn$n.c3)
  15. table(cdn$n.c3)
  16. summary(cdn[,varssocio])
  17.  
  18. table( cdn$c5 )
  19. nivells<-c("Empresario+ de 6 empleados","Empresario - de 6 empleados","Agricultor sin empleados",
  20.   "Profesión liberal/ autónomos","Director y cuadros superiores","Técnicos y especialistas en la Industria",
  21.   "Técnicos y especialistas en Servicios en oficina","Técnicos y especialistas en Servicios fuera de oficina",
  22.   "Funcionarios","Jornaleros y Subalternos","Parados","Jubilados y pensionistas","Otros")
  23. # Els no sabe i no contesta són NA per mi               
  24. cdn$n.c5<-factor(cdn$c5,levels=1:13,labels=nivells)
  25.  
  26. table( cdn$n.c5 )
  27. summary( cdn$n.c5 )
  28.  
  29.  
  30. nivells <- c("Alta","Media alta","Media media","Media baja","Baja")
  31. cdn$n.c6<-factor(cdn$c6,levels=1:5,labels=nivells)
  32. table( cdn$n.c6 )
  33. summary( cdn$n.c6 )
  34.  
  35.  
  36. sel<-which(cdn$c5<1);sel
  37. is.na(cdn$c5)<-sel
  38. cdn$c5[ sel ]<-NA
  39.  
  40. sel<-which(cdn$c6<1);sel
  41. is.na(cdn$c6)<-sel
  42. cdn$c6[ sel ]<-NA
  43.  
  44.  
  45.  
  46.  
  47. prueba<-factor(cdn$c2,levels=6,)
  48. is.na(cdn$c2)<-prueba
  49. cdn$c2[ prueba ]<-NA
  50. prueba<-factor(cdn$c4,levels=6:7,)
  51. is.na(cdn$c4)<-prueba
  52. cdn$c4[ prueba ]<-NA