docs/Self_test.R

In shimlab/HMTree: Hidden Markov tree prior for analysis of multiple functional data

#Instructions: There are two parts, use part two to compare
#Part One
# Create a data set according to old(also known as correct) package
# library("devtools", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
# library("HMTree", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
# library("Rcpp", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
# library("roxygen2", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
# home_dir<-"/Users/haosicheng/Desktop/HMTree"
# setwd(home_dir)
# set.seed(666)
# normal.data<-simu.curves(curve.length = 1024,model.mode = "Normal",num.samples = 1)
# nout<-HMTree.denoise_1D(normal.data)
# set.seed(666)
# pois.data<-simu.curves(curve.length = 1024,model.mode = "Poisson",num.samples = 1)
# pout<-HMTree.denoise_1D(pois.data)

#save the object
# save(pout,file = "pout.RData")
# save(nout,file = "nout.RData")

#Part Two:
##After change the old file:
home_dir<-"/Users/haosicheng/Desktop/HMTree/docs"
#Reload the object 
setwd(home_dir)
load(file = "pout.RData")
load(file = "nout.RData")


#generate a new object
library("devtools", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
library("HMTree", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
library("Rcpp", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
library("roxygen2", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")

setwd(home_dir)
set.seed(666)
normal.data<-simu.curves(curve.length = 1024,model.mode = "Normal",num.samples = 1)
nout.new<-HMTree.denoise_1D(normal.data)
set.seed(666)
pois.data<-simu.curves(curve.length = 1024,model.mode = "Poisson",num.samples = 1)
pout.new<-HMTree.denoise_1D(pois.data)
#compare
#Alpha_vec
old<-unlist(pout$alpha_vec)
new<-unlist(pout.new$alpha_vec)
if(sum(old-new)!=0){
  print("first object wrong!!!")
}
#wavelet_coef
old<-unlist(pout$wavelet_coef)
new<-unlist(pout.new$wavelet_coef)
if(sum(old-new)!=0){
  print("second object wrong!!!")
}
#denoised_wavelet_coef
old<-unlist(pout$denoised_wavelet_coef)
new<-unlist(pout.new$denoised_wavelet_coef)
if(sum(old-new)!=0){
  print("third object wrong!!!")
}
#denoised
old<-unlist(pout$denoised)
new<-unlist(pout.new$denoised)
if(sum(old-new)!=0){
  print("fourth object wrong!!!")
}