22
Le logiciel R
+ fichier,sep=""),paste(getwd(),"/",fichier,sep=""),mode="wb")
> recupere("Mond4D.nii")
> recupere("Mondanat.hdr")
> recupere("Mondanat.img")
La paire de fichiers Mondanat.img/Mondanat.hdr contient une image anatomique du cerveau du sujet. Puis, on installe le package permettant la lecture
de ces donn´ ees.
> install.packages("AnalyzeFMRI") # Choisissez un miroir.
> # Noms des fichiers.
> file.func <- paste(getwd(),"/data/","Mond4D.nii",sep="")
> file.anat <- paste(getwd(),"/data/","Mondanat.img",sep="")
> slice <- 10 # Numéro de la coupe cérébrale.
Les instructions ci-dessous permettent la lecture des donn´ ees (une tranche,
slice en anglais, du cerveau en fait).
> anat.slice <- f.read.nifti.slice(file.anat,slice,1)
> class(anat.slice)
[1] "matrix"
> dim(anat.slice)
[1] 128 128
> func.slice <- f.read.nifti.slice.at.all.timepoints(
+
file.func,slice)
> class(func.slice)
[1] "array"
> dim(func.slice)
[1] 128 128 125
Nous cr´ eons `
a pr´ esent le codage du signal du stimulus visuel (1 = couleur,
0 = absence de couleur).
> stimulus <- c(rep(c(1,1,1,1,1,0,0,0,0,0,0,0,0,0,0),8),1,1,1,1,1)
> plot(stimulus,type="l")
0
2 0
4 0
6 0
8 0
1 0 0
1 2 0
0.0
0.2
0.4
0.6
0.8
1.0
Index
stimulus