R/GDM_swHM_figure.R

In benliemory/adaptiveHM: Adaptive Hierarchical Model (adaptiveHM)

GDM_swHM_figure <-
    function(Control, IPBT.prior=FALSE,winSizeRanges = seq(from = 5,to = 200,by = 5),
             history=NA,IPBT.id=NA)
    {
        if ( (IPBT.prior==FALSE & all(is.na(history)==1)  ) | 
                 (IPBT.prior==TRUE &  all(is.na(IPBT.id)==1) )  ) 
            stop("Historical information is missing!\nPlease provide historical data or use IPBT prior!") 
        
        if(IPBT.prior==FALSE)
        {
            hist_var = apply(history,1,var)
        }
        
        if(IPBT.prior==TRUE)
        {   
            data(IPBT3digits)
            data(SampleSize)
            hist_var = IPBT3digits[,IPBT.id]^2
        }  
        


    
        GDM = sapply(winSizeRanges,function(size){GDM_swHM(hist_var,size,Control) } )
               
        plot(winSizeRanges,GDM,xlab="Window Size",ylab = "GDM",col="blue",lwd = 2 )
        
        GDM
        
}