User Tools

Site Tools


neuroimagen:centiloid

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
neuroimagen:centiloid [2019/04/16 09:18]
osotolongo [Usando el cluster]
neuroimagen:centiloid [2019/04/23 09:45]
osotolongo [Procesando GAAIN]
Line 144: Line 144:
 </code> </code>
  
-Ahora, las imagenes PEt vieen en formato distinto al usado habitualmente (1 solo slice de 20 min) por lo que hay que retocar el script //fbbtemp_reg.sh// un poco para añadirle la posibilidad de proesar esto correctamente.+Ahora, las imagenes PET vienen en formato distinto al usado habitualmente (1 solo slice de 20 min) por lo que hay que retocar el script //fbbtemp_reg.sh// un poco para añadirle la posibilidad de proesar esto correctamente.
  
 <code bash fbbtemp_reg.sh> <code bash fbbtemp_reg.sh>
Line 275: Line 275:
 Tengo un ajuste con $r^2 = 0.97$ y una pendiente de $0.85 \pm 0.02$. La diferencia en los procedimientos esta en que en nuestra pipeline, siguiendo el metodo de ADNI, se toma como region de referencia la materia gris del cerebelo y en el caso de GAAIN se toma todo el cerebelo. No obstante la correspondencia entre las mediciones es buena. Tengo un ajuste con $r^2 = 0.97$ y una pendiente de $0.85 \pm 0.02$. La diferencia en los procedimientos esta en que en nuestra pipeline, siguiendo el metodo de ADNI, se toma como region de referencia la materia gris del cerebelo y en el caso de GAAIN se toma todo el cerebelo. No obstante la correspondencia entre las mediciones es buena.
  
- +Vamos a comprobarlo. Cambio la normalizacion para que normalice por el cerebelo completo y de paso incluyo todos los sujetos, 
 + 
 +<code bash> 
 +[osotolongo@detritus centiloid]$ awk -F ";" '{print $1,$3}' centiloid_fbb_fs_suvr_roisi_wcb.csv > calcs_wcb.dat 
 +[osotolongo@detritus centiloid]$ awk -F ";" '{print $1,$4}' FBB_suvr_centiloid.csv | sed 's/Y1/7/' > reference.dat 
 +[osotolongo@detritus centiloid]$ join calcs_wcb.dat reference.dat  
 +<code> 
 + 
 +Reviso el modelo en R, 
 + 
 +<code> 
 +> d <- read.csv("toreview.dat", sep = " ", header = TRUE) 
 +> m <- lm(d$Global ~ d$WC_suvr) 
 +> summary(m) 
 + 
 +Call: 
 +lm(formula = d$Global ~ d$WC_suvr) 
 + 
 +Residuals: 
 +      Min        1Q    Median        3Q       Max  
 +-0.073212 -0.028196  0.000229  0.026218  0.079672  
 + 
 +Coefficients: 
 +            Estimate Std. Error t value Pr(>|t|)     
 +(Intercept)  0.04626    0.02573   1.798   0.0817 .   
 +d$WC_suvr    0.98001    0.01894  51.731   <2e-16 *** 
 +--- 
 +Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
 + 
 +Residual standard error: 0.04017 on 32 degrees of freedom 
 +Multiple R-squared:  0.9882,    Adjusted R-squared:  0.9878  
 +F-statistic:  2676 on 1 and 32 DF,  p-value: < 2.2e-16 
 +</code> 
 + y miro el ajuste entre los datos en gnuplot, 
 +<code> 
 +gnuplot> f(x) = m*x +n 
 +gnuplot> fit f(x) "toreview.dat" u 2:3 via m,n 
 + 
 +After 4 iterations the fit converged. 
 +final sum of squares of residuals : 0.0531383 
 +rel. change during last iteration : -3.75201e-12 
 + 
 +degrees of freedom    (FIT_NDF)                        : 32 
 +rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0407501 
 +variance of residuals (reduced chisquare) = WSSR/ndf   : 0.00166057 
 + 
 +Final set of parameters            Asymptotic Standard Error 
 +=======================            ========================== 
 + 
 +m               = 1.00834          +/- 0.01949      (1.933%) 
 +n               = -0.0311848       +/- 0.02683      (86.03%) 
 + 
 + 
 +correlation matrix of the fit parameters: 
 + 
 +                    n       
 +m               1.000  
 +n              -0.965  1.000  
 +</code> 
 + 
 +Tengo una pendiente de $1.01 \pm 0.02$ y un ajuste con $R^2 = 0.9878$, lo cual esta muy bien. Este resultado valida el metodo de obtencion de SUVR de nuestra pipeline.  
 + 
 +{{:neuroimagen:wcbl_fit.png|}} 
 + 
 +**Nota:** El SUVR se debe calcular tomando como referencia la materia gris del cerebelo si los valores se van a comparar con datos de ADNI ya que esta es la metodologia que usan pero si se va a calcular centiloides debe tomarse como referencia el cerebelo completo.
 ==== Implementando pipeline Centiloid ===== ==== Implementando pipeline Centiloid =====
  
neuroimagen/centiloid.txt · Last modified: 2020/08/04 10:58 (external edit)