R/lme_mass_LR.R

In Deep-MI/fslmer: Univariate And Mass-Univariate Linear Mixed-Effects Analysis For FreeSurfer Imaging Data

lme_mass_LR<-function(statsfull,statsred,q)
{
    # --------------------------------------------------------------------------
    
    if (exists("nargin") == TRUE) 
    {
        if (nargin < 3) 
        {
            stop(cat("Error: Too Few Inputs"))
        } 
    } 
    
    # --------------------------------------------------------------------------
    
    nv<-length(statsfull)
    pval<-matrix(0,nrow=1,ncol=nv)
    
    for (i in c(1:nv))
    {
        lrstats<-lme_LR(statsfull[i]$lreml,statsred[i]$lreml,q)
        pval[i]<-lrstats$pval
    }
    
    return(pval)    
}