function [ROC GKR RES]=NRL_ROC( file, subject, pov, print )

    if( nargin == 3 )
        print = 0;
    end
    
    load Labels2.mat;

    % This would be more efficient
%    labels = AU_1050_1( :, [1 2 4 5 9 10 12 14 15 17 20] );
%    labels = eval( [ 'load Labels2.mat AU_' subject '_' pov ] );

%    labels = labels( :, [1 2 4 5 9 10 12 14 15 17 20] );
    labels = eval( [ 'AU_' subject '_' pov '( :, [1 2 4 5 9 10 12 14 15 17 20] )' ] );
    labels(:,12) = labels(:,1) & ~labels(:,2);
    labels(:,13) = labels(:,1) & labels(:,2) & labels(:,3);
    labels=logical(labels);
%    row=labels(1,:);

    allAUs = load(file);
    AUvalues = allAUs( :, [ 2:12 19:22 ] );
    GKRvalues = allAUs( :, 49:59 );

    doFirstLast = 1;
    if( doFirstLast )
        [first last] = getFirstLastFrames( strread(subject), strread(pov) );
        AUvalues = AUvalues(first:last,:);
        GKRvalues = GKRvalues(first:last,:);
        labels = labels(first:last,:);
    end

    skipMissedFaces = 0;

    if( skipMissedFaces == 1 )
        K = setdiff( 1:size(AUvalues,1), find(isnan(AUvalues(:,1))) );
    end

%    for i=1:10

    if( skipMissedFaces == 1 )
        ROC=ROC_area( AUvalues(K,:), labels(K,[1:11 12 12 13 13]) );
    else
        ROC=ROC_area( AUvalues, labels(:,[1:11 12 12 13 13]) );
    end

%        disp(nanmean(ROC));
%        labels=[labels(2:3324,:); row];
%    end

    if( skipMissedFaces == 1 )
        GKR=ROC_area( GKRvalues(K,:), labels(K,[1:11]) );
    else
        GKR=ROC_area( GKRvalues, labels );
    end

    tab=sprintf( '\t' );
    nl=sprintf( '\n' );
    str1=sprintf( '%0.4f\t', ROC );
    str2=sprintf( '%0.4f\t', GKR );

    RES = [file nl tab 'No CKF:' tab str1 nl tab 'CKF:' tab str2 nl nl];

    if( print )
        disp( RES );
    end
end
