* Encoding: UTF-8.

*Nacteni dat.
GET DATA 
  /TYPE=XLSX 
  /FILE='C:\Users\koritakova\Documents\Vyuka\DSAN01\2020 jaro\Odpovedniky_IS\Prakticke\neurodata_schizofrenici.xlsx' 
  /SHEET=name 'data' 
  /CELLRANGE=FULL 
  /READNAMES=ON 
  /DATATYPEMIN PERCENTAGE=95.0 
  /HIDDEN IGNORE=YES. 
EXECUTE. 



*** Otazka 1 ***.
* Zadání: K ověření vztahu tloušťky šedé hmoty levého posteriorního cingula a levého postcentrálního gyru 
u všech subjektů dohromady použijeme:.

* Bodový graf s odlišením skupin subjektů.
GRAPH
  /SCATTERPLOT(BIVAR)=lh_posteriorcingulate WITH lh_postcentral BY skupina
  /MISSING=LISTWISE.
* Závěr: v datech jsou 3 skupiny subjektů, které se neprolínají, korelační koeficient počítaný ze všech subjektů
současně nemá smysl počítat, protože by byl zkreslený.



*** Otazka 2 ***.
* Zadání: K otestování vztahu DSS skóre a Stroopova testu u kontrol použijeme:.

*Histogramy.
temp.
sel if skupina=1.
GRAPH
  /HISTOGRAM(NORMAL)=DSStest.

temp.
sel if skupina=1.
GRAPH
  /HISTOGRAM(NORMAL)=Strooptest.

* Bodový graf. 
temp.
sel if skupina=1.
GRAPH
  /SCATTERPLOT(BIVAR)=DSStest WITH Strooptest
  /MISSING=LISTWISE.
*Závěr: už podle histogramu vidíme, že u Stroopova testu je odlehlá hodnota, takže vhodnější bude neparametrický
Spearmanův korelační koeficient. Navíc p-hodnota Pearsonova chí-kvadrát testu by byla zkreslená, protože Stroopův
test nemá normální rozdělení.



*** Otazka 3 ***.
*Zadání: K ověření vztahu tloušťky šedé hmoty levého posteriorního cingula a levého postcentrálního gyru u pacientů 
s chronickou schizofrenií použijeme:.

temp.
sel if skupina=3.
GRAPH
  /HISTOGRAM(NORMAL)=lh_posteriorcingulate.

temp.
sel if skupina=3.
GRAPH
  /HISTOGRAM(NORMAL)=lh_postcentral.

* Bodový graf. 
temp.
sel if skupina=3.
GRAPH
  /SCATTERPLOT(BIVAR)=lh_posteriorcingulate WITH lh_postcentral
  /MISSING=LISTWISE.
*Závěr: V datech se nevyskytují odlehlé hodnoty a proměnné navíc nemají ani odchylky od normálního rozdělení -
můžeme použít Pearsonův korelační koeficient včetně výpočtu p-hodnoty.



*** Otazka 4 ***.
*Zadání: U pacientů s chronickou schizofrenií hodnoty DSS testu s rostoucím věkem:.

*Histogramy.
temp.
sel if skupina=3.
GRAPH
  /HISTOGRAM(NORMAL)=vek.

temp.
sel if skupina=3.
GRAPH
  /HISTOGRAM(NORMAL)=DSStest.

* Bodový graf. 
temp.
sel if skupina=3.
GRAPH
  /SCATTERPLOT(BIVAR)=vek WITH DSStest
  /MISSING=LISTWISE.
*Dle histogramu pro věk a bodového grafu jsou v datech odlehlé hodnoty -> Spearmanův korelační koeficient.

* Spearmanův korelační koeficient. 
temp.
sel if skupina=3.
NONPAR CORR
  /VARIABLES=vek DSStest
  /PRINT=SPEARMAN TWOTAIL NOSIG
  /MISSING=PAIRWISE.
*rho= -0,196; p=0,298
*Závěr: hodnoty DSS testu s rostoucím věkem klesají (korelační koeficient je záporný); vztah je statisticky nevýznamný (p=0,298).



*** Otazka 5 ***. 
*Zadání:Je vztah Stroopova testu a tloušťky šedé hmoty levého entorhinálního kortexu u kontrol statisticky významný?.

