Hello. I have a question about interpreting results from a survival analysis with a circular variable (ang) which runs from -180 to +180 deg.
I generated the sin and cos of the angle and ran a regression (this is cox regression with time updated variables)
Cox regression with Breslow method for ties
No. of subjects = 411 Number of obs = 3,297
No. of failures = 62
Time at risk = 337,274
LR chi2(8) = 23.53
Log likelihood = -321.75789 Prob > chi2 = 0.0027
--------------------------------------------------------------------------------------------
_t | Haz. ratio Std. err. z P>|z| [95% conf. interval]
---------------------------+----------------------------------------------------------------
cos_ang | .3659565 .129723 -2.84 0.005 .1826831 .7330954
sin_ang | 1.712673 .583607 1.58 0.114 .8782552 3.339859
male | 1.161865 .3123606 0.56 0.577 .6859862 1.967869
age | 1.001317 .0093995 0.14 0.889 .9830623 1.01991
cad | .8453366 .3206804 -0.44 0.658 .4019059 1.778013
dm | 1.017383 .007039 2.49 0.013 1.00368 1.031273
dose | 1.000743 .0012915 0.58 0.565 .9982152 1.003278
htn | 1.00317 .0084401 0.38 0.707 .9867635 1.019849
--------------------------------------------------------------------------------------------
The cosine term is significant but the sin term is not.
Based on some other reading and forum posts on circular statistics I understand that since sin and cos are interrelated, its best to consider them together rather than individually.
Is the best way to do this to run a F-test?
test cos_ang sin_ang
( 1) cos_ang = 0
( 2) sin_ang = 0
chi2( 2) = 16.12
Prob > chi2 = 0.0003
The F-test is very significant. Does this mean that I can "ignore" that non significant P value for the sine term and conclude that there is indeed a circular variation in HR based on the angle?
Is this case any different than if both sin and cos had significant individual P values in addition to a significant joint test?
My goal is to graph the HR as a function of the angle using the regression coefficients
thanks for your insight!
hpw
I generated the sin and cos of the angle and ran a regression (this is cox regression with time updated variables)
Cox regression with Breslow method for ties
No. of subjects = 411 Number of obs = 3,297
No. of failures = 62
Time at risk = 337,274
LR chi2(8) = 23.53
Log likelihood = -321.75789 Prob > chi2 = 0.0027
--------------------------------------------------------------------------------------------
_t | Haz. ratio Std. err. z P>|z| [95% conf. interval]
---------------------------+----------------------------------------------------------------
cos_ang | .3659565 .129723 -2.84 0.005 .1826831 .7330954
sin_ang | 1.712673 .583607 1.58 0.114 .8782552 3.339859
male | 1.161865 .3123606 0.56 0.577 .6859862 1.967869
age | 1.001317 .0093995 0.14 0.889 .9830623 1.01991
cad | .8453366 .3206804 -0.44 0.658 .4019059 1.778013
dm | 1.017383 .007039 2.49 0.013 1.00368 1.031273
dose | 1.000743 .0012915 0.58 0.565 .9982152 1.003278
htn | 1.00317 .0084401 0.38 0.707 .9867635 1.019849
--------------------------------------------------------------------------------------------
The cosine term is significant but the sin term is not.
Based on some other reading and forum posts on circular statistics I understand that since sin and cos are interrelated, its best to consider them together rather than individually.
Is the best way to do this to run a F-test?
test cos_ang sin_ang
( 1) cos_ang = 0
( 2) sin_ang = 0
chi2( 2) = 16.12
Prob > chi2 = 0.0003
The F-test is very significant. Does this mean that I can "ignore" that non significant P value for the sine term and conclude that there is indeed a circular variation in HR based on the angle?
Is this case any different than if both sin and cos had significant individual P values in addition to a significant joint test?
My goal is to graph the HR as a function of the angle using the regression coefficients
thanks for your insight!
hpw
Comment