menses_test.m

beh_type = 'BOTH'
beh = false
dirstring = 'menses_test'

% Select variables to use for the model by prefix or by specific name
    pats =["Weather_TMAX", "Moon Phase", "Allergen_Total",...
        "Sleep", "BM", "Menses"];
    % set number of lagged days
    lag =  1;
        mkdir(dirstring);
    % Read in the data
    [Xc, yc, pc, ynan, p, text_vars] = read_data(beh_type, lag, beh,...
        dirstring, pats);
    
       
    females = load('females.csv');