FONIS_Paula.R

rm(list=ls())
estudio <-read.csv("http://www.udec.cl/~jbustosm/consolidado_bases.csv", sep=";", h=T)
estudio0<-estudio[,c(1:51)]
attach(estudio0)
estudio0$cansancio_emocional<-(MBI1.1+MBI2.1+MBI3.1+MBI6.1+MBI8.1+MBI13.1+MBI14.1+
MBI16.1+MBI20.1)
detach(estudio0)
##-------------------------HISTOGRAMA DE FRECUENCIA PARA BURNOUT: CANSANCION EMOCIONAL--------------------------------
hist(estudio0$cansancio_emocional, main="Cansancio Emocional", ylab="Frecuencia", xlab="Puntaje", col = gray.colors(11),
breaks=11)
abline(v=19, col="red")
abline(v=26, col="red")
X11()


##-------------------------HISTOGRAMA DE FRECUENCIA PARA BURNOUT: DESPERSONALIZACIÓN--------------------------------

attach(estudio0)
estudio0$despersonalización<-(MBI5.1+MBI10.1+MBI11.1+MBI15.1+MBI22.1)
detach(estudio0)
hist(estudio0$despersonalización, main="Despersonalización", ylab="Frecuencia", xlab="Puntaje", col = gray.colors(11),
breaks=11)
abline(v=6, col="red")
abline(v=9, col="red")

##-------------------------HISTOGRAMA DE FRECUENCIA PARA BURNOUT: REALIZACION PERSONAL--------------------------------
X11()
attach(estudio0)
estudio0$realizacionpersonal<-(MBI4.1+MBI7.1+MBI9.1+MBI12.1+MBI17.1+MBI18.1+MBI19.1+MBI21.1)
detach(estudio0)
hist(estudio0$realizacionpersonal, main="Realizacion personal", ylab="Frecuencia", xlab="Puntaje", col = gray.colors(11),
breaks=11)
abline(v=34, col="red")
abline(v=39, col="red")

X11()
# Added to the  plot:
par(mar=c(3.1, 3.1, 1.1, 2.1))
hist(estudio0$realizacionpersonal,xlim=c(10,50), col = gray.colors(8), main="Realizacion personal", ylab="Frecuencia", xlab="Puntaje")
boxplot(estudio0$realizacionpersonal, horizontal=TRUE,  outline=TRUE,  ylim=c(10,50), frame=F, col = "white", add = TRUE,width=1)
abline(v=34, col="red")
abline(v=39, col="red")


X11()
##-----------------LAS DIMENSNIONE SEGUN SEXO, EDAD, UNIDAD Y TITULO
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(cansancio_emocional~ SEXO, data = estudio0, col = "lightgray", names=c("Mujeres","Hombres"), ylim=c(0,50),ylab="Puntaje Cansancio Emocional")
     boxplot(despersonalización~ SEXO, data = estudio0, col = "lightgray", names=c("Mujeres","Hombres"), ylim=c(0,50), ylab="Puntaje Despersonalización")
    boxplot(realizacionpersonal~ SEXO, data = estudio0, col = "lightgray", names=c("Mujeres","Hombres"), ylim=c(0,50), ylab="Punataje Realización Personal")



#---------------------------UCI VS ONCO
attach(estudio0)
estudio0$UNIDAD2[UNIDAD=="1"] <- "ONCO"
estudio0$UNIDAD2[UNIDAD=="2"] <- "ONCO"
estudio0$UNIDAD2[UNIDAD=="5"] <- "ONCO"
estudio0$UNIDAD2[UNIDAD=="3"] <- "UCI"
estudio0$UNIDAD2[UNIDAD=="4"] <- "UCI"
detach(estudio0)

X11()
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(cansancio_emocional~ UNIDAD2, data = estudio0, col = "lightgray", ylim=c(0,50),ylab="Puntaje Cansancio Emocional")
     boxplot(despersonalización~ UNIDAD2, data = estudio0, col = "lightgray", ylim=c(0,50), ylab="Puntaje Despersonalización")
    boxplot(realizacionpersonal~ UNIDAD2, data = estudio0, col = "lightgray", ylim=c(0,50),ylab="Punataje Realización Personal")