*Histogramy.
temp.
sel if skupina=1.
GRAPH
  /HISTOGRAM(NORMAL)=Strooptest.

temp.
sel if skupina=1.
GRAPH
  /HISTOGRAM(NORMAL)=lh_entorhinal.

* Bodový graf. 
temp.
sel if skupina=1.
GRAPH
  /SCATTERPLOT(BIVAR)=Strooptest WITH lh_entorhinal
  /MISSING=LISTWISE.
*Dle bodového grafu je v datech minimálně jedna odlehlá hodnota a navíc hodnoty Stroopova testu nemají normální rozdělení
-> Spearmanův korelační koeficient.

temp.
sel if skupina=1.
NONPAR CORR
  /VARIABLES=Strooptest lh_entorhinal
  /PRINT=SPEARMAN TWOTAIL NOSIG
  /MISSING=PAIRWISE.
*r=0,364; p=0,048.
*Závěr: vztah je statisticky významný.



*** Otazka 6 ***. 
*Zadání: Je statisticky významný rozdíl v hodnotách Pearsonova korelačního koeficientu tloušťky šedé hmoty posteriorního  
cingula a levého entorhinálního kortexu u pacientů s první epizodou schizofrenie a pacientů s chronickou schizofrenií?.

temp.
sel if skupina=2.
CORRELATIONS
  /VARIABLES=lh_posteriorcingulate lh_entorhinal
  /PRINT=TWOTAIL NOSIG
  /MISSING=PAIRWISE.
*r=0,437; p=0,016; n=30.

temp.
sel if skupina=3.
CORRELATIONS
  /VARIABLES=lh_posteriorcingulate lh_entorhinal
  /PRINT=TWOTAIL NOSIG
  /MISSING=PAIRWISE.
*r=0,032; p=868; n=30.

*Použijeme webový kalkulátor (http://vassarstats.net/rdiff.html) a zjistíme, že oboustranná p=0,1096 (zaokrouhleně 0,110).
*Závěr: Rozdíl mezi korelačními koeficienty není statisticky významný (tzn. korelační koeficienty se statisticky významně
neliší).



*** Otazka 7 ***.
*Zadání: Pokud by koeficient korelace tloušťky šedé hmoty dané mozkové struktury s věkem u pacientů s chronickou schizofrenií byl
roven -0,48, je tento korelační koeficient statisticky významně nižší než hodnota -0,20, která byla zjištěna při populačním průzkumu?.

*Použijeme webový kalkulátor (http://vassarstats.net/rpop.html) a zjistíme, že oboustranná p=0,097. My se však ptáme, zda je NIŽŠÍ
a skutečně je NIŽŠÍ, p-hodnota tedy bude 0,097/2 = 0,048 (tzn. je to jednostranná p-hodnota vypočítaná kalkulátorem).
*Poznámka: Kdyby nám to ale vyšlo opačně, než jsme čekali (tzn. korelační koeficient by byl VYŠŠÍ), museli bychom počítat 1-p/2.
*Závěr: Korelační koeficient je statisticky významně nižší než hodnota -0,20 zjištěná při populačním průzkumu.



*** Otazka 8 ***.
*Zadání: Pokud chceme srovnat korelační koeficient DSS skóre s věkem a korelační koeficient Stroopova testu s věkem u kontrolních 
subjektů, použijeme:.

*Závěr: nelze srovnat korelační koeficienty pocházející ze stejného vzorku.



*** Otazka 9 ***.
*Zadání: Kolik procent variability DSS skóre je predikováno pomocí věku, pokud r=-0,356?.

*Ručně vypočítáme r^2 = (-0,356)^2 = 0,127 -> 12,7%.



*** Otazka 10 ***. 
*Zadání: Po odstranění vlivu věku z proměnné DSS skóre, je průměrná hodnota nestandardizovaných reziduí u kontrolních subjektů rovna:.

REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN 
  /DEPENDENT DSStest
  /METHOD=ENTER vek
  /SAVE PRED RESID.

temp.
sel if skupina=1.
DESCRIPTIVES VARIABLES=RES_1
  /STATISTICS=MEAN STDDEV MIN MAX.
*Závěr: průměrná hodnota je rovna 4,34.