#~New session: Version: S Tue Jul 16 10:52:22 PDT 1996, Time: Tue Aug 25 21:57:39 1998
#~get "/usr/local/Splus3.4/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/splus/.Functions/.laenv" -1 ...
#~
source("First")
#~error: can't open source file "First"
#~put ".Data/last.dump" 904096607 ...
#~
source(".First")
#~put ".Data/.First" 904096621 ...
#~put ".Data/.Last.value" 904096621 ...
#~
()
#~error: statement with syntax error
#~
q()
#~End session: Time: 904096749; Process: 10976
#~New session: Version: S Tue Jul 16 10:52:22 PDT 1996, Time: Tue Aug 25 21:59:13 1998
#~get "/usr/local/Splus3.4/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/splus/.Functions/.laenv" -1 ...
/home/chelsea2/bruno/Swave_new/Crc_Irrec.S#~error: statement with syntax error
#~put ".Data/last.dump" 904096667 ...
#~
First
#~error: Object "First" not found
#~
.First
#~
dyn.load.shared("./pack.so")
#~put ".Data/.Last.value" 904096813 ...
#~
source("Crazy_Climbers.S")
/home/chelsea2/bruno/Swave_new/Crc_Irrec.S#~error: statement with syntax error
#~put ".Data/last.dump" 904096819 ...
#~
source("Crazy_Climbers.S")
#########################################################################
#      $Log: Crc_Rec.S,v $
# Revision 1.2  1995/04/05  18:56:55  bruno
# *** empty log message ***
#
# Revision 1.1  1995/04/02  01:04:16  bruno
# Initial revision
#
#
#               (c) Copyright  1997                             
#                          by                                   
#      Author: Rene Carmona, Bruno Torresani, Wen-Liang Hwang   
#                  Princeton University
#                  All right reserved                           
#########################################################################


#########################################################################
#
#	Functions to reconstruct a signal from the output of
#	 the crazy climber algorithm.
#
#########################################################################




crcirrec <- function(siginput, inputwt, beemap, noct, nvoice, compr,
	minnbnodes = 2, w0 = 2*pi, bstep = 5,ptile = .01, prob = 0.8,
	epsilon = .5, fast = F, para = 0, real = F, plot=1)
