Hi,
I am currently using a 3dICC.R and something just isn't working right. We are trying to look at the reliability of resting state networks, by looking at the networks of subjects that return for a second FMRI session about 6 weeks after the first session.
My model.txt file is:
Output:test
MASK:
Clusters:1
Subj Session InputFile
s78257 one ../../TBI_RST/ICA_RESULTS/ICA_GB10_20comp/selected_components/NC_RSNs_sub001_comp18+tlrc
s78257 two ../../TBI_RST_visit2/ICA_RESULTS/20_comps/selected_components/NC_RSNs_sub001_comp18+tlrc
s38477 one ../../TBI_RST/ICA_RESULTS/ICA_GB10_20comp/selected_components/NC_RSNs_sub002_comp18+tlrc
s38477 two ../../TBI_RST_visit2/ICA_RESULTS/20_comps/selected_components/NC_RSNs_sub002_comp18+tlrc
s80929 one ../../TBI_RST/ICA_RESULTS/ICA_GB10_20comp/selected_components/NC_RSNs_sub003_comp18+tlrc
s80929 two ../../TBI_RST_visit2/ICA_RESULTS/20_comps/selected_components/NC_RSNs_sub003_comp18+tlrc
...
However, as an output, I only get test+orig.HEAD, without the BRIK file. The output is also in +orig space when my input files are in +tlrc.
From the log file, the error appears to be in this section:
...
Called from: FUN(newX[, i], ...)
Browse[1]> print(sprintf("Analysis finished: %s", format(Sys.time(), "%D %H:%M:%OS3")))
[1] "Analysis finished: 07/07/09 13:45:12.067"
Browse[1]> rm(IData) # retrieve some memory
Warning message:
In rm(IData) : variable "IData" was not found
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]> #outData <- band(outData, 0, 1)
Browse[1]> #dd<-aperm(dd, c(2,3,4,1))
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]> MyLabel <- paste(dimnames(fmAOV)[1][[1]])
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]> write.AFNI(OutFile, outData, MyLabel, note=Data$header$HISTORY_NOTE, origin=Data$origin,
+ delta=Data$delta, idcode="whatever")
Error in paste(label, collapse = "~") : object "MyLabel" not found
Calls: aperm ... write.AFNI -> writeChar -> AFNIheaderpart -> paste -> paste
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]> statpar <- "3drefit"
Browse[1]> statpar <- paste(statpar, " -newid -view ", outView, OutFile)
Browse[1]> system(statpar)
++ 3drefit: AFNI version=AFNI_2008_07_18_1710 (Oct 28 2008) [64-bit]
++ Authored by: RW Cox
** Can't read ANY data from file 'test+orig.HEAD'
** ERROR: Can't open dataset test+orig
++ 3drefit processed 0 datasets
Browse[1]> print(sprintf("Congratulations! You've got output %s+%s.*", Out, outView))
[1] "Congratulations! You've got output test+tlrc.*"
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]>
Called from: FUN(newX[, i], ...)
Browse[1]> # set save defaults using option:
Browse[1]> #options(save.defaults=list(ascii=TRUE, safe=FALSE))
Browse[1]> #save.image()
Browse[1]> #unlink(".RData")
Browse[1]> Called from: FUN(newX[, i], ...)
Browse[1]> Called from: FUN(newX[, i], ...)
...
Any help will be greatly appreciated.
Thank you very much,
Alex