#----------------------------TITULO PROFESIONAL-------------------------------------
attach(estudio0)
estudio0$TITULO2[TITULO=="1"] <- "TÉCNICOS"
estudio0$TITULO2[TITULO=="2"] <- "TÉCNICOS"
estudio0$TITULO2[TITULO=="3"] <- "PROFESIONALES"
estudio0$TITULO2[TITULO=="4"] <- "PROFESIONALES"
estudio0$TITULO2[TITULO=="5"] <- "PROFESIONALES"
estudio0$TITULO2[TITULO=="6"] <- "PROFESIONALES"
estudio0$TITULO2[TITULO=="7"] <- "PROFESIONALES"
estudio0$TITULO2[TITULO=="8"] <- "PROFESIONALES"
detach(estudio0)
X11()
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(cansancio_emocional~ TITULO2, data = estudio0, col = "lightgray", ylim=c(0,50),ylab="Puntaje Cansancio Emocional")
     boxplot(despersonalización~ TITULO2, data = estudio0, col = "lightgray", ylim=c(0,50), ylab="Puntaje Despersonalización")
    boxplot(realizacionpersonal~ TITULO2, data = estudio0, col = "lightgray", ylim=c(0,50),ylab="Punataje Realización Personal")


#----------------------------EDAD-------------------------------------

X11()
plot(estudio0$EDAD, estudio0$cansancio_emocional, xlab="Edad (años)", ylab="Cansancio Emocional")
X11()
plot(estudio0$EDAD, estudio0$despersonalización,  xlab="Edad (años)", ylab="Despersonalización")
X11()
plot(estudio0$EDAD, estudio0$realizacionpersonal, xlab="Edad (años)", ylab="Realizacion Personal")

cor.test(estudio0$despersonalización, estudio0$cansancio_emocional)
cor.test(estudio0$cansancio_emocional, estudio0$realizacionpersonal)
cor.test(estudio0$despersonalización, estudio0$realizacionpersonal)

#----------------------------APOYO EN DUELO-------------------------------------
attach(estudio0)
estudio0$adTotal<-(ED1.1+ED2.1+ ED3.1+ED4.1+ED5.1+ED6.1 +ED7.1+ED8.1+ED9.1+ED10.1+ED11.1+ED12.1+ED13.1+ED14.1+ED15.1)
detach(estudio0)

attach(estudio0)
estudio0$adTotal<-(ED1.1+ED2.1+ ED3.1+ED4.1+ED5.1+ED6.1 +ED7.1+ED8.1+ED9.1+ED10.1+ED11.1+ED12.1+ED13.1+ED14.1+ED15.1)
detach(estudio0)

attach(estudio0)
estudio0$reconocimientodelvinculo<-(ED1.1+ED2.1+ ED3.1+ED4.1+ED5.1)
estudio0$reconocimientodelaperdida<-(ED6.1 +ED7.1+ED8.1+ED9.1+ED10.1)
estudio0$incorritualesfunerarios<-(ED11.1+ED12.1+ED13.1+ED14.1+ED15.1)
detach(estudio0)

X11()
hist(estudio0$adTotal)
X11()
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(reconocimientodelvinculo~ UNIDAD2, data = estudio0, col = "lightgray", ylab="Reconocimiento del vinculo")
t.test(reconocimientodelvinculo~ UNIDAD2, data = estudio0, col = "lightgray")
   
boxplot(reconocimientodelaperdida~ UNIDAD2, data = estudio0, col = "lightgray",ylab="Reconocimiento de la perdida")
   t.test(reconocimientodelaperdida~ UNIDAD2, data = estudio0)

 boxplot(incorritualesfunerarios~ UNIDAD2, data = estudio0, col = "lightgray", ylab="Rituales funerarios")
   t.test(incorritualesfunerarios~ UNIDAD2, data = estudio0, col = "lightgray")



X11()
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(reconocimientodelvinculo~ SEXO, data = estudio0, col = "lightgray",names=c("Mujer","Hombre"), ylab="Puntaje Reconocimiento del vinculo")
t.test(reconocimientodelvinculo~ SEXO, data = estudio0, col = "lightgray")
   
boxplot(reconocimientodelaperdida~ SEXO, data = estudio0, col = "lightgray", names=c("Mujer","Hombre"), ylab="Puntaje Reconocimiento de la perdida")
   t.test(reconocimientodelaperdida~ SEXO, data = estudio0, col = "lightgray")
  

 boxplot(incorritualesfunerarios~ SEXO, data = estudio0, col = "lightgray",names=c("Mujer","Hombre"),  ylab="Puntaje Incorporacion a funerales")
   t.test(incorritualesfunerarios~ SEXO, data = estudio0, col = "lightgray")

