| mni.write.vertex.stats {mni.cortical.statistics} | R Documentation |
Takes the output from mni.vertex.statistics and writes it to a vertstats file using either the new (with headers) or old version of the file format.
mni.write.vertex.stats(vertex.stats, "filename.vertstats", headers=TRUE, mean.stats=ms, glim.matrix=gm)
vertex.stats |
The results from mni.vertex.statistics |
filename |
The filename to write to. |
headers |
Whether to write out header information. If TRUE, then a new style file is created with a hierarchical header; this is the default. If false, then a simple file without headers is created. |
mean.stats |
The results from mni.mean.statistics run on the same dataset. This information is written to the header if headers=TRUE. mean.stats is NULL by default. |
glim.matrix |
The glim matrix used with mni.read.glim.file. If passed to this argument then the matrix is added to the header. Default is NULL |
mni.mean.statistics.
mni.vertex.statistics.
mni.write.vertex.stats.
# A complete statistical run:
# read the design matrix file
gf <- mni.read.glim.file("filename.glim", header=TRUE)
# see if the mean is different
ms <- mni.mean.statistics(gf, 'y ~ 1 + V2')
# look at the results of the mean stats run
summary(ms)
# run a test at each vertex
vs <- mni.vertex.statistics(gf, 'y ~ 1 + V2')
# write the results out to file.
mni.write.vertex.stats(vs, "filename.vertstats", headers=TRUE,
mean.stats=ms, glim.matrix=gf)