Hi,

I apply 3dICC on my datasets with the following model but I couldn't get it through with the ICC maps.

Output:deqi_icc_li4
Mask:deqi_mask+tlrc.BRIK
Clusters:4
Subj InputFile
acup002 acup_002_li4+tlrc
acup005 acup_005_li4+tlrc
acup007 acup_007_li4+tlrc
acup008 acup_008_li4+tlrc
acup010 acup_010_li4+tlrc
acup011 acup_011_li4+tlrc
acup012 acup_012_li4+tlrc
...

Unfortunately I couldn't understand what actually is the problem from the log file (below). Could you help?

Thank you.

Phoebe

****************
...
> print(sprintf("Start to run analysis on %i Z slices: %s", dimz, format(Sys.time(), "%D %H:%M:%OS3")))
[1] "Start to run analysis on 151 Z slices: 05/28/09 16:13:31.615"
> print("You can monitor the progress and estimate the total run time by opening this file from time to time.")
[1] "You can monitor the progress and estimate the total run time by opening this file from time to time."
>
> if (nFact ==2) {
+ nBrick <- 3;
+ outData <- array(0, dim=c(dimx, dimy, dimz, nBrick))
+
+ if (nNodes==1) for (kk in 1:dimz) {
+ outData[,,kk,] <-aperm(apply(IData[,,kk,], c(1,2), runAna2, Model=Model, ModelForm=ModelForm, dof=DF, tag=0), c(2,3,1))
+ cat("Z slice #", kk, "done: ", format(Sys.time(), "%D %H:%M:%OS3"), "\n")
+ }
+
+ if (nNodes>1) {
+ library(snow)
+ #cl <- makeCluster(rep('locahost', nNodes), type = "SOCK")
+ cl <- makeCluster(nNodes, type = "SOCK")
+ for (kk in 1:dimz) {
+ outData[,,kk,] <-aperm(parApply(cl, IData[,,kk,], c(1,2), runAna2, Model=Model, ModelForm=ModelForm, dof=DF, tag=0), c(2,3,1))
+ cat("Z slice #", kk, "done: ", format(Sys.time(), "%D %H:%M:%OS3"), "\n")
+ }
+ stopCluster(cl)
+ }
+ }
>
> if (nFact==3) {
+ nBrick <- 7
+ outData <- array(0, dim=c(dimx, dimy, dimz, nBrick))
+
+ if (nNodes==1) for (kk in 1:dimz) {
+ outData[,,kk,] <-aperm(apply(IData[,,kk,], c(1,2), runAna3, Model=Model, ModelForm=ModelForm, dof=DF, tag=0), c(2,3,1))
+ cat("Z slice #", kk, "done: ", format(Sys.time(), "%D %H:%M:%OS3"), "\n")
+ }
+
+ if (nNodes>1) {
+ library(snow)
+ #cl <- makeCluster(rep('locahost', nNodes), type = "SOCK")
+ cl <- makeCluster(nNodes, type = "SOCK")
+ for (kk in 1:dimz) {
+ outData[,,kk,] <-aperm(parApply(cl, IData[,,kk,], c(1,2), runAna3, Model=Model, ModelForm=ModelForm, dof=DF, tag=0), c(2,3,1))
+ cat("Z slice #", kk, "done: ", format(Sys.time(), "%D %H:%M:%OS3"), "\n")
+ }
+ stopCluster(cl)
+ }
+ }
>
> print(sprintf("Analysis finished: %s", format(Sys.time(), "%D %H:%M:%OS3")))
[1] "Analysis finished: 05/28/09 16:13:32.029"
> rm(IData) # retrieve some memory
>
> #outData <- band(outData, 0, 1)
> #dd<-aperm(dd, c(2,3,4,1))
>
> MyLabel <- paste(dimnames(fmAOV)[1][[1]])
>
> write.AFNI(OutFile, outData, MyLabel, note=Data$header$HISTORY_NOTE, origin=Data$origin,
+ delta=Data$delta, idcode="whatever")
Error in minmax(ttt) : object "outData" not found
Calls: write.AFNI -> minmax
Execution halted