#########################################################################
#     crcirrec:
#     -------
#      Reconstruction of a real valued signal from ridges found by 
#      crazy climbers on a wavelet transform (sampled irregularly).
#
#      input:
#      ------
#      siginput: input signal
#      inputwt: continuous wavelet transform (output of cwt)
#      beemap: output of crazy climber algorithm
#      noct: number of octaves (powers of 2)
#      nvoice: number of different scales per octave
#      compr: subsampling rate for the ridge
#      bstep: used for the chaining
#      ptile: 
#      epsilon: coeff of the Q2 part of reconstruction kernel
#      fast: if set to TRUE, computes the Q2 kernel by Riemann sums
#            if not, uses a Romberg adaptive step quadrature.
#      para:
#      plot: plot the original and reconstructed signal in display
#      prob: the percentile to keep 
#
#      output:
#      -------
#      rec: reconstructed signal
#      ordered: image of the ridges (with different colors)
#      comp: 2D array containing the signals reconstructed from ridges
#
#########################################################################
{

   tmp <- cfamily(beemap,bstep,ptile=ptile);   image(tmp$ordered);   chain <- tmp$chain;   nbchain <- tmp$nbchain;   ordered <- tmp$ordered;   sigsize <- length(siginput);   rec <- numeric(sigsize);   plnb <- 0;
   par(mfrow=c(1,1));   tsplot(siginput);   title("Original signal");
   tmp <- matrix(0,nbchain,length(siginput));
   totnbnodes <- 0;   idx <- numeric(nbchain)	;   p <- 0;
   for (j in 1:nbchain){
      phi.x.min <- 2 * 2^(chain[j,3]/nvoice);      if (chain[j,2] > (para*phi.x.min) ){ 

         cat("Chain number",j);
         phi.x.max <- 2 * 2^(chain[j,(2+chain[j,2])]/nvoice);         x.min <- chain[j,1];         x.max <- chain[j,1] + chain[j,2] - 1;         x.min <- x.min - round(para * phi.x.min);         x.max <- x.max + round(para * phi.x.max);
         tmp2 <- irregrec(siginput[chain[j,1]:(chain[j,1]+chain[j,2]-1)],
            inputwt[chain[j,1]:(chain[j,1]+chain[j,2]-1),],
            chain[j,3:(chain[j,2]+2)], compr,noct,nvoice,
            epsilon, w0 = w0 , prob = prob, fast = fast, para = para,
            minnbnodes = minnbnodes, real = real);


         totnbnodes <- totnbnodes + tmp2$nbnodes;
         np <- length(tmp2$sol);         start <- max(1,x.min);         end <- min(sigsize,x.min+np-1);         start1 <- max(1,2-x.min);         end1 <- min(np, sigsize +1 - x.min);         end <- end1 - start1 + start;         rec[start:end] <- rec[start:end]+tmp2$sol[start1:end1];
         plnb <- plnb + 1;         p <- p + 1;         idx[p] <- j;  
       };   };
   if(plot == 1){
      par(mfrow=c(2,1));      par(cex=1.1);      tsplot(siginput);      title("Original signal");      title("Reconstructed signal");   };   else if (plot == 2){
      par(mfrow=c(plnb+2,1));      par(mar=c(2,4,4,4));      par(cex=1.1);      par(err=-1);      tsplot(siginput);      title("Original signal");
      for (j in 1:p)
         tsplot(tmp[idx[j],]);

      tsplot(Re(rec))	;      title("Reconstructed signal");

   };
   cat("Total number of ridge samples used: ",totnbnodes,"\n");
   par(mfrow=c(1,1));
   list(rec=rec, ordered=ordered,chain = chain, comp=tmp);}



crcirgrec <- function(siginput, inputgt, beemap, nvoice, freqstep, scale,
	compr, prob = 0.8, bstep = 5, ptile = .01, epsilon = .5, fast = T, para = 0,
	minnbnodes=3, hflag = F, plot = 2, real = F)
