我需要找到手握力量在年龄,性别,体重和身高方面标准化的最低三分位的受访者比例。
首先,我尝试将grip
,age
,male
和height
的{{1}}强度归一化。其次,我尝试使用此标准化weight
强度变量的tertiles
创建变量。最后,如果一个人处于较低的三分位数,我会尝试创建一个等于grip
的假人。
到目前为止,我已经构建了这段代码,但它并没有真正起作用:
1
数据可能如下所示:
df <- within(df, normal <- ave(grip, male, age, weight, height ,FUN=function(x) (x-min(x))/diff(range(x))))
df$tertile <- ave(df$normal,
FUN=function(x){cut(x, labels=1:3, breaks=quantile(x, probs = 0:3/3, na.rm = TRUE), include.lowest=TRUE)})
df$lowgrip <- ifelse(df$tertile==1, 1, 0)
然而,我的真实数据有大约7000个观察结果。