$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Diseño del método LBA considerando las perspectivas MPFR y emocionales
Se construye un sistema de evaluación multiperspectiva para analizar el Comportamiento de Aprendizaje en Línea (SOLB) de los estudiantes, y se diseña un modelo MPFR. Para analizar más a fondo los sentimientos subjetivos de los estudiantes durante el aprendizaje, este estudio adopta el algoritmo XGBoost y diseña el algoritmo IGWO para optimizar hiperparámetros y mejorar la precisión de la clasificación.
Construcción del modelo MPFR para LBA
Para analizar SOLB, este estudio parte principalmente de dos perspectivas para formar un plan de análisis completo. En primer lugar, en términos de rendimiento de aprendizaje, este estudio considera tres perspectivas: currículo, individual y clase, para formar el modelo MPFR. En segundo lugar, en términos de análisis de sentimiento de los comentarios de los estudiantes, este estudio diseña un algoritmo XGBoost mejorado. A través del análisis del rendimiento en el aprendizaje y las emociones de retroalimentación de los estudiantes, este estudio puede analizar mejor la SOLB. En cuanto a detalles técnicos específicos del análisis del rendimiento del aprendizaje, este estudio utiliza primero datos de una plataforma de aprendizaje online y los preprocesa. En segundo lugar, este estudio selecciona indicadores clave de evaluación del rendimiento del aprendizaje desde diferentes perspectivas para construir un sistema de evaluación integral. Posteriormente, se construye un modelo MPFR correspondiente para evaluar el comportamiento de aprendizaje.
Este estudio selecciona datos de la plataforma Superstar (fuente: https://www.chaoxing.com/). Contiene información de múltiples dimensiones, como actividades en el aula y evaluación de calificaciones, que puede proporcionar una buena base para futuras asignaturas de LBA. Tras obtener los datos, debe ser preprocesado, incluyendo principalmente la eliminación de campos irrelevantes, el filtrado de datos y las comprobaciones de integridad de los datos. Por ejemplo, en los datos del profesor, es necesario eliminar información irrelevante como el número y el departamento al que pertenece el curso. Posteriormente, este estudio construye la ingeniería de características de LBA basada en los datos de comportamiento de aprendizaje de la plataforma Superstar. El índice LBA de estudiantes proporciona una base de ingeniería de características para el modelo MPFR posterior, y el contenido de este índice se muestra en la Figura 1.

Figura 1: Indicadores para LBA estudiantil. La figura aclara las tres dimensiones principales (currículo, individual, clase) del LBA del estudiante y los indicadores específicos de evaluación bajo cada dimensión, proporcionando soporte de características para el modelo MPFR. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
En la Figura 1, los indicadores del LBA para estudiantes giran principalmente en torno a tres perspectivas: curso, individual y clase, con un claro enfoque en cada dimensión, proporcionando un soporte preciso de características para la construcción de modelos MPFR. Entre ellas, en cuanto a indicadores de dimensión del curso, este estudio selecciona la tasa de finalización de tareas del curso, la duración del aprendizaje y la frecuencia de interacción con el curso. La preocupación central de esta dimensión es la calidad general de finalización de las tareas del curso y la puntualidad de la participación en el aprendizaje. Por ejemplo, la tasa de finalización de las tareas del curso se utiliza para seguir la diferencia en las tasas de finalización antes y después del plazo; La duración del estudio del curso se utiliza para distinguir diferentes intervalos de intensidad de aprendizaje; La frecuencia de interacción se utiliza para filtrar comportamientos interactivos efectivos. A nivel personal, el progreso en el aprendizaje personal, la calidad de la finalización de deberes y las puntuaciones de los exámenes se seleccionan como indicadores de evaluación. Esta dimensión se centra en el grado de coincidencia entre el progreso del aprendizaje personal y los resultados en el dominio del conocimiento. Ejemplos incluyen comparar el progreso personal del aprendizaje con los planes de curso, evaluar la precisión y eficiencia de los trabajos completados y clasificar los niveles de dominio del conocimiento basándose en las puntuaciones de los exámenes. Además, en cuanto a la dimensión de la clase, los indicadores seleccionados incluyen la puntuación media de la clase, la actividad de interacción en clase y el ambiente de aprendizaje en clase. Esta dimensión se centra principalmente en el impacto del entorno de aprendizaje en general en el comportamiento individual. Por ejemplo, la nota media se utiliza para localizar el nivel relativo de rendimiento individual en el grupo; la actividad de interacción de clase se utiliza para reflejar el grado de participación colectiva; El ambiente de aprendizaje se utiliza para analizar el efecto impulsor del entorno grupal en el comportamiento de aprendizaje. Para analizar el comportamiento de aprendizaje de los estudiantes, se construye un modelo de razonamiento difuso y se obtienen las características del comportamiento de aprendizaje desde diferentes perspectivas. El razonamiento difuso es un método de razonamiento que utiliza lógica difusa, que procesa información imprecisa o incierta mediante conjuntos difusos y reglas difusas, y cuenta con las ventajas de una gran flexibilidad y fácil comprensión11,12. El razonamiento difuso, como mecanismo central de computación de los modelos MPFR, se aplica después de la ingeniería de características. El proceso principal de razonamiento difuso se muestra en la Figura 2.

Figura 2: El diagrama principal del razonamiento difuso. La figura muestra el proceso central del razonamiento difuso, incluyendo cuatro pasos clave: fuzzificación, construcción de la base de reglas difusas, razonamiento difuso y desenfoque, y aclara la lógica del modelo MPFR en el manejo de datos de comportamiento de aprendizaje inciertos. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
El razonamiento difuso implica principalmente fuzzificación, establecimiento de una Biblioteca de Reglas Difusas (FRL), razonamiento difuso y desenfoque. Entre ellas, la fuzzificación requiere convertir datos de entrada precisos en conjuntos difusos e implica definir funciones de pertenencia. El FRL contiene una serie de reglas difusas y se utiliza principalmente para describir la relación entre entrada y salida. En la construcción de la base de reglas difusas, el estudio invita a tres profesores asociados de tecnología educativa con ≥ de 8 años de experiencia docente para desarrollar conjuntamente 28 reglas, y determina las reglas finales mediante tres rondas de iteración utilizando el "método Delphi". Por ejemplo, Regla 1: SI la tasa de finalización de la tarea del curso es inferior al 30% y el progreso individual está 2 semanas retrasado respecto al calendario → ENTONCES el estado de alto riesgo. El razonamiento difuso es el proceso de inferir datos de entrada difusos construidos sobre un FRL y obtener resultados de salida difusos. Finalmente, la eliminación de difuminado convertirá el resultado de salida borrosa en un valor de salida preciso. El estudio selecciona una función de pertenencia triangular, que es una función ampliamente utilizada en sistemas lógicos difusos y que presenta ventajas como gran flexibilidad y alta eficiencia computacional. La expresión de la función μA(x) se muestra en la ecuación (1).
(1)
En la ecuación (1), x es el valor de entrada específico. A, B y C son los tres vértices de un triángulo. En el desenfoque, este artículo utiliza el método del centroide para convertir los resultados de salida. El método del centroide es una técnica de desenfoque comúnmente utilizada en lógica difusa, que presenta ventajas como una fuerte intuición, cálculo sencillo y estabilidad13. La expresión del método del centroide se muestra en la ecuación (2)14.
(2)
En la ecuación (2), f* es el valor de salida tras el desenfoque, que es el centroide del conjunto difuso. μ(x) es la función de membresía.
es la suma ponderada de todos los puntos en un conjunto difuso.
es la integral de la función de pertenencia sobre todos los valores posibles de x. La racionalidad de usar reglas y funciones de pertenencia definidas manualmente en el modelo MPFR se refleja en tres aspectos. En primer lugar, garantía de cualificación experta: los responsables de la normativa son tres profesores asociados de tecnología educativa y han sido certificados por el "Analista de Comportamiento de Aprendizaje en Línea" del "Centro de Investigación en Educación Online del Ministerio de Educación" para garantizar que las normas cumplan con las leyes de educación y enseñanza. En segundo lugar, ventaja de eficiencia: El tiempo de inferencia de las reglas manuales es mucho menor que el de los métodos basados en datos, lo que lo hace más adecuado para las necesidades de "inferencia en tiempo real" de las plataformas educativas online y no requiere una gran cantidad de datos anotados, reduciendo así los costes de recogidade datos 15. En tercer lugar, adaptándose al requisito fundamental de "interpretabilidad" en escenarios de educación online, las reglas definidas manualmente y las funciones de pertenencia triangular (Ecuación 1) pueden corresponder directamente a la cognición intuitiva en escenarios de enseñanza. Los educadores no necesitan una formación técnica compleja para entender la razón por la que un estudiante se considera en riesgo. Métodos basados en datos, como los sistemas neuro-fuzzy, pueden optimizar automáticamente las reglas. Sin embargo, la mayoría de las reglas generadas son expresiones matemáticas complejas, que resultan difíciles de interpretar para los educadores, lo que reduce la aceptabilidad del modelo en la intervención docente real.
Diseño de un modelo mejorado de análisis de sentimiento de comentarios estudiantiles para XGBoost
El modelo MPFR diseñado puede analizar SOLB desde tres aspectos: curso, individual y clase. Sin embargo, el análisis externo de datos conductuales tiene limitaciones a la hora de expresar los sentimientos subjetivos de los estudiantes. Por ello, para analizar más a fondo los sentimientos subjetivos de los estudiantes sobre el aprendizaje, se recopilan y preprocesan datos adicionales de las plataformas de aprendizaje MOOC que contienen información de comentarios de los estudiantes. Posteriormente, el XGBoost se utiliza para construir el modelo de clasificación de sentimiento. El IGWO está diseñado para optimizar sus parámetros y mejorar la precisión del modelo de clasificación. La mejora de XGBoost se refleja en su optimización de parámetros mediante el algoritmo IGWO. El preprocesamiento de datos es un primer paso crucial antes de la construcción y el análisis del modelo. El proceso de preprocesamiento de datos se muestra en la Figura 3.

Figura 3: El proceso de preprocesamiento de datos. El proceso de preprocesamiento de datos de nueve pasos incluye limpieza de texto, segmentación de palabras, eliminación de palabras de parada, extracción de tallos y reconocimiento de palabras de sentimiento, proporcionando datos de alta calidad para el análisis de comportamiento y la clasificación de sentimientos posteriores. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
Los pasos de preprocesamiento de datos son los siguientes: El primer paso es realizar una limpieza de texto para eliminar campos irrelevantes en los datos de la plataforma Chaoxing y filtrar muestras inválidas. Al mismo tiempo, los campos no textuales y las reseñas breves se eliminan de los datos de reseñas en la plataforma MOOC. El segundo paso es realizar el procesamiento de segmentación. Entre ellos, los comentarios chinos utilizan el "modo preciso" de Jieba para el procesamiento de segmentación de texto, mientras que los comentarios en inglés emplean la función de segmentación de palabras de la biblioteca NLTK para su procesamiento. El tercer paso es eliminar palabras clave comunes como ''de'', ''yes'', ''in'', etc. Entre ellas, las palabras clave chinas utilizan la lista de palabras clave (https://github.com/goto456/stopwords/blob/master/hit_stopwords.txt) del Instituto de Tecnología de Harbin. Las palabras de paso se eliminan mediante coincidencias de palabras. Las palabras de oclusiva en inglés utilizan la lista universal de palabras de parada integrada en la biblioteca NLTK, que también se elimina mediante coincidencia palabra por palabra. El cuarto paso es utilizar el stemmer de la biblioteca NLTK para extraer stems del texto y restaurar los verbos en los datos ingleses a su forma stem. El quinto paso es identificar las palabras emocionales positivas y negativas en los comentarios basándose en el diccionario emocional de CNKI para proporcionar una base previa a la anotación para la posterior clasificación emocional. El sexto paso es realizar la extracción de características. Entre ellos, los datos estructurados en la plataforma Chaoxing están estandarizados, mientras que los indicadores de clasificación se codifican por separado. Al mismo tiempo, estos datos de comentarios de la plataforma MOOC utilizan un modelo de lenguaje preentrenado (Bidirectional Encoder Representation from Transformers, BERT) para extraer características de texto. La principal ventaja del modelo BERT en la extracción de características de texto es que puede generar dinámicamente vectores de palabra conscientes del contexto mediante una arquitectura Transformer bidireccional profunda, abordando efectivamente la ambigüedad que los modelos estáticos tradicionales de incrustación de palabras no pueden manejar. El séptimo paso es abordar el problema de las categorías de sentimiento desequilibradas en los datos de comentarios. Se emplea la Técnica de Sobremuestreo de Minorías Sintéticas (SMOTE) para procesarlo. Se seleccionan cinco muestras de vecinos más cercanos para interpolación y así generar muestras sintéticas de la clase minoritaria, con una semilla aleatoria fija de 42 para asegurar la reproducibilidad del experimento. El octavo paso es anotar y particionar los datos. El noveno paso es rellenar los valores que faltan en los datos de superestrellas y eliminar los valores atípicos. Al realizar el procesamiento de balanceo de datos, SMOTE crea nuevas muestras compuestas interpolando entre muestras de clases minoritarias, aumentando así el número de muestras de clases minoritarias y haciendo que la distribución por clases del conjunto de datos sea másequilibrada 16,17. La expresión de SMOTE se muestra en la ecuación (3).
Bmn = dm + rand(0,1) x (dmn -d m) (3)
En la ecuación (3), Bmn es una muestra de clase minoritaria construida artificialmente, dm es la i-ésima muestra de clase minoritaria, y dmn es la n-ésima muestra entre los k vecinos más cercanos de dm. rand(0,1) es un número aleatorio entre 0 y 1. XGBoost mejora el rendimiento predictivo añadiendo iterativamente árboles de predicción. Tiene ventajas como alta precisión, gran flexibilidad y facilidad de uso18,19. Los pasos para construir iterativamente un modelo de árbol usando XGBoost incluyen: inicializar el modelo, construir iterativamente el árbol, la función objetivo y el término de regularización. La salida
predicha en la t-ésima iteración se muestra en la ecuación (4).
(4)
En la ecuación (4),
es el valor de predicción del modelo tras la t - 1ª iteración, ft(h) es la función de predicción del modelo en árbol añadida en la t-ésima iteración, y h es el vector de características de entrada. La función objetivo para minimizar XGBoost se muestra en la ecuación (5).
(5)
En la ecuación (5), i es el número de muestra. I y J son el número total de muestras y árboles, y j es el número de árboles.
es la función de pérdida, y gi es la verdadera etiqueta de la i -ésima muestra.
es el valor predicho del modelo para la i -ésima muestra tras la t-ésima iteración. Ω(f t) es el término de regularización, y (ft) es la función de predicción del j -ésimo árbol. La expresión general Ω(f) para la regularización se muestra en la ecuación (6).
(6)
En la ecuación (6), γ significa el coeficiente de penalización para la cantidad de nodos hoja, λ es el coeficiente de regularización L2 para el peso de los nodos hoja, y w es el peso de los nodos hoja. Sin embargo, la elección de los hiperparámetros XGBoost tiene un efecto directo y significativo en el rendimiento. Los hiperparámetros comunes de XGBoost incluyen la tasa de aprendizaje, la profundidad máxima del árbol y los parámetros de regularización. Debido al potencial gran espacio de hiperparámetros, ajustar manualmente estos parámetros no solo consume mucho tiempo, sino que también es difícil encontrar la combinación óptima de parámetros. Por ello, este estudio diseña IGWO para optimizar los hiperparámetros de XGBoost. GWO busca soluciones óptimas simulando la jerarquía social y las estrategias de caza de los lobos grises, con ventajas como menos parámetros y una fuerteadaptabilidad 20,21. En GWO, se genera la posición inicial de la población como se muestra en la ecuación (7).
(7)
En la ecuación (7), Puv es la posición del u-ésimo individuo en la v-ésima dimensión.
y
son los límites superior e inferior del v-ésimo espacio de búsqueda. rand() es un número aleatorio entre [0,1]. Sin embargo, el algoritmo GWO puede encontrarse con problemas como la lenta velocidad de convergencia y quedarse atascado en óptimos locales en las etapas intermedias y posteriores. Esto también significa que, en este momento, GWO puede no ser capaz de explorar eficazmente todo el espacio de soluciones, lo que dificulta encontrar la solución óptima global. Para abordar este problema, este estudio mejora GWO desde dos aspectos: la introducción del Mapa Caótico de la Tienda (TCM) y el Factor de Convergencia No Lineal (NCF). La MTC es un mapeo caótico simple que garantiza que las regiones óptimas de solución potenciales no se pierdan durante el proceso de búsqueda y evita búsquedas repetidas en la misma región22. Por tanto, mediante la MTC, se puede generar una población inicial más uniforme y aleatoria, y se puede mejorar la capacidad del algoritmo para superar óptimos locales. El cálculo de la MTC se muestra en la ecuación (8).
(8)
En la ecuación (8), y es el parámetro de control. Rt y Rt+1 son las variables de estado del sistema en las t-ésima y t+1-ésima iteraciones. La fórmula para la NCF z se muestra en la ecuación (9).
(9)
En la ecuación (9), zmax es el factor máximo de convergencia y tmax es el número máximo de iteraciones. Este método no lineal de disminución permite a GWO mantener un gran tamaño de paso para la búsqueda global en la etapa inicial. En la etapa intermedia de la búsqueda, la velocidad de convergencia se acelera, mientras que en la etapa posterior la búsqueda local se realiza con pasos más pequeños, equilibrando eficazmente las capacidades de búsqueda global y local y mejorando el rendimiento de optimización. El algoritmo IGWO se utiliza específicamente para la etapa de optimización de hiperparámetros del clasificador de sentimiento XGBoost, y su proceso principal se muestra en la Figura 4.

Figura 4: El proceso principal de IGWO. La figura describe el proceso de ejecución de IGWO, incluyendo la inicialización de poblaciones basada en mapeo caótico de tendas, actualización de posición de factores de convergencia no lineales y cribado individual óptimo, y aclara la lógica para optimizar los hiperparámetros XGBoost. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
En la Figura 4, el proceso IGWO incluye: inicializar el algoritmo; el uso de la MTC para inicializar la población; calcular el valor inicial de aptitud de cada individuo; usando NCF para actualizar las posiciones de los lobos grises; seleccionar la solución con la mejor forma física como nuevo individuo óptimo; y la salida de la persona óptima. La combinación óptima de hiperparámetros XGBoost puede enviarse a través de IGWO. El proceso general de clasificación de sentimiento integra el preprocesamiento de datos, la optimización IGWO y el modelo final XGBoost. El proceso del modelo de clasificación basado en IGWO-XGBoost se muestra en la Figura 5.

Figura 5: El proceso del modelo de clasificación basado en IGWO-XGBoost. La figura muestra el proceso completo del modelo de clasificación de sentimiento IGWO-XGBoost, desde la entrada y preprocesamiento de datos hasta la partición de conjuntos de datos, la optimización de hiperparámetros de IGWO, la construcción del modelo XGBoost y la salida de resultados de clasificación, presentando claramente la cadena de operaciones del modelo. Por favor, haz clic aquí para ver una versión ampliada de esta figura.
El modelo basado en IGWO-XGBoost incluye procesos como datos de entrada, preprocesamiento de datos, partición de conjuntos de datos, optimización de hiperparámetros del algoritmo IGWO, combinación óptima de hiperparámetros, construcción de XGBoost basada en la combinación óptima de hiperparámetros y salida de resultados de clasificación. A través de este modelo de clasificación, las emociones en los comentarios de los estudiantes pueden clasificarse, lo que permite una comprensión más intuitiva de los sentimientos subjetivos de los estudiantes sobre el aprendizaje. Las versiones específicas del software, semillas aleatorias, especificaciones de hardware, especificaciones ambientales y fuentes de conjuntos de datos utilizadas en el estudio se muestran en la Tabla 1.
| Tipo de decorado | Modelo y contenido específicos |
| Software básico | Python 3.9.7 |
| Biblioteca central | XGBoost 1.7.5、Scikit-learn 1.2.2、Jieba 0.42.1、NLTK 3.8.1、Pandas 1.5.3、NumPy 1.24.3、Matplotlib 3.7.1、Imbalanced-learn 0.10.1 |
| Semilla aleatoria | Establece la semilla aleatoria global en 42 |
| Especificaciones de hardware/entorno | 1. Sistema operativo: Windows 10 Edición Profesional (64 bits, número de versión 22H2) |
| 2. Procesador: Intel Core i5-12600KF (con frecuencia principal de 3,7 GHz y frecuencia de aceleración dinámica de 4,9 GHz) |
| 3. Memoria: 64GB DDR4 (frecuencia 3200MHz, soporta hasta 128GB de memoria) |
| 4. Almacenamiento: SSD de 1TB (Samsung 980 Pro, velocidad de lectura 7450MB/s) |
| 5. Tarjeta gráfica: NVIDIA GeForce RTX 3060 (12GB de memoria de vídeo) |
| Fuentes de conjuntos de datos y detalles de acceso | 1. Conjunto de datos de plataformas superstar: |
| -Localizador uniforme de recursos de la fuente (URL): https://www.chaoxing.com/ |
| -Método de adquisición: solicitar a través de la plataforma abierta de big data de Chaoxing Education (ruta de aplicación: plataforma web oficial → centro de desarrollo → interfaz de datos → conjunto de datos de comportamiento de aprendizaje) |
| 2. Conjunto de datos de comentarios de la plataforma MOOC: |
| -URL de la fuente: https://www.icourse163.org/.cn |
| -Método de adquisición: Solicitar a través del canal de "soporte de investigación de datos" de la plataforma MOOC |
| Scripts o modelos personalizados | 1. Script de preprocesamiento de datos |
| 2. Script modelo MPFR |
3. Script modelo IGWO-XGBoost
|
Tabla 1: Versiones específicas de software, semillas aleatorias, especificaciones de hardware, especificaciones ambientales y fuentes de conjuntos de datos utilizadas en la investigación. Proporciona una lista detallada del entorno de software y hardware requerido, configuraciones semillada aleatorias, fuentes de conjuntos de datos y el rango de búsqueda de hiperparámetros XGBoost para el estudio para garantizar la reproducibilidad y estandarización del experimento.
La Tabla 1 muestra que el estudio adopta XGBoost 1.7.5 como marco central del modelo de clasificación de sentimiento e introduce Scikit learn 1.2.2 para preprocesamiento de datos, extracción de características y evaluación de modelos. Además, el estudio establece uniformemente la semilla aleatoria en 42 para asegurar la reproducibilidad de la partición de datos, el sobremuestreo SMOTE, la optimización de hiperparámetros IGWO y los resultados del entrenamiento del modelo IGWO-XGBoost. Además, el algoritmo IGWO establece el espacio de búsqueda para los hiperparámetros XGBoost. El rango de búsqueda para la tasa de aprendizaje es [0,001, 0,3] en la escala logarítmica, y la profundidad máxima del árbol se busca en el conjunto entero {3, 4,..., 15}. El rango de optimización del término de regularización L2 es [0,1, 5,0], y la proporción de muestreo de cada subconjunto de características del árbol está optimizada dentro del rango de [0,6, 1,0]. El rango de sintonización para el peso mínimo de los nodos hoja es [1, 10].
El conjunto de datos y el código de preprocesamiento y análisis de datos relacionados fueron generados y analizados por el propio equipo. Los scripts básicos para el preprocesamiento de datos se muestran en la Tabla 2.
| Pandas importados como PD, Jieba, Re, Numpy, como NP |
| de la word_tokenize de importación de nltk.tokenize; de la importación nltk.stem PorterStemmer |
| Definitivamente limpio (texto, l): |
| return re.sub(r"[^\u4e00-\u9fa5]" si l=="zh" si no r"[^a-zA-Z]", " ", text).strip() |
| Def Process (texto, l): |
| t = jieba.lcut(texto) si l=="zh" si no word_tokenize(texto) |
| return " ".join([PorterStemmer().stem(w) si l=="en" si no w por w en t si w no está en open("hit_stopwords.txt").read().split()]) |
| def main(c,m,l): |
| cx=pd.read_csv(c).query("duración del curso>=1 & resultados del examen.notna()"); mc=pd.read_csv(m).query("text de comentario.str.len()>=5") |
| mc["t"]=mc["texto de comentario"].apply(clean,args=("zh",).apply(proceso,args=("zh",)) |
| np.savez("out.npz",**{k:v for k,v in locals().items()}) |
Tabla 2: Script central para el preprocesamiento de datos. Presentar la información central del script para el preprocesamiento de datos, aclarar los pasos de preprocesamiento correspondientes al script y proporcionar referencias técnicas para la replicación de métodos de investigación.