#########################################################################
#     crcirgrec:
#     ----------
#      Reconstruction of a real valued signal from ridges found by 
#      crazy climbers on the gabor transform (irregularly sampled).
#
#      input:
#      ------
#      siginput: input signal
#      inputgt: continuous gabor transform (output of cgt)
#      beemap: output of crazy climber algorithm
#      nvoice: number of different frequencies
#      freqstep: difference between two consecutive frequencies
#      scale: scale of the window
#      compr: subsampling rate for the ridge
#      bstep: used for the chaining
#      ptile: 
#      epsilon: coeff of the Q2 part of reconstruction kernel
#      fast: if set to TRUE, computes the Q2 kernel by Riemann sums
#            if not, uses a Romberg adaptive step quadrature.
#      para:
#      prob: percentile to keep 
#      minnbnodes: minimal number of nodes for a sampled ridge.
#      hflag: if set to FALSE, uses the identity as first term
#             in the reconstruction kernel. If not, uses Q1 instead. 
#      plot: if set to 1, displays the signal, the components, and
#            signal and reconstruction one after another. If set to
#            2, displays the signal, the components, and the
#            reconstruction on the same page. Else, no plot.
#      real: if set to true, uses only real constraints.
#
#      output:
#      -------
#      rec: reconstructed signal
#      ordered: image of the ridges (with different colors)
#      comp: 2D array containing the signals reconstructed from ridges
#
#########################################################################
{
   tmp <- cfamily(beemap,bstep,ptile=ptile);   image(tmp$ordered);   chain <- tmp$chain;   nbchain <- tmp$nbchain;   ordered <- tmp$ordered;   sigsize <- length(siginput);   rec <- numeric(sigsize);   plnb <- 0;
   par(mfrow=c(1,1));   tsplot(siginput);   title("Original signal");
   totnbnodes <- 0;   tmp <- matrix(0,nbchain,length(siginput));   for (j in 1:nbchain){
      if (chain[j,2] > scale){
         nbnodes <- round(chain[j,2]/compr);
         if(nbnodes < minnbnodes)
            nbnodes <- minnbnodes;
         totnbnodes <- totnbnodes + nbnodes;
         cat("Chain number",j);
         phi.x.min <- scale;         phi.x.max <- scale;         x.min <- chain[j,1];         x.max <- chain[j,1] + chain[j,2] - 1;         x.min <- x.min - round(para * phi.x.min);         x.max <- x.max + round(para * phi.x.max);         x.inc <- 1;         np <- as.integer((x.max-x.min)/x.inc) + 1;
         tmp2 <- girregrec(siginput[chain[j,1]:(chain[j,1]+chain[j,2]-1)],
            inputgt[chain[j,1]:(chain[j,1]+chain[j,2]-1),],
            chain[j,3:(chain[j,2]+2)], nbnodes,nvoice,freqstep,scale,
            epsilon,fast,prob=prob,para=para, hflag = hflag, real = real);


         start <- max(1,x.min);         end <- min(sigsize,x.min+np-1);         start1 <- max(1,2-x.min);         end1 <- min(np, sigsize +1 - x.min);         end <- end1 - start1 + start;         rec[start:end] <- rec[start:end]+tmp2$sol[start1:end1];         tmp[j,start:end] <- tmp2$sol[start1:end1];
         plnb <- plnb + 1;         tsplot(tmp[j,]);
      };
   };
   if(plot == 1){
      par(mfrow=c(2,1));      par(cex=1.1);      tsplot(siginput);      title("Original signal");      tsplot(rec);      title("Reconstructed signal");   };   else if (plot == 2){
      par(mfrow=c(plnb+2,1));      par(mar=c(2,4,4,4));      par(cex=1.1);      par(err=-1);      tsplot(siginput);      title("Original signal");
      for (j in 1:nbchain)
         tsplot(tmp[j,]);
        
      tsplot(rec);      title("Reconstructed signal");   };
   cat("Total number of ridge samples used: ",totnbnodes,"\n");
   par(mfrow=c(1,1));   list(rec=rec, ordered=ordered,chain = chain, comp=tmp);}















      tsplot(rec)
      title("Reconstructed signal")
   }
#~error: statement with syntax error
#~put ".Data/last.dump" 904097046 ...
#~
.System("more Crazy_Cimbers.S")
#~put ".Data/.Last.value" 904097510 ...
#~
.System("more Crazy_Climbers.S")
#~put ".Data/.Last.value" 904098026 ...
#~
.System("more Crazy_Climbers.S")
#~put ".Data/.Last.value" 904098052 ...
#~
#~New session: Version: S Tue Jul 16 10:52:22 PDT 1996, Time: Tue Aug 25 23:13:24 1998
#~get "/usr/local/Splus3.4/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/splus/.Functions/.laenv" -1 ...
#########################################################################
#      $Log: Crc_Rec.S,v $
# Revision 1.2  1995/04/05  18:56:55  bruno
# *** empty log message ***
#
# Revision 1.1  1995/04/02  01:04:16  bruno
# Initial revision
#
#
#               (c) Copyright  1997                             
#                          by                                   
#      Author: Rene Carmona, Bruno Torresani, Wen-Liang Hwang   
#                  Princeton University
#                  All right reserved                           
#########################################################################


#########################################################################
#
#	Functions to reconstruct a signal from the output of
#	 the crazy climber algorithm.
#
#########################################################################




crcirrec <- function(siginput, inputwt, beemap, noct, nvoice, compr,
	minnbnodes = 2, w0 = 2*pi, bstep = 5,ptile = .01, prob = 0.8,
	epsilon = .5, fast = F, para = 0, real = F, plot=1)
