FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;

FILENAME url URL "http://192.38.117.59/~linearpredictors/datafiles/readVitaminD.sas";
%include url;

/*
All women
*/ 

DATA women;
	SET vitamind (WHERE = (category=2));

        logvitd=log10(vitd);
        logvitdintake=log10(vitdintake);

RUN;

title 'medians for adjusting estimates';
proc means N median DATA=women;
    var bmi logvitdintake;
run;

options ls=85 ps=300;

title 'Model with two interactions';
title2 'Specific estimates at the buttom';
PROC GLM DATA=women;
    CLASS sunexp country;
    MODEL logvitd = country bmi
                    country*logvitdintake country*sunexp
                  / noint;
    format country countryf.;
    format sunexp sunf.;
ESTIMATE 'Effect of log(vitamin D intake), DK' country*logvitdintake 1 0 0 0;
ESTIMATE 'Effect of log(vitamin D intake), SF' country*logvitdintake 0 1 0 0;
ESTIMATE 'Effect of log(vitamin D intake), EI' country*logvitdintake 0 0 1 0;
ESTIMATE 'Effect of log(vitamin D intake), PL' country*logvitdintake 0 0 0 1;

ESTIMATE 'sometimes, adjusted to medians, DK'
    country 1 0 0 0 bmi 27.153 country*logvitdintake 0.8961 0 0 0
    country*sunexp 0 0 0 0 0 0 0 0 1 0 0 0;
ESTIMATE 'sometimes, adjusted to medians, SF'
    country 0 1 0 0 bmi 27.153 country*logvitdintake 0 0.8961 0 0 
    country*sunexp 0 0 0 0 0 0 0 0 0 1 0 0;
ESTIMATE 'sometimes, adjusted to medians, EI'
    country 0 0 1 0 bmi 27.153 country*logvitdintake 0 0 0.8961 0 
    country*sunexp 0 0 0 0 0 0 0 0 0 0 1 0;
ESTIMATE 'sometimes, adjusted to medians, PL'
    country 0 0 0 1 bmi 27.153 country*logvitdintake 0 0 0 0.8961
    country*sunexp 0 0 0 0 0 0 0 0 0 0 0 1;

ESTIMATE 'Prefer vs. sometimes, DK' country*sunexp 0 0 0 0 1 0 0 0 -1 0 0 0;
ESTIMATE 'Prefer vs. sometimes, SF' country*sunexp 0 0 0 0 0 1 0 0 0 -1 0 0;
ESTIMATE 'Prefer vs. sometimes, EI' country*sunexp 0 0 0 0 0 0 1 0 0 0 -1 0;
ESTIMATE 'Prefer vs. sometimes, PL' country*sunexp 0 0 0 0 0 0 0 1 0 0 0 -1;

ESTIMATE 'Avoid vs. sometimes, DK' country*sunexp 1 0 0 0 0 0 0 0 -1 0 0 0;
ESTIMATE 'Avoid vs. sometimes, SF' country*sunexp 0 1 0 0 0 0 0 0 0 -1 0 0;
ESTIMATE 'Avoid vs. sometimes, EI' country*sunexp 0 0 1 0 0 0 0 0 0 0 -1 0;
ESTIMATE 'Avoid vs. sometimes, PL' country*sunexp 0 0 0 1 0 0 0 0 0 0 0 -1;

run;