Hi Everyone,
I need your help to be able to use the appropriate statistical procedure for my dataset below
Here is the detail of the data.
sub_id = subjects, measurement is taken on subject at different time point.
position: where measurement is taken (anterior or posterior)
time_in_seconds: time to apply treatment
resistant:
force: force of treatment application
The y variable is whether there is a block after treatment has been applied i.e. y =1 or there is a gap y= 0 and it also contains missing data.
Question:
I want to know how Time+Force+resistance affect the outcome i.e y =gap/block
Sub_ID | Position | Time_in_Seconds | Resistance_ ohms | Force_N | y_variable |
1 | Ant | 3.82 | -23 | 37 | 1 |
1 | Ant | 4.94 | -24 | 4.4 | 0 |
1 | Ant | 3.49 | -41.5 | 15.8 | 1 |
1 | Ant | 3.07 | -39.5 | 21.15 | 1 |
1 | Post | 3.43 | -39 | 29 | 1 |
1 | Post | 3.53 | -45.5 | 14.15 | 1 |
1 | Post | 4.44 | -46 | 9.55 | 1 |
1 | Post | 3.37 | -19 | 12.9 | 1 |
1 | Ant | 3.35 | -46.5 | 7.2 | 1 |
1 | Ant | 3.19 | -43 | 14.2 | 1 |
1 | Ant | 3.61 | -41 | 24.55 | 1 |
1 | Ant | 4.24 | -48 | 23.15 | 1 |
1 | Post | 2.09 | -33 | 27.25 | 1 |
1 | Post | 2.83 | -32 | 21.2 | 1 |
1 | Post | 3.26 | -28 | 29.7 | 0 |
1 | Post | 3.34 | -41.5 | 10.15 | 1 |
2 | Ant | 5.29 | -30 | 10.9 | 1 |
2 | Ant | 4.22 | -29 | 11.3 | 1 |
2 | Ant | 2.40 | -15 | 10.2 | 1 |
2 | Ant | 3.32 | -18.5 | 8.75 | 1 |
2 | Post | 1.85 | -27 | 9.3 | 1 |
2 | Post | 4.31 | -37 | 8.75 | 1 |
2 | Post | 1.82 | -29.5 | 13.95 | 1 |
2 | Post | 3.98 | -24.5 | 9.8 | 0 |
2 | Ant | 4.06 | -39.5 | 21.45 | 1 |
2 | Ant | 2.90 | -35.5 | 16.5 | 1 |
2 | Ant | 3.23 | -35.5 | 18.2 | 1 |
2 | Ant | 4.24 | -31 | 14.3 | 1 |
2 | Post | 2.45 | -31 | 9.6 | 1 |
2 | Post | 3.51 | -20 | 6 | 1 |
2 | Post | 4.27 | -17.5 | 8.4 | 1 |
2 | Post | 2.67 | -25.5 | 25 | 0 |
3 | Ant | 3.065092996 | -40 | 10.7 | 1 |
3 | Ant | 3.74 | -38 | 17.8 | . |
3 | Ant | 3.61 | -27 | 10.1 | 0 |
3 | Ant | 2.08 | -26.5 | 6.45 | . |
3 | Post | 2.12 | -35 | 20.4 | 1 |
3 | Post | 3.244 | -39 | 27.5 | 1 |
3 | Post | 4.02 | -42 | 19.9 | 1 |
3 | Post | 1.94 | -19 | 16.6 | 1 |
3 | Ant | 4.37 | -14 | 4.2 | 0 |
3 | Ant | 4.68 | -33 | 6.9 | 0 |
3 | Ant | 3.35 | -30.5 | 8.65 | 1 |
3 | Ant | 1.72 | -33 | 14.1 | 1 |
3 | Post | 0.81 | -27 | 12.2 | 1 |
3 | Post | 3.90 | -26 | 18.35 | 1 |
3 | Post | 4.19 | -29 | 9.4 | 1 |
3 | Post | 4.46 | -19 | 10.2 | 1 |
4 | Ant | 2.89 | -42 | 11.3 | 1 |
4 | Ant | 2.20 | -28 | 12.45 | 1 |
4 | Ant | 2.97 | -31 | 19.5 | 1 |
4 | Ant | 2.06 | -31 | 22.3 | 1 |
4 | Post | 3.35 | -44.5 | 32.9 | 1 |
4 | Post | 2.10 | -35 | 15.3 | 1 |
4 | Post | 3.42 | -35 | 8.35 | 1 |
4 | Post | 4.16 | -33 | 20.9 | 1 |
4 | Ant | 4.06 | -15.5 | 6 | 1 |
4 | Ant | 5.00 | -25 | 21.5 | 1 |
4 | Ant | 4.13 | -33.5 | 24.25 | 1 |
4 | Ant | 5.56 | -34 | 16.7 | 1 |
4 | Post | 4.14 | -35 | 31.75 | 1 |
4 | Post | 4.49 | -33.5 | 25.4 | 1 |
4 | Post | 4.17 | -29 | 41.9 | 1 |
4 | Post | 3.85 | -28 | 28.8 | 1 |
5 | Ant | 2.67 | -23 | 28.2 | 0 |
5 | Ant | 1.68 | -23 | 10.3 | 1 |
5 | Ant | 2.07 | -19.5 | 9.85 | 1 |
5 | Ant | 1.06 | -25 | 12.7 | 1 |
5 | Post | 5.02 | -31 | 10.4 | 0 |
5 | Post | 2.53 | -23 | 11.7 | 1 |
5 | Post | 3.40 | -71.5 | 27.15 | 1 |
5 | Post | 4.78 | -41.5 | 31.85 | 1 |
5 | Ant | 2.15 | -42 | 15.95 | 1 |
5 | Ant | 2.89 | -26.5 | 16.9 | 1 |
5 | Ant | 2.25 | -33 | 9.8 | 1 |
5 | Ant | 2.76 | -28 | 12.1 | 0 |
5 | Post | 3.04 | -22 | 9 | 1 |
5 | Post | 4.45 | -37.5 | 9.25 | 1 |
5 | Post | 4.10 | -20.5 | 15.3 | 1 |
5 | Post | 4.41 | -31.5 | 18.05 | 1 |
6 | Ant | 1.61 | -26 | 18.7 | 1 |
6 | Ant | 1.68 | -26 | 7.4 | 0 |
6 | Ant | 3.93 | -29 | 7.4 | 0 |
6 | Ant | 4.45 | -21.5 | 9.15 | . |
6 | Post | 5.48 | -28 | 25.05 | 1 |
6 | Post | 4.11 | -48 | 30.7 | 1 |
6 | Post | 3.20 | -23 | 22.3 | 1 |
6 | Post | 2.77 | -28 | 15.3 | 1 |
6 | Ant | 2.34 | -22 | 10.95 | 1 |
6 | Ant | 2.25 | -30 | 6.2 | 1 |
6 | Ant | 4.16 | -24.5 | 6.6 | |
6 | Ant | 4.62 | -33 | 12 | 0 |
6 | Post | 2.32 | -31 | 16.65 | 1 |
6 | Post | 4.03 | -31 | 19 | 1 |
6 | Post | 3.43 | -17 | 12 | . |
6 | Post | 3.51 | -14 | 11.1 | 0 |
7 | Ant | 2.99 | -30 | 7.65 | 1 |
7 | Ant | 1.80517419 | -25 | 15.95 | 1 |
7 | Ant | 2.106053494 | -32 | 13.8 | 1 |
7 | Ant | 3.096114016 | -29 | 17.55 | 1 |
7 | Post | 3.167542074 | -24 | 11.45 | 1 |
7 | Post | 3.338268984 | -24 | 22.7 | 1 |
7 | Post | 2.659685183 | -32.5 | 21.95 | 1 |
7 | Post | 3.751749917 | -17 | 8.25 | 1 |
7 | Ant | 2.197529839 | -25 | 8.8 | 0 |
7 | Ant | 3.664137015 | -35 | 10.1 | 0 |
7 | Ant | 3.545702335 | -39 | 9.4 | 0 |
7 | Ant | 2.023625001 | -35 | 10.1 | 1 |
7 | Post | 2.372086883 | -33 | 24.25 | 1 |
7 | Post | 3.582104056 | -37.5 | 17.3 | 1 |
7 | Post | 3.45055168 | -38 | 12.6 | 1 |
7 | Post | 3.841677068 | -23 | 9.8 | 1 |
8 | Ant | 3.432008272 | -25 | 18.6 | 0 |
8 | Ant | 2.136605495 | -32.5 | 14.1 | 1 |
8 | Ant | 1.841190586 | -18 | 7 | 0 |
8 | Ant | 2.15865836 | -25.5 | 6.8 | 0 |
8 | Post | 3.359805409 | -30 | 15.5 | 1 |
8 | Post | 3.631259765 | -31 | 28 | 1 |
8 | Post | 4.674356585 | -32.5 | 20.9 | 1 |
8 | Post | 4.044977037 | -25 | 12.9 | 0 |
8 | Ant | 3.346860731 | -28 | 14.6 | 1 |
8 | Ant | 3.850582629 | -46 | 24.5 | 1 |
8 | Ant | 5.340021635 | -31.5 | 8.5 | 1 |
8 | Ant | 3.980653721 | -26 | 11.6 | 1 |
8 | Post | 3.704121331 | -26 | 14.8 | 1 |
8 | Post | 3.848852913 | -31.5 | 15.8 | 1 |
8 | Post | 4.939191479 | -18 | 9.8 | 1 |
8 | Post | 3.134066196 | -15 | 9.1 | 1 |
9 | Ant | 1.309024248 | -25 | 9 | 1 |
9 | Ant | 3.369404446 | -27 | 11.55 | 1 |
9 | Ant | 1.841284373 | -26 | 13.15 | 1 |
9 | Ant | 3.675231524 | -21 | 12.4 | 1 |
9 | Post | 3.06826061 | -22 | 24.2 | 1 |
9 | Post | 3.59966626 | -19 | 17.3 | 1 |
9 | Post | 4.466268907 | -16.5 | 15.7 | 1 |
9 | Post | 1.882204503 | -27 | 11.3 | . |
9 | Ant | 3.89896461 | -25 | 21.1 | 1 |
9 | Ant | 2.295202494 | -26 | 9 | 1 |
9 | Ant | 3.272687389 | -24.5 | 4.95 | 0 |
9 | Ant | 4.201883396 | -13 | 2.7 | 0 |
9 | Post | 4.374845784 | -25 | 10.8 | 0 |
9 | Post | 4.459715564 | -15.5 | 11.8 | 1 |
9 | Post | 3.227763303 | -19.5 | 21.05 | 0 |
9 | Post | 2.517233031 | -20.5 | 31.6 | 1 |
This may serve as a starting point for anyone wishing to tackle this question:
data have;
infile datalines truncover;
input ID Position $ Time Resistance Force y;
datalines;
1 Ant 3.82 -23 37 1
1 Ant 4.94 -24 4.4 0
1 Ant 3.49 -41.5 15.8 1
1 Ant 3.07 -39.5 21.15 1
1 Post 3.43 -39 29 1
1 Post 3.53 -45.5 14.15 1
1 Post 4.44 -46 9.55 1
1 Post 3.37 -19 12.9 1
1 Ant 3.35 -46.5 7.2 1
1 Ant 3.19 -43 14.2 1
1 Ant 3.61 -41 24.55 1
1 Ant 4.24 -48 23.15 1
1 Post 2.09 -33 27.25 1
1 Post 2.83 -32 21.2 1
1 Post 3.26 -28 29.7 0
1 Post 3.34 -41.5 10.15 1
2 Ant 5.29 -30 10.9 1
2 Ant 4.22 -29 11.3 1
2 Ant 2.40 -15 10.2 1
2 Ant 3.32 -18.5 8.75 1
2 Post 1.85 -27 9.3 1
2 Post 4.31 -37 8.75 1
2 Post 1.82 -29.5 13.95 1
2 Post 3.98 -24.5 9.8 0
2 Ant 4.06 -39.5 21.45 1
2 Ant 2.90 -35.5 16.5 1
2 Ant 3.23 -35.5 18.2 1
2 Ant 4.24 -31 14.3 1
2 Post 2.45 -31 9.6 1
2 Post 3.51 -20 6 1
2 Post 4.27 -17.5 8.4 1
2 Post 2.67 -25.5 25 0
3 Ant 3.065092996 -40 10.7 1
3 Ant 3.74 -38 17.8 .
3 Ant 3.61 -27 10.1 0
3 Ant 2.08 -26.5 6.45 .
3 Post 2.12 -35 20.4 1
3 Post 3.244 -39 27.5 1
3 Post 4.02 -42 19.9 1
3 Post 1.94 -19 16.6 1
3 Ant 4.37 -14 4.2 0
3 Ant 4.68 -33 6.9 0
3 Ant 3.35 -30.5 8.65 1
3 Ant 1.72 -33 14.1 1
3 Post 0.81 -27 12.2 1
3 Post 3.90 -26 18.35 1
3 Post 4.19 -29 9.4 1
3 Post 4.46 -19 10.2 1
4 Ant 2.89 -42 11.3 1
4 Ant 2.20 -28 12.45 1
4 Ant 2.97 -31 19.5 1
4 Ant 2.06 -31 22.3 1
4 Post 3.35 -44.5 32.9 1
4 Post 2.10 -35 15.3 1
4 Post 3.42 -35 8.35 1
4 Post 4.16 -33 20.9 1
4 Ant 4.06 -15.5 6 1
4 Ant 5.00 -25 21.5 1
4 Ant 4.13 -33.5 24.25 1
4 Ant 5.56 -34 16.7 1
4 Post 4.14 -35 31.75 1
4 Post 4.49 -33.5 25.4 1
4 Post 4.17 -29 41.9 1
4 Post 3.85 -28 28.8 1
5 Ant 2.67 -23 28.2 0
5 Ant 1.68 -23 10.3 1
5 Ant 2.07 -19.5 9.85 1
5 Ant 1.06 -25 12.7 1
5 Post 5.02 -31 10.4 0
5 Post 2.53 -23 11.7 1
5 Post 3.40 -71.5 27.15 1
5 Post 4.78 -41.5 31.85 1
5 Ant 2.15 -42 15.95 1
5 Ant 2.89 -26.5 16.9 1
5 Ant 2.25 -33 9.8 1
5 Ant 2.76 -28 12.1 0
5 Post 3.04 -22 9 1
5 Post 4.45 -37.5 9.25 1
5 Post 4.10 -20.5 15.3 1
5 Post 4.41 -31.5 18.05 1
6 Ant 1.61 -26 18.7 1
6 Ant 1.68 -26 7.4 0
6 Ant 3.93 -29 7.4 0
6 Ant 4.45 -21.5 9.15 .
6 Post 5.48 -28 25.05 1
6 Post 4.11 -48 30.7 1
6 Post 3.20 -23 22.3 1
6 Post 2.77 -28 15.3 1
6 Ant 2.34 -22 10.95 1
6 Ant 2.25 -30 6.2 1
6 Ant 4.16 -24.5 6.6
6 Ant 4.62 -33 12 0
6 Post 2.32 -31 16.65 1
6 Post 4.03 -31 19 1
6 Post 3.43 -17 12 .
6 Post 3.51 -14 11.1 0
7 Ant 2.99 -30 7.65 1
7 Ant 1.80517419 -25 15.95 1
7 Ant 2.106053494 -32 13.8 1
7 Ant 3.096114016 -29 17.55 1
7 Post 3.167542074 -24 11.45 1
7 Post 3.338268984 -24 22.7 1
7 Post 2.659685183 -32.5 21.95 1
7 Post 3.751749917 -17 8.25 1
7 Ant 2.197529839 -25 8.8 0
7 Ant 3.664137015 -35 10.1 0
7 Ant 3.545702335 -39 9.4 0
7 Ant 2.023625001 -35 10.1 1
7 Post 2.372086883 -33 24.25 1
7 Post 3.582104056 -37.5 17.3 1
7 Post 3.45055168 -38 12.6 1
7 Post 3.841677068 -23 9.8 1
8 Ant 3.432008272 -25 18.6 0
8 Ant 2.136605495 -32.5 14.1 1
8 Ant 1.841190586 -18 7 0
8 Ant 2.15865836 -25.5 6.8 0
8 Post 3.359805409 -30 15.5 1
8 Post 3.631259765 -31 28 1
8 Post 4.674356585 -32.5 20.9 1
8 Post 4.044977037 -25 12.9 0
8 Ant 3.346860731 -28 14.6 1
8 Ant 3.850582629 -46 24.5 1
8 Ant 5.340021635 -31.5 8.5 1
8 Ant 3.980653721 -26 11.6 1
8 Post 3.704121331 -26 14.8 1
8 Post 3.848852913 -31.5 15.8 1
8 Post 4.939191479 -18 9.8 1
8 Post 3.134066196 -15 9.1 1
9 Ant 1.309024248 -25 9 1
9 Ant 3.369404446 -27 11.55 1
9 Ant 1.841284373 -26 13.15 1
9 Ant 3.675231524 -21 12.4 1
9 Post 3.06826061 -22 24.2 1
9 Post 3.59966626 -19 17.3 1
9 Post 4.466268907 -16.5 15.7 1
9 Post 1.882204503 -27 11.3 .
9 Ant 3.89896461 -25 21.1 1
9 Ant 2.295202494 -26 9 1
9 Ant 3.272687389 -24.5 4.95 0
9 Ant 4.201883396 -13 2.7 0
9 Post 4.374845784 -25 10.8 0
9 Post 4.459715564 -15.5 11.8 1
9 Post 3.227763303 -19.5 21.05 0
9 Post 2.517233031 -20.5 31.6 1
;
proc sort data=have; by position id; run;
proc sgscatter data=have;
where resistance gt -50; /* remove outlier */
by position;
matrix force time resistance / group=y;
run;
I would start by performing a logistic regression. (PROC LOGISTIC)
It is repeated measured model, therefore you need Mixed Model of Logistic.
@SteveDenham @lvm @StatDave could give you right sas code .
But you better to post it at Stat forum.
https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.