#########################################################################
#     crcirrec:
#     -------
#      Reconstruction of a real valued signal from ridges found by 
#      crazy climbers on a wavelet transform (sampled irregularly).
#
#      input:
#      ------
#      siginput: input signal
#      inputwt: continuous wavelet transform (output of cwt)
#      beemap: output of crazy climber algorithm
#      noct: number of octaves (powers of 2)
#      nvoice: number of different scales per octave
#      compr: subsampling rate for the ridge
#      bstep: used for the chaining
#      ptile: 
#      epsilon: coeff of the Q2 part of reconstruction kernel
#      fast: if set to TRUE, computes the Q2 kernel by Riemann sums
#            if not, uses a Romberg adaptive step quadrature.
#      para:
#      plot: plot the original and reconstructed signal in display
#      prob: the percentile to keep 
#
#      output:
#      -------
#      rec: reconstructed signal
#      ordered: image of the ridges (with different colors)
#      comp: 2D array containing the signals reconstructed from ridges
#
#########################################################################
{

   tmp <- cfamily(beemap,bstep,ptile=ptile);   image(tmp$ordered);   chain <- tmp$chain;   nbchain <- tmp$nbchain;   ordered <- tmp$ordered;   sigsize <- length(siginput);   rec <- numeric(sigsize);   plnb <- 0;
   par(mfrow=c(1,1));   tsplot(siginput);   title("Original signal");
   tmp <- matrix(0,nbchain,length(siginput));
   totnbnodes <- 0;   idx <- numeric(nbchain)	;   p <- 0;
   for (j in 1:nbchain){
      phi.x.min <- 2 * 2^(chain[j,3]/nvoice);      if (chain[j,2] > (para*phi.x.min) ){ 

         cat("Chain number",j);
         phi.x.max <- 2 * 2^(chain[j,(2+chain[j,2])]/nvoice);         x.min <- chain[j,1];         x.max <- chain[j,1] + chain[j,2] - 1;         x.min <- x.min - round(para * phi.x.min);         x.max <- x.max + round(para * phi.x.max);
         tmp2 <- irregrec(siginput[chain[j,1]:(chain[j,1]+chain[j,2]-1)],
            inputwt[chain[j,1]:(chain[j,1]+chain[j,2]-1),],
            chain[j,3:(chain[j,2]+2)], compr,noct,nvoice,
            epsilon, w0 = w0 , prob = prob, fast = fast, para = para,
            minnbnodes = minnbnodes, real = real);


         totnbnodes <- totnbnodes + tmp2$nbnodes;
         np <- length(tmp2$sol);         start <- max(1,x.min);         end <- min(sigsize,x.min+np-1);         start1 <- max(1,2-x.min);         end1 <- min(np, sigsize +1 - x.min);         end <- end1 - start1 + start;         rec[start:end] <- rec[start:end]+tmp2$sol[start1:end1];
         plnb <- plnb + 1;         p <- p + 1;         idx[p] <- j;  
       };   };
   if(plot == 1){
      par(mfrow=c(2,1));      par(cex=1.1);      tsplot(siginput);      title("Original signal");      title("Reconstructed signal");   };   else if (plot == 2){
      par(mfrow=c(plnb+2,1));      par(mar=c(2,4,4,4));      par(cex=1.1);      par(err=-1);      tsplot(siginput);      title("Original signal");
      for (j in 1:p)
         tsplot(tmp[idx[j],]);

      tsplot(Re(rec))	;      title("Reconstructed signal");

   };
   cat("Total number of ridge samples used: ",totnbnodes,"\n");
   par(mfrow=c(1,1));
   list(rec=rec, ordered=ordered,chain = chain, comp=tmp);}



crcirgrec <- function(siginput, inputgt, beemap, nvoice, freqstep, scale,
	compr, prob = 0.8, bstep = 5, ptile = .01, epsilon = .5, fast = T, para = 0,
	minnbnodes=3, hflag = F, plot = 2, real = F)
