Saltar al contenido principal

# 7.5: Fondo: Interpolación de canales defectuosos

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$

( \newcommand{\kernel}{\mathrm{null}\,}\) $$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\id}{\mathrm{id}}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\kernel}{\mathrm{null}\,}$$

$$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$

$$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$

$$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

$$\newcommand{\vectorA}[1]{\vec{#1}} % arrow$$

$$\newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow$$

$$\newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vectorC}[1]{\textbf{#1}}$$

$$\newcommand{\vectorD}[1]{\overrightarrow{#1}}$$

$$\newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}}$$

$$\newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}}$$

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\avec}{\mathbf a}$$ $$\newcommand{\bvec}{\mathbf b}$$ $$\newcommand{\cvec}{\mathbf c}$$ $$\newcommand{\dvec}{\mathbf d}$$ $$\newcommand{\dtil}{\widetilde{\mathbf d}}$$ $$\newcommand{\evec}{\mathbf e}$$ $$\newcommand{\fvec}{\mathbf f}$$ $$\newcommand{\nvec}{\mathbf n}$$ $$\newcommand{\pvec}{\mathbf p}$$ $$\newcommand{\qvec}{\mathbf q}$$ $$\newcommand{\svec}{\mathbf s}$$ $$\newcommand{\tvec}{\mathbf t}$$ $$\newcommand{\uvec}{\mathbf u}$$ $$\newcommand{\vvec}{\mathbf v}$$ $$\newcommand{\wvec}{\mathbf w}$$ $$\newcommand{\xvec}{\mathbf x}$$ $$\newcommand{\yvec}{\mathbf y}$$ $$\newcommand{\zvec}{\mathbf z}$$ $$\newcommand{\rvec}{\mathbf r}$$ $$\newcommand{\mvec}{\mathbf m}$$ $$\newcommand{\zerovec}{\mathbf 0}$$ $$\newcommand{\onevec}{\mathbf 1}$$ $$\newcommand{\real}{\mathbb R}$$ $$\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}$$ $$\newcommand{\laspan}[1]{\text{Span}\{#1\}}$$ $$\newcommand{\bcal}{\cal B}$$ $$\newcommand{\ccal}{\cal C}$$ $$\newcommand{\scal}{\cal S}$$ $$\newcommand{\wcal}{\cal W}$$ $$\newcommand{\ecal}{\cal E}$$ $$\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}$$ $$\newcommand{\gray}[1]{\color{gray}{#1}}$$ $$\newcommand{\lgray}[1]{\color{lightgray}{#1}}$$ $$\newcommand{\rank}{\operatorname{rank}}$$ $$\newcommand{\row}{\text{Row}}$$ $$\newcommand{\col}{\text{Col}}$$ $$\renewcommand{\row}{\text{Row}}$$ $$\newcommand{\nul}{\text{Nul}}$$ $$\newcommand{\var}{\text{Var}}$$ $$\newcommand{\corr}{\text{corr}}$$ $$\newcommand{\len}[1]{\left|#1\right|}$$ $$\newcommand{\bbar}{\overline{\bvec}}$$ $$\newcommand{\bhat}{\widehat{\bvec}}$$ $$\newcommand{\bperp}{\bvec^\perp}$$ $$\newcommand{\xhat}{\widehat{\xvec}}$$ $$\newcommand{\vhat}{\widehat{\vvec}}$$ $$\newcommand{\uhat}{\widehat{\uvec}}$$ $$\newcommand{\what}{\widehat{\wvec}}$$ $$\newcommand{\Sighat}{\widehat{\Sigma}}$$ $$\newcommand{\lt}{<}$$ $$\newcommand{\gt}{>}$$ $$\newcommand{\amp}{&}$$ $$\definecolor{fillinmathshade}{gray}{0.9}$$

En la demostración de video, el canal C5 parecía bastante problemático, y mencioné que podríamos querer interpolar ese canal. Además, el canal F8 tenía mucho ruido de alta frecuencia, por lo que también deberíamos considerar interpolar ese canal. En este ejercicio, discutiremos cómo decidir si un canal debe ser interpolado y luego veremos el proceso de interpolación real.

La interpolación se realiza por separado en cada punto de tiempo, utilizando los valores de voltaje de los otros canales en ese punto de tiempo. Existen varios algoritmos que se pueden utilizar para estimar un valor razonable para un canal sobre la base de los otros canales, y todos funcionan razonablemente bien. Me gusta usar el algoritmo esférico proporcionado por EEGLAB. Trata los electrodos como si estuvieran ubicados en una cabeza esférica, y ajusta una función polinómica a la distribución de voltajes, dejando fuera el canal a ser interpolado. Entonces, el valor de la función polinómica en la ubicación a ser interpolada se usa como el voltaje estimado en esa ubicación. Este proceso se repite de forma independiente en cada momento de la muestra. Obtienes una forma de onda bastante razonable en la ubicación interpolada, pero ten en cuenta que es solo una estimación imperfecta de la forma de onda verdadera.

La decisión sobre si interpolar un canal dado en última instancia se reduce a nuestro objetivo fundamental de responder con precisión a la pregunta científica que el experimento fue diseñado para abordar. Voy a proporcionar alguna orientación general, pero en última instancia hay que pensar si la interpolación sirve a ese objetivo. ¿Te acerca a la verdad o más lejos de la verdad? Imagina que has enviado un manuscrito a una revista, y uno de los revisores visita tu laboratorio para ver cómo procesaste realmente tus datos (lo que nunca sucedería en la realidad, por supuesto). ¿Estarías feliz de explicarle a este revisor cómo decidiste si interpolar? ¿O te sentirías un poco avergonzado?

Empecemos con un caso extremo. Imagínese que se rompió un electrodo y la señal de ese electrodo era puro ruido. E imagina que el electrodo roto no estaba siendo utilizado en ninguno de los análisis clave. Por ejemplo, medimos la amplitud MMN en FCz para los análisis principales en el artículo ERP CORE (Kappenman et al., 2021), por lo que los electrodos C5 y F8 no jugaron un papel importante en nuestros análisis (aunque sí hicieron una contribución menor cuando trazamos mapas de cuero cabelludo). Una forma de lidiar con el electrodo roto sería descartar completamente los datos de este participante. Esto reduciría el tamaño de la muestra en nuestros análisis, lo que tiende a disminuir nuestra capacidad de sacar conclusiones precisas sobre la población, por lo que no es una gran opción. La otra opción principal sería interpolar los datos del electrodo roto. Dado que el electrodo roto no contribuye a los análisis principales, poder incluir al participante interpolando y con ello aumentando el tamaño de tu muestra parece que sirve a la verdad mucho más que excluir al participante de todos los análisis.

Ahora consideremos un caso extremo en sentido contrario. Imagina que estás analizando los datos recopilados en otro laboratorio con métodos de grabación realmente pobres (o un participante realmente desafiante), y los datos de FCz y los 8 electrodos más cercanos se ven terribles. No sólo es malo el canal clave para el análisis, sino que también lo son los canales circundantes, lo que dificultará la interpolación precisa. En este caso, incluir al participante en los análisis finales parece que no agregará información real y hará más daño que bien.