topic Re: Piecewise linear regression in SAS Programming

Piecewise linear regression

PunkinSAS08 — Sat, 09 Dec 2023 06:10:03 GMT

Hi all,

I am completing an assignment that involves constructing new variables in my code I am supposed to:

1. create a formula for BMI

2. compute a Categorical BMI variable in which 1=normal (<25) , 2=overweight (25 to <30) and 3=obese

3. three piecewise variables for the BMI intervals as defined in the Categorical BMI.

My question is, does my code satisfy the piecewise variables? In our lecture, there are a lot of long and drawn out ways to create piecewise variables but I want to ensure this is being done correctly for analysis. Here is my code:

proc format;

value BMIcatf 1='normal' 2='overweight' 3= 'obese';

run;

data final;

set totalcohort;

weight_kg = (weight / 2.205);

height_m = (height * 0.0254);

BMI = (weight_kg) / (height_m **2); *should we have to square the value for correct BMI?*;

if BMI >=19;

if BMI <25 then BMIcat = 1;

else if BMI >= 25 and BMI < 30 then BMIcat = 2;

else BMIcat = 3;

BMI1 = BMI < 25;

BMI2 = BMI >= 25 and BMI <30;

BMI3 = BMI >=30;

waist_hip_ratio = (waist / hip);

meanSBP = (bp_1s + bp_2s) / 2;

log_glyhb = log(glyhb);

if gender = 'male' then male =1;

else male = 0;

format BMIcat BMIcatf.;

run;

Re: Piecewise linear regression

Ksharp — Sat, 09 Dec 2023 09:39:04 GMT

Check @Rick_SAS 's blogs:
https://blogs.sas.com/content/iml/2017/04/05/nonsmooth-models-spline-effects.html
https://blogs.sas.com/content/iml/2020/12/14/segmented-regression-sas.html

Re: Piecewise linear regression

Rick_SAS — Sun, 10 Dec 2023 10:39:20 GMT

> 3. three piecewise variables for the BMI intervals

You have constructed one categorical variable (BMIcat), which has the values 1,2,3.

You have created three binary indicator variables (often called "dummy variables"), which you have named BMI1, BMI2, and BMI3.

I don't know how your instructor defines a "piecewise variable." If you use that term to mean a binary indicator variable, then this looks fine. If it means something else, we need to know the definition.

I would delete the line

if BMI >=19;

unless you purposely want to exclude underweight people.