#########################################################################
#     crcirgrec:
#     ----------
#      Reconstruction of a real valued signal from ridges found by 
#      crazy climbers on the gabor transform (irregularly sampled).
#
#      input:
#      ------
#      siginput: input signal
#      inputgt: continuous gabor transform (output of cgt)
#      beemap: output of crazy climber algorithm
#      nvoice: number of different frequencies
#      freqstep: difference between two consecutive frequencies
#      scale: scale of the window
#      compr: subsampling rate for the ridge
#      bstep: used for the chaining
#      ptile: 
#      epsilon: coeff of the Q2 part of reconstruction kernel
#      fast: if set to TRUE, computes the Q2 kernel by Riemann sums
#            if not, uses a Romberg adaptive step quadrature.
#      para:
#      prob: percentile to keep 
#      minnbnodes: minimal number of nodes for a sampled ridge.
#      hflag: if set to FALSE, uses the identity as first term
#             in the reconstruction kernel. If not, uses Q1 instead. 
#      plot: if set to 1, displays the signal, the components, and
#            signal and reconstruction one after another. If set to
#            2, displays the signal, the components, and the
#            reconstruction on the same page. Else, no plot.
#      real: if set to true, uses only real constraints.
#
#      output:
#      -------
#      rec: reconstructed signal
#      ordered: image of the ridges (with different colors)
#      comp: 2D array containing the signals reconstructed from ridges
#
#########################################################################
{
   tmp <- cfamily(beemap,bstep,ptile=ptile);   image(tmp$ordered);   chain <- tmp$chain;   nbchain <- tmp$nbchain;   ordered <- tmp$ordered;   sigsize <- length(siginput);   rec <- numeric(sigsize);   plnb <- 0;
   par(mfrow=c(1,1));   tsplot(siginput);   title("Original signal");
   totnbnodes <- 0;   tmp <- matrix(0,nbchain,length(siginput));   for (j in 1:nbchain){
      if (chain[j,2] > scale){
         nbnodes <- round(chain[j,2]/compr);
         if(nbnodes < minnbnodes)
            nbnodes <- minnbnodes;
         totnbnodes <- totnbnodes + nbnodes;
         cat("Chain number",j);
         phi.x.min <- scale;         phi.x.max <- scale;         x.min <- chain[j,1];         x.max <- chain[j,1] + chain[j,2] - 1;         x.min <- x.min - round(para * phi.x.min);         x.max <- x.max + round(para * phi.x.max);         x.inc <- 1;         np <- as.integer((x.max-x.min)/x.inc) + 1;
         tmp2 <- girregrec(siginput[chain[j,1]:(chain[j,1]+chain[j,2]-1)],
            inputgt[chain[j,1]:(chain[j,1]+chain[j,2]-1),],
            chain[j,3:(chain[j,2]+2)], nbnodes,nvoice,freqstep,scale,
            epsilon,fast,prob=prob,para=para, hflag = hflag, real = real);


         start <- max(1,x.min);         end <- min(sigsize,x.min+np-1);         start1 <- max(1,2-x.min);         end1 <- min(np, sigsize +1 - x.min);         end <- end1 - start1 + start;         rec[start:end] <- rec[start:end]+tmp2$sol[start1:end1];         tmp[j,start:end] <- tmp2$sol[start1:end1];
         plnb <- plnb + 1;         tsplot(tmp[j,]);
      };
   };
   if(plot == 1){
      par(mfrow=c(2,1));      par(cex=1.1);      tsplot(siginput);      title("Original signal");      tsplot(rec);      title("Reconstructed signal");   };   else if (plot == 2){
      par(mfrow=c(plnb+2,1));      par(mar=c(2,4,4,4));      par(cex=1.1);      par(err=-1);      tsplot(siginput);      title("Original signal");
      for (j in 1:nbchain)
         tsplot(tmp[j,]);
        
      tsplot(rec);      title("Reconstructed signal");   };
   cat("Total number of ridge samples used: ",totnbnodes,"\n");
   par(mfrow=c(1,1));   list(rec=rec, ordered=ordered,chain = chain, comp=tmp);}















      tsplot(rec)
      title("Reconstructed signal")
   }