X11()
par(mfrow = c(1, 3))
 ## boxplot on a formula:
     boxplot(reconocimientodelvinculo~ TITULO2, data = estudio0, col = "lightgray", ylab="Puntaje Reconocimiento del vinculo")
t.test(reconocimientodelvinculo~ TITULO2, data = estudio0, col = "lightgray")
   
boxplot(reconocimientodelaperdida~ TITULO2, data = estudio0, col = "lightgray", ylab="Puntaje Reconocimiento de la perdida")
   t.test(reconocimientodelaperdida~ TITULO2, data = estudio0, col = "lightgray")
  

 boxplot(incorritualesfunerarios~ TITULO2, data = estudio0, col = "lightgray",  ylab="Puntaje Incorporacion a funerales")
   t.test(incorritualesfunerarios~ TITULO2, data = estudio0, col = "lightgray")



X11()
plot(estudio0$EDAD, estudio0$reconocimientodelvinculo, xlab="Edad (años)", ylab="Puntaje Reconocimiento del vinculo")
X11()
plot(estudio0$EDAD, estudio0$reconocimientodelaperdida, xlab="Edad (años)", ylab="Puntaje Reconocimiento de la perdida")
X11()
plot(estudio0$EDAD, estudio0$incorritualesfunerarios, xlab="Edad (años)", ylab="Puntaje Rituales Funerarios")

X11()
 boxplot(reconocimientodelvinculo~ T.Prof, data = estudio0,ylab="Puntaje Reconocimiento del vinculo", main="Tiempo en la Profesión",names=c("0 a 1 año",  
"1 a 5 años","6 a 10 años","11 a 15 años","15 años o más"),pos=3 )
X11()
 boxplot(reconocimientodelaperdida~ T.Prof, data = estudio0,ylab="Puntaje Reconocimiento de la perdida", main="Tiempo en la Profesión",names=c("0 a 1 año",  
"1 a 5 años","6 a 10 años","11 a 15 años","15 años o más"),pos=3 )

 boxplot(incorritualesfunerarios~ T.Prof, data = estudio0,ylab="Puntaje Rituales Funerarios", main="Tiempo en la Profesión",names=c("0 a 1 año",  
"1 a 5 años","6 a 10 años","11 a 15 años","15 años o más"),pos=3 )


X11()
par(mfrow = c(1, 3))
 boxplot(reconocimientodelvinculo~ CURSOS, data = estudio0, ylab="Puntaje Reconocimiento del vinculo", ylab="Cursos", names=c("si","no"))
 boxplot(reconocimientodelaperdida~ CURSOS, data = estudio0, ylab="Puntaje Reconocimiento de la perdida",ylab="Cursos", names=c("si","no"))
 boxplot(incorritualesfunerarios~ CURSOS, data = estudio0, ylab="Puntaje Rituales Funerarios",ylab="Cursos", names=c("si","no"))
title("Asistencia a Cursos de Duelo")

#-----------Segunda parte: reunión del viernes 5 de febrero en clinica puc
#comparación de las escalas (sentir, perdida y apoyo)del duelo
#con variables sociodemograficas en UCI y ONCO

attach(estudio0)
estudio0$DueloPuntajeSentir<-ED1.1+ED2.1+ED3.1+ED4.1+ED5.1
estudio0$DueloPuntajePerdida<-ED6.1+ED7.1+ED8.1+ED9.1+ED10.1
estudio0$DueloPuntajeApoyo<-ED11.1+ED12.1+ED13.1+ED14.1+ED15.1
detach(estudio0)

X11()
par(mfrow = c(1, 3))
 boxplot(DueloPuntajeSentir~ UNIDAD2, data = estudio0, ylab="Puntaje Escala Duelo (Sentir)")
 boxplot(DueloPuntajePerdida~ UNIDAD2, data = estudio0, ylab="Puntaje Escala Duelo (Perdida)")
 boxplot(DueloPuntajeApoyo~ UNIDAD2, data = estudio0, ylab="Puntaje Escala Duelo (Apoyo)")





#B de burnout vs sentir, perdida y apoyo.
Tabla de frecuencia de 

Sentir	Leve	Moderado	Alto
Perdida
Apoyo

En t1, t2,..., tn