1. Introducción una aproximación hacia el sentido de pertenencia en estudiantes chilenos

La preocupación por el sentido de pertenencia proviene de su importancia en cuanto ésta puede generar impactos positivos tanto como para estudiantes como para la comunidad educativa en su conjunto (Juvonen, 2006). De hecho, existen estudios que ponen el acento en cómo el sentido de pertenencia en la escuela puede ayudar a los jóvenes a integrarse de cierta forma al interior de la sociedad, sobre todo cuando se habla de estudiantes en contextos vulnerables (Riley, 2013). Frente a este respecto, debe entenderse el sentido de pertenencia como “el nivel de apego de un estudiante a su escuela (inversión personal en su escuela), el compromiso (cumplimiento de las reglas y expectativas de la escuela), la participación, (compromiso académicos de la escuela y actividades extracurriculares), y la convicción en su escuela (la fe en sus valores y su significado)” (Wehlage, Rutter, Smith, Lesko, y Fernández, 1989 en Kia-Keating M. y Ellis, 2007). En tanto, son varios los factores que investigaciones empíricas han indicado como importantes para comprender cómo se desarrolla en sí el sentido de pertenencia en la escuela, como lo es la relación intracomunidad (profesor-estudiante; estudiante-estudiante) (Hopenhayn et al. (2011), también factores como la motivación personal que tenga el estudiante a la hora de participar en la comunidad escolar. Aun cuando se tengan tales precendentes con respecto al sentido de pertenencia, asimismo se debe rescatar un rasgo cultural presente en los espacios de socialización secundaria que se manifiesta en un proceso de homosocialización por parte de la población masculina presente entre la población masculina que se relaciona al interior de los espacios escolares (Saavedra, 2016).A la luz de lo anterior, los espacios de interacción escolares resultarían ser ambientes favorables para cristalizar rasgos elementales en el proceso identitario de los jóvenes varones que se ven reflejados en ciertas disposiciones corporales y apropiaciones espaciales al interior de los establecimientos de educación. Ejemplo de esto son las apropiaciones espaciales de los patios por parte de varones, que, en los recreos tenderían a acaparar este lugar para practicar deportes.

Asumiendo la diferencia entre escuelas en el contexto nacional, resulta relevante establecer una distinción apropiada entre escuelas más allá de poder levantar un análisis que refiera solamente a un caso concreto sin tener presente su especificidad contextual, por ello se ha optado por pesquisar la base de datos de la prueba PISA que en una de sus mediciones especifica en un índice construido el sentido de pertenencia de los distintos establecimientos educacionales a lo que también se suma un índice construido por la misma institución internacional referida al nivel sociocultural y económico de las escuelas contempladas en la medición.

A partir de las anteriores afirmaciones que se tratará de responder la siguiente pregunta de investigación:

Pregunta 1 ¿Qué parte de la varianza de la variable BELONG es explicada por el contexto de las escuelas (tamaño de la escuela y promedio de Estatus socioeconómico y cultural)? H1: Se estima que las varianzas entre escuelas son altas, en tanto una escuela pequeña habría de tener mayor sentido de pertenencia y una escuela con un promedio de Estatus socioeconómico y cultural mayor también. Pregunta 2 ¿Cuáles son los efectos de las variables ST004D01T (género); MOTIVAT (motivación al éxito); COOPERATE (disposición a la cooperación) y unfairteacher (Parcialidad docente) sobre el sentido de pertenencia de los estudiantes? H1 (género): Debido a que las escuelas pueden presentarse como espacios donde se produce una apropiación de los espacios por parte de hombres, se cree que el sentido de pertenencia en mujeres sería menor, por lo que habría una relación entre variables. H2 (motivación al éxito): Producto de que teóricamente se pondera la motivación como uno de los factores fundamentales del sentido de pertenencia, se cree que existiría una relación significativa y positiva entre motivación y sentido de pertenencia. H3 (disposición a la cooperación): Debido a que el clima entre compañeros se considera fundamental para hablar de sentido de pertenencia en la teoría, se cree que las variables estarían relacionadas en un sentido positivo. H4: (parcialidad docente): Producto de que teóricamente la relación profesor-estudiante es vista como uno de los factores del sentido de pertenencia en la escuela, se cree que existiría una relación negativa entre las variables, es decir, mientras más parcial el profesor menos sentido de pertenencia tendría el estudiante. Pregunta 3 ¿Cuál es el efecto de la variación entre de las variables ST004D01T (género); MOTIVAT (motivación al éxito); COOPERATE (disposición a la cooperación) y unfairteacher (Parcialidad docente) y el sentido de pertenencia de los estudiantes entre contextos escolares (tamaño de la escuela y promedio de Estatus socioeconómico y cultural)? H5: Probablemente los efectos de estas variables para el modelo multinivel tengan que ver con efectos de interrelación entre las mismas variables. Es por esto que creemos que primero debemos profundizar en contenidos acerca de interacción entre niveles para elaborar una hipótesis más apropiada para esta pregunta de investigación.

2. Metodología 2.1 Descripción de la base Como ya se ha hecho mención, se ha utilizado la prueba PISA del año 2015 . La prueba en cuestión consiste en un trabajo colaborativo de los países miembros de la OCDE para medir el rendimiento de los estudiantes de 15 años que están ad portas de terminar el ciclo escolar obligatorio con el fin de observar si los y las estudiantes cuentan con las competencias para desenvolverse en los posibles desafíos dela vida cotidiana una vez terminado el periodo escolar. La encuesta PISA se realiza cada 3 años y, asimismo, para obtener una visión global, la metodología de la encuesta comprende la aplicación de cuestionarios hacia el mismo establecimiento, los profesores y los apoderados (OCDE, 2015).

El tamaño muestral para la aplicación del año 2015 cuenta con 6640 observaciones a nivel general. Por otra parte, las variables de nivel dos compuestas por el “promedios estatus socioeconómico y género por escuela” son de: (colocar)

#Librerias
pacman::p_load(ggplot2, foreign, lme4, lmtest, lattice, stargazer, nortest, ggplot2,dplyr,car,corrplot,reghelper,haven,texreg)

#abrir dataset
dataset = read.csv("CY6_MS_CMB_MERGE_F.csv")

#Selección de variables
dataset=dataset %>% select(CNTSCHID,BELONG,ST004D01T,MOTIVAT,COOPERATE,unfairteacher,ESCS,SCHSIZE) %>% 
  as.data.frame()

#converción de variables a numeric
dataset$BELONG<- as.numeric(dataset$BELONG)
dataset$MOTIVAT<-as.numeric(dataset$MOTIVAT)
dataset$SCHSIZE<- as.numeric(dataset$SCHSIZE)
dataset$COOPERATE<-as.numeric(dataset$COOPERATE)
dataset$unfairteacher<-as.numeric(dataset$unfairteacher)
dataset$ESCS<-as.numeric(dataset$ESCS)

dataset$unfairteacher[dataset$unfairteacher==99] <- NA
dataset$SCHSIZE[dataset$SCHSIZE==99999999] <- NA
dataset = na.omit(dataset)

2.2 Metodología Variable dependiente Sentido de pertenencia (BELONG): Dicha variable refiere a un índice creado al interior de la encuesta PISA, este indicador se compone de seis items ya preguntados en PISA 2012. El formato de respuesta corresponde a una escala likert de 4 puntos con las categorías de respuesta que van de “muy de acuerdo” a “muy en desacuerdo” (OCDE, 2015). Los Items ST034Q02TA,ST034Q03TA y ST034Q05TA fueron recodificados a la inversa, por lo que el número mayor indica un mayor grado de sentido de pertenencia.

Variables independientes Motivación al éxito(MOTIVAT) :dice construido a partir de una batería de preguntas tipo escala likert de 4(muy de acuerdo a muy en desacuerdo) categorías por parte de la división técnica de PISA 2015. Esta variable pretende dar cuenta de preferencias y pensamientos del o la estudiante.

Disposición a la cooperación (COOPERATE) : PISA 2015 incluyó preguntas referidas la cooperación y al trabajo en equipo para medir el dominio de resolución de problemas de forma colaborativa. Esta variable corresponde a un índice elaborado de un conjunto de preguntas tipo escala likert de 4 categorías “Muy de acuerdo” a “Muy en desacuerdo”.

Parcialidad docente (unfairteacher): PISA 2015 incluyó esta pregunta como una variable proxy para medir clima escolar. Dicha variable tiene como propósito observar comportamiento “injusto” por parte del profesorado. durante los últimos 12 meses. La variable se construyó a partir de una escala de 4 puntos cuyas categorías de respuesta son de “Nunca o casi nunca”, “algunas veces durante el año”, “algunas veces al mes”, “una vez a la semana o más”. Para efectos de nuestro análisis la variable fue recodificada para cambiar la direccionalidad del efecto. De esta manera, a mayor magnitud, más “justo” el o la docente.

Variables independientes de nivel 2 Tamaño de la escuela (SCHSIZE): variable recuperada desde el cuestionario aplicado a nivel de escuela que contiene el total de matriculados. La información fue suministrada por los directores de los establecimientos.

Promedio de Estatus socioeconómico y cultural: El índice de Estatus socioeconómico y cultural toma en consideración la educación de los padres, la ocupación de mayor estatus de alguno de los padres y la posesión de libros en el hogar, entre otras variables como ingreso, etc. Se considerará el promedio del estatus socioeconómico y cultural por escuela.

Promedio estudiantes por género: Para examinar el comportamiento de la variable dependiente teniendo en cuenta el género de las y los estudiantes, se decide calcular el promedio de la variable género por escuela.

#Recodificar variable unfairteacher para cambiar su direccionalidad 
dataset$fairteacher <- 24 - dataset$unfairteacher

#cálculo promedios estatus socioeconómico y género por escuela
dataset<-dataset %>%
  group_by(CNTSCHID) %>%
  mutate(m_gen=mean(ST004D01T))
summary(dataset$m_gen)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.000   1.423   1.500   1.503   1.579   2.000
dataset<-dataset %>%
  group_by(CNTSCHID) %>%
  mutate(m_escs=mean(ESCS))
# (dataset$m_escs)
#descriptivos por variable 
stargazer(dataset, title = "Descriptivos generales", type='text')
## 
## Descriptivos generales
## ===================================================
## Statistic N Mean St. Dev. Min Pctl(25) Pctl(75) Max
## ===================================================
#exploración de datos
hist(dataset$BELONG, xlim = c(0,1300)) #distribución variable dependiente a nivel estudiante

#histograma de BELONG por escuela
dataset <- dataset %>%
  group_by(CNTSCHID) %>%
  mutate(meanbelong = mean(BELONG))
hist(dataset$meanbelong)

#matriz de correlaciones 
cormat=dataset %>% select(BELONG,ST004D01T,MOTIVAT,COOPERATE,fairteacher,ESCS,SCHSIZE,m_gen,m_escs) %>% cor()
## Adding missing grouping variables: `CNTSCHID`
round(cormat, digits=2)
##             CNTSCHID BELONG ST004D01T MOTIVAT COOPERATE fairteacher  ESCS
## CNTSCHID        1.00  -0.02     -0.02    0.00      0.01        0.01 -0.01
## BELONG         -0.02   1.00     -0.03    0.02      0.03        0.04 -0.01
## ST004D01T      -0.02  -0.03      1.00    0.02      0.00       -0.10  0.00
## MOTIVAT         0.00   0.02      0.02    1.00      0.01        0.02 -0.02
## COOPERATE       0.01   0.03      0.00    0.01      1.00        0.04  0.00
## fairteacher     0.01   0.04     -0.10    0.02      0.04        1.00  0.00
## ESCS           -0.01  -0.01      0.00   -0.02      0.00        0.00  1.00
## SCHSIZE         0.13  -0.01      0.01   -0.01      0.00        0.02  0.00
## m_gen          -0.05  -0.02      0.41   -0.01      0.00       -0.04  0.00
## m_escs         -0.02  -0.01      0.00    0.01      0.01        0.03  0.22
##             SCHSIZE m_gen m_escs
## CNTSCHID       0.13 -0.05  -0.02
## BELONG        -0.01 -0.02  -0.01
## ST004D01T      0.01  0.41   0.00
## MOTIVAT       -0.01 -0.01   0.01
## COOPERATE      0.00  0.00   0.01
## fairteacher    0.02 -0.04   0.03
## ESCS           0.00  0.00   0.22
## SCHSIZE        1.00  0.03  -0.01
## m_gen          0.03  1.00   0.00
## m_escs        -0.01  0.00   1.00
corrplot.mixed(cormat)

Al examinar la distribución de la variable BELONG con un histograma nos damos cuenta que en los niveles más bajos de Sentido de Pertenencia hay una concentración de casos que va decayendo al avanzar el índice. Luego, al llegar a los valores superiores a 1100, existe una mayor cantidad de casos. La distribución termina por tener un sesgo hacia la derecha.

Al examinar las correlaciones entre las variables que se tendrán en cuenta, gran parte de estas son bastante pequeñas, ninguna pasando de 0,18. Esto, con excepción de las variables de estatus socioeconómico y cultural y género, con sus respectivos promedios (correlación de 0,69 y 0,41 respectivamente). Esto es esperable dado que, como se mencionó, son los promedios de las variables.

3. Resultados: Aplicación de modelos 3.1 Modelos de regresión de mínimos cuadrados de nivel uno

reg<- lm(BELONG~ST004D01T+MOTIVAT+COOPERATE+fairteacher+ESCS+SCHSIZE+m_gen+m_escs, data=dataset)
stargazer(reg, title = "Regresión datos individuales", type='text')
## 
## Regresión datos individuales
## ===============================================
##                         Dependent variable:    
##                     ---------------------------
##                               BELONG           
## -----------------------------------------------
## ST004D01T                    -17.220**         
##                               (8.622)          
##                                                
## MOTIVAT                        0.036           
##                               (0.028)          
##                                                
## COOPERATE                     0.150**          
##                               (0.060)          
##                                                
## fairteacher                  3.728***          
##                               (1.248)          
##                                                
## ESCS                          -0.001           
##                               (0.002)          
##                                                
## SCHSIZE                       -0.006           
##                               (0.006)          
##                                                
## m_gen                         -9.438           
##                              (20.937)          
##                                                
## m_escs                        -0.005           
##                               (0.010)          
##                                                
## Constant                    633.766***         
##                              (47.332)          
##                                                
## -----------------------------------------------
## Observations                   6,640           
## R2                             0.004           
## Adjusted R2                    0.003           
## Residual Std. Error     318.870 (df = 6631)    
## F Statistic           3.303*** (df = 8; 6631)  
## ===============================================
## Note:               *p<0.1; **p<0.05; ***p<0.01
#Evaluación de residuos
influenceIndexPlot(reg, vars=c("Cook", "Studentized", "hat"), id.n=5)

influencePlot(reg, id.n=3)