#~error: statement with syntax error
#~put ".Data/last.dump" 904101118 ...
#~
.First
#~
source("Util.S")
#~put ".Data/smoothwt" 904101140 ...
#~put ".Data/smoothts" 904101140 ...
#~put ".Data/adjust.length" 904101140 ...
#~put ".Data/npl" 904101140 ...
#~put ".Data/SpecGen" 904101140 ...
#~put ".Data/SampleGen" 904101140 ...
#~put ".Data/hurst.est" 904101140 ...
#~put ".Data/wspec.pl" 904101140 ...
#~put ".Data/.Last.value" 904101140 ...
#~
source("extrema.S")
#~put ".Data/ext" 904101144 ...
#~put ".Data/.Last.value" 904101144 ...
#~
source("svd.S")
#~put ".Data/SVD" 904101148 ...
#~put ".Data/.Last.value" 904101148 ...
#~
source("mw.S")
#~put ".Data/check.maxresoln" 904101151 ...
#~put ".Data/mw" 904101151 ...
#~put ".Data/.Last.value" 904101151 ...
#~
source("mreconst.S")
#~put ".Data/mrecons" 904101157 ...
#~put ".Data/dwinverse" 904101157 ...
#~put ".Data/.Last.value" 904101157 ...
#~
source("plot.S")
#~put ".Data/plotwt" 904101161 ...
#~put ".Data/plot.result" 904101161 ...
#~put ".Data/wpl" 904101161 ...
#~put ".Data/epl" 904101161 ...
#~put ".Data/.Last.value" 904101161 ...
#~
trellis.device(motif)
#~get "/usr/local/Splus3.4/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.paper.regions" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/library/trellis/.Data/xcm.trellis.motif" -1 ...
#~get "/usr/local/Splus3.4/library/trellis/.Data/trellis.settings.motif" -1 ...
#~put "Session_DB$.Device" 0 ...
#~put "Session_DB$.Devices" 0 ...
#~put ".Data/.Last.value" 904101167 ...
#~
dyn.load.shared("./pack.so")
#~put ".Data/.Last.value" 904101180 ...
#~
attach("/home/chelsea/rcarmona/data/siemens/blood/.Data")
#~put ".Data/.Last.value" 904101209 ...
#~
ls()
#~put ".Data/.Last.value" 904101211 ...
#~
.System(" ls /home/chelsea/rcarmona/data/siemens/blood/.Data")
#~put ".Data/.Last.value" 904101242 ...
#~
tsplotF2N)
#~error: statement with syntax error
#~
tsplot(F2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/F2N" 904079705 ...
#~
wpl(mF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/mF2N" 904081504 ...
#~get "/usr/local/Splus3.4/s/.Datasets/gr.pars" -1 ...
#~
epl(eF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/eF2N" 904080323 ...
#~get "/usr/local/Splus3.4/s/.Datasets/gr.pars" -1 ...
#~
wpl(mF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/mF2N" 904081504 ...
#~get "/usr/local/Splus3.4/s/.Datasets/gr.pars" -1 ...
#~
wpl(mF2N, scale = T)
#~put ".Data/last.dump" 904101584 ...
#~
source("plot.S")
#~put ".Data/plotwt" 904101699 ...
#~put ".Data/plot.result" 904101699 ...
#~put ".Data/wpl" 904101700 ...
#~put ".Data/epl" 904101700 ...
#~put ".Data/.Last.value" 904101700 ...
#~
wpl(mF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/mF2N" 904081504 ...
#~get "/usr/local/Splus3.4/s/.Datasets/gr.pars" -1 ...
#~
source("plot.S")
#~put ".Data/plotwt" 904101798 ...
#~put ".Data/plot.result" 904101798 ...
#~put ".Data/wpl" 904101798 ...
#~put ".Data/epl" 904101798 ...
#~put ".Data/.Last.value" 904101798 ...
#~
TST <- mreconst(teF2N)
#~error: couldn't find function "mreconst"
#~
mreconst
#~error: Object "mreconst" not found
#~
source("mreconst.S")
#~put ".Data/mrecons" 904101856 ...
#~put ".Data/dwinverse" 904101856 ...
#~put ".Data/.Last.value" 904101856 ...
#~
TST <- mrecons(teF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/teF2N" 904080753 ...
#~put ".Data/TST" 904101940 ...
#~
mrecons
#~
tsplot(TST$g)
#~get ".Data/TST" 904101940 ...
#~
wpl(mF2N)
#~get "/home/chelsea/rcarmona/data/siemens/blood/.Data/mF2N" 904081504 ...
#~get "/usr/local/Splus3.4/s/.Datasets/gr.pars" -1 ...
#~
erecF2N <- TST
#~get ".Data/TST" 904101940 ...
#~put ".Data/erecF2N" 904102145 ...
#~
mrecF2N <- dwinverse(mF2N)
#~error: "Sinverse_wavelet_transform" is not a symbol in the load table
#~put ".Data/last.dump" 904102226 ...
#~
source("mw.S")
#~put ".Data/check.maxresoln" 904102247 ...
#~put ".Data/mw" 904102247 ...
#~put ".Data/.Last.value" 904102247 ...
#~
source("Util.S")
#~put ".Data/smoothwt" 904102256 ...
#~put ".Data/smoothts" 904102256 ...
#~put ".Data/adjust.length" 904102256 ...
#~put ".Data/npl" 904102256 ...
#~put ".Data/SpecGen" 904102256 ...
#~put ".Data/SampleGen" 904102256 ...
#~put ".Data/hurst.est" 904102256 ...
#~put ".Data/wspec.pl" 904102256 ...
#~put ".Data/.Last.value" 904102257 ...
#~
mrecF2N <- dwinverse(mF2N)
#~error: "Sinverse_wavelet_transform" is not a symbol in the load table
#~put ".Data/last.dump" 904102286 ...
#~
dyn.load.shared("./pack.so")
#~put ".Data/.Last.value" 904102306 ...
#~
mrecF2N <- dwinverse(mF2N)
#~error: "Sinverse_wavelet_transform" is not a symbol in the load table
#~put ".Data/last.dump" 904102311 ...
#~
q()
#~End session: Time: 904102442; Process: 10985
.System("pwd")
#~put ".Data/.Last.value" 904102366 ...
#~
mrecF2N <- dwinverse(mF2N)
#~error: "Sinverse_wavelet_transform" is not a symbol in the load table
#~put ".Data/last.dump" 904102383 ...
#~
.System("ls *.c")
#~put ".Data/.Last.value" 904102407 ...
#~
q()
#~End session: Time: 904102719; Process: 11222
#~New session: Version: S Tue Jul 16 10:52:22 PDT 1996, Time: Tue Aug 25 23:41:55 1998
#~get "/usr/local/Splus3.4/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/splus/.Functions/.laenv" -1 ...
#~error: Can't load (dlopen) library ./pack.so: 12332:/usr/local/Splus3.4/cmd/Sqpe: rld: Fatal Error: unresolvable symbol in ./pack.so: signal_W_tilda
#~put ".Data/last.dump" 904102830 ...
#~
.First
#~
source("Util.S")
#~put ".Data/smoothwt" 904102845 ...
#~put ".Data/smoothts" 904102845 ...
#~put ".Data/adjust.length" 904102845 ...
#~put ".Data/npl" 904102845 ...
#~put ".Data/SpecGen" 904102846 ...
#~put ".Data/SampleGen" 904102846 ...
#~put ".Data/hurst.est" 904102846 ...
#~put ".Data/wspec.pl" 904102846 ...
#~put ".Data/.Last.value" 904102846 ...
#~
source("extrema.S")
#~put ".Data/ext" 904102849 ...
#~put ".Data/.Last.value" 904102849 ...
#~
dyn.load.shared("./pack.so")
#~error: Can't load (dlopen) library ./pack.so: 12332:/usr/local/Splus3.4/cmd/Sqpe: rld: Fatal Error: unresolvable symbol in ./pack.so: signal_W_tilda
#~put ".Data/last.dump" 904102856 ...
#~
q()
#~End session: Time: 904103009; Process: 12332
#~New session: Version: S Tue Jul 16 10:52:22 PDT 1996, Time: Tue Aug 25 23:45:38 1998
#~get "/usr/local/Splus3.4/splus/.Datasets/..min.script.id" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/.Copyright" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/splus/.Functions/.laenv" -1 ...
#~error: Can't load (dlopen) library ./pack.so: 13382:/usr/local/Splus3.4/cmd/Sqpe: rld: Fatal Error: unresolvable symbol in ./pack.so: signal_W_tilda
#~put ".Data/last.dump" 904103053 ...
#~
.First
#~
options(object.size = 1e+21)
#~put ".Data/.Last.value" 904103077 ...
#~
dyn.load.shared("./pack.so")
#~error: Can't load (dlopen) library ./pack.so: 13382:/usr/local/Splus3.4/cmd/Sqpe: rld: Fatal Error: unresolvable symbol in ./pack.so: signal_W_tilda
#~put ".Data/last.dump" 904103080 ...
#~
source("plot.S")
#~put ".Data/plotwt" 904103087 ...
#~put ".Data/plot.result" 904103087 ...
#~put ".Data/wpl" 904103088 ...
#~put ".Data/epl" 904103088 ...
#~put ".Data/.Last.value" 904103088 ...
#~
source("mreconst.S")
#~put ".Data/mrecons" 904103091 ...
#~put ".Data/dwinverse" 904103091 ...
#~put ".Data/.Last.value" 904103091 ...
#~
source("mw.S")
#~put ".Data/check.maxresoln" 904103094 ...
#~put ".Data/mw" 904103094 ...
#~put ".Data/.Last.value" 904103094 ...
#~
source("svd.S")
#~put ".Data/SVD" 904103097 ...
#~put ".Data/.Last.value" 904103097 ...
#~
mw
#~
mwinverse
#~error: Object "mwinverse" not found
#~
dwinverse
#~
mrecF2N <- dwinverse(mF2N)
#~error: Object "mF2N" not found
#~put ".Data/last.dump" 904103160 ...
#~
attach("/home/chelsea/data/siemens/blood/.Data")
#~error: Can't find "/home/chelsea/data/siemens/blood/.Data" to use as a database
#~put ".Data/last.dump" 904103184 ...
#~
attach("/home/chelsea/rcarmona/data/siemens/blood/.Data")
#~put ".Data/.Last.value" 904103195 ...
#~
mrecF2N <- dwinverse(mF2N)
#~error: "Sinverse_wavelet_transform" is not a symbol in the load table
#~put ".Data/last.dump" 904103201 ...
#~
.System("more dwinverse.c")
#~put ".Data/.Last.value" 904103254 ...
#~
dyn.load.shared("./pack.so")
#~error: Can't load (dlopen) library ./pack.so: 13382:/usr/local/Splus3.4/cmd/Sqpe: rld: Fatal Error: unresolvable symbol in ./pack.so: signal_W_tilda
#~put ".Data/last.dump" 904103342 ...
#~
tsplot(erecF2N)
#~put ".Data/last.dump" 904104398 ...
#~
ls()
#~put ".Data/.Last.value" 904104411 ...
#~
tsplot(erecF2N)
#~put ".Data/last.dump" 904104468 ...
#~
trellis.device(motif)
#~get "/usr/local/Splus3.4/splus/.Datasets/.PostScript.Options" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.preamble.std" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.setcolor.rgb" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.setfont.std" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/ps.paper.regions" -1 ...
#~get "/usr/local/Splus3.4/splus/.Datasets/version" -1 ...
#~get "/usr/local/Splus3.4/library/trellis/.Data/xcm.trellis.motif" -1 ...
#~get "/usr/local/Splus3.4/library/trellis/.Data/trellis.settings.motif" -1 ...
#~put "Session_DB$.Device" 0 ...
#~put "Session_DB$.Devices" 0 ...
#~put ".Data/.Last.value" 904104526 ...
#~
tsplot(erecF2N)
#~put ".Data/last.dump" 904104531 ...
#~
tsplot(erecF2N$g)
#~get ".Data/erecF2N" 904102145 ...
#~
.System("mv ps.out.0001.ps tess.rec.ps")
#~put ".Data/.Last.value" 904104583 ...
#~
q()
#~End session: Time: 904104673; Process: 13382
