zhangyanhui-web
commented
3 years ago hello，Sorry to be a bother. I found that author code calculates s_ r is missing kjkj s_r = ((1-(rchar^2))/(rchar^2))(((age.range[2]-age.range[1])^2)/(12*nrow(agev))) I don't know if I'm right. If possible, can you give me a reply？ Thank you very much
Hello, I met with a problem in the results of kdm_calc function.
Under the kdmfun.R, s2 represents the varience of the difference between estimated biological age(Be) and chronological age(C), and s_r represents the varience of the difference between estimated biological age(Be) and ideal biological age(B). The difference between s2 and s_r is used as the sba2(the varience of the ideal biological age(B)) by default to calculate Bec(another biological age adjusted by C based on Be).
However, there is a problem, when s_r is bigger than s2, sba2 gets a negative value. Then, the Bec will be calculated in a different result compared with Be. In detail, ( Bec-C ) has a opposite sign with ( Be-C ). For example, people, who is 45 years old in reality, his Be is 35 years old, but his Bec would be 55 years old! His ba.ca( Be-C )<0 but his baaccel_diff( Bec-C)>0.
Could you give me some ideas to solve this problem?