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", headers=TRUE, mean.stats=NULL, glim.matrix=NULL)
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 ## Not run: gf <- mni.read.glim.file("filename.glim", header=TRUE) # see if the mean is different ## Not run: ms <- mni.mean.statistics(gf, 'y ~ 1 + V2') # look at the results of the mean stats run ## Not run: summary(ms) # run a test at each vertex ## Not run: vs <- mni.vertex.statistics(gf, 'y ~ 1 + V2') # write the results out to file. ## Not run: mni.write.vertex.stats(vs, "filename.vertstats", headers=TRUE, mean.stats=ms, glim.matrix=gf) ## End(Not run)